Energetic characterization of short helical polyalanine peptides in water: analysis of 13C=O chemical shift data

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

Capping parallel β-sheets of acetyl(Ala)6NH2 with an acetyl(Ala)5ProNH2 can arrest the growth of the sheet, suggesting a potential for curtailing amyloid growth. An ONIOM and density functional theory study

We present ONIOM calculations using B3LYP/d95(d,p) as the high level and AM1 as the medium level on parallel β-sheets containing four strands of Ac-AAAAAA-NH2 capped with either Ac-AAPAAA-NH2 or Ac-AAAPAA-NH2. Because Pro can form H-bonds from only one side of the peptide linkage (that containing the C═O H-bond acceptor), only one of the two Pro-containing strands can favorably add to the sheet on each side. Surprisingly, when the sheet is capped with AAPAAA-NH2 at one edge, the interaction between the cap and sheet is slightly more stabilizing than that of another all Ala strand. Breaking down the interaction enthalpies into H-bonding and distortion energies shows the favorable interaction to be due to lower distortion energies in both the strand and the four-stranded sheet. Because another strand would be inhibited for attachment to the other side of the capping (Pro-containing) strand, we suggest the possible use of Pro residues in peptides designed to arrest the growth of many amyloids.

Comparison of some dispersion-corrected and traditional functionals as applied to peptides and conformations of cyclohexane derivatives

We compare the energetic and structural properties of fully optimized α-helical and antiparallel β-sheet polyalanines and the energetic differences between axial and equatorial conformations of three cyclohexane derivatives (methyl, fluoro, and chloro) as calculated using several functionals designed to treat dispersion (B97-D, ωB97x-D, M06, M06L, and M06-2X) with other traditional functionals not specifically parametrized to treat dispersion (B3LYP, X3LYP, and PBE1PBE) and with experimental results. Those functionals developed to treat dispersion significantly overestimate interaction enthalpies of folding for the α-helix and predict unreasonable structures that contain Ramachandran φ and ψ and C = O...N H-bonding angles that are out of the bounds of databases compiled the β-sheets. These structures are consistent with overestimation of the interaction energies. For the cyclohexanes, these functionals overestimate the stabilities of the axial conformation, especially when used with smaller basis sets. Their performance improves when the basis set is improved from D95∗∗ to aug-cc-pVTZ (which would not be possible with systems as large as the peptides).

Distinctions in early stage unwinding mechanisms of zwitterionic, capped, and neutral forms of different α-helical homopolymeric peptides

Molecular dynamics simulations of α-helical polyalanine, polyleucine, polylysine, and poly(glutamic acid) with different forms of terminal groups in water at 300 K showed sharp distinctions in their unwinding mechanisms. Zwitterionic, capped, and neutral forms of polyalanine, polyleucine, and polylysine have been explored to elucidate their unwinding mechanism at very early stage, e.g., initial time window. Role of water in the unwinding mechanisms of the various helices has been envisaged. Also, it is evident from our calculations that the short- and long-range nonbonded interactions among the side chains is an important factor determining the unwinding mechanisms of the various homopolymeric α-helices. These findings can be helpful in constructing predictive models for understanding of the unwinding of α-helical proteins and peptides.

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.

Helix-coil energetics for helix formers and breakers reflect context and temperature: mutants of helically robust, guest-sensitive homopeptide hosts

The natural amino acids are primarily helix breakers at the low assignment temperatures characteristic of many studies, but recent genomic analyses of thermophilic proteins suggest that at high temperatures, some breakers may become strong helix formers. Moreover, the breaker/former inventory has not been previously characterized at the physiologically relevant temperature of 37 degrees C. The versatility of 13C==O NMR chemical shifts as helicity reporters allows construction of two mutant peptide series, tailored to expand the range of temperature assignments for helical propensities and derived from the core hosts tL-Ala9XxxAla9-tL and tL-AlaNva4XxxNva4Ala9-tL, Nva=norvaline. For three limiting guests Xxx, the helix former Nva and the breakers Gly and Pro, we report wXxx[T] assignments at seven temperatures from 2 to 80 degrees C, validating our reasoning and paving the way for assignment of a definitive wXxx[T] data-base.

Context-independent, temperature-dependent helical propensities for amino acid residues

Assigned from data sets measured in water at 2, 25, and 60 degrees C containing (13)C=O NMR chemical shifts and [theta](222) ellipticities, helical propensities are reported for the 20 genetically coded amino acids, as well as for norvaline and norleucine. These have been introduced by chemical synthesis at central sites within length-optimized, spaced, solubilized Ala(19) hosts. The resulting polyalanine-derived, quantitative propensity sets express for each residue its temperature-dependent but context-independent tendency to forego a coil state and join a preexisting helical conformation. At 2 degrees C their rank ordering is: P << G < H < C, T, N < S < Y, F, W < V, D < K < Q < I < R, M < L < E < A; at 60 degrees C the rank becomes: H, P < G < C < R, K < T, Y, F < N, V < S < Q < W, D < I, M < E < A < L. The DeltaDeltaG values, kcal/mol, relative to alanine, for the cluster T, N, S, Y, F, W, V, D, Q, imply that at 2 degrees C all are strong breakers: DeltaDeltaG(mean) = +0.63 +/- 0.11, but at 60 degrees C their breaking tendencies are dramatically attenuated and converge toward the mean: DeltaDeltaG(mean) = +0.25 +/- 0.07. Accurate modeling of helix-rich proteins found in thermophiles, mesophiles, and organisms that flourish near 0 degrees C thus requires appropriately matched propensity sets. Comparisons are offered between the temperature-dependent propensity assignments of this study and those previously assigned by the Scheraga group; the special problems that attend propensity assignments for charged residues are illustrated by lysine guest data; and comparisons of errors in helicity assignments from shifts and ellipticity data show that the former provide superior precision and accuracy.

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.

Temperature- and length-dependent energetics of formation for polyalanine helices in water: assignment of w(Ala)(n,T) and temperature-dependent CD ellipticity standards

Length-dependent helical propensities w(Ala)(n,T) at T = 10, 25, and 60 degrees C are assigned from t/c values and NMR 13C chemical shifts for series 1 peptides TrpLys(m)Inp2(t)Leu-Ala(n)(t)LeuInp2Lys(m)NH2, n = 15, 19, and 25, m = 5, in water. Van't Hoff analysis of w(Ala)(n,T) show that alpha-helix formation is primarily enthalpy-driven. For series 2 peptides Ac-Trp Lys5Inp2(t)Leu-(beta)AspHel-Ala(n)-beta-(t)LeuInp2Lys5NH2, n = 12 and 22, which contain exceptionally helical Ala(n) cores, protection factor-derived fractional helicities FH are assigned in the range 10-30 degrees C in water and used to calibrate temperature-dependent CD ellipticities [theta](lambda,H,n,T). These are applied to CD data for series 1 peptides, 12 < or = n < or = 45, to confirm the w(Ala)(n,T) assignments at T = 25 and 60 degrees C. The [theta](lambda,H,n,T) are temperature dependent within the wavelength region, 222 +/- 12 nm, and yield a temperature correction for calculation of FH from experimental values of [theta](222,n,T,Exp).

Finding the transition state without initial guess: The growing string method for Newton trajectory to isomerization and enantiomerization reaction of alanine dipeptide and poly(15)alanine

We report a new, high-dimensional application of a method for finding a transition state (TS) between a reactant and a product on the potential energy surface: the search of a growing string along a reaction path defined by any Newton trajectory in combination with the Berny method (Quapp, J Chem Phys (2005), 122, 174106; we have provided this algorithm on a web page). Two given minima are connected by a one-dimensional, but usually curvilinear reaction coordinate. It leads to the TS region. The application of the method to alanine dipeptide finds the TS of the isomerisation C(7 ax) --> C(5), some TSs of the enantiomerisation of C(7 ax) from L-form to quasi-D-form, and it finds the TS region of a transition of a partly unfolded, bent structure which turns back into a mainly alpha-helix in the Ac(Ala)(15)NHMe polyalanine (all at the quantum mechanical level B3LYP/6-31G: the growing string calculation is interfaced with the Gaussian03 package). The formation or dissolvation of some alpha- or 3(10)-hydrogen bonds of the helix are discussed along the TS pathway, as well as the case of an enantiomer at the central residue of the helix.

UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting

We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.