New stochastic strategy to analyze helix folding

We propose an alternative stochastic strategy to search secondary structures based on the generalized simulated annealing (GSA) algorithm, by using conformational preferences based on the Ramachandran map. We optimize the search for polypeptide conformational space and apply to peptides considered to be good alpha-helix promoters above a critical number of residues. Our strategy to obtain conformational energies consist in coupling a classical force field (THOR package) with the GSA procedure, biasing the Phi x Psi backbone angles to the allowed regions in the Ramachandran map. For polyalanines we obtained stable alpha-helix structures when the number of residues were equal or exceeded 13 amino acids residues. We also observed that the energy gap between the global minimum and the first local minimum tends to increase with the polypeptide size. These conformations were generated by performing 2880 stochastic molecular optimizations with a continuum medium approach. When compared with molecular dynamics or Monte Carlo methods, GSA can be considered the fastest.     

Calibration of effective van der Waals atomic contact radii for proteins and peptides

Effective van der Waals radii were calibrated in such a way that molecular models built from standard bond lengths and bond angles reproduced the amino acid conformations observed by crystallography in proteins and peptides. The calibrations were based on the comparison of the Ramachandran plots prepared from high-resolution X-ray data of proteins and peptides with the allowed phi, psi torsional angle space for the dipeptide molecular models. The calibrated radii are useful as criteria with which to filter energetically improbable conformations in molecular modeling studies of proteins and peptides.     

Peptides and peptoids—A quantum chemical structure comparison

Peptoids represent a very interesting structure alternative to peptides. Based on ab initio MO theory employing the 6-31G* and 3-21G basis sets and considering correlation energy, a systematic structure comparison between the basic structure units of peptoids and peptides is performed. The calculations show three minimum conformations denoted as C_7β, α_D, and α that do not correspond to conformers on the peptide potential energy hypersurface. The possibility of cis peptide bonds in the peptoids was examined. The solvent influence on the structure was estimated by means of various quantum chemical continuum models. © 1996 John Wiley & Sons, Inc.     

Design and conformational analysis of peptoids containing N-hydroxy amides reveals a unique sheet-like secondary structure

N-hydroxy amides can be found in many naturally occurring and synthetic compounds and are known to act as both strong proton donors and chelators of metal cations. We have initiated studies of peptoids, or N-substituted glycines which contain N-hydroxy amide side chains to investigate the potential effects of these functional groups on peptoid backbone amide rotamer equilibria and local conformations. We reasoned that the propensity of these functional groups to participate in hydrogen bonding could be exploited to enforce intramolecular or intermolecular interactions that yield new peptoid structures. Here, we report the design, synthesis, and detailed conformational analysis of a series of model N-hydroxy peptoids. These peptoids were readily synthesized, and their structures were analyzed in solution by 1D and 2D NMR and in the solid-state by X-ray crystallography. The N-hydroxy amides were found to strongly favor trans conformations with respect to the peptoid backbone in chloroform. More notably, unique sheet-like structures held together via intermolecular hydrogen bonds were observed in the X-ray crystal structures of an N-hydroxy amide peptoid dimer, which to our knowledge represent the first structure of this type reported for peptoids. These results suggest that the N-hydroxy amide can be utilized to control both local backbone geometries and longer-range intermolecular interactions in peptoids, and represents a new functional group in the peptoid design toolbox.     

Molecular mechanics evaluation of the proposed mechanisms for the degradation of urea by urease

A thorough conformational search of all the conformations available to oxygen-bound urea within wild-type urease was carried out. Identical low energy urea conformations were obtained by a Ramachandran type plot for the NHis272-Ni1-O-Curea, and Ni1-O-Curea-Nurea dihedral angles. Ramachandran plots, with active sites and protonation states modified to model the different urease mechanisms, were used to evaluate the different mechanisms. Based upon the low energy conformations available to urea in the active site of wild-type urease one can conclude that the traditional "His320 acts as a base" mechanism is unlikely. while the N,O urea bridged and the reverse protonation mechanisms cannot be ruled out. A consensus hydrogen-bonding network that does not favor any of the mechanisms has been reconfirmed by the extensive conformational search.     

Multiple solution conformations and internal rotations of the decapeptide gramicidin S.

The conformations of every C alpha H-C beta H2 moiety of the peptide gramicidin S are reported. Internal rotation occurs, but distinct preferences for one side chain rotamer, greater than 80%, are found for the D-phenylalanine and ornithine residues. Leucine and valine exhibit more extensive averaging while proline is shown to be at least 90% in the Ramachandran B conformation. The data are consistent with the coexistence of many tertiary conformations of gramicidin S; the statistical weights of the twelve major tertiary conformations consistent with the rotamer populations are reported. The relative statistical weights of the tertiary conformers depend upon temperature and solvent. A comparison of the conclusions from this publication and conformations derived by energy minimization procedures is made. Partial agreement was found, but the calculations have not yet predicted the wealth of coexisting tertiary conformations nor accounted for the subtle effects of solvent. It is proposed that a more complete picture of the conformational dynamics of gramicidin S and other peptides will result from calculations which use as a basis the extensive data reported here.     

Molecular Mechanics Evaluation of the Proposed Mechanisms for the Degradation of Urea by Urease

Abstract

A thorough conformational search of all the conformations available to oxygen-bound urea within wild-type urease was carried out. Identical low energy urea conformations were obtained by a Ramachandran type plot for the N_His272-Ni1-O-C_urea and Ni1-O-Curea-N_urea dihedral angles. Ramachandran plots, with active sites and protonation states modified to model the different urease mechanisms, were used to evaluate the different mechanisms. Based upon the low energy conformations available to urea in the active site of wild-type urease one can conclude that the traditional “His320 acts as a base” mechanism is unlikely, while the N,O urea bridged and the reverse protonation mechanisms cannot be ruled out. A consensus hydrogen-bonding network that does not favor any of the mechanisms has been reconfirmed by the extensive conformational search.     

Rings and ribbons in protein structures: Characterization using helical parameters and Ramachandran plots for repeating dipeptides

Helical parameters displayed on a Ramachandran plot allow peptide structures with successive residues having identical main chain conformations to be studied. We investigate repeating dipeptide main chain conformations and present Ramachandran plots encompassing the range of possible structures. Repeating dipeptides fall into the categories: rings, ribbons, and helices. Partial rings occur in the form of “nests” and “catgrips”; many nests are bridged by an oxygen atom hydrogen bonding to the main chain NH groups of alternate residues, an interaction optimized by the ring structure of the nest. A novel recurring feature is identified that we name unpleated β, often situated at the ends of a β‐sheet strand. Some are partial rings causing the polypeptide to curve gently away from the sheet; some are straight. They lack β‐pleat and almost all incorporate a glycine. An example is the first glycine in the GxxxxGK motif of P‐loop proteins. Ribbons in repeating dipeptides can be either flat, as seen in repeated type II and type II′ β‐turns, or twisted, as in multiple type I and type I′ β‐turns. Hexa‐ and octa‐peptides in such twisted ribbons occur frequently in proteins, predominantly with type I β‐turns, and are the same as the “β‐bend ribbons” hitherto identified only in short peptides. One is seen in the GTPase‐activating protein for Rho in the active, but not the inactive, form of the enzyme. It forms a β‐bend ribbon, which incorporates the catalytic arginine, allowing its side chain guanidino group to approach the active site and enhance enzyme activity. Proteins 2014; 82:230–239. © 2013 Wiley Periodicals, Inc.     

Stability and amino acid preferences of type VIII reverse turn: the most common turn in peptides?

Free energies of the alpha(r)beta and betabeta conformations of 14 tetrapeptides, based on the sequence SALN and protein X-ray structures, were calculated using molecular dynamics simulations and MM-PBSA calculations. The alphaalpha conformations of five of the tetrapeptides were also studied. SALN has been earlier shown by molecular dynamics simulations and NMR spectroscopy to have a tendency to form an alpha(r)beta turn. The gas-phase energy of the molecular mechanical force field (CHARMM), the electrostatic and non-polar solvation free energies and solute entropies were used to explain the free energy differences of the alphaalpha, betabeta and alpha(r)beta conformations of the peptides. The alpha(r)beta conformation of SALN and SATN was predicted to be slightly more stable than the extended conformation (betabeta), in agreement with experimental results. The SALN mutants SAIN, SAVN, SATN, SSIN and MSHV, were also predicted to be potential alpha(r)beta turn-forming peptides. We report also revised positional potentials for the type VIII turn, based on a non-homologous set of protein structures. This protein databank analysis confirms the main results of the earlier analyses and reveals several new amino acid residues with a significant positional preference. The results of this work led us to suggest that the alpha(r)beta turn may be the most common turn type in peptides. Such turns may be readily formed in aqueous solution and thereby play important roles in the protein folding process by serving as an initiation point for structure formation.     

Carbonyl-carbonyl interactions stabilize the partially allowed Ramachandran conformations of asparagine and aspartic acid

Asparagine and aspartate are known to adopt conformations in the left-handed alpha-helical region and other partially allowed regions of the Ramachandran plot more readily than any other non-glycyl amino acids. The reason for this preference has not been established. An examination of the local environments of asparagine and aspartic acid in protein structures with a resolution better than 1.5 A revealed that their side-chain carbonyls are frequently within 4 A of their own backbone carbonyl or the backbone carbonyl of the previous residue. Calculations using protein structures with a resolution better than 1.8 A reveal that this close contact occurs in more than 80% of cases. This carbonyl-carbonyl interaction offers an energetic sabilization for the partially allowed conformations of asparagine and aspartic acid with respect to all other non-glycyl amino acids. The non-covalent attractive interactions between the dipoles of two carbonyls has recently been calculated to have an energy comparable to that of a hydrogen bond. The preponderance of asparagine in the left-handed alpha-helical region, and in general of aspartic acid and asparagine in the partially allowed regions of the Ramachandran plot, may be a consequence of this carbonyl-carbonyl stacking interaction.     

Molecular mechanics studies on the conformations of 2',3'-dideoxy-2',3'-didehydroguanine nucleoside, D4G.

Conformational energy calculations have been presented on guanine nucleoside in which the furanose ring is replaced by 2',3'-dideoxy-2',3'-didehydrofuran using molecular mechanics and conformational analysis. Conformational energies have been evaluated using the MM2 and AMBER94 force field parameters at two different dielectric constants. The results are presented in terms of isoenergy contours in the conformational space of the glycosidic (chi) and C4'-C5' (gamma) bonds torsions. In general, the chi-gamma interrelationships differ from the corresponding plots for unmodified nucleosides and nucleotides, reported previously. Consistency of the calculated preferred conformations with the x-ray data is sensitive to the force field employed.     

Experimental conformational energy maps of proteins and peptides

We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X‐ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (? , ψ ) plots in which global minima are predominantly observed either in the right‐handed α‐helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (? ,ψ ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol^?1) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X‐ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979–1001. © 2017 Wiley Periodicals, Inc.     

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I′ band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.     

Comparison of multiple Amber force fields and development of improved protein backbone parameters

The ff94 force field that is commonly associated with the Amber simulation package is one of the most widely used parameter sets for biomolecular simulation. After a decade of extensive use and testing, limitations in this force field, such as over-stabilization of alpha-helices, were reported by us and other researchers. This led to a number of attempts to improve these parameters, resulting in a variety of "Amber" force fields and significant difficulty in determining which should be used for a particular application. We show that several of these continue to suffer from inadequate balance between different secondary structure elements. In addition, the approach used in most of these studies neglected to account for the existence in Amber of two sets of backbone phi/psi dihedral terms. This led to parameter sets that provide unreasonable conformational preferences for glycine. We report here an effort to improve the phi/psi dihedral terms in the ff99 energy function. Dihedral term parameters are based on fitting the energies of multiple conformations of glycine and alanine tetrapeptides from high level ab initio quantum mechanical calculations. The new parameters for backbone dihedrals replace those in the existing ff99 force field. This parameter set, which we denote ff99SB, achieves a better balance of secondary structure elements as judged by improved distribution of backbone dihedrals for glycine and alanine with respect to PDB survey data. It also accomplishes improved agreement with published experimental data for conformational preferences of short alanine peptides and better accord with experimental NMR relaxation data of test protein systems.     

Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

Computational and experimental determination of the alpha-helix unfolding reaction coordinate

We demonstrate a calculated alpha-helix peptide folding energy landscape which accurately simulates the first experimentally measured alpha-helix melting energy landscape. We examine a 21-amino acid, mainly polyalanine peptide and calculate the free energy along the Psi Ramachandran angle secondary folding coordinate. The experimental free energy landscape was determined using UV resonance Raman spectroscopy. The relative free energy values are very close as are the equilibrium peptide conformations. We find 2.3 kcal/mol activation barriers between the alpha-helix-like and PPII-like basins. We also find that the alpha-helix-like conformations are quite defective and the alpha-helix-like structure dynamically samples 310-helix and pi-bulges.     

Controlling the secondary-structure-forming tendencies of proteins by a backbone torsion-energy term

We examined a new backbone torsion-energy term proposed by us in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, namely the backbone dihedral angles φ and ψ. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both φ and ψ can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase or decrease the α-helix-forming-tendencies by lowering or raising the torsion-energy surface in the α-helix region and likewise increase or decrease the β-sheet-forming tendencies by lowering or raising the surface in the β-sheet region in the Ramachandran space. We applied this torsion-energy modification method to six force fields, AMBER parm94, AMBER parm96, AMBER parm99, CHARMM27, OPLS-AA and OPLS-AA/L, and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming tendencies by performing folding simulations of α-helical and β-hairpin peptides.     

Peptides and Peptoids - A Systematic Structure Comparison

A systematic analysis of the conformational space of the basic structure unit of peptoids in comparison to the corresponding peptide unit was performed based on ab initio MO theory and complemented by molecular mechanics (MM) and molecular dynamics (MD) calculations both in the gas phase and in aqueous solution.The calculations show three minimum conformations denoted as C_7ß, a_D and a that do not correspond to conformers on the gas phase peptide potential energy hypersurface. The influence of aqueous solvation was estimated by means of continuum models. The MD simulations indicate the a_D form as the preferred conformation in solution both in cis and trans peptide bond orientations.     

Hydrophobic interactions modulate antimicrobial peptoid selectivity towards anionic lipid membranes

Hydrophobic interactions govern specificity for natural antimicrobial peptides. No such relationship has been established for synthetic peptoids that mimic antimicrobial peptides. Peptoid macrocycles synthesized with five different aromatic groups are investigated by minimum inhibitory and hemolytic concentration assays, epifluorescence microscopy, atomic force microscopy, and X-ray reflectivity. Peptoid hydrophobicity is determined using high performance liquid chromatography. Disruption of bacterial but not eukaryotic lipid membranes is demonstrated on the solid supported lipid bilayers and Langmuir monolayers. X-ray reflectivity studies demonstrate that intercalation of peptoids with zwitterionic or negatively charged lipid membranes is found to be regulated by hydrophobicity. Critical levels of peptoid selectivity are demonstrated and found to be modulated by their hydrophobic groups. It is suggested that peptoids may follow different optimization schemes as compared to their natural analogues.