Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Dipole Cooperativity and Polarization Frustration Determine the Secondary Structure Distribution of Short Alanine Peptides in Water

The physical driving forces for secondary structure preferences of hydrated alanine peptide are investigated with B3LYP-D3(BJ) and the adaptive force matching (AFM) method. The AFM fit to the DFT surface, ALA2022, provides excellent agreement with the nuclear magnetic resonance scalar coupling constants from experiments. In turn, the model is used to gain insight into the physical driving forces behind secondary structure preferences of hydrated peptides. DFT calculations with and without the Conductor-like Screening Model (COSMO) show that the α helix is stabilized by solvent polarization due to dipole cooperativity. The two adjacent amide groups in β strand form a near-planar trapezoid that is not much larger than the size of water molecules. When the finite size of a water molecule is considered, the stabilization from solvent polarization for such a trapezoid is frustrated. Water molecules cannot find orientations to properly stabilize all four polar regions close to each other with such an awkward arrangement. This leads to quite substantial reduction in polarization stabilization. Although the polyproline II (PP-II) conformation is very similar to the β strand, the small twist in the backbone angles allowed much improved polarization stabilization. The improved polarization, when combined with favorable intrapeptide interactions, leads to the PP-II to be lowest in free energy. Other factors, such as the entropic TΔS and the ϕ, ψ coupling terms, are also studied but are found to play only a minor role. The insight shown in this work helps to better understand the structure of globular and intrinsic disordered proteins and facilitate future force field development.

Water model for hydrophobic cavities: structure and energy from quantum-chemical calculations

This ab initio study aims to design a series of large water clusters having a hollow clathrate-like cage able to host hydrophobic solutes of various sizes. Starting from the (H₂O)_n (n = 18, 20, 24 and 28) hollow cages, water layers have been added in a stepwise manner in order to model the configuration of water molecules beyond the primary shell. The large (H₂O)₁₀₀, (H₂O)₁₂₀ and (H₂O)₁₄₀ clusters complete the hydrogen bonding network of the cage with optimal and regular tiling of the do-, tetra-decahedron and hexa-decahedron, respectively. This study is corroborated by an investigation of dense water clusters up to the (H₂O)₁₂₃ one, being highly consistent with experimental data on ice concerning the electronic and zero-point energies for aggregate formation at 0 K and enthalpy and entropy at 273 K. The cavity creation profoundly alters the orientation of water molecules compared with those found in dense clusters. Nevertheless, such a large reorganization is necessary to maximize the water-water attraction by making it similar to the one found in dense clusters. The cage formation is an endothermic process; however, the computed values are large compared with previous reports for hydrocarbon aqueous solutions. Larger clusters are required for a more fruitful comparison.

Forbidden Secondary Structures Found in Gel-Forming Fibrils of Glycylphenylalanylglycine

The zwitterionic l-tripeptide glycylphenylalanylglycine self-assembles into very long crystalline fibrils in an aqueous solution, which causes the formation of an exceptionally strong gel phase (G' ∼ 5 × 10⁶ Pa). The Rietveld refinement analysis of its powder X-ray diffraction (PXRD) pattern reveals a unit cell with four peptides forming a P2₁2₁2₁ space group and adopting an inverse polyproline II conformation, that is, a right-handed helical structure that occupies the "forbidden" region of the Ramachandran plot. This unusual structure is stabilized by a plethora of intermolecular interactions facilitated by the large number of different functional groups of the unblocked tripeptide. Comparisons of simulated and experimental Fourier transform infrared and vibrational circular dichroism (VCD) amide I' profiles corroborate the PXRD structure. Our experimental setup reduces the sample to a quasi-two-dimensional network of fibrils. We exploited the influence of this reduced dimensionality on the amide I VCD to identify the main fibril axis. We demonstrate that PXRD, vibrational spectroscopy, and amide I simulations provide a powerful toolset for secondary structure and fibril axis determination.

On the size, shape and energetics of the hydration shell around alkanes

A large number of clathrate-like cages have been proposed as the very first hydration shell of alkanes. The cages include canonical structures commonly found in clathrate hydrates and many others, not previously reported, derived from the carbon fullerene cavities. These structures have a rich and variegated form, which can adapt to the shape and conformation of the solute. They avoid "wasting" hydrogen bonds, while minimizing the volume cage and maximizing the solute-solvent van der Waals interactions. DFT/M06-2X and MP2 ab initio calculations give comparable structural and energetic results although the latter predicts slightly larger cages for a given solute. It is shown that the van der Waals interactions are substantial and the large exoenergetic values found for isobutane and cyclopentane provide an explanation for the surprising high melting points of related hydrates at room pressure. The encaging enthalpy for various hydrocarbons is similar to the enthalpy of solution measured at a temperature just above the melting point of aqueous hydrocarbon solutions, thus indicating that water molecules should not deviate too much from the configuration with O-H bonds tangentially oriented with respect to the solute surface. The computed trend differs from the enthalpy of solution measured at room temperature, thus the very first hydration shell departs, up to a certain degree, from the clathrate-like structures.

Water-Mediated Electronic Structure of Oligopeptides Probed by Their UV Circular Dichroism, Absorption Spectra, and Time-Dependent DFT Calculations

We investigate the UV absorption spectra of a series of cationic GxG peptides (where x denotes a guest residue) in aqueous solution and find that only a subset of these spectra show a strong dependence with temperature. To explore whether or not this observation reflects conformational dependencies, we carry out time-dependent density functional calculations for the polyproline II (pPII) and β-strand conformations in implicit and explicit water. We find that the calculated CD spectra for pPII can qualitatively account for the experimental spectra irrespective of the water model. The β-strand UV-CD spectra, however, require the explicit consideration of water. Contrary to conventional wisdom, we find that both the NV₁ and NV₂ band are the envelopes of contributions from multiple transitions that involve more than just the HOMOs and LUMOs of the peptide groups. A natural transition orbital analysis reveals that some of the transitions have a charge-transfer character. The overall manifold of transitions depends on the peptide's backbone conformation, peptide hydration, and side chain of the guest residue. Our results reveal that peptide groups, side chains, and hydration shells must be considered as an entity for a physically valid characterization of UV absorbance and circular dichroism.

Quantum Mechanics Study on Hydrophilic and Hydrophobic Interactions in the Trivaline-Water System

With the aim to elucidate hydrophobic effects in the unfolded state of peptides, DFT-M062X computations on the Val₃H⁺· nH₂O ( n up to 22) clusters have been accomplished. As far as the main chain is concerned, four conformers with β-strand and/or polyproline type II conformations, PPII (indicated as β-β, β-PPII, PPII-β, and PPII-PPII), have been found by changing the ϕ and ψ angles. For bare peptide, the side chain (isopropyl) of each residue can independently take on three different orientations with negligible effects on energetics. The great isopropyl spatial separations in β-β and β-PPII conformers allow for the construction of synergic and extensive water-water and water-peptide H-bonding in the minimal hydration Val₃H⁺·22H₂O models without significant steric encumbrance. Conversely, due to the proximity of the isopropyl of the central residue with the other two, some restrictions in the water shell construction around the peptide become evident for the PPII-PPII conformer and the number of energetically accessible structures decreases. This is indicative of correlated motion involving isopropyls and backbone mediated by water molecules, the origin of the nearest neighbor effects. Comparing the thermodynamic data of Ala₃H⁺·22H₂O and Val₃H⁺·22H₂O, what emerges is that both hydration enthalpy and entropy drive the β-strand stability of the latter.

Effect of the solvent on the conformational behavior of the alanine dipeptide deduced from MD simulations

In general, peptides do not exhibit a well-defined conformational profile in solution. However, despite the experimental blurred picture associated with their structure, compelling spectroscopic evidence shows that peptides exhibit local order. The conformational profile of a peptide is the result of a balance between intramolecular interactions between different atoms of the molecule and intermolecular interactions between atoms of the molecule and the solvent. Accordingly, the conformational profile of a peptide will change upon the properties of the solvent it is soaked. To get insight into the balance between intra- and intermolecular interactions on the conformational preferences of the peptide backbone we have studied the conformational profile of the alanine dipeptide in diverse solvents using molecular dynamics as sampling technique. Solvents studied include chloroform, methanol, dimethyl sulfoxide, water and N-methylacetamide. Different treatments of the solvent have been studied in the present work including explicit solvent molecules, a generalized Born model and using the bulk dielectric constant of the solvent. The diverse calculations identify four major conformations with different populations in the diverse solvents: the C₇^eq only sampled in chloroform; the C₅ or extended conformation; the polyproline (PII) conformation and the right-handed α-helix conformation (α_R). The results of present calculations permit to analyze how the balance between intra- and intermolecular interactions explains the populations of the diverse conformations observed.

Quantum Mechanics Approach to Hydration Energies and Structures of Alanine and Dialanine

A systematic approach to the phenomena related to hydration of biomolecules is reported at the state of the art of electronic-structure methods. Large-scale CCSD(T), MP4-SDQ, MP2, and DFT(M06-2X) calculations for some hydrated complexes of alanine and dialanine (Ala⋅13 H₂ O, Ala₂ H⁺ ⋅18 H₂ O, and Ala₂ ⋅18 H₂ O) are compared with experimental data and other elaborate modeling to assess the reliability of a simple bottom-up approach. The inclusion of a minimal number of water molecules for microhydration of the polar groups together with the polarizable continuum model is sufficient to reproduce the relative bulk thermodynamic functions of the considered biomolecules. These quantities depend on the adopted electronic-structure method, which should be chosen with great care. Nevertheless, the computationally feasible MP2 and M06-2X functionals with the aug-cc-pVTZ basis set satisfactorily reproduce values derived by high-level CCSD(T) and MP4-SDQ methods, and thus they are suitable for future developments of more elaborate and hence more biochemically significant peptides.

Conformational analysis of short polar side-chain amino-acids through umbrella sampling and DFT calculations

Molecular and quantum mechanics calculations were carried out in a series of tripeptides (GXG, where X = D, N and C) as models of the unfolded states of proteins. The selected central amino acids, especially aspartic acid (D) and asparagine (N) are known to present significant average conformations in partially allowed areas of the Ramachandran plot, which have been suggested to be important in unfolded protein regions. In this report, we present the calculation of the propensity values through an umbrella sampling procedure in combination with the calculation of the NMR J-coupling constants obtained by a DFT model. The experimental NMR observations can be reasonably explained in terms of a conformational distribution where PP_II and β basins sum up propensities above 0.9. The conformational analysis of the side chain dihedral angle (χ₁), along with the computation of ³J(H^αH^β), revealed a preference for the g _- and g ₊ rotamers. These may be connected with the presence of intermolecular H-bonding and carbonyl-carbonyl interactions sampled in the PP_II and β basins. Taking into account all those results, it can be established that these residues show a similar behavior to other amino acids in short peptides regarding backbone φ,ψ dihedral angle distribution, in agreement with some experimental analysis of capped dipeptides.

Assessing backbone solvation effects in the conformational propensities of amino acid residues in unfolded peptides

Conformational ensembles of individual amino acid residues within model GxG peptides (x representing different amino acid residues) are dominated by a mixture of polyproline II (pPII) and β-strand like conformations. We recently discovered rather substantial differences between the enthalpic and entropic contributions to this equilibrium for different amino acid residues. Isoleucine and valine exceed all other amino acid residues in terms of their rather large enthalpic stabilization and entropic destabilization of polyproline II. In order to shed light on these underlying physical mechanisms, we performed high-level DFT calculations to explore the energetics of four representative GxG peptides where x = alanine (A), leucine (L), valine (V), and isoleucine (I) in explicit water (10 H2O molecules with a polarizable continuum water model) and in vacuo. We found that the large energetic contributions to the stabilization of pPII result, to a major extent, from peptide-water, water-water interactions, and changes of the solvent self-energy. Differences between the peptide-solvent interaction energies of hydration in pPII and β-strand peptides are particularly important for the pPII ⇌ β equilibria of the more aliphatic peptides GIG and GLG. Furthermore, we performed a vibrational analysis of the four peptides in both conformations and discovered a rather substantial mixing between water motions and peptide vibrations below 700 cm(-1). We found that the respective vibrational entropies are substantially different for the considered conformations, and their contributions to the Gibbs/Helmholtz energy stabilize β-strand conformations. Taken together, our results underscore the notion of the solvent being the predominant determinant of peptide (and protein) conformations in the unfolded state.

Water-Centered Interpretation of Intrinsic pPII Propensities of Amino Acid Residues: In Vitro-Driven Molecular Dynamics Study

Amino acid residues of unfolded peptides in water sample only a few basins in the Ramachandran plot, including prominent polyproline II-like (pPII) conformations. Dynamics of guest residues, X, in GXG peptides in water were recently reported to be dominated by pPII and β-strand-like (β) conformations, resulting in an enthalpy-entropy compensation at ∼300 K. Using molecular dynamics (MD) in explicit solvent, we here examine pPII and β conformational ensembles of 15 guest residues in GXG peptides, quantify local orientation of water around their side chains through novel water orientation plots, and study their hydration and hydrogen bonding properties. We show that pPII and β ensembles are characterized by distinct water orientations: pPII ensembles are associated with an increased population of water oriented in parallel to the side chain surface whereas β ensembles exhibit more heterogeneous water orientations. The backbone hydration is significantly higher in pPII than in β ensembles. Importantly, pPII to β hydration differences and the solvent accessible surface area of Cβ hydrogens both correlate with experimental pPII propensities. We propose that pPII conformations are stabilized by a local, hydrogen-bonded clathrate-like water structure and that residue-specific intrinsic pPII propensities reflect distinct abilities of side chains to template this water structure.

Interfacial water at the trialanine hydrophilic surface: a DFT electronic structure and bottom-up investigation

A model describing a network of hydrogen bonded water-trialanine has been developed to estimate hydration effects on various conformers of the peptide.