Characterizing aqueous solution conformations of a peptide backbone using Raman optical activity computations

Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase

Investigating vibrational behavior in solution is crucial for understanding molecular dynamics within a solvent environment. Notably, the analysis of Raman spectra for molecules in solution is important owing to its ability to unveil intricate solute-solvent interactions. Previous studies have effectively employed frequency calculations utilizing the reference interaction site model self-consistent field method in conjunction with constrained spatial electron density distribution (RISM-SCF-cSED) to understand molecular vibrations in solution, primarily focusing on fundamental vibrational modes. However, the oversight of overtones and combination tones in these studies prompted us to combine the vibrational self-consistent field (VSCF) and vibrational second-order Mo̷ller-Plesset perturbation (VMP2) methods with RISM-SCF-cSED to address these aspects theoretically. Illustrating the efficacy of this integrated approach, we computed the Raman spectra of sodium formate (NaHCOO) in water, revealing the necessity of accounting for molecular anharmonicity in solution vibrational analysis. Our findings underscore the potency of VSCF and VMP2 in conjunction with RISM-SCF-cSED as a robust theoretical framework for such calculations.

The effect of protein backbone hydration on the amide vibrations in Raman and Raman optical activity spectra

Raman and specifically Raman optical activity (ROA) spectroscopy are very sensitive to the solution structure and conformation of biomolecules. Because of this strong conformational sensitivity, density functional theory (DFT) calculations are often used to get a better understanding of the experimentally observed spectral patterns. While e.g. for carbohydrate structure the water molecules that surround the solute have been demonstrated to be of vital importance to get accurate modelled ROA spectra, the effect of explicit water molecules on the calculated ROA patterns of peptides and proteins is less well studied. Here, the effect of protein backbone hydration was studied using DFT calculations of HCO-(l-Ala)5-NH2 in specific secondary structure conformations with different treatments of the solvation. The effect of the explicit water molecules on the calculated spectra mainly arises from the formation of hydrogen bonds with the amide C[double bond, length as m-dash]O and N-H groups. Hydrogen bonding of water with the C[double bond, length as m-dash]O group determines the shape and position of the amide I band. The C[double bond, length as m-dash]O bond length increases upon formation of C[double bond, length as m-dash]OH2O hydrogen bonds. The effect of the explicit water molecules on the amide III vibrations arises from hydrogen bonding of the solvent with both the C[double bond, length as m-dash]O and N-H group, but their contributions to this spectral region differ: geometrically, the formation of a C[double bond, length as m-dash]OH2O bond decreases the C-N bond length, while upon forming a N-HH2O hydrogen bond, the N-H bond length increases.

Ramachandran mapping of peptide conformation using a large database of computed Raman and Raman optical activity spectra

A novel ROA database is reported that assigns peptide structures in detail by pattern recognition of the experimental spectrum.

Conformational dynamics of carbohydrates: Raman optical activity of D-glucuronic acid and N-acetyl-D-glucosamine using a combined molecular dynamics and quantum chemical approach

As two biologically and medically relevant monosaccharides, the constituents of hyaluronic acid, d-glucuronic acid and N-acetyl-d-glucosamine, constitute perfect test cases for the development of carbohydrate-specific structural methods. These two molecules have been analysed by Raman optical activity (ROA), a spectroscopic technique exhibiting exquisite sensitivity to stereochemistry. We show that it is possible to support the experiment with a simulation approach combining density functional theory (DFT) and molecular dynamics (MD), both using explicit solvation. Thus, we have gained new insight into the crucial hydration effects that contribute to the conformational dynamics of carbohydrates and managed to characterize in detail the poorly understood vibrational nature of this class of biomolecules. Experimental and calculated ROA spectra of these two molecules are reported and excellent agreement has been found. More specifically, comparison has been made with the more commonly used gas phase and implicitly solvated calculation approaches, which offer poor or zero modelling of solvent interactions. The calculated spectra have been used to resolve the structural origins of the observed bands, a current challenge in the study of carbohydrates due to a lack of definitive vibrational assignments. We report and analyse major features in the fingerprint region of the ROA spectra, with recurrent structural and spectral features between the two monosaccharides observed.

The Raman optical activity of β-D-xylose: where experiment and computation meet

Besides its applications in bioenergy and biosynthesis, β-d-xylose is a very simple monosaccharide that exhibits relatively high rigidity. As such, it provides the best basis to study the impact of different solvation shell radii on the computation of its Raman optical activity (ROA) spectrum. Indeed, this chiroptical spectroscopic technique provides exquisite sensitivity to stereochemistry, and benefits much from theoretical support for interpretation. Our simulation approach combines density functional theory (DFT) and molecular dynamics (MD) in order to efficiently account for the crucial hydration effects in the simulation of carbohydrates and their spectroscopic response predictions. Excellent agreement between the simulated spectrum and the experiment was obtained with a solvation radius of 10 Å. Vibrational bands have been resolved from the computed ROA data, and compared with previous results on different monosaccharides in order to identify specific structure-spectrum relationships and to investigate the effect of the solvation environment on the conformational dynamics of small sugars. From the comparison with ROA analytical results, a shortcoming of the classical force field used for the MD simulations has been identified and overcome, again highlighting the complementary role of experiment and theory in the structural characterisation of complex biomolecules. Indeed, due to unphysical puckering, a spurious ring conformation initially led to erroneous conformer ratios, which are used as weights for the averaging of the spectral average, and only by removing this contribution was near perfect comparison between theory and experiment achieved.

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.

The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins

Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.

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.

PPII propensity of multiple-guest amino acids in a proline-rich environment

There has been considerable debate about the intrinsic PPII propensity of amino acid residues in denatured polypeptides. Experimentally, this scale is based on the behavior of guest amino acid residues placed in the middle of proline-based hosts. We have used classical molecular dynamics simulations combined with replica-exchange methods to carry out a comprehensive analysis of the conformational equilibria of proline-based host oligopeptides with multiple guest amino acids including alanine, glutamine, valine, and asparagine. The tracked structural characteristics include the secondary structural motifs based on the Ramachandran angles and the cis/trans isomerization of the prolyl bonds. In agreement with our recent study of single amino acid guests, we did not observe an intrinsic PPII propensity in any of the guest amino acids in a multiple-guest setting. Instead, the experimental results can be explained in terms of (i) the steric restrictions imposed on the C-terminal guest amino acid that is immediately followed by a proline residue and (ii) an increase in the trans content of the prolyl bonds due to the presence of guest residues. In terms of the latter, we found that the more guests added to the system, the larger the increase in the trans content of the prolyl bonds, which results in an effective increase in the PPII content of the peptide.

A statistical analysis of the PPII propensity of amino acid guests in proline-rich peptides

There has been considerable debate about the intrinsic PPII propensity of amino-acid residues in denatured polypeptides. Experimentally, the propensity scale is based on the behavior of guest amino-acid residues placed in the middle of polyproline hosts. We have used classical molecular dynamics simulations, with state-of-the-art force fields to carry out a comprehensive analysis of the conformational equilibria of the proline-based host oligopeptides with single guests. The tracked structural characteristics include the PPII content, the cis/trans isomerization of the prolyl bonds, the puckering of the pyrrolidine rings of the proline residues, and the secondary structural motifs. We find no evidence for an intrinsic PPII propensity in any of the guest amino acids other than proline. Instead, the PPII content as derived from experiments may be explained in terms of: 1), a local correlation between the dihedral angles of the guest amino acid and the proline residue immediately preceding it; and 2), a nonlocal correlation between the cis/trans states of the peptide bonds. In terms of the latter, we find that the presence of a guest (other than proline, tyrosine, or tryptophan) increases the trans content of most of the prolyl bonds, which results in an effective increase of the peptide PPII content. With respect to the local dihedral correlations, we find that these are well described in terms of the so-called odds-ratio statistic. Expressed in terms of free energy language, the PPII content based on the odds-ratio of the relevant residues correlate well with the experimentally measured PPII content.

Quantitative Assessment of Force Fields on Both Low-Energy Conformational Basins and Transition-State Regions of the (ϕ-ψ) Space

The free energy surfaces (FESs) of alanine dipeptide are studied to illustrate a new strategy to assess the performance of classical molecular mechanics force field on the full range of the (ϕ-ψ) conformational space. The FES is obtained from metadynamics simulations with five commonly used force fields and from ab initio density functional theory calculations in both gas phase and aqueous solution. The FESs obtained at the B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d,p) level of theory are validated by comparison with previously reported MP2 and LMP2 results as well as with experimentally obtained probability distribution between the C5-β (or β-PPII) and αR states. A quantitative assessment is made for each force field in three conformational basins, LeRI (C5-β-C7eq), LeRII (β2-αR), and LeRIII(αL-C7ax-αD) as well as three transition-state regions linking the above conformational basins. The performance of each force field is evaluated in terms of the average free energy of each region in comparison with that of the ab initio results. We quantify how well a force field FES matches the ab initio FES through the calculation of the standard deviation of a free energy difference map between the two FESs. The results indicate that the performance varies largely from region to region or from force field to force field. Although not one force field is able to outperform all others in all conformational areas, the OPLSAA/L force field gives the best performance overall, followed by OPLSAA and AMBER03. For the three top performers, the average free energies differ from the corresponding ab initio values from within the error range (<0.4 kcal/mol) to ∼1.5 kcal/mol for the low-energy regions and up to ∼2.0 kcal/mol for the transition-state regions. The strategy presented and the results obtained here should be useful for improving the parametrization of force fields targeting both accuracy in the energies of conformers and the transition-state barriers.

Theoretical study on conformational preferences of ribose in 2-thiouridine--the role of the 2'OH group

Conformational changes in ribose are well-known to play a significant role in biomolecular identification. The mechanism of selectivity towards C3'-endo conformation (conformer b) in ribose of 2-thiouridine has been studied using DFT (B3LYP) and MP2 methodology, together with 6-31+G(d,p) basis set. The polarity of the C2S2 bond is enhanced due to the orientation of H2' towards the S2 atoms, which leads to a difference in the corresponding bond lengths, the atomic charges and the vO2'H2' stretch vibrations in all the conformers. NBO analysis shows that charge transfer mainly occurs in the C2N3 and C2S2 orbitals. The higher stability of conformer b is attributed to its larger orbital interaction energies within the 2-thiouracil base, and total orbital interaction energies of conformer b. Our conclusion is that the distant electrostatic rather than hydrogen bonding effects between 2'OH and the S2 atoms play the dominant role in the orbital interaction, and enhance the selectivity towards the C3'-endo conformation of ribose.

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

Optical signatures of molecular dissymmetry: combining theory with experiments to address stereochemical puzzles

Modern chemistry emerged from the quest to describe the three-dimensional structure of molecules: van't Hoff's tetravalent carbon placed symmetry and dissymmetry at the heart of chemistry. In this Account, we explore how modern theory, synthesis, and spectroscopy can be used in concert to elucidate the symmetry and dissymmetry of molecules and their assemblies. Chiroptical spectroscopy, including optical rotatory dispersion (ORD), electronic circular dichroism (ECD), vibrational circular dichroism (VCD), and Raman optical activity (ROA), measures the response of dissymmetric structures to electromagnetic radiation. This response can in turn reveal the arrangement of atoms in space, but deciphering the molecular information encoded in chiroptical spectra requires an effective theoretical approach. Although important correlations between ECD and molecular stereochemistry have existed for some time, a battery of accurate new theoretical methods that link a much wider range of chiroptical spectroscopies to structure have emerged over the past decade. The promise of this field is considerable: theory and spectroscopy can assist in assigning the relative and absolute configurations of complex products, revealing the structure of noncovalent aggregates, defining metrics for molecular diversity based on polarization response, and designing chirally imprinted nanomaterials. The physical organic chemistry of chirality is fascinating in its own right: defining atomic and group contributions to optical rotation (OR) is now possible. Although the common expectation is that chiroptical response is determined solely by a chiral solute's electronic structure in a given environment, chiral imprinting effects on the surrounding medium and molecular assembly can, in fact, dominate the chiroptical signatures. The theoretical interpretation of chiroptical markers is challenging because the optical properties are subtle, resulting from the strong electric dipole and the weaker electric quadrupole and magnetic dipole perturbations by the electromagnetic field. Moreover, OR arises from a combination of nearly canceling contributions to the electronic response. Indeed, the challenge posed by the chiroptical properties delayed the advent of even qualitatively accurate descriptions for some chiroptical signatures until the past decade when, for example, prediction of the observed sign of experimental OR became accessible to theory. The computation of chiroptical signatures, in close coordination with synthesis and spectroscopy, provides a powerful framework to diagnose and interpret the dissymmetry of chemical structures and molecular assemblies. Chiroptical theory now produces new schemes to elucidate structure, to describe the specific molecular sources of chiroptical signatures, and to assist in our understanding of how dissymmetry is templated and propagated in the condensed phase.