A conformational comparison of N- and C-protected methionine and N- and C-protected homocysteine

An improved two-rotor function for conformational potential energy surfaces of 20 amino acid diamides

Predicting the three-dimensional structure of a protein from its amino acid sequence requires a complete understanding of the molecular forces that influences the protein folding process. Each possible conformation has its corresponding potential energy, which characterizes its thermodynamic stability. This is needed to identify the primary intra- and inter-molecular interactions, so that we can reduce the dimensionality of the problem, and create a relatively simple representation of the system. Investigating this problem using quantum chemical methods produces accurate results; however, this also entails large computational resources. In this study, an improved two-rotor potential energy function is proposed to represent the backbone interactions in amino acids through a linear combination of a Fourier series and a mixture of Gaussian functions. This function is applied to approximate the 20 amino acid diamide Ramachandran-type PESs, and results yielded an average RMSE of 2.36 kJ mol−1, which suggest that the mathematical model precisely captures the general topology of the conformational potential energy surface. Furthermore, this paper provides insights on the conformational preferences of amino acid diamides through local minima geometries and energy ranges, using the improved mathematical model. The proposed mathematical model presents a simpler representation that attempts to provide a framework on building polypeptide models from individual amino acid functions, and consequently, a novel method for rapid but accurate evaluation of potential energies for biomolecular simulations.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

Peptide models XLV: Conformational properties of N-formyl-L-methioninamide and its relevance to methionine in proteins

The conformational space of the most biologically significant backbone folds of a suitable methionine peptide model was explored by density functional computational method. Using a medium [6-31G(d)] and a larger basis set [6-311++G(2d,2p)], the systematic exploration of low-energy backbone structures restricted for the "L-region" in the Ramachandran map of N-formyl-L-methioninamide results in conformers corresponding to the building units of an extended backbone structure (betaL), an inverse gamma-turn (gammaL), or a right-handed helical structure (alphaL). However, no poly-proline II type (epsilonL) fold was found, indicating that this conformer has no intrinsic stability, and highlighting the effect of molecular environment in stabilizing this backbone structure. This is in agreement with the abundance of the epsilonL-type backbone conformation of methionine found in proteins. Stability properties (DeltaE) and distinct backbone-side-chain interactions support the idea that specific intramolecular contacts are operative in the selection of the lowest energy conformers. Apart from the number of different folds, all stable conformers are within a 10 kcal x mol(-1) energy range, indicating the highly flexible behavior of methionine. This conformational feature can be important in supporting catalytic processes, facilitating protein folding and dimerization via metal ion binding. In both of the biological examples discussed (HIV-1 reverse transcriptase and PcoC copper-resistant protein), the conformational properties of Met residues were found to be of key importance. Spatial proximity to other types of residues or the same type of residue seems to be crucial for the structural integrity of a protein, whether Met is buried or exposed.