Deciphering the Structural Enigma of HLA Class-II Binding Peptides for Enhanced Immunoinformatics-based Prediction of Vaccine Epitopes.

Vaccines remain the most efficacious means to avoid and eliminate morbid diseases associated with high morbidity and mortality. Clinical trials indicate the gaining impetus of peptide vaccines against diseases for which an effective treatment still remains obscure. CD4 T-cell-based peptide vaccines involve immunization with antigenic determinants from pathogens or neoplastic cells that possess the ability to elicit a robust T helper cell response, which subsequently activates other arms of the immune system. The available in silico predictors of human leukocyte antigen II (HLA-II) binding peptides are sequence-based techniques, which ostensibly have balanced sensitivity and specificity. Structural analysis and understanding of the cognate peptide and HLA-II interactions are essential to empirically derive a successful peptide vaccine. However, the availability of structure-based epitope prediction algorithms is inadequate compared with sequence-based prediction methods. The present study is an attempt to understand the structural aspects of HLA-II binders by analyzing the Protein Data Bank (PDB) complexes of pHLA-II. Furthermore, we mimic the peptide exchange mechanism and demonstrate the structural implication of an acidic environment on HLA-II binders. Finally, we discuss a structure-guided approach to decipher potential HLA-II binders within an antigenic protein. This strategy may accurately predict the peptide epitopes and thus aid in designing successful peptide vaccines.

Sugar Puckering Drives G-Quadruplex Refolding: Implications for V-Shaped Loops.

A DNA G-quadruplex adopting a (3+1) hybrid structure was modified in two adjacent syn positions of the antiparallel strand with anti-favoring 2'-deoxy-2'-fluoro-riboguanosine (F rG) analogues. The two substitutions promoted a structural rearrangement to a topology with the 5'-terminal G residue located in the central tetrad and the two modified residues linked by a V-shaped zero-nucleotide loop. Strikingly, whereas a sugar pucker in the preferred north domain is found for both modified nucleotides, the F rG analogue preceding the V-loop is forced to adopt the unfavored syn conformation in the new quadruplex fold. Apparently, a preferred C3'-endo sugar pucker within the V-loop architecture outweighs the propensity of the F rG analogue to adopt an anti glycosidic conformation. Refolding into a V-loop topology is likewise observed for a sequence modified at corresponding positions with two riboguanosine substitutions. In contrast, 2'-F-arabinoguanosine analogues with their favored south-east sugar conformation do not support formation of the V-loop topology. Examination of known G-quadruplexes with a V-shaped loop highlights the critical role of the sugar conformation for this distinct structural motif.

LCS-TA to identify similar fragments in RNA 3D structures.

BACKGROUND: In modern structural bioinformatics, comparison of molecular structures aimed to identify and assess similarities and differences between them is one of the most commonly performed procedures. It gives the basis for evaluation of in silico predicted models. It constitutes the preliminary step in searching for structural motifs. In particular, it supports tracing the molecular evolution. Faced with an ever-increasing amount of available structural data, researchers need a range of methods enabling comparative analysis of the structures from either global or local perspective. RESULTS: Herein, we present a new, superposition-independent method which processes pairs of RNA 3D structures to identify their local similarities. The similarity is considered in the context of structure bending and bonds' rotation which are described by torsion angles. In the analyzed RNA structures, the method finds the longest continuous segments that show similar torsion within a user-defined threshold. The length of the segment is provided as local similarity measure. The method has been implemented as LCS-TA algorithm (Longest Continuous Segments in Torsion Angle space) and is incorporated into our MCQ4Structures application, freely available for download from http://www.cs.put.poznan.pl/tzok/mcq/ . CONCLUSIONS: The presented approach ties torsion-angle-based method of structure analysis with the idea of local similarity identification by handling continuous 3D structure segments. The first method, implemented in MCQ4Structures, has been successfully utilized in RNA-Puzzles initiative. The second one, originally applied in Euclidean space, is a component of LGA (Local-Global Alignment) algorithm commonly used in assessing protein models submitted to CASP. This unique combination of concepts implemented in LCS-TA provides a new perspective on structure quality assessment in local and quantitative aspect. A series of computational experiments show the first results of applying our method to comparison of RNA 3D models. LCS-TA can be used for identifying strengths and weaknesses in the prediction of RNA tertiary structures.

C9/12 Ribbon-Like Structures in Hybrid Peptides Alternating α- and Thiazole-Based γ-Amino Acids.

According to their restricted conformational freedom, heterocyclic γ-amino acids are usually considered to be related to Z-vinylogous γ-amino acids. In this context, oligomers alternating α-amino acids and thiazole-based γ-amino acids (ATCs) were expected to fold into a canonical 12-helical shape as described for α/γ-hybrid peptides composed of cis-α/ß-unsaturated γ-amino acids. However, through a combination of X-ray crystallography, NMR spectroscopy, FTIR experiments, and DFT calculations, it was determined that the folding behavior of ATC-containing hybrid peptides is much more complex. The homochiral α/(S)-ATC sequences were unable to adopt a stable conformation, whereas the heterochiral α/(R)-ATC peptides displayed novel ribbon structures stabilized by unusual C9/12 -bifurcated hydrogen bonds. These ribbon structures could be considered as a succession of pre-organized γ/α dipeptides and may provide the basis for designing original α-helix mimics.

Lowering the reduction potential of a boron compound by means of the substituent effect of the boryl group: one-electron reduction of an unsymmetrical diborane(4).

We have clarified and observed the high electron affinity of pinB-BMes2 (1; Mes = mesityl, pin = pinacolato). By using electrochemistry, it was shown that 1 has a higher electron affinity than those of B2pin2 and Mes3B. One-electron reduction of 1 gave the corresponding radical anion. The ESR spectroscopy and DFT calculation revealed the unsymmetrical distribution of electron density over the B-B bond. UV/Vis spectroscopy showed that the SOMO-related absorption supports the deep purple color of the radical anion. DFT studies on the torsion angle dependency of the LUMO levels and relative energies revealed the reason why 1 has high electron affinity as a result of the substituent effect of the Bpin group.

Mechanosensitive membrane probes.

This article assembles pertinent insights behind the concept of planarizable push-pull probes. As a response to the planarization of their polarized ground state, a red shift of their excitation maximum is expected to report on either the disorder, the tension, or the potential of biomembranes. The combination of chromophore planarization and polarization contributes to various, usually more complex processes in nature. Examples include the color change of crabs or lobsters during cooking or the chemistry of vision, particularly color vision. The summary of lessons from nature is followed by an overview of mechanosensitive organic materials. Although often twisted and sometimes also polarized, their change of color under pressure usually originates from changes in their crystal packing. Intriguing exceptions include the planarization of several elegantly twisted phenylethynyl oligomers and polymers. Also mechanosensitive probes in plastics usually respond to stretching by disassembly. True ground-state planarization in response to molecular recognition is best exemplified with the binding of thoughtfully twisted cationic polythiophenes to single- and double-stranded oligonucleotides. Molecular rotors, en vogue as viscosity sensors in cells, operate by deplanarization of the first excited state. Pertinent recent examples are described, focusing on λ-ratiometry and intracellular targeting. Complementary to planarization of the ground state with twisted push-pull probes, molecular rotors report on environmental changes with quenching or shifts in emission rather than absorption. The labeling of mechanosensitive channels is discussed as a bioengineering approach to bypass the challenge to create molecular mechanosensitivity and use biological systems instead to sense membrane tension. With planarizable push-pull probes, this challenge is met not with twistome screening, but with "fluorescent flippers," a new concept to insert large and bright monomers into oligomeric probes to really feel the environment and also shine when twisted out of conjugation.

High accuracy of Karplus equations for relating three-bond J couplings to protein backbone torsion angles.

(3) JC'C' and (3) JHNHα couplings are related to the intervening backbone torsion angle ${\varphi }$ by standard Karplus equations. Although these couplings are known to be affected by parameters other than ${\varphi }$, including H-bonding, valence angles and residue type, experimental results and quantum calculations indicate that the impact of these latter parameters is typically very small. The solution NMR structure of protein GB3, newly refined by using extensive sets of residual dipolar couplings, yields 50-60 % better Karplus equation agreement between ${\varphi }$ angles and experimental (3) JC'C' and (3) JHNHα values than does the high-resolution X-ray structure. In intrinsically disordered proteins, (3) JC'C' and (3) JHNHα couplings can be measured at even higher accuracy, and the impact of factors other than the intervening torsion angle on (3) J will be smaller than in folded proteins, making these couplings exceptionally valuable reporters on the ensemble of ${\varphi }$ angles sampled by each residue.

Emerging role of carbonyl-carbonyl interactions in the classification of beta turns.

Carbonyl-carbonyl interactions in peptides and proteins attracted considerable interest in recent years. Here, we report a survey of carbonyl-carbonyl interactions in cyclic peptides, depsipeptides, peptoids and discuss the relationship between backbone torsion angles and CO∙∙∙CO distances. In general, φ values in the range between -40° and -90° and between 40° and 90° correspond to CO∙∙∙CO distances below 3.22 Å. By extending the analysis of carbonyl-carbonyl interactions in different types of beta turns in proteins, we also highlight the role of direct or reciprocal carbonyl-carbonyl interactions in stabilizing the beta turn conformation for each specific type. We confirmed the new type II beta turn, detected by Dunbrack and coworkers, and named Pa, and detect the presence of a direct carbonyl-carbonyl interaction between the second and third residues of the turn. We also evidenced the existence of another new type II beta turn, named pA (following Dunbrack's notation), which represents the alternative conformation of Pa with opposite φ and ψ values and is characterized by a direct carbonyl-carbonyl interaction between the second and third residues of the turn. Finally, we show that the occurrence of CO∙∙∙CO interactions could be also advocated to explain from a chemical point of view the diversity of turn types.

Analysis and validation of overall N-glycan conformation in Privateer.

The oligosaccharides in N-glycosylation provide key structural and functional contributions to a glycoprotein. These contributions are dependent on the composition and overall conformation of the glycans. The Privateer software allows structural biologists to evaluate and improve the atomic structures of carbohydrates, including N-glycans; this software has recently been extended to check glycan composition through the use of glycomics data. Here, a broadening of the scope of the software to analyse and validate the overall conformation of N-glycans is presented, focusing on a newly compiled set of glycosidic linkage torsional preferences harvested from a curated set of glycoprotein models.

A new circle method for measuring humeral torsion on MRI-scans less sensitive to Hill-Sachs lesions.

Objectives: The literature on humeral torsion angles (retrotorsion) reveals great inconsistencies between methodology and values. Decreased retrotorsion was suspected to correlate with instability, but evidence is contradictory. The measurement according to the gold standard method of Bernageau and Godefroy (B&G) can be challenging especially in the presence of Hill-Sachs-lesions. Therefore, we have developed and evaluated a new measurement method for the humeral torsion angle on MRI-scans. Materials and Methods: Three investigators have measured 67 patients (35 with shoulder instability, 32 healthy) on axial MRIs with 603 measurements used for reliability calculation. The new Circle-method determines the retrotorsion by overlaying two circles on the transversal section of the humeral head. The first circle is adjusted congruent with the margin of the humeral head, whereas the second circle is adjusted to the greater tubercle. The line bisecting the centres of these circles is defined as the humeral head axis. This method was compared to B&G. Results: The mean retrotorsion angle of all patients was 25°± 25° (mean ± SD) with B&G, and 24° ± 27° with the Circle-method. Neither method revealed a significant difference between stable and unstable shoulders (p = 0.47). Of the 35 patients with unstable shoulders 21 (60%) presented Hill-Sachs lesions. No significant differences between patients with or without Hill-Sachs lesions (Circle-method: p = 0.61; B&G: p = 0.67). The reliability parameters for both methods were similar. Conclusions: The new Circle-method is as precise as the method of B&G. It may yield more consistent values in cases with substantial Hill-Sachs-lesions. Our data do not suggest retrotorsion as a predictor of instability.

Annotating the protein-RNA interaction sites in proteins using evolutionary information and protein backbone structure.

RNA-protein interactions play important roles in various biological processes. The precise detection of RNA-protein interaction sites is very important for understanding essential biological processes and annotating the function of the proteins. In this study, based on various features from amino acid sequence and structure, including evolutionary information, solvent accessible surface area and torsion angles (φ, ψ) in the backbone structure of the polypeptide chain, a computational method for predicting RNA-binding sites in proteins is proposed. When the method is applied to predict RNA-binding sites in three datasets: RBP86 containing 86 protein chains, RBP107 containing 107 proteins chains and RBP109 containing 109 proteins chains, better sensitivities and specificities are obtained compared to previously published methods in five-fold cross-validation tests. In order to make further examination for the efficiency of our method, the RBP107 dataset is used as training set, RBP86 and RBP109 datasets are used as the independent test sets. In addition, as examples of our prediction, RNA-binding sites in a few proteins are presented. The annotated results are consistent with the PDB annotation. These results show that our method is useful for annotating RNA binding sites of novel proteins.

Local Backbone Geometry Plays a Critical Role in Determining Conformational Preferences of Amino Acid Residues in Proteins.

The definition of the structural basis of the conformational preferences of the genetically encoded amino acid residues is an important yet unresolved issue of structural biology. In order to gain insights into this intricate topic, we here determined and compared the amino acid propensity scales for different (φ, ψ) regions of the Ramachandran plot and for different secondary structure elements. These propensities were calculated using the Chou-Fasman approach on a database of non-redundant protein chains retrieved from the Protein Data Bank. Similarities between propensity scales were evaluated by linear regression analyses. One of the most striking and unexpected findings is that distant regions of the Ramachandran plot may exhibit significantly similar propensity scales. On the other hand, contiguous regions of the Ramachandran plot may present anticorrelated propensities. In order to provide an interpretative background to these results, we evaluated the role that the local variability of protein backbone geometry plays in this context. Our analysis indicates that (dis)similarities of propensity scales between different regions of the Ramachandran plot are coupled with (dis)similarities in the local geometry. The concept that similarities of the propensity scales are dictated by the similarity of the NCαC angle and not necessarily by the similarity of the (φ, ψ) conformation may have far-reaching implications in the field.

An Effective Cumulative Torsion Angles Model for Prediction of Protein Folding Rates.

BACKGROUND: Protein folding rate is mainly determined by the size of the conformational space to search, which in turn is dictated by factors such as size, structure and amino-acid sequence in a protein. It is important to integrate these factors effectively to form a more precisely description of conformation space. But there is no general paradigm to answer this question except some intuitions and empirical rules. Therefore, at the present stage, predictions of the folding rate can be improved through finding new factors, and some insights are given to the above question. OBJECTIVE: Its purpose is to propose a new parameter that can describe the size of the conformational space to improve the prediction accuracy of protein folding rate. METHODS: Based on the optimal set of amino acids in a protein, an effective cumulative backbone torsion angles (CBTAeff) was proposed to describe the size of the conformational space. Linear regression model was used to predict protein folding rate with CBTAeff as a parameter. The degree of correlation was described by the coefficient of determination and the mean absolute error MAE between the predicted folding rates and experimental observations. RESULTS: It achieved a high correlation (with the coefficient of determination of 0.70 and MAE of 1.88) between the logarithm of folding rates and the (CBTAeff)0.5 with experimental over 112 twoand multi-state folding proteins. CONCLUSION: The remarkable performance of our simplistic model demonstrates that CBTA based on optimal set was the major determinants of the conformation space of natural proteins.

SPIDER2: A Package to Predict Secondary Structure, Accessible Surface Area, and Main-Chain Torsional Angles by Deep Neural Networks.

Predicting one-dimensional structure properties has played an important role to improve prediction of protein three-dimensional structures and functions. The most commonly predicted properties are secondary structure and accessible surface area (ASA) representing local and nonlocal structural characteristics, respectively. Secondary structure prediction is further complemented by prediction of continuous main-chain torsional angles. Here we describe a newly developed method SPIDER2 that utilizes three iterations of deep learning neural networks to improve the prediction accuracy of several structural properties simultaneously. For an independent test set of 1199 proteins SPIDER2 achieves 82 % accuracy for secondary structure prediction, 0.76 for the correlation coefficient between predicted and actual solvent accessible surface area, 19° and 30° for mean absolute errors of backbone φ and ψ angles, respectively, and 8° and 32° for mean absolute errors of Cα-based θ and τ angles, respectively. The method provides state-of-the-art, all-in-one accurate prediction of local structure and solvent accessible surface area. The method is implemented, as a webserver along with a standalone package that are available in our website: http://sparks-lab.org .

Factors affecting the amplitude of the τ angle in proteins: a revisitation.

The protein folded state is the result of the fine balance of a variety of different forces. Even minor structural perturbations may have a significant impact on the stability of these macromolecules. Studies carried out in recent decades have led to the convergent view that proteins are endowed with a flexible spine. One of the open issues related to protein local backbone geometry is the identification of the factors that influence the amplitude of the τ (N-Cα-C) angle. Here, statistical analyses performed on an updated ensemble of X-ray protein structures by dissecting the contribution of the major factors that can potentially influence the local backbone geometry of proteins are reported. The data clearly indicate that the local backbone conformation has a prominent impact on the modulation of the τ angle. Therefore, a proper assessment of the impact of the other potential factors can only be appropriately evaluated when small (ϕ, ψ) regions are considered. Here, it is shown that when the contribution of the backbone conformation is removed by considering small (ϕ, ψ) areas, an impact of secondary structure, as defined by DSSP, and/or the residue type on τ is still detectable, although to a limited extent. Indeed, distinct τ-value distributions are detected for Pro/Gly and ß-branched (Ile/Val) residues. The key role of the local backbone conformation highlighted here supports the use of variable local backbone geometry in protein refinement protocols.

A user-friendly web portal for analyzing conformational changes in structures of Mycobacterium tuberculosis.

Initiation of the Tuberculosis Structural Consortium has resulted in the expansion of the Mycobacterium tuberculosis (MTB) protein structural database. Currently, 969 experimentally solved structures are available for 354 MTB proteins. This includes multiple crystal structures for a given protein under different functional conditions, such as the presence of different ligands or mutations. In depth analysis of the multiple structures reveal that subtle differences exist in conformations of a given protein under varied conditions. Therefore, it is immensely important to understand the conformational differences between the multiple structures of a given protein in order to select the most suitable structure for molecular docking and structure-based drug designing. Here, we introduce a web portal ( http://bmi.icmr.org.in/mtbsd/torsion.php ) that we developed to provide comparative data on the ensemble of available structures of MTB proteins, such as Cα root means square deviation (RMSD), sequence identity, presence of mutations and torsion angles. Additionally, torsion angles were used to perform principal component analysis (PCA) to identify the conformational differences between the structures. Additionally, we present a few case studies to demonstrate this database. Graphical Abstract Conformational changes seen in the structures of the enoyl-ACP reductase protein encoded by the Mycobacterial gene inhA.

Design and Synthesis of Mixed Oligomers with Thiophenes, Dithienothiophene S,S-Dioxides, Thieno[3,4]pyrazines and 2,1,3-Benzothiadiazoles: Flipper Screening for Mechanosensitive Systems.

Monomers with large surface area and high quantum yield, that is fluorescent flippers, have been engineered into twisted push-pull oligomers to create membrane probes with high mechanosensitivity and long fluorescence lifetime. Here, the synthesis and characterization of thieno[3,4]pyrazines and 2,1,3-benzothiadiazoles are described in comparison with the original dithienothiophene S,S-dioxides. Dithienothiophene S,S-dioxide flippers are confirmed as the best reported so far, and poor results with single flipper probes support that two flippers are needed for the probe to really "swim", that is, for high mechanosensitivity.

The gramicidin channel ion permeation free-energy profile: direct and indirect effects of CHARMM force field improvements.

A revised CHARMM force field for tryptophan residues is studied as well as a new grid-based correction algorithm, called CMAP, using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. The conformational stability of the interfacial side chains is studied, which shows good stability on the 10 ns time scale. The revised force field for the tryptophan side chain produces, in the decomposition, a Na(+) PMF(Trp) profile that is consonant with the prediction from the experimental results, analyzed with rate theory by Durrant et al. (2006), but in stark contrast to the prediction of the original CHARMM force field, version 22. However, the effect is diluted in the PMF profile due to indirect effects mediated by other components of the system (polypeptide, lipid molecules, ions, and water molecules). CMAP corrections to the L-amino acids help reduce the excessive translocation barrier. Decomposition demonstrates that this effect is due to effects on the K(+) PMF(H(2)O) profile rather than on the K(+) PMF(gA) profile. The results have been confirmed to be robust using an alternative umbrella-potential method. Further force field balancing efforts (direct and indirect) are required for future studies to evaluate whether these effects give rise to predictions that are consistent with those observables extracted from real experiments.