Enthalpies of hydration of N-methylacetamide by one, two, and three waters and the effect upon the C=O stretching frequency. An Ab initio DFT study

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.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

N-Methylacetamide Aqueous Solutions: A Neutron Diffraction Study

The hydration of N-methylacetamide (NMA) in solution has been determined by neutron diffraction with isotopic Hydrogen/Deuterium substitution (NDIS), augmented by Monte Carlo simulation. This study is representative of the hydration of the peptide bonds characteristic of proteins and might shed light on aggregation phenomena in intrinsically disordered proteins. It is found that NMA forms hydrogen bonds with water at both O and H peptide sites, although of different lengths and strengths. The comparison with the case of tripeptide glutathione evidences differences in both hydration and propensity for aggregation.

Insights into the self-assembly steps of cyanuric acid toward rosette motifs: a DFT study

The nature of non-covalent interactions in self-assembling systems is a topic that has aroused great attention in literature. In this field, the 1,3,5-triazinane-2,4,6-trione or cyanuric acid (CA) is one of the most widely used molecules to formulate self-assembled materials or monolayers. In the present work, a variety of molecular aggregates of CA are examined using three different DFT functionals (B3LYP, B3LYP-D3, and ω-B97XD) in the framework of the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Herein, a step by step aggregation path is proposed and the origin of cooperative effects is also examined. It is shown that a greater cooperativity is not always associated with a greater binding energy, and the greatest cooperative effect occurs with highly directional hydrogen bonds. The intramolecular charge transfers play a key role in this effect. Graphical abstract The noncovalent interactions in cyanuric acid supramolecules were analyzed. The calculations provide insights into the self-assembly steps from dimers to rosette-like motif. The complexes with collinear hydrogen bonds show positive cooperativity, while in the arrangement with double hydrogen bonds the cooperative effect is essentially zero.

Anharmonic vibrations of N-H in Cl(-)(N-methylacetamide)1(H2O)(0-2)Ar2 cluster ions. Combined IRPD experiments and BOMD simulations

Infrared Predissociation (IRPD) spectra of Cl(-)(NMA)1(H2O)0-2Ar2 combined with Born-Oppenheimer Molecular Dynamics (BOMD) IR spectra have been acquired, providing the structure and dynamics of these systems. We show that the chloride ion is bound to the hydrogen of the amide N-H group, forming a strong ionic hydrogen bond, weakening the N-H stretch, and shifting it to lower frequency. The presence of water molecules enhances the ionic hydrogen bond by binding to the amide carbonyl oxygen of NMA and shifts the N-H stretch further to lower frequency. The BOMD IR spectra can recapture all, but about 100 cm(-1), of the 600 to 700 cm(-1) shifts due to the strong N-H stretch anharmonicities observed in experiments. This residual error was found to be due to the lack of zero point energy in the classical treatment of motion in the BOMD method.

"Additive" cooperativity of hydrogen bonds in complexes of catechol with proton acceptors in the gas phase: FTIR spectroscopy and quantum chemical calculations

Experimental study of hydrogen bond cooperativity in hetero-complexes in the gas phase was carried out by IR-spectroscopy method. Stretching vibration frequencies of O-H groups in phenol and catechol molecules as well as of their complexes with nitriles and ethers were determined in the gas phase using a specially designed cell. O-H groups experimental frequency shifts in the complexes of catechol induced by the formation of intermolecular hydrogen bonds are significantly higher than in the complexes of phenol due to the hydrogen bond cooperativity. It was shown that the cooperativity factors of hydrogen bonds in the complexes of catechol with nitriles and ethers in the gas phase are approximately the same. Quantum chemical calculations of the studied systems have been performed using density functional theory (DFT) methods. It was shown, that theoretically obtained cooperativity factors of hydrogen bonds in the complexes of catechol with proton acceptors are in good agreement with experimental values. Cooperative effects lead to a strengthening of intermolecular hydrogen bonds in the complexes of catechol on about 30%, despite the significant difference in the proton acceptor ability of the bases. The analysis within quantum theory of atoms in molecules was carried out for the explanation of this fact.

Aqueous solvation of polyalanine α-helices with specific water molecules and with the CPCM and SM5.2 aqueous continuum models using density functional theory

We present density functional theory (DFT) calculations at the X3LYP/D95(d,p) level on the solvation of polyalanine α-helices in water. The study includes the effects of discrete water molecules and the CPCM and AMSOL SM5.2 solvent continuum model both separately and in combination. We find that individual water molecules cooperatively hydrogen-bond to both the C- and N-termini of the helix, which results in increases in the dipole moment of the helix/water complex to more than the vector sum of their individual dipole moments. These waters are found to be more stable than in bulk solvent. On the other hand, individual water molecules that interact with the backbone lower the dipole moment of the helix/water complex to below that of the helix itself. Small clusters of waters at the termini increase the dipole moments of the helix/water aggregates, but the effect diminishes as more waters are added. We discuss the somewhat complex behavior of the helix with the discrete waters in the continuum models.

A Density Functional Theory Evaluation of Hydrophobic Solvation: Ne, Ar and Kr in a 50-Water Cluster. Implications for the Hydrophobic Effect

The physical explanation for the hydrophobic effect has been the subject of disagreement. Physical organic chemists tend to use a explanation related to pressure, while many biochemists prefer an explanation that involves decreased entropy of the aqueous solvent. We present DFT calculations at the B3LYP/6-31G(d,p) and X3LYP/6-31G(d,p) levels on the solvation of three noble gases (Ne, Ar, and Kr) in clusters of 50 waters. Vibrational analyses show no substantial decreases in the vibrational entropies of the waters in any of the three clusters. The observed positive free energies of transfer from the gas phase or from nonpolar solvents to water appear to be due to the work needed to make a suitable hole in the aqueous solvent. We distinguish between hydrophobic solvations (explicitly studied here) and the hydrophobic effect that occurs when a solute (or transition state) can decrease its volume through conformational change (which is not possible for the noble gases).

The effect of polarization on multiple hydrogen-bond formation in models of self-assembling materials

We report density functional theory calculations at the B3LYP/D95(d,p) level on several different cyclic H-bonding dimers, where the monomers of each are connected by a pair of N-H···O=C H-bonding interactions, and the H-bonding donors and acceptors on each monomer are separated by polarizable spacers. Depending on the structures, the individual H-bonds vary in strength (enthalpy) by over a factor of four, from 2.41 to 10.99 kcal/mol. We attribute most of the variation in interaction energies to differences in the extent of polarization due to each of the H-bonds, which can either combine constructively or destructively. The dipole-dipole interactions between the pair of H-bonds also contribute somewhat to the relative stabilities. The relevance of these results to the design of self-assembling materials is discussed.

FTIR study of H-bonds cooperativity in complexes of 1,2-dihydroxybenzene with proton acceptors in aprotic solvents: influence of the intramolecular hydrogen bond

FTIR spectroscopic study of hydrogen bonding of 1,2-dihydroxybenzene (catechol) with proton acceptors has been carried out. The influence of intramolecular and intermolecular hydrogen bonds on the strengths of each other in complexes of 1,2-dihydroxybenzene with various proton acceptors has been analyzed. It was shown that intramolecular hydrogen bond is strengthened when 1,2-dihydroxybenzene interacts with bases (ethers, amines, nitriles, etc.) in inert solvents. The contribution of the cooperativity of intramolecular hydrogen bonds in the frequency of stretching vibrations of O-H groups linearly depends on the proton acceptor ability of the bases. The solvent effect on hydrogen bond cooperativity in 1,2-dihydroxybenzene-base complexes has been studied. The approach to determine the influence of cooperative effects on the formation of intermolecular complexes with 1,2-dihydroxybenzene is proposed. It was shown that the strength of intramolecular hydrogen bonds in the complexes of 1,2-dihydroxybenzene with bases due to cooperativity of interactions increases by 30-70%, and the strength of intermolecular hydrogen bond by 7-22%.

Estimation of N-H...O=C intramolecular hydrogen bond energy in polypeptides

The previously proposed molecular tailoring approached (MTA) [Deshmukh, M. M.; Gadre, S. R.; Bartolotti, L. J. J. Phys. Chem. A 2006, 110, 12519] for the estimation of intramolecular O-H...O hydrogen bond energy is extended to that for the N-H...O=C bond within polypeptides. The methodology is initially tested on a tetrapeptide containing two types of N-H...O=C hydrogen bonds and is found to distinguish between them. The estimated values are in good agreement with the trends predicted by the geometrical parameters. Furthermore, this methodology is applied to partially as well as fully substituted, capped polyglycines that contain five glycine residues (acetyl-(gly)(5)-NH(2)) to check the effect of substituents on the energetics of hydrogen bonds. The estimated N-H...O=C bond energy values lie in the range of 4-6 kcal/mol. These estimated values are not only in concurrence with the geometric parameters but also able to reflect the subtle effects of substituents for the substituted polypeptides studied in the present work.

Ramachandran revisited. DFT energy surfaces of diastereomeric trialanine peptides in the gas phase and aqueous solution

We report DFT calculations at the B3LYP/D95(d,p) level on the gas phase, aqueous solvation and solvated energies as functions of the central psi and phi dihedral angles (in steps of 5 degrees each) of acetyl-(L)Ala-(L)Ala-(L)Ala-NH(2) (3AL) and its diastereomer, acetyl-(L)Ala-(D)Ala-(L)Ala-NH(2) (3AD). In addition to structures without internal H-bonds (C(5) interactions are neglected), many (95) structures containing internal H-bonds were completely optimized. The only minima for non-H-bonding structures in the gas phase correspond to extended beta-strands for both diastereomers. Some (but not all) structures with internal H-bonds are more stable than those without them. The energy landscapes for the solvated species show multiple minima for the non-H-bonding species and a single minimum for the H-bonding species (3(10)-helix), suggesting that the equilibrium conformational mixture in water be composed of the extended beta-strand, polyproline II, 3(10)-helix, and alpha-helix-like (with no H-bonds) conformations which are all within about 1 kcal/mol of each other. Most H-bonding structures are destabilized relative to the non-H-bonding structures in aqueous solution, but some with large dipole moments are not. The large dipole moment of the alpha-helix-like conformation leads to its increased stability in water (vs the gas phase). Significant qualitative and quantitative differences are reported for the energy landscapes of the two diastereomers when one is compared with the mirror image of the other landscape (particularly in the beta-turn region), suggesting that the differences in the energies of the unfolded peptides need to be considered when considering the stabilities of folded peptides and proteins with single amino acid mutations.

Solvent effect on H-bond cooperativity factors in ternary complexes of methanol, octan-1-ol, 2,2,2-trifluoroethanol with some bases

Cooperative hydrogen bonds in ternary complexes (ROH)(2)...B (ROH-alcohols; B-bases) formed in pure bases (B) and solutions in n-hexane, carbon tetrachloride, benzene and 1,2-dichloroethane were studied by FTIR spectroscopy. Based on the observations, the authors were able to propose an original method of evaluating solvent effects on cooperativity factors in the complexes. Frequencies of cooperative hydrogen bonds OH...B (nu(b)) were determined for ternary complexes of pyridine with aliphatic alcohols (methanol, octan-1-ol) and for 2,2,2-trifluoroethanol with three different bases (acetonitrile, diethyl ether, tetrahydrofuran). The solvent shifts of nu(b) were found to correlate with an empirical thermochemical parameter of the solvent, S(VW). The cooperativity factors were determined for the complexes (ROH)(2)...B in all studied media. It has been found that the cooperativity factors are almost independent of the solvent. In addition, a method was proposed of estimating the frequencies and cooperativity factors for ternary complexes (ROH)(2)...B in the gas phase. It has been found that in gas phase the cooperativity factors are practically the same as in condensed media.

Fluoroolefins as peptide mimetics: a computational study of structure, charge distribution, hydration, and hydrogen bonding

The design of peptide mimetic compounds is greatly facilitated by the identification of functionalities that can act as peptide replacements. The fluoroalkene moiety has recently been employed for that purpose. The purpose of this work is to characterize prototypical fluoroalkenes (fluoroethylene and 2-fluoro-2-butene) with respect to key properties of peptides (amides) including structure, charge distribution, hydration, and hydrogen bonding. The results are compared to those obtained for model peptides (formamide, N-methylacetamide). Calculations have been carried out at the MP2 and B3LYP levels of theory with the 6-311++G(2d,p) and 6-311++G(2d,2p) basis sets. The results suggest that the fluoroalkene is similar in steric requirements to a peptide bond but that there is less charge separation. Calculations of the hydration free energies with the PCM bulk continuum solvent model indicate that the fluoroalkene has much smaller hydration free energies than an amide but that the difference in solvation free energy for cis and trans isomers is comparable. In studies of complexes with water molecules, the fluoroalkene is found to engage in interactions that are analogous to backbone hydrogen-bonding interactions that govern many properties of natural peptides and proteins but with smaller interaction energies. In addition, key structural differences are noted when the fluoroalkene is playing the role of hydrogen-bond acceptor which may have implications in binding, aggregation, and conformational preferences in fluoroalkene peptidomimetics. The issue of cooperativity in hydrogen-bonding interactions in complexes with multiple waters has also been investigated. The fluoroalkene is found to exhibit cooperative effects that mirror those of the peptide but are smaller in magnitude. Thus, pairwise addivitity of interactions appears to more adequately describe the fluoroalkenes than the peptides they are intended to mimic.

Modeling competitive interactions in proteins: vibrational spectroscopy of M+(n-methylacetamide)1(H2O)n=0-3, M=Na and K, in the 3 microm region

To properly understand the preferred structures and biological properties of proteins, it is important to understand how they are influenced by their immediate environment. Competitive intrapeptide, peptide...water, ion...water, and ion...peptide interactions, such as hydrogen bonding, play a key role in determining the structures, properties, and functionality of proteins. The primary types of hydrogen bonding involving proteins are intramolecular amide...amide (N-H...O=C) and intermolecular amide...water (O-H...O=C and H-O...H-N). n-Methylacetamide (NMA) is a convenient model for investigating these competitive interactions. An analysis of the IR photodissociation (IRPD) spectra of M+(n-methylacetamide)1(H2O)n=0-3 (M=Na and K) in the O-H and N-H spectral regions is presented. Ab initio calculations (MP2/cc-pVDZ) are used as a guide in identifying both the type and location of hydrogen bonds present. In larger clusters, where several structural isomers may be present in the molecular beam, ab initio calculations are also used to suggest assignments for the observed spectral features. The results presented offer insight to the nature of ion...NMA interactions in an aqueous environment and reveal how different ion...ligand pairwise interactions direct the extent of water...water and water...NMA hydrogen bonding observed.

Surface change of Ras enabling effector binding monitored in real time at atomic resolution

Ras, the prototype of the Ras superfamily, acts as a molecular switch for cell growth. External growth signals induce a GDP-to-GTP exchange. This modifies the Ras surface (Ras(on)GTP) and enables effector binding, which then activates signal-transduction pathways. GTP hydrolysis, catalysed by Ras and GAP, returns the signal to "off" (Ras(off)GDP). Oncogenic mutations in Ras prevent this hydrolysis, and thereby cause uncontrolled cell growth. In the Ras(off)-to-Ras(on) transition, the Ras surface is changed by a movement of the switch I loop that controls effector binding. We monitored this surface change at atomic resolution in real time by time-resolved FTIR (trFTIR) spectroscopy. In the transition from Ras(off) to Ras(on) a GTP-bound intermediate is now identified, in which effector binding is still prevented (Ras(off)GTP). The loop movement from Ras(off)GTP to Ras(on)GTP was directly monitored by the C=O vibration of Thr35. The structural change creates a binding site with a rate constant of 5 s(-1) at 260 K. A small molecule that shifted the equilibrium from the Ras(on)GTP state towards the Ras(off)GTP state would prevent effector binding, even if hydrolysis were blocked by oncogenic mutations. We present a spectroscopic fingerprint of both states that can be used as an assay in drug screening for such small molecules.

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.

The Effect of Aqueous Solvation upon α-Helix Formation for Polyalanines

The incremental free energies of aqueous solution for acetyl(ala)NNH2 in its extended unfolded and alpha-helical conformations are compared using the SM5.2 solvation method of Cramer and Truhlar. A combination of density functional theory (DFT) at the B3LYP/D95(d,p) and AM1 has been employed using the ONIOM method. The incremental solvation energies of alpha-helical structures are very similar for both ONIOM and AM1 optimized structures as these structures do not significantly change upon solution. However, the conformations of the unfolded peptides change from extended beta-strand to polyproline II conformations upon aqueous solution. The incremental solvation free energy per residue of the polyproline II structure is about 2 kcal/mol/residue greater than that for the alpha-helix, representing an upper limit for the difference between the solvation energies. However, most of this difference disappears when the energy required to distort the optimized gas-phase extended beta-strand structure to the optimized polyproline II solution structure is included in the analysis, leaving an estimated difference in incremental solvation free energy of 0.3-0.5 kcal/mol favoring the unfolded structure. The solution structure sacrifices the stability derived from the intramolecular C5 H-bonds for more favorable interactions with the aqueous solvent.