Main-chain Torsion Research Articles

Molecular Mechanic Studies on the Zwitterionic and Neutral Conformation Stability of Biological Active Peptide

In this current work the conformational properties of the L-Alanyl-L-Glutamine dipeptide which was biological active dipeptide has been investigated by theoretical conformational analysis method using program which is based on Molecular Mechanic calculations in order to determine the structure function relation. The calculations on zwitterionic and neutral form of Ala-Gln which is formed side and main chain torsion angles let us to define their preferable energetically conformers. The side and main chain status of the stable conformations of Ala-Gln were acquired based on the results of the conformational analysis. By using Ramachandran maps, the global conformation, which has the lowest energy of the dipeptide, has been determined. The goal of this work was to explain structure-activity relationship by investigating the bioactive conformation of biological active peptide.

The role of the backbone torsion in protein folding

BackgroundThe set of forces and sequence of events that govern the transition from an unfolded polypeptide chain to a functional protein with correct spatial structure remain incompletely known, despite the importance of the problem and decades of theory development, computer simulations, and laboratory experiments. Information about the correctly folded state of most proteins is likely to be present in their sequences, and yet many proteins fail to attain native structure after overexpression in a non-native environment or upon experimental denaturation and refolding.Presentation of the hypothesisWe hypothesize that correct protein folding in vivo is an active, energy-dependent process that most likely applies torque force co-translationally to all proteins and possibly also post-translationally to many proteins in every cell. When a site on an unfolded polypeptide is rotationally constrained, torsion applied at another site would induce twisting of the main chain, which would initiate the formation of a local secondary structure, such as an alpha-helical turn or a beta-turn/beta-hairpin. The nucleation of structural elements is a rate-limiting, energetically unfavorable step in the process of protein folding, and energy-dependent chain torsion is likely to help overcome this barrier in vivo. Several molecular machines in a cell, primarily ribosomes, but also possibly signal recognition particles and chaperone systems, may play a role in applying torque to an unfolded protein chain, using the energy of GTP or ATP hydrolysis. Lack of such force in the in vitro systems may be the main reason of the failure of many longer proteins to attain the correct functional conformation.Testing the hypothesisThe hypothesis can be tested using single-molecule approaches, by measuring directly the forces applied to polypeptide chains under controlled conditions in vitro, and in bulk, by assessing folding rates and extent of misfolding in proteins that are engineered to experience transient spatial constraint during their synthesis.Implications of the hypothesisLearning about the role of main chain torsion in protein folding will improve our understanding of folding mechanisms and may lead to bioengineering solutions that would enhance the yields of correctly folded proteins in heterologous expression systems.ReviewersThis article was reviewed by Frank Eisenhaber, Igor Berezovsky and Michael Gromiha.

Correct folding of an α-helix and a β-hairpin using a polarized 2D torsional potential.

A new modification to the AMBER force field that incorporates the coupled two-dimensional main chain torsion energy has been evaluated for the balanced representation of secondary structures. In this modified AMBER force field (AMBER032D), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by high-level quantum mechanical calculations of small peptides in solution. Molecular dynamics simulations are performed to study the folding of two model peptides adopting either α-helix or β-hairpin structures. Both peptides are successfully folded into their native structures using an AMBER032D force field with the implementation of a polarization scheme (AMBER032Dp). For comparison, simulations using a standard AMBER03 force field with and without polarization, as well as AMBER032D without polarization, fail to fold both peptides successfully. The correction to secondary structure propensity in the AMBER03 force field and the polarization effect are critical to folding Trpzip2; without these factors, a helical structure is obtained. This study strongly suggests that this new force field is capable of providing a more balanced preference for helical and extended conformations. The electrostatic polarization effect is shown to be indispensable to the growth of secondary structures.

Coupled two-dimensional main-chain torsional potential for protein dynamics II: performance and validation.

The accuracy of force fields is of utmost importance in molecular modeling of proteins. Despite successful applications of force fields for about the past 30 years, some inherent flaws lying in force fields, such as biased secondary propensities and fixed atomic charges, have been observed in different aspects of biomolecular research; hence, a correction to current force fields is desirable. Because of the simplified functional form and the limited number of parameters for main chain torsion (MCT) in traditional force fields, it is not easy to propose an exquisite force field that is well-balanced among various conformations. Recently, AMBER-compatible force fields with coupled MCT term have been proposed, which show some improvement over AMBER03 and AMBER99SB force fields. In this work, further calibration of the torsional parameters has been conducted by changing the solvation model in quantum mechanical calculation and minimizing the deviation from the nuclear magnetic resonance experiments for some benchmark model systems and a folded protein. The results show that the revised force fields give excellent agreement with experiments in J coupling, chemical shifts, and secondary structure populations. In addition, the polarization effect is found to be crucial for the systems with ordered secondary structures.

A Coupled Two-Dimensional Main Chain Torsional Potential for Protein Dynamics

A new AMBER compatible force field is proposed for balanced representation of secondary structures. In this modified AMBER force field (AMBER2D), the main chain torsion energy is represented by 2-dimensional Fourier expansions with parameters fitted to the potential energy surface generated by quantum mechanical calculations of small peptides in solution at M06 2X/aug-cc-pvtz//HF/6-31G∗∗ level. Solvation model based on solute electron density (SMD) developed in Truhlar's group was considered. The benchmark systems used in the validation of this force field include capped dipeptides (Ace-X-NME, XP), tripeptides (XXX, XA, G, V}; GYG, Y {A, V, F, L, S, E, K, M}), alanine tetrapeptide, Ac-(AAQAA)3 -NH2, and ubiquitin. Besides, we also investigate the folding of two representative proteins (PDB ID 2I9M and 1LE1). The results demonstrated that this 2D main chain torsion is effective in delineating the energy variation associated with main chain torsions. Furthermore, the electrostatic polarization effect is very important for long peptides or proteins. This work also serves as an implication for the necessity of more accurate functional forms for main chain torsions in the future development of force field.

Conformational analysis of a polyconjugated protein-binding ligand by joint quantum chemistry and polarizable molecular mechanics. Addressing the issues of anisotropy, conjugation, polarization, and multipole transferability.

We investigate the conformational properties of a potent inhibitor of neuropilin-1, a protein involved in cancer processes and macular degeneration. This inhibitor consists of four aromatic/conjugated fragments: a benzimidazole, a methylbenzene, a carboxythiourea, and a benzene-linker dioxane, and these fragments are all linked together by conjugated bonds. The calculations use the SIBFA polarizable molecular mechanics procedure. Prior to docking simulations, it is essential to ensure that variations in the ligand conformational energy upon rotations around its six main-chain torsional bonds are correctly represented (as compared to high-level ab initio quantum chemistry, QC). This is done in two successive calibration stages and one validation stage. In the latter, the minima identified following independent stepwise variations of each of the six main-chain torsion angles are used as starting points for energy minimization of all the torsion angles simultaneously. Single-point QC calculations of the minimized structures are then done to compare their relative energies ΔE conf to the SIBFA ones. We compare three different methods of deriving the multipoles and polarizabilities of the central, most critical moiety of the inhibitor: carboxythiourea (CTU). The representation that gives the best agreement with QC is the one that includes the effects of the mutual polarization energy E pol between the amide and thioamide moieties. This again highlights the critical role of this contribution. The implications and perspectives of these findings are discussed.

Folding simulation of Trp-cage utilizing a new AMBER compatible force field with coupled main chain torsions

A newly developed AMBER compatible force field with coupled backbone torsion potential terms (AMBER032D) is utilized in a folding simulation of a mini-protein Trp-cage. Through replica exchange and direct molecular dynamics (MD) simulations, a multi-step folding mechanism with a synergetic folding of the hydrophobic core (HPC) and the α-helix in the final stage is suggested. The native structure has the lowest free energy and the melting temperature predicted from the specific heat capacity Cvis only 12 K higher than the experimental measurement. This study, together with our previous study, shows that AMBER032Dis an accurate force field that can be used for protein folding simulations.

Correlation as a Determinant of Configurational Entropy in Supramolecular and Protein Systems

For biomolecules in solution, changes in configurational entropy are thought to contribute substantially to the free energies of processes like binding and conformational change. In principle, the configurational entropy can be strongly affected by pairwise and higher-order correlations among conformational degrees of freedom. However, the literature offers mixed perspectives regarding the contributions that changes in correlations make to changes in configurational entropy for such processes. Here we take advantage of powerful techniques for simulation and entropy analysis to carry out rigorous in silico studies of correlation in binding and conformational changes. In particular, we apply information-theoretic expansions of the configurational entropy to well-sampled molecular dynamics simulations of a model host–guest system and the protein bovine pancreatic trypsin inhibitor. The results bear on the interpretation of NMR data, as they indicate that changes in correlation are important determinants of entropy changes for biologically relevant processes and that changes in correlation may either balance or reinforce changes in first-order entropy. The results also highlight the importance of main-chain torsions as contributors to changes in protein configurational entropy. As simulation techniques grow in power, the mathematical techniques used here will offer new opportunities to answer challenging questions about complex molecular systems.

Electronic polarization stabilizes tertiary structure prediction of HP-36

Molecular dynamic (MD) simulations with both implicit and explicit solvent models have been carried out to study the folding dynamics of HP-36 protein. Starting from the extended conformation, the secondary structure of all three helices in HP-36 was formed in about 50 ns and remained stable in the remaining simulation. However, the formation of the tertiary structure was difficult. Although some intermediates were close to the native structure, the overall conformation was not stable. Further analysis revealed that the large structure fluctuation of loop and hydrophobic core regions was devoted mostly to the instability of the structure during MD simulation. The backbone root-mean-square deviation (RMSD) of the loop and hydrophobic core regions showed strong correlation with the backbone RMSD of the whole protein. The free energy landscape indicated that the distribution of main chain torsions in loop and turn regions was far away from the native state. Starting from an intermediate structure extracted from the initial AMBER simulation, HP-36 was found to generally fold to the native state under the dynamically adjusted polarized protein-specific charge (DPPC) simulation, while the peptide did not fold into the native structure when AMBER force filed was used. The two best folded structures were extracted and taken into further simulations in water employing AMBER03 charge and DPPC for 25 ns. Result showed that introducing polarization effect into interacting potential could stabilize the near-native protein structure.

Structural basis of conformational transitions in the active site and 80′s loop in the FK506-binding protein FKBP12

The extensive set of NMR doublings exhibited by the immunophilin FKBP12 (FK506-binding protein 12) arose from a slow transition to the cis-peptide configuration at Gly89 near the tip of the 80′s loop, the site for numerous protein-recognition interactions for both FKBP12 and other FKBP domain proteins. The 80′s loop also exhibited linebroadening, indicative of microsecond to millisecond conformational dynamics, but only in the trans-peptide state. The G89A variant shifted the trans–cis peptide equilibrium from 88:12 to 33:67, whereas a proline residue substitution induced fully the cis-peptide configuration. The 80′s loop conformation in the G89P crystal structure at 1.50 Å resolution differed from wild-type FKBP12 primarily at residues 88, 89 and 90, and it closely resembled that reported for FKBP52. Structure-based chemical-shift predictions indicated that the microsecond to millisecond dynamics in the 80′s loop probably arose from a concerted main chain (ψ88 and ϕ89) torsion angle transition. The indole side chain of Trp59 at the base of the active-site cleft was reoriented ~90o and the adjacent backbone was shifted in the G89P crystal structure. NOE analysis of wild-type FKBP12 demonstrated that this indole populates the perpendicular orientation at 20%. The 15N relaxation analysis was consistent with the indole reorientation occurring in the nanosecond timeframe. Recollection of the G89P crystal data at 1.20 Å resolution revealed a weaker wild-type-like orientation for the indole ring. Differences in the residues that underlie the Trp59 indole ring and altered interactions linking the 50′s loop to the active site suggested that reorientation of this ring may be disfavoured in the other six members of the FKBP domain family that bear this active-site tryptophan residue.

Conformational Entropy of the RNA Phosphate-Sugar Backbone

While major contributors to the free energy of RNA tertiary structures such as base pairing, base stacking and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the mainchain torsion angle conformation space between every pair of nucleobases along any sequence and the corresponding backbone entropy. We compare the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several thousand backbone conformations between nearest-neighbor nucleobases as well as against computer simulations. Using this calculation we have analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized. This suggests the folding and stability of the RNA structures are critically dependent on these local entropic strain domains and leave the rest of the sequence relatively flexible. The total entropic penalty in the backbone for the native fold can be as high as 0.7 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 40 kcal/mol for a 200-nucleotide structure. We also look at the correlation of nucleotide conformations along several sample loop sequences to determine the effect of nucleotides beyond the nearest-neighbor.

A coupled two-dimensional main chain torsional potential for protein dynamics: generation and implementation

Main chain torsions of alanine dipeptide are parameterized into coupled 2-dimensional Fourier expansions based on quantum mechanical (QM) calculations at M06 2X/aug-cc-pvtz//HF/6-31G** level. Solvation effect is considered by employing polarizable continuum model. Utilization of the M06 2X functional leads to precise potential energy surface that is comparable to or even better than MP2 level, but with much less computational demand. Parameterization of the 2D expansions is against the full main chain torsion space instead of just a few low energy conformations. This procedure is similar to that for the development of AMBER03 force field, except unique weighting factor was assigned to all the grid points. To avoid inconsistency between quantum mechanical calculations and molecular modeling, the model peptide is further optimized at molecular mechanics level with main chain dihedral angles fixed before the calculation of the conformational energy on molecular mechanical level at each grid point, during which generalized Born model is employed. Difference in solvation models at quantum mechanics and molecular mechanics levels makes this parameterization procedure less straightforward. All force field parameters other than main chain torsions are taken from existing AMBER force field. With this new main chain torsion terms, we have studied the main chain dihedral distributions of ALA dipeptide and pentapeptide in aqueous solution. The results demonstrate that 2D main chain torsion is effective in delineating the energy variation associated with rotations along main chain dihedrals. This work is an implication for the necessity of more accurate description of main chain torsions in the future development of ab initio force field and it also raises a challenge to the development of quantum mechanical methods, especially the quantum mechanical solvation models.

Sub-similarity matching based on data mining with dihedral angles

Protein sub-similarity matching remains largely unknown even though it is becoming one of the most important open problems in bioinformatics for drug and vaccine design. Variations in human immune responses to vaccines are, and thus responses, fail. We propose a new matching and protein alignment method based on clustering and Longest Common Subsequence (LCS) techniques. After clustering, we found LCS between a candidate protein and meningitis outer membrane antigen for each candidate. Each similarity was scored, and closest similarities were determined with statistical methods. We located three closely matching proteins among a total of 50 human immune system proteins. Moreover, we selected a HIV-1 related protein from one of scenarios, because it revealed a relationship between HIV and meningitis patients. We also found that Ω main chain torsion angle for atoms CA, C and N is the best angle for determining sub-similarities between meningitis antigen and immune proteins.

A proteomic Ramachandran plot (PRplot)

Each protein structure can be characterized by the average values of the main chain torsion angles ϕ and ψ and, as a consequence, be plotted on a bidimensional diagram, which resembles the Ramachandran plot. Here, we describe a proteomic ϕ-ψ plot (PRplot) where each protein structure is associated with one point, allowing in this way to represent the entire protein structure universe. It was verified that the PRplot is a robust tool since it does not depend on the dimension of the proteins, on the crystallographic resolution of the structures, nor on the biological source; moreover, it is little affected by disordered and structurally uncharacterized residues. The proteins mapped on the PRplot tend to cluster in three regions that correspond to the structures rich in alpha-helices, in beta-strands, and in both helices and strands, and are distributed along a sigmoidal curve that connect these three highly populated regions. PRplots are a unique instrument to project all protein structures on a single bidimensional plane where the entire structural complexity is reduced to a striking simplicity, with the sigmoid curve clearly delineating the space fraction accessible to a stable protein.

Invariant Gly Residue Is Important for α-Defensin Folding, Dimerization, and Function: A CASE STUDY OF THE HUMAN NEUTROPHIL α-DEFENSIN HNP1

The human α-defensins (HNP) are synthesized in vivo as inactive prodefensins, and contain a conserved glycine, Gly(17), which is part of a β-bulge structure. It had previously been shown that the glycine main chain torsion angles are in a D-configuration, and that d-amino acids but not L-alanine could be substituted at that position to yield correctly folded peptides without the help of a prodomain. In this study, the glycine to L-alanine mutant defensin was synthesized in the form of a prodefensin using native chemical ligation. The ligation product folded correctly and yielded an active peptide upon CNBr cleavage. The L-Ala(17)-HNP1 crystal structure depicted a β-bulge identical to wild-type HNP1. However, dimerization was perturbed, causing one monomer to tilt with respect to the other in a dimerization model. Inhibitory activity against the anthrax lethal factor showed a 2-fold reduction relative to wild-type HNP1 as measured by the inhibitory concentration IC(50). Self-association was slightly reduced, as detected by surface plasmon resonance measurements. According to the results of the virtual colony count assay, the antibacterial activity against Escherichia coli, Staphylococcus aureus, and Bacillus cereus exhibited a less than 2-fold reduction in virtual lethal dose values. Prodefensins with two other L-amino acid substitutions, Arg and Phe, at the same position did not fold, indicating that only small side chains are tolerable. These results further elucidate the factors governing the region of the β-bulge structure that includes Gly(17), illuminating why glycine is conserved in all mammalian α-defensins.

The aberrance of the 4S diastereomer of 4-hydroxyproline.

Prolyl 4-hydroxylases install a hydroxyl group in the 4R configuration on the gamma-carbon atom of certain (2S)-proline (Pro) residues in tropocollagen, elastin, and other proteins to form (2S,4R)-4-hydroxyproline (Hyp). The gauche effect arising from this prevalent post-translational modification enforces a C(gamma)-exo ring pucker and stabilizes the collagen triple helix. The Hyp diastereomer (2S,4S)-4-hydroxyproline (hyp) has not been observed in a protein, despite the ability of electronegative 4S substituents to enforce the more common C(gamma)-endo ring pucker of Pro. Here, we use density functional theory, spectroscopy, crystallography, and calorimetry to explore the consequences of hyp incorporation on protein stability using a collagen model system. We find that the 4S-hydroxylation of Pro to form hyp does indeed enforce a C(gamma)-endo ring pucker but a transannular hydrogen bond between the hydroxyl moiety and the carbonyl of hyp distorts the main-chain torsion angles that typically accompany a C(gamma)-endo ring pucker. This same transannular hydrogen bond enhances an n-->pi* interaction that stabilizes the trans conformation of the peptide bond preceding hyp, endowing hyp with the unusual combination of a C(gamma)-endo ring pucker and high trans/cis ratio. O-Methylation of hyp to form (2S,4S)-4-methoxyproline (mop) eliminates the transannular hydrogen bond and restores a prototypical C(gamma)-endo pucker. mop residues endow the collagen triple helix with much more conformational stability than do hyp residues. These findings highlight the critical importance of the configuration of the hydroxyl group installed on C(gamma) of proline residues.

Ethyl 2-acetyl-3-(4-bromoanilino)butanoate

The title compound, C14H18BrNO3, adopts an extended conformation, with all of the main-chain torsion angles associated with the ester and amino groups close to trans. In the crystal, inversion dimers linked by pairs of N—H⋯O hydrogen bonds are observed.

Stereoelectronic and steric effects in side chains preorganize a protein main chain.

Preorganization is shown to endow a protein with extraordinary conformational stability. This preorganization is achieved by installing side-chain substituents that impose stereoelectronic and steric effects that restrict main-chain torsion angles. Replacing proline residues in (ProProGly)(7) collagen strands with 4-fluoroproline and 4-methylproline leads to the most stable known triple helices, having T ( m ) values that are increased by > 50 degrees C. Differential scanning calorimetry data indicate an entropic basis to the hyperstability, as expected from an origin in preorganization. Structural data at a resolution of 1.21 A reveal a prototypical triple helix with insignificant deviations to its main chain, even though 2/3 of the residues are nonnatural. Thus, preorganization of a main chain by subtle changes to side chains can confer extraordinary conformational stability upon a protein without perturbing its structure.

Ethyl 2-acetyl-3-(4-chloroanilino)butanoate

The title compound, C14H18ClNO3, adopts an extended conformation, with all of the main chain torsion angles associated with the ester and amino groups trans. In the crystal, inversion dimers linked by pairs of N—H⋯O hydrogen bonds are observed.

Citrullination-dependent Differential Presentation of a Self-peptide by HLA-B27 Subtypes

Inflammatory processes are accompanied by the posttranslational modification of certain arginine residues within proteins to yield citrulline, although it is largely unknown how this modification influences antigen presentation. We employed crystallographic and functional studies to investigate whether the exchange of arginine to citrulline affects the display of a peptide by two human major histocompatibility antigen class I subtypes, HLA-B(*)2705 and HLA-B(*)2709. Both differ only in residue 116 within the peptide binding groove despite their differential association with ankylosing spondylitis, an inflammatory rheumatic disorder. The crystal structures described here show that a modified self-peptide, pVIPR-U5 (RRKWURWHL; U = citrulline), is presented by the two HLA-B27 molecules in distinct conformations. These binding modes differ not only drastically from each other but also from the conformations exhibited by the non-citrullinated peptide in a given subtype. The differential reactivity of HLA-B27-restricted cytotoxic T cells with modified or unmodified pVIPR supports the structural findings and shows that the presentation of citrullinated peptides has the potential to influence immune responses.