Conservation of cis prolyl bonds in proteins during evolution

Insights into the solution structure and transcriptional regulation of the MazE9 antitoxin in Mycobacterium tuberculosis

The present study endeavors to decode the details of the transcriptional autoregulation effected by the MazE9 antitoxin of the Mycobacterium tuberculosis MazEF9 toxin-antitoxin system. Regulation of this bicistronic operon at the level of transcription is a critical biochemical process that is key for the organism's stress adaptation and virulence. Here, we have reported the solution structure of the DNA binding domain of MazE9 and scrutinized the thermodynamic and kinetic parameters operational in its interaction with the promoter/operator region, specific to the mazEF9 operon. A HADDOCK model of MazE9 bound to its operator DNA has been calculated based on the information on interacting residues obtained from these studies. The thermodynamics and kinetics of the interaction of MazE9 with the functionally related mazEF6 operon indicate that the potential for intracellular cross-regulation is unlikely. An interesting feature of MazE9 is the cis ⇌ trans conformational isomerization of proline residues in the intrinsically disordered C-terminal domain of this antitoxin.

Distance-based global analysis of consistent cis-bonds in protein backbones

Biological polypeptides are known to contain cis-linkage in their main chain as a minor but important feature. Such anomalous connection of amino acids has different structural and functional effects on proteins. Experimental evidence of cis-bonds in proteins is mainly obtained using X-ray crystallography and other methods in the field of structural biology. To date, extensive analyses have been carried out on the experimentally found cis-bonds using the Protein Data Bank (PDB) entry-wise or residue-wise; however, their consistency in each protein has not been examined on a global scale. Data accumulation and advances in computational methodology enable the use of new approaches from a proteomic point of view. Here, we sought to carry out protein-wise analysis and describe a simple procedure for the detection and confirmation of cis-bonds from a set of experimental PDB chains for a protein to discriminate this type of bond from isomerizable and/or misassigned bonds. The resulting set of consistent cis bonds (found at identical positions in multiple chains) provides unprecedented insights into the trend of "high cis content" proteins and the upper limit of consistent cis bonds per polypeptide length. Recognizing such limit would not only be important for a practical check of upcoming structures, but also for the design of novel protein folds beyond the evolutionally-acquired repertoire.

Synthesis and Conformational Analysis of N-Aromatic Acetamides Bearing Thiophene: Effect of Intramolecular Chalcogen-Chalcogen Interaction on Amide Conformational Stability

The conformations of aromatic amides bearing an N-(2-thienyl) or N-(3-thienyl) group were investigated in solution and in the crystal state. NMR spectral data indicate that the conformational preferences of these amides in solution are dependent not only on the relative π-electron densities of the N-aromatic moieties, but also on the three-dimensional relationship between carbonyl oxygen and the N-aromatic moieties. A comparison of the conformational preferences of N-(2-thienyl)amides and N-(3-thienyl)amides revealed that the Z-conformers of N-(2-thienyl)acetamides are stabilized by 1,5-type intramolecular S···O═C interactions between amide carbonyl and thiophene sulfur. The crystal structures of these compounds were similar to the solution structures. The stabilization energy due to 1,5-type intramolecular S···O═C interaction in N-aryl-N-(2-thienyl)acetamides and N-methyl-N-(2-thienyl)acetamide was estimated to be ca. 0.74 and 0.93 kcal/mol, respectively.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Which DFT levels of theory are appropriate in predicting the prolyl cis–trans isomerization in solution?

DFTs were assessed for the conformational preferences of the peptides containing Pro and its derivatives in chloroform and water.

Structures of two fimbrial adhesins, AtfE and UcaD, from the uropathogenProteus mirabilis

The important uropathogenProteus mirabilisencodes a record number of chaperone/usher-pathway adhesive fimbriae. Such fimbriae, which are used for adhesion to cell surfaces/tissues and for biofilm formation, are typically important virulence factors in bacterial pathogenesis. Here, the structures of the receptor-binding domains of the tip-located two-domain adhesins UcaD (1.5 Å resolution) and AtfE (1.58 Å resolution) from twoP. mirabilisfimbriae (UCA/NAF and ATF) are presented. The structures of UcaD and AtfE are both similar to the F17G type of tip-located fimbrial receptor-binding domains, and the structures are very similar despite having only limited sequence similarity. These structures represent an important step towards a molecular-level understanding ofP. mirabilisfimbrial adhesins and their roles in the complex pathogenesis of urinary-tract infections.

Detecting Proline and Non-Proline Cis Isomers in Protein Structures from Sequences Using Deep Residual Ensemble Learning

It has been long established that cis conformations of amino acid residues play many biologically important roles despite their rare occurrence in protein structure. Because of this rarity, few methods have been developed for predicting cis isomers from protein sequences, most of which are based on outdated datasets and lack the means for independent testing. In this work, using a database of >10000 high-resolution protein structures, we update the statistics of cis isomers and develop a sequence-based prediction technique using an ensemble of residual convolutional and long short-term memory bidirectional recurrent neural networks that allow learning from the whole protein sequence. We show that ensembling eight neural network models yields maximum Matthews correlation coefficient values of approximately 0.35 for cis-Pro isomers and 0.1 for cis-nonPro residues. The method should be useful for prioritizing functionally important residues in cis isomers for experimental validations and improving the sampling of rare protein conformations for ab initio protein structure prediction.

Molecular Dynamics Information Improves cis-Peptide-Based Function Annotation of Proteins

cis-Peptide bonds, whose occurrence in proteins is rare but evolutionarily conserved, are implicated to play an important role in protein function. This has led to their previous use in a homology-independent, fragment-match-based protein function annotation method. However, proteins are not static molecules; dynamics is integral to their activity. This is nicely epitomized by the geometric isomerization of cis-peptide to trans form for molecular activity. Hence we have incorporated both static (cis-peptide) and dynamics information to improve the prediction of protein molecular function. Our results show that cis-peptide information alone cannot detect functional matches in cases where cis-trans isomerization exists but 3D coordinates have been obtained for only the trans isomer or when the cis-peptide bond is incorrectly assigned as trans. On the contrary, use of dynamics information alone includes false-positive matches for cases where fragments with similar secondary structure show similar dynamics, but the proteins do not share a common function. Combining the two methods reduces errors while detecting the true matches, thereby enhancing the utility of our method in function annotation. A combined approach, therefore, opens up new avenues of improving existing automated function annotation methodologies.

Peptide Bond Isomerization in High-Temperature Simulations

Force fields for molecular simulation are generally optimized to model macromolecules such as proteins at ambient temperature and pressure. Nevertheless, elevated temperatures are frequently used to enhance conformational sampling, either during system setup or as a component of an advanced sampling technique such as temperature replica exchange. Because macromolecular force fields are now put upon to simulate temperatures and time scales that greatly exceed their original design specifications, it is appropriate to re-evaluate whether these force fields are up to the task. Here, we quantify the rates of peptide bond isomerization in high-temperature simulations of three octameric peptides and a small fast-folding protein. We show that peptide octamers with and without proline residues undergo cis/trans isomerization every 1-5 ns at 800 K with three classical atomistic force fields (AMBER99SB-ILDN, CHARMM22/CMAP, and OPLS-AA/L). On the low microsecond time scale, these force fields permit isomerization of nonprolyl peptide bonds at temperatures ≥500 K, and the CHARMM22/CMAP force field permits isomerization of prolyl peptide bonds ≥400 K. Moreover, the OPLS-AA/L force field allows chiral inversion about the Cα atom at 800 K. Finally, we show that temperature replica exchange permits cis peptide bonds developed at 540 K to subsequently migrate back to the 300 K ensemble, where cis peptide bonds are present in 2 ± 1% of the population of Trp-cage TC5b, including up to 4% of its folded state. Further work is required to assess the accuracy of cis/trans isomerization in the current generation of protein force fields.

The structure of a prophenoloxidase (PPO) from Anopheles gambiae provides new insights into the mechanism of PPO activation

BACKGROUND

Phenoloxidase (PO)-catalyzed melanization is a universal defense mechanism of insects against pathogenic and parasitic infections. In mosquitos such as Anopheles gambiae, melanotic encapsulation is a resistance mechanism against certain parasites that cause malaria and filariasis. PO is initially synthesized by hemocytes and released into hemolymph as inactive prophenoloxidase (PPO), which is activated by a serine protease cascade upon recognition of foreign invaders. The mechanisms of PPO activation and PO catalysis have been elusive.

RESULTS

Herein, we report the crystal structure of PPO8 from A. gambiae at 2.6 Å resolution. PPO8 forms a homodimer with each subunit displaying a classical type III di-copper active center. Our molecular docking and mutagenesis studies revealed a new substrate-binding site with Glu364 as the catalytic residue responsible for the deprotonation of mono- and di-phenolic substrates. Mutation of Glu364 severely impaired both the monophenol hydroxylase and diphenoloxidase activities of AgPPO8. Our data suggested that the newly identified substrate-binding pocket is the actual site for catalysis, and PPO activation could be achieved without withdrawing the conserved phenylalanine residue that was previously deemed as the substrate 'placeholder'.

CONCLUSIONS

We present the structural and functional data from a mosquito PPO. Our results revealed a novel substrate-binding site with Glu364 identified as the key catalytic residue for PO enzymatic activities. Our data offered a new model for PPO activation at the molecular level, which differs from the canonical mechanism that demands withdrawing a blocking phenylalanine residue from the previously deemed substrate-binding site. This study provides new insights into the mechanisms of PPO activation and enzymatic catalysis of PO.

Identifying functionally important cis-peptide containing segments in proteins and their utility in molecular function annotation

Cis-peptide embedded segments are rare in proteins but often highlight their important role in molecular function when they do occur. The high evolutionary conservation of these segments illustrates this observation almost universally, although no attempt has been made to systematically use this information for the purpose of function annotation. In the present study, we demonstrate how geometric clustering and level-specific Gene Ontology molecular-function terms (also known as annotations) can be used in a statistically significant manner to identify cis-embedded segments in a protein linked to its molecular function. The present study identifies novel cis-peptide fragments, which are subsequently used for fragment-based function annotation. Annotation recall benchmarks interpreted using the receiver-operator characteristic plot returned an area-under-curve > 0.9, corroborating the utility of the annotation method. In addition, we identified cis-peptide fragments occurring in conjunction with functionally important trans-peptide fragments, providing additional insights into molecular function. We further illustrate the applicability of our method in function annotation where homology-based annotation transfer is not possible. The findings of the present study add to the repertoire of function annotation approaches and also facilitate engineering, design and allied studies around the cis-peptide neighborhood of proteins.

Force-dependent isomerization kinetics of a highly conserved proline switch modulates the mechanosensing region of filamin

Proline switches, controlled by cis-trans isomerization, have emerged as a particularly effective regulatory mechanism in a wide range of biological processes. In this study, we use single-molecule mechanical measurements to develop a full kinetic and energetic description of a highly conserved proline switch in the force-sensing domain 20 of human filamin and how prolyl isomerization modulates the force-sensing mechanism. Proline isomerization toggles domain 20 between two conformations. A stable cis conformation with slow unfolding, favoring the autoinhibited closed conformation of filamin's force-sensing domain pair 20-21, and a less stable, uninhibited conformation promoted by the trans form. The data provide detailed insight into the folding mechanisms that underpin the functionality of this binary switch and elucidate its remarkable efficiency in modulating force-sensing, thus combining two previously unconnected regulatory mechanisms, proline switches and mechanosensing.

Tracking molecular recognition at the atomic level with a new protein scaffold based on the OB-fold

The OB-fold is a small, versatile single-domain protein binding module that occurs in all forms of life, where it binds protein, carbohydrate, nucleic acid and small-molecule ligands. We have exploited this natural plasticity to engineer a new class of non-immunoglobulin alternatives to antibodies with unique structural and biophysical characteristics. We present here the engineering of the OB-fold anticodon recognition domain from aspartyl tRNA synthetase taken from the thermophile Pyrobaculum aerophilum. For this single-domain scaffold we have coined the term OBody. Starting from a naïve combinatorial library, we engineered an OBody with 3 nM affinity for hen egg-white lysozyme, by optimising the affinity of a naïve OBody 11,700-fold over several affinity maturation steps, using phage display. At each maturation step a crystal structure of the engineered OBody in complex with hen egg-white lysozyme was determined, showing binding elements in atomic detail. These structures have given us an unprecedented insight into the directed evolution of affinity for a single antigen on the molecular scale. The engineered OBodies retain the high thermal stability of the parental OB-fold despite mutation of up to 22% of their residues. They can be expressed in soluble form and also purified from bacteria at high yields. They also lack disulfide bonds. These data demonstrate the potential of OBodies as a new scaffold for the engineering of specific binding reagents and provide a platform for further development of future OBody-based applications.

Intelligent Consensus Modeling for Proline Cis-Trans Isomerization Prediction

Proline cis-trans isomerization (CTI) plays a key role in the rate-determining steps of protein folding. Accurate prediction of proline CTI is of great importance for the understanding of protein folding, splicing, cell signaling, and transmembrane active transport in both the human body and animals. Our goal is to develop a state-of-the-art proline CTI predictor based on a biophysically motivated intelligent consensus modeling through the use of sequence information only (i.e., position specific scores generated by PSI-BLAST). The current computational proline CTI predictors reach about 70-73 percent Q2 accuracies and about 0.40 Matthew correlation coefficient (Mcc) through the use of sequence-based evolutionary information as well as predicted protein secondary structure information. However, our approach that utilizes a novel decision tree-based consensus model with a powerful randomized-metal earning technique has achieved 86.58 percent Q2 accuracy and 0.74 Mcc, on the same proline CTI data set, which is a better result than those of any existing computational proline CTI predictors reported in the literature.

Local control of cis-peptidyl-prolyl bonds mediated by CH···π interactions: the Xaa-Pro-Tyr motif

Compared to generic peptide bonds, the peptidyl-prolyl bond shows a strong propensity for the cis conformer. The presence of a sequence-contiguous aromatic (Aro) residue can further stabilize the cis conformer, as observed for the Aro-Pro motif. The cis propensity of the reverse sequence motif, Pro-Aro, is not so well understood, especially the effect of N-capping the Pro-Aro motif with different amino acid residues. From a comparative nuclear magnetic resonance study of two peptide series with the general sequences Ac-Xaa-Pro-Tyr-NH2 and Ac-Xaa-Pro-Ala-NH2, we present a relative thermodynamic scale that reflects how the nature of the Xaa side chain influences the cis propensity of the Xaa-Pro-Tyr motif, with Gly, Pro, and Ala at position Xaa giving the greatest enhancement of the cis-peptidyl-prolyl population. We also show that CH···π interaction between Xaa and Tyr is responsible for the enhanced cis population. However, the mere presence of the CH···π interaction does not guarantee that the peptidyl-prolyl bond will have a higher cis content in Xaa-Pro-Tyr than in Xaa-Pro-Ala. Xaa-dependent intramolecular interactions present in Xaa-trans-Pro-Tyr can nullify favorable CH···π interactions in Xaa-cis-Pro-Tyr. The relative cis-peptidyl-prolyl stabilizing propensities of Xaa (Xaa-Pro-Tyr) in proteins and in our peptide series show strong linear correlation except when Xaa is aromatic. We also explore the Xaa-Pro-Gly-Tyr sequence motif and show that mediated by a Pro-Tyr CH···π interaction, the cis-peptidyl-prolyl bond in the motif is stabilized when Xaa is Pro.

Cis-trans isomerization of omega dihedrals in proteins

Peptide bonds in protein structures are mainly found in trans conformation with a torsion angle ω close to 180°. Only a very low proportion is observed in cis conformation with ω angle around 0°. Cis-trans isomerization leads to local conformation changes which play an important role in many biological processes. In this paper, we reviewed the recent discoveries and research achievements in this field. First, we presented some interesting cases of biological processes in which cis-trans isomerization is directly implicated. It is involved in protein folding and various aspect of protein function like dimerization interfaces, autoinhibition control, channel gating, membrane binding. Then we reviewed conservation studies of cis peptide bonds which emphasized evolution constraints in term of sequence and local conformation. Finally we made an overview of the numerous molecular dynamics studies and prediction methodologies already developed to take into account this structural feature in the research area of protein modeling. Many cis peptide bonds have not been recognized as such due to the limited resolution of the data and to the refinement protocol used. Cis-trans proline isomerization reactions represents a vast and promising research area that still needs to be further explored for a better understanding of isomerization mechanism and improvement of cis peptide bond predictions.

Proline editing: a general and practical approach to the synthesis of functionally and structurally diverse peptides. Analysis of steric versus stereoelectronic effects of 4-substituted prolines on conformation within peptides

Functionalized proline residues have diverse applications. Herein we describe a practical approach, proline editing, for the synthesis of peptides with stereospecifically modified proline residues. Peptides are synthesized by standard solid-phase peptide synthesis to incorporate Fmoc-hydroxyproline (4R-Hyp). In an automated manner, the Hyp hydroxyl is protected and the remainder of the peptide synthesized. After peptide synthesis, the Hyp protecting group is orthogonally removed and Hyp selectively modified to generate substituted proline amino acids, with the peptide main chain functioning to "protect" the proline amino and carboxyl groups. In a model tetrapeptide (Ac-TYPN-NH2), 4R-Hyp was stereospecifically converted to 122 different 4-substituted prolyl amino acids, with 4R or 4S stereochemistry, via Mitsunobu, oxidation, reduction, acylation, and substitution reactions. 4-Substituted prolines synthesized via proline editing include incorporated structured amino acid mimetics (Cys, Asp/Glu, Phe, Lys, Arg, pSer/pThr), recognition motifs (biotin, RGD), electron-withdrawing groups to induce stereoelectronic effects (fluoro, nitrobenzoate), handles for heteronuclear NMR ((19)F:fluoro; pentafluorophenyl or perfluoro-tert-butyl ether; 4,4-difluoro; (77)SePh) and other spectroscopies (fluorescence, IR: cyanophenyl ether), leaving groups (sulfonate, halide, NHS, bromoacetate), and other reactive handles (amine, thiol, thioester, ketone, hydroxylamine, maleimide, acrylate, azide, alkene, alkyne, aryl halide, tetrazine, 1,2-aminothiol). Proline editing provides access to these proline derivatives with no solution-phase synthesis. All peptides were analyzed by NMR to identify stereoelectronic and steric effects on conformation. Proline derivatives were synthesized to permit bioorthogonal conjugation reactions, including azide-alkyne, tetrazine-trans-cyclooctene, oxime, reductive amination, native chemical ligation, Suzuki, Sonogashira, cross-metathesis, and Diels-Alder reactions. These proline derivatives allowed three parallel bioorthogonal reactions to be conducted in one solution.

Secreted trypanosome cyclophilin inactivates lytic insect defense peptides and induces parasite calcineurin activation and infectivity

The mechanisms by which Trypanosoma cruzi survives antimicrobial peptides and differentiates during its transit through the gastrointestinal tract of the reduviid vector are unknown. We show that cyclophilin, a peptidyl-prolyl isomerase secreted from T. cruzi epimastigotes, binds to and neutralizes the reduviid antimicrobial peptide trialysin promoting parasite survival. This is dependent on a singular proline residue in trialysin and is inhibited by the cyclophilin inhibitor cyclosporine A. In addition, cyclophilin-trialysin complexes enhance the production of ATP and reductase responses of parasites, which are inhibited by both calcineurin-specific inhibitors cyclosporine A and FK506. Calcineurin phosphatase activity of cyclophilin-trialysin-treated parasites was higher than in controls and was inhibited by preincubation by either inhibitor. Parasites exposed to cyclophilin-trialysin have enhanced binding and invasion of host cells leading to higher infectivity. Leishmanial cyclophilin also mediates trialysin protection and metabolic stimulation by T. cruzi, indicating that extracellular cyclophilin may be critical to adaptation in other insect-borne protozoa. This work demonstrates that cyclophilin serves as molecular sensor leading to the evasion and adaptive metabolic response to insect defense peptides.

Characterization of secondary amide peptide bond isomerization: thermodynamics and kinetics from 2D NMR spectroscopy

Secondary amide cis peptide bonds are of even lower abundance than the cis tertiary amide bonds of prolines, yet they are of biochemical importance. Using 2D NMR exchange spectroscopy (EXSY) we investigated the formation of cis peptide bonds in several oligopeptides: Ac-G-G-G-NH(2) , Ac-I-G-G-NH(2) , Ac-I-G-G-N-NH(2) and its cyclic form: I-G-G-N in dimethylsulfoxide (DMSO). From the NMR studies, using the amide protons as monitors, an occurrence of 0.13-0.23% of cis bonds was obtained at 296 K. The rate constants for the trans to cis conversion determined from 2D EXSY spectroscopy were 4-9 × 10(-3) s(-1) . Multiple minor conformations were detected for most peptide bonds. From their thermodynamic and kinetic properties the cis isomers are distinguished from minor trans isomers that appear because of an adjacent cis peptide bond. Solvent and sequence effects were investigated utilizing N-methylacetamide (NMA) and various peptides, which revealed a unique enthalpy profile in DMSO. The cyclization of a tetrapeptide resulted in greatly lowered cis populations and slower isomerization rates compared to its linear counterpart, further highlighting the impact of structural constraints.

Unraveling the role of peptidyl-prolyl isomerases in neurodegeneration

Immunophilins are a family of highly conserved proteins with a peptidyl-prolyl isomerase activity that binds immunosuppressive drugs such as FK506, cyclosporin A, and rapamycin. Immunophilins can be divided into two subfamilies, the cyclophilins, and the FK506 binding proteins (FKBPs). Next to the immunophilins, a third group of peptidyl-prolyl isomerases exist, the parvulins, which do not influence the immune system. The beneficial role of immunophilin ligands in neurodegenerative disease models has been known for more than a decade but remains largely unexplained in terms of molecular mechanisms. In this review, we summarize reported effects of parvulins, immunophilins, and their ligands in the context of neurodegeneration. We focus on the role of FKBP12 in Parkinson's disease and propose it as a novel drug target for therapy of Parkinson's disease.

Human GLTP: Three distinct functions for the three tryptophans in a novel peripheral amphitropic fold

Human glycolipid transfer protein (GLTP) serves as the GLTP-fold prototype, a novel, to our knowledge, peripheral amphitropic fold and structurally unique lipid binding motif that defines the GLTP superfamily. Despite conservation of all three intrinsic Trps in vertebrate GLTPs, the Trp functional role(s) remains unclear. Herein, the issue is addressed using circular dichroism and fluorescence spectroscopy along with an atypical Trp point mutation strategy. Far-ultraviolet and near-ultraviolet circular dichroism spectroscopic analyses showed that W96F-W142Y-GLTP and W96Y-GLTP retain their native conformation and stability, whereas W85Y-W96F-GLTP is slightly altered, in agreement with relative glycolipid transfer activities of >90%, ∼85%, and ∼45%, respectively. In silico three-dimensional modeling and acrylamide quenching of Trp fluorescence supported a nativelike folding conformation. With the Trp⁹⁶-less mutants, changes in emission intensity, wavelength maximum, lifetime, and time-resolved anisotropy decay induced by phosphoglyceride membranes lacking or containing glycolipid and by excitation at different wavelengths along the absorption-spectrum red edge indicated differing functions for W142 and W85. The data suggest that W142 acts as a shallow-penetration anchor during docking with membrane interfaces, whereas the buried W85 indole helps maintain proper folding and possibly regulates membrane-induced transitioning to a glycolipid-acquiring conformation. The findings illustrate remarkable versatility for Trp, providing three distinct intramolecular functions in the novel amphitropic GLTP fold.

Structural determinants of protein folding

The last several decades have seen an explosion of knowledge in the field of structural biology. With critical advances in spectroscopic techniques in examining structures of biomacromolecules, in maturation of molecular biology techniques, as well as vast improvements in computation prowess, protein structures are now being elucidated at an unprecedented rate. In spite of all the recent advances, the protein folding puzzle remains as one of the fundamental biochemical challenges. A facet to this empiric problem is the structural determinants of protein folding. What are the driving forces that pivot a polypeptide chain to a specific conformation amongst the vast conformation space? In this review, we shall discuss some of the structural determinants to protein folding that have been identified in the recent decades.

Detection of discriminative sequence patterns in the neighborhood of proline cis peptide bonds and their functional annotation

Background

Polypeptides are composed of amino acids covalently bonded via a peptide bond. The majority of peptide bonds in proteins is found to occur in the trans conformation. In spite of their infrequent occurrence, cis peptide bonds play a key role in the protein structure and function, as well as in many significant biological processes.

Results

We perform a systematic analysis of regions in protein sequences that contain a proline cis peptide bond in order to discover non-random associations between the primary sequence and the nature of proline cis/trans isomerization. For this purpose an efficient pattern discovery algorithm is employed which discovers regular expression-type patterns that are overrepresented (i.e. appear frequently repeated) in a set of sequences. Four types of pattern discovery are performed: i) exact pattern discovery, ii) pattern discovery using a chemical equivalency set, iii) pattern discovery using a structural equivalency set and iv) pattern discovery using certain amino acids' physicochemical properties. The extracted patterns are carefully validated using a specially implemented scoring function and a significance measure (i.e. log-probability estimate) indicative of their specificity. The score threshold for the first three types of pattern discovery is 0.90 while for the last type of pattern discovery 0.80. Regarding the significance measure, all patterns yielded values in the range [-9, -31] which ensure that the derived patterns are highly unlikely to have emerged by chance. Among the highest scoring patterns, most of them are consistent with previous investigations concerning the neighborhood of cis proline peptide bonds, and many new ones are identified. Finally, the extracted patterns are systematically compared against the PROSITE database, in order to gain insight into the functional implications of cis prolyl bonds.

Conclusion

Cis patterns with matches in the PROSITE database fell mostly into two main functional clusters: family signatures and protein signatures. However considerable propensity was also observed for targeting signals, active and phosphorylation sites as well as domain signatures.

Prediction of cis/trans isomerization using feature selection and support vector machines

In protein structures the peptide bond is found to be in trans conformation in the majority of the cases. Only a small fraction of peptide bonds in proteins is reported to be in cis conformation. Most of these instances (>90%) occur when the peptide bond is an imide (X-Pro) rather than an amide bond (X-nonPro). Due to the implication of cis/trans isomerization in many biologically significant processes, the accurate prediction of the peptide bond conformation is of high interest. In this study, we evaluate the effect of a wide range of features, towards the reliable prediction of both proline and non-proline cis/trans isomerization. We use evolutionary profiles, secondary structure information, real-valued solvent accessibility predictions for each amino acid and the physicochemical properties of the surrounding residues. We also explore the predictive impact of a modified feature vector, which consists of condensed position-specific scoring matrices (PSSMX), secondary structure and solvent accessibility. The best discriminating ability is achieved using the first feature vector combined with a wrapper feature selection algorithm and a support vector machine (SVM). The proposed method results in 70% accuracy, 75% sensitivity and 71% positive predictive value (PPV) in the prediction of the peptide bond conformation between any two amino acids. The output of the feature selection stage is investigated in order to identify discriminatory features as well as the contribution of each neighboring residue in the formation of the peptide bond, thus, advancing our knowledge towards cis/trans isomerization.

Predicting peptide bond conformation using feature selection and the Naïve Bayes approach

Distinguishing cis peptide bonds from trans isomers in protein sequences facilitates the exploration of protein structures and functions. In this study, we evaluated the effect of a large and informative feature vector, towards the reliable prediction of peptide bond conformation between any two amino acids. We used multiple sequence alignment, secondary structure information, real valued solvent accessibility predictions for each amino acid and physicochemical properties of the surrounding residues. A three stage schema was developed, comprising of feature extraction, feature selection and peptide bond classification between any two amino acids. We also explored the performance achieved when using the full feature vector without performing feature selection. The best discriminating ability was achieved using a Naïve Bayes classifier, combined with wrapper feature selection. The proposed approach yielded prediction accuracy 86%, sensitivity 82% and specificity 90% in discriminating cis and trans peptide bond conformations.

First-Principle Computational Study on the Full Conformational Space of l-Threonine Diamide, the Energetic Stability of Cis and Trans Isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.

Effects of i and i+3 residue identity on cis-trans isomerism of the aromatic(i+1)-prolyl(i+2) amide bond: implications for type VI beta-turn formation

Cis-trans isomerization of amide bonds plays critical roles in protein molecular recognition, protein folding, protein misfolding, and disease. Aromatic-proline sequences are particularly prone to exhibit cis amide bonds. The roles of residues adjacent to a tyrosine-proline residue pair on cis-trans isomerism were examined. A short series of peptides XYPZ was synthesized and cis-trans isomerism was analyzed. Based on these initial studies, a series of peptides XYPN, X = all 20 canonical amino acids, was synthesized and analyzed by NMR for i residue effects on cis-trans isomerization. The following effects were observed: (a) aromatic residues immediately preceding Tyr-Pro disfavor cis amide bonds, with K(trans/cis)= 5.7-8.0, W > Y > F; (b) proline residues preceding Tyr-Pro lead to multiple species, exhibiting cis-trans isomerization of either or both X-Pro amide bonds; and (c) other residues exhibit similar values of K(trans/cis) (= 2.9-4.2), with Thr and protonated His exhibiting the highest fraction cis. beta-Branched and short polar residues were somewhat more favorable in stabilizing the cis conformation. Phosphorylation of serine at the i position modestly increases the stability of the cis conformer. In addition, the effect of the i+3 residue was examined in a limited series of peptides TYPZ. NMR data indicated that aromatic residues, Pro, Asn, Ala, and Val at the i+3 residue all favor cis amide bonds, with aromatic residues and Asn favoring more compact phi at Tyr(cis) and Ala and Pro favoring more extended phi at Tyr(cis). D-Alanine at the i+3 position particularly disfavors cis amide bonds.

Electronic control of amide cis-trans isomerism via the aromatic-prolyl interaction

The cis-trans isomerization of prolyl amide bonds results in large structural and functional changes in proteins and is a rate-determining step in protein folding. We describe a novel electronic strategy to control cis-trans isomerization, based on the demonstration that interactions between aromatic residues and proline are tunable by aromatic electronics. A series of peptides of sequence TXPN, X = Trp, pyridylalanine, pentafluorophenylalanine, or 4-Z-phenylalanine derivatives (Z = electron-donating, electron-withdrawing, or electron-neutral substituents), was synthesized and Ktrans/cis analyzed by NMR. Electron-rich aromatic residues stabilized cis amide bond formation, while electron-poor aromatics relatively favored trans amide bond formation. A Hammett correlation between aromatic electronics and cis-trans isomerization was observed. These results indicate that the interaction between aromatic residues and proline, which is observed to stabilize cis amide bonds and is also a general stabilizing interaction ubiquitous in proteins and protein-protein complexes, is not stabilized exclusively by a classical hydrophobic effect. To a large extent, the aromatic-prolyl interaction is driven and controllable by an electronic effect between the aromatic ring pi-electrons and the proline ring, consistent with a C-H-pi interaction as the key stabilizing force. The aromatic-prolyl interaction is electronically tunable by 0.9 kcal/mol and is enthalpic in nature. In addition, by combining aromatic ring electronics and stereoelectronic effects using 4-fluoroprolines, we demonstrate broad tuning (2.0 kcal/mol) of cis-trans isomerism in tetrapeptides. We demonstrate a simple tetrapeptide, TWflpN, that exhibits 60% cis amide bond and adopts a type VIa1 beta-turn conformation.

Analysis and prediction of helix–helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.

Prediction of cis/trans isomerization in proteins using PSI-BLAST profiles and secondary structure information

Background

The majority of peptide bonds in proteins are found to occur in the trans conformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt the cis form. Proline cis/trans isomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of proline cis/trans isomerization in proteins would have many important applications towards the understanding of protein structure and function.

Results

In this paper, we propose a new approach to predict the proline cis/trans isomerization in proteins using support vector machine (SVM). The preliminary results indicated that using Radial Basis Function (RBF) kernels could lead to better prediction performance than that of polynomial and linear kernel functions. We used single sequence information of different local window sizes, amino acid compositions of different local sequences, multiple sequence alignment obtained from PSI-BLAST and the secondary structure information predicted by PSIPRED. We explored these different sequence encoding schemes in order to investigate their effects on the prediction performance. The training and testing of this approach was performed on a newly enlarged dataset of 2424 non-homologous proteins determined by X-Ray diffraction method using 5-fold cross-validation. Selecting the window size 11 provided the best performance for determining the proline cis/trans isomerization based on the single amino acid sequence. It was found that using multiple sequence alignments in the form of PSI-BLAST profiles could significantly improve the prediction performance, the prediction accuracy increased from 62.8% with single sequence to 69.8% and Matthews Correlation Coefficient (MCC) improved from 0.26 with single local sequence to 0.40. Furthermore, if coupled with the predicted secondary structure information by PSIPRED, our method yielded a prediction accuracy of 71.5% and MCC of 0.43, 9% and 0.17 higher than the accuracy achieved based on the singe sequence information, respectively.

Conclusion

A new method has been developed to predict the proline cis/trans isomerization in proteins based on support vector machine, which used the single amino acid sequence with different local window sizes, the amino acid compositions of local sequence flanking centered proline residues, the position-specific scoring matrices (PSSMs) extracted by PSI-BLAST and the predicted secondary structures generated by PSIPRED. The successful application of SVM approach in this study reinforced that SVM is a powerful tool in predicting proline cis/trans isomerization in proteins and biological sequence analysis.

Bacillus circulans MH-K1 Chitosanase: Amino Acid Residues Responsible for Substrate Binding

To identify the amino acids responsible for the substrate binding of chitosanase from Bacillus circulans MH-K1 (MH-K1 chitosanase), Tyr148 and Lys218 of the chitosanase were mutated to serine and proline, respectively, and the mutated chitosanases were characterized. The enzymatic activities of Y148S and K218P were found to be 12.5% and 0.16% of the wild type, respectively. When the (GlcN)3 binding ability to the chitosanase was evaluated by fluorescence spectroscopy and thermal unfolding experiments, the binding abilities of both mutant enzymes were markedly reduced as compared with the wild type enzyme. The affinity of the enzyme for the trisaccharide decreased by 1.0 kcal/mol of binding free energy for Y148S, and 3.7 kcal/mol for K218P. The crystal structure of K218P revealed that Pro218 forms a cis-peptide bond and that the state of the flexible loop containing the 218th residue is considerably affected by the mutation. Thus, we conclude that the flexible loop containing Lys218 plays an important role in substrate binding, and that the role of Tyr148 is less critical, but still important, due to a stacking interaction or hydrogen bond.