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 b-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 b-strand peptides are particularly important for the pPII " b 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 . We found that the respective vibrational entropies are substantially different for the considered conformations, and their contributions to the Gibbs/Helmholtz energy stabilize b-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.


Introduction
Our current understanding of protein folding processes relies on a detailed knowledge of the conformations they can adopt in their folded and unfolded states.2][3][4] This predominant view is largely based on Flory's classical independent site model, which allows the various protein residues to sample the entire allowed region of the Ramachandran space. 2 However, over the last 15 years, considerable evidence has accumulated suggesting that significantly more restricted conformational ensembles of individual amino acid residues are present in unfolded peptides [5][6][7][8][9][10][11][12] and restricted coil libraries. 13,14Specifically, amino acid residues of short peptides in solution as well as in truncated coil libraries were found to exhibit much more restricted conformational sampling and to differ in terms of their conformational propensities.A recent conformational analysis of GxG peptides in water based on vibrational spectroscopic and NMR data revealed that most residues do predominantly sample the upper left quadrant of the Ramachandran plot, which can be subdivided into two sub-conformational ensembles associated with the b-strand and less extended polyproline II (pPII) conformations. 12,15,16ndividual amino acid residues differ in terms of the equilibrium between these two conformations -while alanine shows a clear preference for pPII, b-strand like conformations are slightly more preferred for valine and isoleucine.The DG pPII-b values This journal is © the Owner Societies 2015 associated with these equilibria range between À3.5 and 2 kJ mol À1 , with most values lying in the RÁT energy range. 17owever, these DG pPII-b values were found to be composed of much larger DH and DS values, which also exhibit much larger variations with changing residues than DG.For two subsets of data, the analysis resulted in isoequilibria, e.g.temperatures at which the DG values of the different members of the subset were practically identical.This was interpreted as an indication that peptide-water interactions are the key determinant of the residue's propensities.In this context it was particular striking that valine and isoleucine were found to exhibit rather large DH and TDS values (in the 40-60 kJ mol À1 range at room temperature) that exceed those of other residues with aromatic or aliphatic side chains.The fact that the DH and TDS values of isoleucine exceed the corresponding values for the isomeric amino acid leucine by a factor of B2.5 was particular astonishing and has thus far been rationalized, in general terms, as an intricate interplay between backbone and side chain solvation. 17However, these results cannot be understood in terms of solvent accessible surface areas that are normally considered to be the decisive parameters in side chain-solvent interactions.Specifically, since isoleucine (I) and leucine (L) have practically the same accessible surface area, which exceeds that of valine (V), 18 this suggests a hierarchy with I, L 4 V, contrary to the observed I 4 V c L hierarchy.
Multiple lines of evidence suggest that interactions between unfolded peptides/proteins and the solvent are pivotal for conformational preferences exhibited by individual amino acid residues.0][21][22][23][24][25][26][27][28] In non-aqueous solvents like DMSO and primary alcohols, the sampling of pPII is generally reduced or even absent. 25,28The role of water in stabilizing the pPII conformation of alanine residues was first proposed by Han et al. based on density functional theory (DFT) calculations of N-acetyl-L-alanine N 0 -methylamide (alanine dipeptide) in complexes with four water molecules. 19These authors attributed the stabilization of pPII to water bridges between CO and NH groups of the two peptide linkages of the molecule.Garcia, by means of molecular dynamics (MD) simulations with a modified AMBER force field, suggested that pPII conformations allowed for an optimal packing of water molecules in the hydration shell. 21Fleming et al. proposed a correlation between pPII propensity and changes in the solvent accessibility area of side chains associated with transitions between pPII and b-strand like conformations. 18However, this notion is at odds with results of MD studies reported by Mezei et al., who found instead that the pPII preference of alanine over b-strand results from stronger backbone-water interactions in the former. 24omputational studies by Drozdov et al. 20 as well as Avbelj and Baldwin 29,30 led these authors to suggest a more indirect role of solvation in the stabilization of pPII, namely the screening of electrostatic interactions between peptide groups, which would produce a preference for more extended b-strand conformations.Amino residues with side chains bulkier than alanine were predicted to perturb the hydration shell and thus destabilize pPII, a notion consistent with experimental results. 31A recent DFT-based study by Lanza and Chiacchio on N-acetyl-L-alanine amide complexes with up to 13 water molecules revealed the existence of water clusters that were hydrogen bonded to the functional groups of the peptide backbone (i.e.CO and NH).They concluded that pPII conformations exhibit an energetic preference that leads to a greater reduction of water mobility, 32 further suggesting that pPII is energetically favored and entropically disfavored, in agreement with experimental data. 31 more recent study by these authors confirmed this view in which the number of water molecules in the hydration shell was substantially increased up to 37. 33 A somewhat lower number of water molecules (22) was found to be necessary for obtaining stabilization energies of pPII relative to the b-strand that are comparable with experimentally obtained values.
The present study is aimed at shedding some light on the influence of side chains on the backbone hydration of GxG peptides.Toal et al. recently reported rather diverse enthalpy and entropy values for the pPII " b-strand equilibria of fifteen GxG peptides in aqueous solution. 31To this end, we performed several DFT calculations on a selected set of the earlier investigated GxG peptides surrounded by 10 explicit water molecules with a polarizable continuum water model.As representative peptides, we selected a set of aliphatic amino acid residues, namely, alanine (A), leucine (L), valine (V), and isoleucine (I) for our studies, as these residues showed the largest variance in thermodynamic values. 179,30,32,[34][35][36][37] The other three aliphatic side residues were selected because the above thermodynamic analysis has yielded large and surprisingly different DH and TDS values for their respective pPII " b equilibria.The experimental DH and DS values for these peptides as reported by Toal et al. are displayed in Fig. 1. 31 In this study, we focus on determining the energetics of the pPII " b equilibria for the investigated peptides rather than on a thorough exploration of the energy landscapes of these peptides.In view of the established preponderance of pPII and b-strand conformations in their Ramachandran plots, this restriction is justified. 12,15,16eaction field.The initial geometries of the water complexes were taken from Lanza and Chiacchio. 32Geometry optimizations and harmonic frequencies at the same level of theory were calculated to verify that these stationary points were local minima.

Energetics of GxG pPII " b equilibria in water
Quantum calculations were performed for the selected cationic GxG peptides in vacuo and in a bath of 10 explicit H 2 O molecules with the entire system embedded in a polarizable continuum solvent model, as described in the Computational methods section.The number of water molecules is substantially less than what has recently been employed by Lanza and Chiacchio. 33owever, as shown in Fig. 1 our approach allowed us to provide a minimal hydration of all functional groups of the peptide backbone, which includes 2 water molecules as hydrogenbonded donors for the two carbonyl groups and one water molecule as a hydrogen-bond acceptor for the two amide protons.The remaining six water molecules were used to hydrate the terminal groups and to connect the different water molecules hydrogen bonded to the peptide.It is important to note that there exists many conformational isomers in each tripeptide, and we have only focused on a single conformer in each of the four GxG tripeptides.However, the goal of this study is the comparison of pPII " b equilibria rather than an accurate calculation of absolute values of energy differences or a thorough sampling of the conformational space of water molecules in the hydration shell.The chosen setup of water molecules does not hydrate the side chains of the chosen x-residues, hence the contributions from side chain solvent interactions are neglected in our study (as they were in the studies of Lanza and Chiacchio 32,33 ).Furthermore, the previous studies of Lanza and Chiacchio investigated both AAA and AdP-type peptides, while our work focuses on the different GxG peptide geometries.While a comparison between our work and Lanza and Chiacchio's study is not entirely straightforward, we were able to confirm that the backbone geometry for our computed GAG structure has a similar topology with their previously published structure for AdP.This comparison, as well as our results described further below, demonstrate that our approach is capable of elucidating the influence of the side chains on backbone hydration.We selected the fully protonated form to directly correspond to the experiments this study is referring to. 12,16,31Since we were solely interested in the pPII " b equilibria of the investigated peptides, we carried out two fully unconstrained geometry optimizations per peptide, one starting in the pPII conformation and the other in the b-strand region.After a full optimization, a normal mode analysis was carried out.These calculations served three purposes: (1) determining the internal energy difference between b and pPII, (2) analyzing, in particular, the manifold of peptide modes that are vibrationally mixed with the hydration shell of 10 water molecules, and (3) checking for the appearance of any imaginary frequencies which would indicate that the observed conformation does not represent a true minimum in the energy landscape of the system.
Tables 1 and 2 list the dihedral angles of the obtained optimized structures, the energy differences between the pPII and b-strand geometries, and the permanent dipole moments for the various GxG systems in vacuo and in explicit water, respectively.Table 2 also lists the dihedral angles of earlier reported maxima of pPII and b-strand sub-distributions obtained from amide I 0 profiles and various J-coupling constants. 12,15,16,47ost of the obtained values are in reasonable and some in excellent agreement with the reported experimentally-derived values.Exceptions are, to a minor extent, the f-value of the b-strand conformation of GLG and, to a major extent, the c-value of the b-strand GVG.The experimentally-obtained Table 1 Dihedral angles f and c, electronic energies, and total permanent dipole moment of the geometry-optimized pPII and b-strand type conformations of GxG peptide + reaction sphere.DU e is the electronic energy difference calculated as DU e (P s b s w n ) described in the main text  values of pPII are generally very well reproduced by our calculations.][50][51][52][53][54] Fig. 1 depicts the conformations resulting from the geometry optimization of the investigated GAG + 10 H 2 O configurations.The optimization process, to a major extent, has maintained the water-mediated hydrogen bonding network and the hydration of the backbone groups, consistent with the initial setup.These conformations resemble what Ben Naim has termed the HfI interaction, which are hydrophilic in character. 55A notable difference is observed for the arrangements of the water molecules above the peptide backbone.In the b-strand conformation, a water-mediated hydrogen bonding loop is maintained between the N-and the C-terminal which involves hydrogen bonding to and from the amide group of the C-terminal and the carbonyl group of the N-terminal peptide.In the pPII conformation, however, this loop is broken owing to the oxygen atom of the C-terminal's water molecule being now hydrogen bonded to the amide proton of the C-terminal peptide.On the contrary, the hydrogen bonding network below the peptide backbone is more extended in the pPII conformation.A distance analysis revealed that the CQOÁ Á ÁH-O-H, NHÁ Á ÁO, and H-O-HÁ Á ÁOH 2 distances are generally shorter in the b-strand than in the corresponding pPII conformations, thus stabilizing the former.Corresponding bending angles are also closer to their ideal 1801 value in the b-strand conformation than in pPII.Such a stabilization of the b-strand conformation by hydrogen bonding was earlier reported by Mezei et al. for a polyalanine peptide in explicit water. 24Lanza and Chiacchio reported a similar finding for their AcANH 2 Á13H 2 O complex. 32ased on the above analyses, one is tempted to suspect an overall energetic stabilization of b-strand over pPII.However, the corresponding electronic energy differences, DU e (pPII,b) = U e (pPII) À U e (b), as plotted in Fig. 2  as internal energies that can generally be decomposed into electronic, vibrational, rotational, and translational contributions.For our system we can certainly ignore translation and rotation and we can further expect that the electronic part exceeds by far the vibrational contribution (vide infra).Hence, it is justified to compare the computational DU e with the experimental DU int values reported by Toal et al. 27 The calculated and experimental DU values are practically identical for GLG, whereas calculated values for GVG and GIG are lower than the experimental values by factors of 2 and 3, respectively.Surprisingly, our calculation leads to an overestimation of DU for GAG; the value we obtained more resembles the experimentally-obtained AAA and alanine dipeptide structures. 27However, our calculated values reproduced the experimental hierarchy of I 4 V 4 L, which is one of the primary goals of this study.The corresponding DU-values for GxG in vacuo are also displayed in Fig. 2. They are all positive, indicating a stabilization of the b-strand conformation and confirming that water is indeed pivotal for pPII stabilization.Interestingly, GLG and GIG in vacuo now exhibit very similar energies, as one would expect from the fact that their side chains are chemical isomers.This suggests that the dissimilarity between pPII/b energy differences obtained (both experimentally and computationally) for hydrated leucine and isoleucine (Fig. 1) is in part a direct result of differences in side-chain and backbone solvation. 31Since the only difference between these two residues is the position of -CH 3 branching on the side-chain (i.e., b-branching in isoleucine and g-branching is leucine), it is likely that the hydration shell about this group is the pivotal determinant.In general, these results strongly corroborate the notion that the unfolded state of peptides in water would be dominated by pPII-like conformations if the favorable hydration energy (Fig. 1) was not compensated and sometimes overcompensated by entropic contributions.They further Table 2 Dihedral angles f and c, electronic energies, and total permanent dipole moment of the geometry-optimized pPII and b-strand type conformations of GxG peptides + 10 H 2 O + polarizable continuum.DU e is the energy difference calculated as DU e (P s b s w n ) described in the main text.Numbers shown in parenthesis reflect the maxima positions of subpopulations observed from experimental data 12,16 f [1]  c  suggest that the pPII " b equilibrium enthalpy/energy of GxG peptides is predominantly determined by peptide-solvent interactions. 28Large differences between DU (and also between corresponding DS values) for different amino acid residues can give rise to large disparities between their respective pPII " b equilibria at conditions significantly above and below their compensation temperatures.This further suggests a sequence dependence of the conformational entropies of unfolded states of peptides and proteins at high (thermal denaturation) and low temperatures (cold denaturation). 56he discrepancy between the calculated energetics and the b-strand stabilization suggested by the hydrogen bond network analysis is surprising, but not unprecedented.Fleming et al., from an MD-based analysis of a 12-residue polyalanine in explicit water, found no evidence for water bridges in pPII at all. 18evertheless, they obtained a strong stabilization of pPII, attributable to peptide-water interactions, in agreement with our findings for all the peptides investigated.The DFT-based optimization of various AcANH 2 ÁnH 2 O complexes by Lanza and Chiacchio yielded water for both conformations, but with more stable ones in b-strands.With regard to the total energy, however, hydration was found to stabilize pPII energetically. 32The underlying reason for this overall stabilization of pPII remains elusive in all these studies.
Taken together, our analysis suggests that it is difficult to pinpoint the reason for the enthalpic (energetic) stabilization of pPII in water.It is clear from our study and previous experimental results that the interaction between the peptide and the first hydration shell is pivotal in this regard.Furthermore, the present calculations confirm the notion of Toal et al. 31 that the hydration energy difference between pPII and b-strand is heavily side chain dependent.However, the net electronic energies that emerged from our calculations are likely to be a superposition of competing contributions, which might involve through-space electrostatic in addition to through-hydrogen-bond interactions.To quantify these various interactions, we carried out a detailed analysis of the configurations for all four peptides using the energy-decomposition analysis of Mirkin and Krimm. 57Within this energy-decomposition scheme, one can rigorously define the total energy of a composite system as the sum of the energies of its individual components plus the interaction energy between them.For example, the energy U e ( j s w n ), of the conformation j = P, b in explicit water is given by U e ( j s w n ) = U e ( j s ) + U e, j (w n ) + U e,i ( j), (1)   where U e ( j s ) is the energy of the peptide geometry adopting the state j in the solvated system, U e, j (w n ) is the self-energy of the n water molecules surrounding a peptide adopting the conformation j, and U e,i ( j) represents the total intermolecular interaction energies in the composite state j (which includes water-water and water-peptide hydrogen bonds as well as all electrostatic interaction effects).Thus, the difference between pPII and b-strand energies of the solvated systems is given by DU e (P s b s w n ) = DU e (P s b s ) + DU e,Pb (w n ) + DU e,i (Pb), where  2) can be directly obtained from DFT total energy calculations of the composite peptide-water system, the peptide geometry, and the water geometries, respectively.As the original energy-decomposition analysis by Mirkin and Krimm 57 does not account for basis set superposition errors, we neglect these contributions which is also justified by our fairly large cc-pVTZ basis set.With all of the energies in eqn (2) determined, the intermolecular interaction energy DU e , i (Pb) can be obtained as well.To further understand how our results would change with other basis sets, we also carried out this analysis with the smaller cc-pVDZ basis (see ESI †).Upon optimization of the various GxG peptide structures with the cc-pVDZ basis, we found that several of the geometries exhibited significant distortions for both the f and c dihedral angles.As such, these findings highlight the importance of using fairly large basis sets, and we only discuss the oB97X-D/ cc-pVTZ results in the main text.Using this decomposition scheme, we calculated the various energy contributions to the obtained DU e (P s b s w n ) values for all 4 of the investigated peptides, which are summarized in Table 3.All of the DU e (P s b s ) values are positive (cf.Fig. 2), indicating that the b conformation in the geometry of the solvated system is more stable (the latter statement is also true for the b conformation in vacuo).This suggests that the consideration of only indirect solvation is insufficient to describe pPII preferences.The intermolecular interaction energy differences, DU e,i (Pb) as well as the solvent selfenergies (DU e , Pb (w n )) are all negative, and both contribute significantly to the stabilization of the pPII structure.The latter are nearly identical for GLG, GVG, and GIG, while the value for GAG is only slightly lower.Apparently, the obtained differences DU e (P s b s w n ) between the total energies, particularly between those of GLG, GVG, and GIG result from the more residue-specific peptide-water and water-water interactions.With regard to GAG, our results are at variance with the findings of Mirkin and Krimm, who identified the solvent energy as the main contributor to the stabilization of the pPII conformation of an alanine dipeptide-like compound. 57

Estimation of vibrational entropies and energies
Generally, DFT energy calculations do not provide any direct information about entropies and, therefore, do not enable one to gauge relative stabilities of conformations.However, we wondered to what  extent the hydration shell could contribute to entropic differences between the considered conformations solely through peptide-water and water-water vibrational mixing.Such vibrational mixing is wellestablished for amide I and the water bending modes. 58,59However, similar interactions between low-frequency modes of the peptide and collective water modes have not been fully investigated on a quantum chemical level.Our vibrational analysis reveals that nearly all modes below 700 cm À1 are involved in peptide-water mixing.As an example, Fig. 3 shows four mixed water-peptide modes of the GAG-water complex, two for pPII and two for b-strand.
The strong mixing between the water and peptide modes clearly indicates that peptide-water interactions constitute a dynamic entity with significant coherent vibrational dynamics.It is obvious that the degree of this mixing and, thus, the vibrational energy density should depend on the backbone conformation.To check this conjecture we calculated the difference between the (harmonic) entropies of pPII and b at room temperature by utilizing the well-known relationship: where N A is the Avogadro constant, h the Planck constant, c the vacuum velocity of light, ñi the wavenumber of the ith peptide-H 2 O vibration, R the gas constant, k B the Boltzmann constant, and T the absolute temperature.To avoid any arbitrary cutoffs, we included the entire vibrational manifold of N modes as obtained from the DFT based normal mode calculations.We are aware that this approach neglects the anharmonicity to be expected for low frequency modes and possible inaccuracies of our calculated frequencies.However, it is reasonable to expect that our calculations should provide us with reliable trends in the magnitude and sign of hydration-induced differences between the vibrational entropies of pPII and b-strand conformations.This is corroborated by recent comparisons between experimentally-obtained and computed entropy differences between pPII and b-strand of AcANH 2 in water. 32,33he vibrational entropy contributions to the Helmholtz energies of the investigated GxG peptides as obtained by eqn (3) are visualized in Fig. 4. Apparently, all vibrational entropies stabilize the b-strand conformation in explicit water by contributing between 5 to 18 kJ mol À1 at room temperature.Interestingly, the obtained side chain dependence of the vibrational entropies correlate with the corresponding DU values.Even though they display the same I 4 V 4 L hierarchy, they Fig. 3 Representative mixed water-peptide modes of the GAG-water complex.For both the pPII and b-strand conformations, there is a strong mixing between water and peptide modes.
underestimate the respective experimental values for these residues and overestimate the entropy of GAG.The underestimated entropy of GLG, GVG, and GIG can again be related to the inability of the considered 10 water molecules to fully hydrate these peptides' side chains.However, our results clearly suggest that the vibrational dynamics of the combined peptidewater system is a major source of the entropic contribution to the Gibbs or Helmholtz energy differences between pPII and the b-strand, in agreement with Lanza and Chiacchio. 33n order to further corroborate the contributions of hydration to the vibrational entropy differences between pPII and b-strand conformations, we calculated the entropy differences between the respective vibrational manifolds of the investigated GxGs without solvent molecules.The respective contributions to the total free energy are also plotted in Fig. 4. For GAG and GIG, these (positive) contributions are negligible.For GLG and GVG, the calculations actually yielded a room temperature stabilization of the pPII conformation by 4.3 and 1.6 kJ mol À1 , respectively.Together, all these calculations demonstrate that peptide solvation inverts the sign of the internal energy and entropy contribution to the Helmholtz energy governing the pPII " b equilibrium.
Finally, for the sake of completeness, we also estimated the vibrational contribution to the internal energy for the above considered vibrational manifold of GxGÁ10H 2 O.This can be done by using the textbook equation: We thus found that the contribution of the vibrational energy to the total internal energy does not vary significantly with x, and the corresponding values lie between À1.4 and À1.7 kJ mol À1 .This justifies our focus on the electronic energies derived from our structure optimization procedure.

Conclusion
Taken together with our earlier reported experimental results, the complementary DFT calculations reveal the dominant role of backbone hydration in stabilizing the pPII conformation of amino acid residues of short peptides in water.Specifically, these interactions stabilize the pPII conformation energetically but destabilizing it entropically.We were able to reproduce an earlier observed I 4 V 4 L hierarchy of the experimentallyobtained DH/DU and TDS values associated with the pPII " b equilibrium of GxG peptides in water on a qualitative level.We decomposed the obtained energies of the investigated pPII " b equilibria into intramolecular (peptide and water) and intermolecular contributions (water-peptide).The results of this procedure revealed that peptide-water interactions and solvent self-energies promote the energetic stabilization of pPII over water.Our results show that the hydration of the peptide backbone critically depends on the backbone structure and on the nature of the side chain.This notion seems to be consistent with the conditional solvation model of Ben Naim, but at variance with a model constructed from context free analogues of side chains. 55It is further consistent with predictions from electrostatic theories. 30Moreover, we provide evidence for the notion that the rather large entropic differences between pPII and b-strand are assignable to the vibrational dynamics of the peptide-water entity.Our calculations with only ten water molecules actually underestimate these contributions for residues with large and bulky side chains.Mezei et al. performed a molecular-dynamics-based analysis of a polyalanine 12-mer in explicit water and reported that the number of water molecules in the first hydration shell of both pPII and b is approximately 120. 24Any further characterization of peptide hydration shells in terms of hydrogen bonding and orientational distribution based on quantum chemical calculations will be the subject of future studies.

Fig. 1
Fig. 1 pPII and b-strand conformation of a cationic GAG + 10 H 2 O complex obtained from geometry optimizations as described in the text.
do not meet this expectation.These DU e (pPII,b) values can be directly compared to the experimental DH(pPII,b) (= H(pPII) À H(b)) values reported by Toal et al. 31 Since no pressure and significant volume changes are involved, one can expect that the enthalpy DH equals the internal energy DU int ; i.e. the Gibbs free energy equals the Helmholtz energy.Therefore, in the following, we interpret the experimentally-obtained DH values reported by Toal et al.

Fig. 2
Fig.2Comparison of experimental enthalpies (black bars)28 and calculated energies (light grey: explicit water; grey: in vacuo) for the pPII " b-strand equilibrium of the indicated amino acid residues in cationic GxG peptides in H 2 O.
DU e (P s b s w n ) = U e (P s w n ) À U e (b s w n ), DU e (P s b s ) = U e (P s ) À U e (b s ), DU e,Pb (w n ) = U e,P (w n ) À U e,b (w n ), and DU e (Pb) = U e (P) À U e (b).The energies DU e (P s b s w n ), DU e (P s b s ), and DU e,Pb (w n ) in eqn (

a
Electronic energy difference between pPII and b solvated systems consisting of the peptide and 10 water molecules.b Electronic energy difference between individual pPII and b structures in the geometry of the solvated system.c Electronic energy difference between 10 water molecules associated with the pPII and b structures.d Electronic interaction energy difference between pPII and b systems.This journal is © the Owner Societies 2015

Fig. 4
Fig. 4 Difference between the Gibbs/Helmholtz energy contributions of the vibrational entropy of the indicated amino acid residue in cationic GxG peptides in explicit water (grey) and in vacuo (light grey) in the pPII and b-strand conformation.Experimentally obtained TDS values of the corresponding the experimentally obtained pPII " b-strand equilibrium are plotted in black.In vacuo values for L and I are so small that they are barely visible on the chosen scale.