My 65 years in protein chemistry

This is a tour of a physical chemist through 65 years of protein chemistry from the time when emphasis was placed on the determination of the size and shape of the protein molecule as a colloidal particle, with an early breakthrough by James Sumner, followed by Linus Pauling and Fred Sanger, that a protein was a real molecule, albeit a macromolecule. It deals with the recognition of the nature and importance of hydrogen bonds and hydrophobic interactions in determining the structure, properties, and biological function of proteins until the present acquisition of an understanding of the structure, thermodynamics, and folding pathways from a linear array of amino acids to a biological entity. Along the way, with a combination of experiment and theoretical interpretation, a mechanism was elucidated for the thrombin-induced conversion of fibrinogen to a fibrin blood clot and for the oxidative-folding pathways of ribonuclease A. Before the atomic structure of a protein molecule was determined by x-ray diffraction or nuclear magnetic resonance spectroscopy, experimental studies of the fundamental interactions underlying protein structure led to several distance constraints which motivated the theoretical approach to determine protein structure, and culminated in the Empirical Conformational Energy Program for Peptides (ECEPP), an all-atom force field, with which the structures of fibrous collagen-like proteins and the 46-residue globular staphylococcal protein A were determined. To undertake the study of larger globular proteins, a physics-based coarse-grained UNited-RESidue (UNRES) force field was developed, and applied to the protein-folding problem in terms of structure, thermodynamics, dynamics, and folding pathways. Initially, single-chain and, ultimately, multiple-chain proteins were examined, and the methodology was extended to protein-protein interactions and to nucleic acids and to protein-nucleic acid interactions. The ultimate results led to an understanding of a variety of biological processes underlying natural and disease phenomena.

Effects of side chains in helix nucleation differ from helix propagation

Helix-coil transition theory connects observable properties of the α-helix to an ensemble of microstates and provides a foundation for analyzing secondary structure formation in proteins. Classical models account for cooperative helix formation in terms of an energetically demanding nucleation event (described by the σ constant) followed by a more facile propagation reaction, with corresponding s constants that are sequence dependent. Extensive studies of folding and unfolding in model peptides have led to the determination of the propagation constants for amino acids. However, the role of individual side chains in helix nucleation has not been separately accessible, so the σ constant is treated as independent of sequence. We describe here a synthetic model that allows the assessment of the role of individual amino acids in helix nucleation. Studies with this model lead to the surprising conclusion that widely accepted scales of helical propensity are not predictive of helix nucleation. Residues known to be helix stabilizers or breakers in propagation have only a tenuous relationship to residues that favor or disfavor helix nucleation.

Side-chain-side-chain interactions and stability of the helical state

Understanding the driving forces that lead to the stability of the secondary motifs found in proteins, namely α-helix and β-sheet, is a major goal in structural biology. The thermodynamic stability of these repetitive units is a result of a delicate balance between many factors, which in addition to the peptide chain involves also the solvent. Despite the fact that the backbones of all amino acids are the same (except of that of proline), there are large differences in the propensity of the different amino acids to promote the helical structure. In this paper, we investigate by explicit-solvent molecular dynamics simulations the role of the side chains (modeled as coarse-grained single sites) in stabilizing α helices in an aqueous solution. Our model systems include four (six-mer-nine-mer) peptide lengths in which the magnitude of the effective attraction between the side chains is systematically increased. We find that these interactions between the side chains can induce (for the nine-mer almost completely) a transition from a coil to a helical state. This transition is found to be characterized by three states in which the intermediate state is a partially folded α-helical conformation. In the absence of any interactions between the side chains the free energy change for helix formation has a small positive value indicating that favorable contributions from the side chains are necessary to stabilize the helical conformation. Thus, the helix-coil transition is controlled by the effective potentials between the side-chain residues and the magnitude of the required attraction per residue, which is on the order of the thermal energy, reduces with the length of the peptide. Surprisingly, the plots of the population of the helical state (or the change in the free energy for helix formation) as a function of the total effective interactions between the side chains in the helical state for all peptide lengths fall on the same curve.

Nucleation effects in peptide foldamers

Oligomers composed of β(3)-amino acid residues and a mixture of α- and β(3)-residues have emerged as proteolytically stable structural mimics of α-helices. An attractive feature of these oligomers is that they adopt defined conformations in short sequences. In this manuscript, we evaluate the impact of β(3)-residues as compared to their α-amino acid analogs in prenucleated helices. Our hydrogen-deuterium exchange results suggest that heterogeneous sequences composed of "αααβ" repeats are conformationally more rigid than the corresponding homogeneous α-peptide helices, with the macrocycle templating the helical conformation having a significant influence.

Empirical protein partition functions

In the present paper, we outline how to construct the partition function for a protein using empirical heat capacity data. The procedure is based on the calculation of a set of energy moments from the temperature dependence of the heat capacity. Given a set of energy moments, one can then use the maximum-entropy method to calculate an approximate energy distribution for the protein; the more energy moments one has, the better the approximation. The energy distribution can then be used to calculate the probability that the molecule is in a given energy level, which, using standard statistical mechanics, gives the degeneracy of the particular energy level. The degeneracy as a function of energy is the central ingredient in the construction of the partition function. Given the partition function, one can calculate all of the thermodynamic functions of the protein (free energy, energy, entropy, heat capacity, and energy probability distribution) as a function of temperature. The three-dimensional plot of the probability that the protein has a given energy at a given temperature tells one graphically (without imposing the assumption) whether or not it is a good approximation to divide the terms in the partition function into two or more groups, reflecting, for example, the presence of distinct native and denatured species.

Respice, adspice, and prospice

The title, "Look to the past, Look to the present, and Look to the future," the motto of City College of New York, expresses how my family life and education led me to an academic career in physical chemistry and ultimately to a study of proteins. The economic depression of the 1930s left a lasting impression on my outlook and career aspirations. With fortunate experiences at several stages in my life, I was able to participate in the great adventure of the last half of the twentieth century: the revolution in biology that advanced the field of protein chemistry to so great an extent. The future is bright and limitless, with greater understanding of biology yet to come.

Circular dichroism and ultraviolet resonance Raman indicate little Arg-Glu side chain α-helix peptide stabilization

Electrostatic interactions between side chains can control the conformation and folding of peptides and proteins. We used circular dichroism (CD) and ultraviolet (UV) resonance Raman spectroscopy (UVRR) to examine the impact of side chain charge on the conformations of two 21 residue mainly polyala peptides with a few Arg and Glu residues. We expected that attractions between Arg-10 and Glu-14 side chains would stabilize the α-helix conformation compared to a peptide with an Arg-14. Surprisingly, CD suggests that the peptide with the Glu-14 is less helical. In contrast, the UVRR show that these two peptides have similar α-helix content. We conclude that the peptide with Glu-14 has the same net α-helix content as the peptide with the Arg but has two α-helices of shorter length. Thus, side chain interactions between Arg-10 and Glu-14 have a minor impact on α-helix stability. The thermal melting of these two peptides is similar. However the Glu-14 peptide pH induced melting forms type III turn structures that form α-helix-turn-α-helix conformations.

Roles of non-native hydrogen-bonding interaction in helix-coil transition of a single polypeptide as revealed by comparison between Gō-like and non-Gō models

To explore the role of non-native interactions in the helix-coil transition, a detailed comparison between a Gō-like model and a non-Gō model has been performed via lattice Monte Carlo simulations. Only native hydrogen bonding interactions occur in the Gō-like model, and the non-native ones with sequence interval more than 4 is also included into the non-Gō model. Some significant differences between the results from those two models have been found. The non-native hydrogen bonds were found most populated at temperature around the helix-coil transition. The rearrangement of non-native hydrogen bonds into native ones in the formation of alpha-helix leads to the increase of susceptibility of chain conformation, and even two peaks of susceptibility of radius of gyration versus temperature exist in the case of non-Gō model for a non-short peptide, while just a single peak exists in the case of Gō model for a single polypeptide chain with various chain lengths. The non-native hydrogen bonds have complicated the temperature-dependence of Zimm-Bragg nucleation constant. The increase of relative probability of non-native hydrogen bonding for long polypeptide chains leads to non-monotonous chain length effect on the transition temperature.

SNARE complex zipping as a driving force in the dilation of proteinaceous fusion pores

The assembly of SNARE proteins into a tight complex has been hypothesized to drive membrane fusion. A model of the initial fusion pore as a proteinaceous channel formed by SNARE proteins places their membrane anchors in separate membranes. This leaves the possibility of a final assembly step that brings the membrane anchors together and drives fusion pore expansion. The present study develops a model for expansion in which the final SNARE complex zipping step drives a transition from a proteinaceous fusion pore to a lipidic fusion pore. An estimate of the energy released upon merger of the helical segments of the SNARE motifs with the helical segments of the membrane anchors indicates that completing the assembly of a few SNARE complexes can overcome the elastic energy that opposes lipid bilayer deformation into a narrow fusion pore. The angle between the helical axes of the membrane anchor and SNARE motif serves as a useful reaction coordinate for this transition. Energy was calculated as a function of this angle, incorporating contributions from membrane bending, SNARE complex assembly, membrane anchor flexing and hydrophobic interactions. The rate of this transition was evaluated as a process of diffusion over the barrier imposed by these combined energies, and the rates estimated were consistent with experimental measurements.

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.

Intrinsic amino acid side-chain hydrophilicity/hydrophobicity coefficients determined by reversed-phase high-performance liquid chromatography of model peptides: comparison with other hydrophilicity/hydrophobicity scales

An accurate determination of the intrinsic hydrophilicity/hydrophobicity of amino acid side-chains in peptides and proteins is fundamental in understanding many area of research, including protein folding and stability, peptide and protein function, protein-protein interactions and peptide/protein oligomerization, as well as the design of protocols for purification and characterization of peptides and proteins. Our definition of intrinsic hydrophilicity/hydrophobicity of side-chains is the maximum possible hydrophilicity/hydrophobicity of side-chains in the absence of any nearest-neighbor effects and/or any conformational effects of the polypeptide chain that prevent full expression of side-chain hydrophilicity/hydrophobicity. In this review, we have compared an experimentally derived intrinsic side-chain hydrophilicity/hydrophobicity scale generated from RP-HPLC retention behavior of de novo designed synthetic model peptides at pH 2 and pH 7 with other RP-HPLC-derived scales, as well as scales generated from classic experimental and calculation-based methods of octanol/water partitioning of Nalpha-acetyl-amino-acid amides or free energy of transfer of free amino acids. Generally poor correlation was found with previous RP-HPLC-derived scales, likely due to the random nature of the peptide mixtures in terms of varying peptide size, conformation and frequency of particular amino acids. In addition, generally poor correlation with the classical approaches served to underline the importance of the presence of a polypeptide backbone when generating intrinsic values. We have shown that the intrinsic scale determined here is in full agreement with the structural characteristics of amino acid side-chains.

Statistical mechanical model for helix-sheet-coil transitions in homopolypeptides

In this paper, we propose a simple statistical mechanical model to study the conformation transition between the alpha helix, beta sheet, and random coil in homopolypeptides. In our model, five parameters are introduced to obtain the partition function. There are two factors for helical propagation and initiation, which are the same as those used in the Zimm-Bragg model, and three newly introduced parameters for beta structures: the strand propagation factor for residues in beta strands and two correction factors for the initiation effect of the beta strand and beta sheet. Our model shows that the variation of these parameters may induce conformation transition from alpha helix or random coil to beta sheet. The sharpness of the transition depends on the initiation factors.

The helix-coil transition revisited

In this article, we perform a dynamic Monte Carlo simulation study of the helix-coil transition by using a bond-fluctuation lattice model. The results of the simulations are compared with those predicted by the Zimm-Bragg statistical thermodynamic theory with propagation and nucleation parameters determined from simulation data. The Zimm-Bragg theory provides a satisfactory description of the helix-coil transition of a homopolypeptide chain of 32 residues (N = 32). For such a medium-length chain, however, the analytical equation based on a widely-used large-N approximation to the Zimm-Bragg theory is not suitable to predict the average length of helical blocks at low temperatures when helicity is high. We propose an analytical large-eigenvalue (lambda) approximation. The new equation yields a significantly improved agreement on the average helix-block length with the original Zimm-Bragg theory for both medium and long chain lengths in the entire temperature range. Nevertheless, even the original Zimm-Bragg theory does not provide an accurate description of helix-coil transition for longer chains. We assume that the single-residue nucleation of helix formation as suggested in the original Zimm-Bragg model might be responsible for this deviation. A mechanism of nucleation by a short helical block is proposed by us and provides a significantly improved agreement with our simulation data.

Transformation of an alpha-helix peptide into a beta-hairpin induced by addition of a fragment results in creation of a coexisting state

Intrinsic rules of determining the tertiary structure of a protein have been unknown partly because physicochemical factors that contribute to stabilization of a protein structure cannot be represented as a linear combination of local interactions. To clarify the rules on the nonlinear term caused by nonlocal interaction in a protein, we tried to transform a peptide that has a fully helical structure (Target Peptide or TP) into a peptide that has a beta-hairpin structure (Designed Peptide or DP) by adding seven residues to the C terminus of TP. According to analyses of nuclear magnetic resonance measurements, while the beta-hairpin structure is stabilized in some DPs, it is evident that the helical structure observed in TP is also persistent and even extended throughout the length of the molecule. As a result, we have produced a peptide molecule that contains both the alpha-helix and beta-hairpin conformation at an almost equally populated level. The helical structures contained in these DPs were more stable than the helix in TP, suggesting that stabilizing one conformation does not result in destabilizing the other conformation. These DPs can thus be regarded as an isolated peptide version of the chameleon sequence, which has the capability of changing the secondary structure depending on the context of the surrounding environment in a protein structure. The fact that the transformation of one secondary structure caused stabilization of both the original and the induced structure would shed light on the mechanism of protein folding.

Equilibrium unfolding of the poly(glutamic acid)20 helix

The equilibrium structural ensemble of a 20-residue polyglutamic acid peptide (E(20)) was studied with FRET, circular dichroism, and molecular dynamics (MD) simulations. A FRET donor, o-aminobenzamide, and acceptor, 3-nitrotyrosine, were introduced at the N- and C-termini, respectively. Circular dichroism, steady state FRET, and time-resolved FRET measurements were employed to characterize the fraction helix and end-to-end distance under different pH conditions: pH 4 (60% alpha-helix), pH 6 (0% alpha-helix), and pH 9 (0% alpha-helix). At pH 4, the end-to-end distance was measured at 24 A and determined to be considerably less than the 31 A predicted for an alpha-helix of the same length. At pH 6 and 9, the end-to-end distance was measured at > 31 and 39 A respectively, both which are determined to be considerably greater than the 27 A predicted for a freely jointed random coil of the same length. To better understand the physical forces underlying the unusual helix-coil transition in this peptide, three theoretical MD models of E(20) were constructed: (1) a pure alpha-helix, (2) an alpha-helix with equivalent attractive intramolecular contacts, and (3) a weak alpha-helix with termini-weighted intramolecular contacts ("sticky ends"). Using MD simulations, the bent helix structure calculated from Model 3 was found to be the closest in agreement with the experimental data.

Alpha-helix stabilization by alanine relative to glycine: roles of polar and apolar solvent exposures and of backbone entropy

The energetics of alpha-helix formation are fairly well understood and the helix content of a given amino acid sequence can be calculated with reasonable accuracy from helix-coil transition theories that assign to the different residues specific effects on helix stability. In internal helical positions, alanine is regarded as the most stabilizing residue, whereas glycine, after proline, is the more destabilizing. The difference in stabilization afforded by alanine and glycine has been explained by invoking various physical reasons, including the hydrophobic effect and the entropy of folding. Herein, the contribution of these two effects and that of hydrophilic area burial is evaluated by analyzing Ala and Gly mutants implemented in three helices of apoflavodoxin. These data, combined with available data for similar mutations in other proteins (22 Ala/Gly mutations in alpha-helices have been considered), allow estimation of the difference in backbone entropy between alanine and glycine and evaluation of its contribution and that of apolar and polar area burial to the helical stabilization typically associated to Gly-->Ala substitutions. Alanine consistently stabilizes the helical conformation relative to glycine because it buries more apolar area upon folding and because its backbone entropy is lower. However, the relative contribution of polar area burial (which is shown to be destabilizing) and of backbone entropy critically depends on the approximation used to model the structure of the denatured state. In this respect, the excised-peptide model of the unfolded state, proposed by Creamer and coworkers (1995), predicts a major contribution of polar area burial, which is in good agreement with recent quantitations of the relative enthalpic contribution of Ala and Gly residues to alpha-helix formation.

Kinetic Studies of Folding of the B-domain of Staphylococcal Protein A with Molecular Dynamics and a United-residue (UNRES) Model of Polypeptide Chains

Langevin dynamics is used with our physics-based united-residue (UNRES) force field to study the folding pathways of the B-domain of staphylococcal protein A (1BDD (alpha; 46 residues)). With 400 trajectories of protein A started from the extended state (to gather meaningful statistics), and simulated for more than 35 ns each, 380 of them folded to the native structure. The simulations were carried out at the optimal folding temperature of protein A with this force field. To the best of our knowledge, this is the first simulation study of protein-folding kinetics with a physics-based force field in which reliable statistics can be gathered. In all the simulations, the C-terminal alpha-helix forms first. The ensemble of the native basin has an average RMSD value of 4 A from the native structure. There is a stable intermediate along the folding pathway, in which the N-terminal alpha-helix is unfolded; this intermediate appears on the way to the native structure in less than one-fourth of the folding pathways, while the remaining ones proceed directly to the native state. Non-native structures persist until the end of the simulations, but the native-like structures dominate. To express the kinetics of protein A folding quantitatively, two observables were used: (i) the average alpha-helix content (averaged over all trajectories within a given time window); and (ii) the fraction of conformations (averaged over all trajectories within a given time window) with Calpha RMSD values from the native structure less than 5 A (fraction of completely folded structures). The alpha-helix content grows quickly with time, and its variation fits well to a single-exponential term, suggesting fast two-state kinetics. On the other hand, the fraction of folded structures changes more slowly with time and fits to a sum of two exponentials, in agreement with the appearance of the intermediate, found when analyzing the folding pathways. This observation demonstrates that different qualitative and quantitative conclusions about folding kinetics can be drawn depending on which observable is monitored.

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

Measured at 2 degrees C in water, NMR chemical shifts of (13)C=O labeled central alanine residues of peptides W-Lys(5)-(t)L(3)-Ala(n)-(t)L(3)-Lys(5)NH(2), n = 9, 11, 13, 15, 19 and W-Lys(5)-(t)L(3)-a-Ala(n)-A-Inp-(t)L(2)-Lys(5)NH(2) (a = D-Ala; (t)L = tert-leucine; Inp = 4-carboxypiperidine) are used to assign jt(L) and ct(L), the N- and C-terminal (t)L capping parameters and length-dependent values for w(Ala)(n), the alanine helical propensity for Ala(n) peptides. These parameters allow Lifson-Roig characterization of the stabilities of Ala(n)() helices in water. To facilitate chemical shift characterization, different (13)C/(12)C ratios are incorporated into specific Ala sites to code up to six residue sites per peptide. Large left/right chemical shift anisotropies are intrinsic to helical polyalanines, and a correcting L-R-based model is introduced. Capping parameters jt(L) = ct(L) lie in the range of 0.3 to 0.5; the (t)L residues are thus moderately helix-destabilizing. For helical conformations of lengths shorter than eight residues, assigned values for w(Ala) approach 1.0 but increase monotonically with length to a value of 1.59 for w(Ala)(19).

Electronic properties of amino acid side chains: quantum mechanics calculation of substituent effects

BACKGROUND

Electronic properties of amino acid side chains such as inductive and field effects have not been characterized in any detail. Quantum mechanics (QM) calculations and fundamental equations that account for substituent effects may provide insight into these important properties. PM3 analysis of electron distribution and polarizability was used to derive quantitative scales that describe steric factors, inductive effects, resonance effects, and field effects of amino acid side chains.

RESULTS

These studies revealed that: (1) different semiempirical QM methods yield similar results for the electronic effects of side chain groups, (2) polarizability, which reflects molecular deformability, represents steric factors in electronic terms, and (3) inductive effects contribute to the propensity of an amino acid for alpha-helices.

CONCLUSION

The data provide initial characterization of the substituent effects of amino acid side chains and suggest that these properties affect electron density along the peptide backbone.

Position dependence of the 13C chemical shifts of alpha-helical model peptides. Fingerprint of the 20 naturally occurring amino acids

The position dependence of the (13)C chemical shifts was investigated at the density functional level for alpha-helical model peptides represented by the sequence Ac-(Ala)(i)-X-(Ala)(j)-NH(2), where X represents any of the 20 naturally occurring amino acids, with 0 < or = i < or = 8 and i + j = 8. Adoption of the locally dense basis approach for the quantum chemical calculations enabled us to reduce the length of the chemical-shift calculations while maintaining good accuracy of the results. For the 20 naturally occurring amino acids in alpha-helices, there is (1) significant variability of the computed (13)C shielding as a function of both the guest residue (X) and the position along the sequence; for example, at the N terminus, the (13)C(alpha) and (13)C(beta) shieldings exhibit a uniform pattern of variation with respect to both the central or the C-terminal positions; (2) good agreement between computed and observed (13)C(alpha) and (13)C(beta) chemical shifts in the interior of the helix, with correlation coefficients of 0.98 and 0.99, respectively; for (13)C(alpha) chemical shifts, computed in the middle of the helix, only five residues, namely Asn, Asp, Ser, Thr, and Leu, exhibit chemical shifts beyond the observed standard deviation; and (3) better agreement for four of these residues (Asn, Asp, Ser, and Thr) only for the computed values of the (13)C(alpha) chemical shifts at the N terminus. The results indicate that (13)C(beta), but not (13)C(beta), chemical shifts are sensitive enough to reflect the propensities of some amino acids for specific positions within an alpha-helix, relative to the N and C termini of peptides and proteins.

Pairwise Coupling in an Arg-Phe-Met Triplet Stabilizes α-Helical Peptide via Shared Rotamer Preferences

The hydrophobic Arg-Phe and Phe-Met side chain interactions stabilize the alpha-helix by -0.29 and -0.59 kcal/mol, respectively, when placed i, i + 4 in an alanine-based peptide. When both interactions are present simultaneously, however, they stabilize the helix by an additional -0.75 kcal/mol, nearly as much as the sum of its parts. We attribute this coupling to a shared rotamer preference, as the central Phe is t in both bonds. The energetic cost of restricting the Phe residue into a t conformation is only paid once in the triplet, rather than twice when the interactions are separate. Coupling is thus demonstrated to have large effects on protein stability.

Metal clips induce folding of a short unstructured peptide into an alpha-helix via turn conformations in water. Kinetic versus thermodynamic products

Short peptides corresponding to two to four alpha-helical turns of proteins are not thermodynamically stable helices in water. Unstructured octapeptide Ac-His1-Ala2-Ala3-His4-His5-Glu6-Leu7-His8-NH(2) (1) reacts with two [Pd((15)NH(2)(CH(2))(2)(15)NH(2))(NO(3))(2)] in water to form a kinetically stable intermediate, [[Pden](2)[(1,4)(5,8)-peptide]](2), in which two 19-membered metallocyclic rings stabilize two peptide turns. Slow subsequent folding to a thermodynamically more stable two-turn alpha-helix drives the equilibrium to [[Pden](2)[(1,5)(4,8)-peptide]] (3), featuring two 22-membered rings. This transformation from unstructured peptide via turns to an alpha-helix suggests that metal clips might be useful probes for investigating peptide folding.

Thermodynamics Of α‐Helix Formation

The alpha-helix was the first proposed and experimentally confirmed secondary structure. The elegant simplicity of the alpha-helical structure, stabilized by hydrogen bonding between the backbone carbonyl oxygen and the peptide amide four residues away, has captivated the scientific community. In proteins, alpha-helices are also stabilized by the so-called capping interactions that occur at both the C- and the N-termini of the helix. This chapter provides a brief historical overview of the thermodynamic studies of the energetics of helix formation, and reviews recent progress in our understanding of the thermodynamics of helix formation.

The helical alanine controversy: an (Ala)6 insertion dramatically increases helicity

Employing chemical shift melts and hydrogen/deuterium exchange NMR techniques, we have determined the stabilization of the Trp-cage miniprotein due to multiple alanine insertions within the N-terminal alpha-helix. Alanine is shown to be uniquely helix-stabilizing and this stabilization is reflected in the global fold stability of the Trp-cage. The associated free energy change per alanine can be utilized to calculate the alanine propagation value. From the Lifson-Roig formulation, the calculated value (wAla = 1.6) is comparable to those obtained for short, solubilized, alanine-rich helices and is much larger than the values obtained by prior host-guest techniques or in N-terminally templated helices and peptides bearing long contiguous strings of alanines with no capping or solubilizing units present.

A structural model of polyglutamine determined from a host-guest method combining experiments and landscape theory

Modeling the structure of natively disordered peptides has proved difficult due to the lack of structural information on these peptides. In this work, we use a novel application of the host-guest method, combining folding theory with experiments, to model the structure of natively disordered polyglutamine peptides. Initially, a minimalist molecular model (C(alpha)C(beta)) of CI2 is developed with a structurally based potential and captures many of the folding properties of CI2 determined from experiments. Next, polyglutamine "guest" inserts of increasing length are introduced into the CI2 "host" model and the polyglutamine is modeled to match the resultant change in CI2 thermodynamic stability between simulations and experiments. The polyglutamine model that best mimics the experimental changes in CI2 thermodynamic stability has 1), a beta-strand dihedral preference and 2), an attractive energy between polyglutamine atoms 0.75-times the attractive energy between the CI2 host Go-contacts. When free-energy differences in the CI2 host-guest system are correctly modeled at varying lengths of polyglutamine guest inserts, the kinetic folding rates and structural perturbation of these CI2 insert mutants are also correctly captured in simulations without any additional parameter adjustment. In agreement with experiments, the residues showing structural perturbation are located in the immediate vicinity of the loop insert. The simulated polyglutamine loop insert predominantly adopts extended random coil conformations, a structural model consistent with low resolution experimental methods. The agreement between simulation and experimental CI2 folding rates, CI2 structural perturbation, and polyglutamine insert structure show that this host-guest method can select a physically realistic model for inserted polyglutamine. If other amyloid peptides can be inserted into stable protein hosts and the stabilities of these host-guest mutants determined, this novel host-guest method may prove useful to determine structural preferences of these intractable but biologically relevant protein fragments.

Water Molecule Adsorption on Short Alanine Peptides: How Short Is the Shortest Gas-Phase Alanine-Based Helix?

Water adsorption measurements have been performed under equilibrium conditions for unsolvated Ac-A(n)K+H(+) and Ac-KA(n)+H(+) peptides with n = 4 - 10. Previous work on larger alanine peptides has shown that two dominant conformations (helices and globules) are present for these peptides and that water adsorbs much more strongly to the globules than to the helices. All the Ac-KA(n)+H(+) peptides studied here (which are expected to be globular) adsorb water strongly, and so do the Ac-A(n)K+H(+) peptides with n < 8. However, for Ac-A(n)K+H(+) with n = 8-10 there is a substantial drop in the propensity to adsorb water. This result suggests that Ac-A(8)K+H(+) is the smallest Ac-A(n)K+H(+) peptide to have a significant helical content in the gas phase. Water adsorption measurements for Ac-V(n)K+H(+) and Ac-L(n)K+H(+) with n = 5-10 suggest that the helix emerges at n = 8 for these peptides as well.

α-Turn Mimetics: Short Peptide α-Helices Composed of Cyclic Metallopentapeptide Modules

Alpha-Helices are key structural components of proteins and important recognition motifs in biology. Short peptides (<or=15 residues) corresponding to these helical sequences are rarely helical away from their stabilizing protein environments. New techniques for stabilizing short peptide helices could be valuable for studying protein folding, modeling proteins, creating artificial proteins, and may aid the design of inhibitors or mimics of protein function. This study reports the facile incorporation of 3- and 4-alpha turns in 10-15 residue peptides through formation in situ of multiple cyclic metallopeptide modules [Pd(en)(H*XXXH*)](2+). The nonhelical peptides Ac-H*ELTH*H*VTDH*-NH(2) (1), Ac-H*ELTH*AVTDYH*ELTH*-NH(2) (2), and Ac-H*AAAH*HELTH*H*VTDH*-NH(2) (3) (H is histidine-methylated at imidazole-N3) react in N,N-dimethylformamide (DMF) or water with 2, 2, and 3 molar equivalents, respectively, of [Pd(en)(NO(3))(2)] to form exclusively [Pd(2)(en)(2)(Ac-H*ELTH*H*VTDH*-NH(2))](4+) (4), [Pd(2)(en)(2)(Ac-H*ELTH*AVTDYH*ELTH*-NH(2))](4+) (5), and [Pd(3)(en)(3)(Ac-H*AAAH*HELTH*H*VTDH*-NH(2))](6+) (6), characterized by mass spectrometry, 1D and 2D (1)H- and 1D (15)N-NMR spectroscopy. Despite the presence of multiple histidines and other possible metal-binding residues in these peptides, 2D (1)H NMR spectra reveal that Pd(en)(2+) is remarkably specific in coordinating to imidazole-N1 of only (i, i + 4) pairs of histidines (i.e., only those separated by three amino acids), resulting in 4-6 made up of cyclic metallopentapeptide modules ([Pd(en)(H*XXXH*)](2+))(n), n = 2, 2, 3, respectively, each cycle being a 22-membered ring. We have previously shown that a single metallopentapeptide can nucleate alpha-helicity (Kelso et al., Angew. Chem., Int. Ed. 2003, 42, 421-424.). We now demonstrate its use as an alpha-turn-mimicking module for the facile conversion of unstructured short peptides into helices of macrocycles and provide 1D and 2D NMR spectroscopic data, structure calculations via XPLOR and NMR analysis of molecular flexibility in solution (NAMFIS), and CD spectra in support of the alpha-helical nature of these monomeric metallopeptides in solution.

A coarse-grained model and associated lattice Monte Carlo simulation of the coil–helix transition of a homopolypeptide

A new coarse-grained lattice model neglecting atomic details is proposed for the coil-helix transition and a new physical parameter is suggested to characterize a helical structure. In our model, each residue is represented by eight lattice sites, and side groups are not considered explicitly. Chirality and hydrogen bonding are taken into consideration in addition to chain connectivity and the excluded volume effect. Through a dynamic Monte Carlo simulation, the physical properties of the coil-helix transition of a single homopolypeptide have been produced successfully within a short computing time on the PC. We also examined the effects of the variation of chain configurations including chain size and chain shape, etc. A spatial correlation function has been introduced in order to characterize periodicity of a helical chain in a simple way. A propagation parameter and a nucleation parameter have also been calculated, which compares favorably with the results of the Zimm-Bragg theory for the coil-helix transition.

Probing helix formation in unsolvated peptides

Ion mobility measurements have been used to examine helix formation in unsolvated glycine-based peptides containing three alanine residues. Nine sequence isomers of Ac-[12G3A]K+H(+) were studied (Ac = acetyl, G = glycine, A = alanine, and K = lysine). The amount of helix present for each peptide was examined using two metrics, and it is strongly dependent on the proximity and the location of the alanine residues. Peptides with three adjacent alanines have the highest helix abundances, and those with well-separated alanines have the lowest. The helix abundances for most of the peptides can be fit reasonably well using a modified Lifson-Roig theory. However, Lifson-Roig theory fails to account for several key features of the experimental results. The most likely explanation for the correlation between helix abundances and the number of adjacent alanines is that neighboring alanines promote helix nucleation.

Electrostatic interactions in a peptide--RNA complex

Parallel experimental measurements and theoretical calculations have been used to investigate the energetics of electrostatic interactions in the complex formed between a 22 residue, alpha-helical peptide from the N protein of phage lambda and its cognate 19 nucleotide box B RNA hairpin. Salt-dependent free energies were measured for both peptide folding from coil to helix and peptide binding to RNA, and from these the salt-dependence of binding pre-folded, helical peptide to RNA was determined ( partial differential (DeltaG degrees (dock))/ partial differential log[KCl]=5.98(+/-0.21)kcal/mol). (A folding transition taking place in the RNA hairpin loop was shown to have a negligible dependence on salt concentration.) The non-linear Poisson-Boltzmann equation was used to calculate the same salt dependence of the binding free energy as 5.87(+/-0.22)kcal/mol, in excellent agreement with the measured value. Close agreement between experimental measurements and calculations was also obtained for two variant peptides in which either a basic or acidic residue was replaced with an uncharged residue, and for an RNA variant with a deletion of a single loop nucleotide. The calculations suggest that the strength of electrostatic interactions between a peptide residue and RNA varies considerably with environment, but that all 12 positive and negative N peptide charges contribute significantly to the electrostatic free energy of RNA binding, even at distances up to 11A from backbone phosphate groups. Calculations also show that the net release of ions that accompanies complex formation originates from rearrangements of both peptide and RNA ion atmospheres, and includes accumulation of ions in some regions of the complex as well as displacement of cations and anions from the ion atmospheres of the RNA and peptide, respectively.

Noncharged amino acid residues at the solvent-exposed positions in the middle and at the C terminus of the alpha-helix have the same helical propensity

It was established previously that helical propensities of different amino acid residues in the middle of alpha-helix in peptides and in proteins are very similar. The statistical analysis of the protein helices from the known three-dimensional structures shows no difference in the frequency of noncharged residues in the middle and at the C terminus. Yet, experimental studies show distinctive differences for the helical propensities of noncharged residues in the middle and in the C terminus in model peptides. Is this a general effect, and is it applicable to protein helices or is it specific to the model alanine-based peptides? To answer this question, the effects of substitutions at positions 28 (middle residue) and 32 (C2 position at the C terminus) of the alpha-helix of ubiquitin on the stability of this protein are measured by using differential scanning calorimetry. The two data sets produce similar values for intrinsic helix propensity, leading to a conclusion that noncharged amino acid residues at the solvent-exposed positions in the middle and at the C terminus of the alpha-helix have the same helical propensity. This conclusion is further supported with an excellent correlation between the helix propensity scale obtained for the two positions in ubiquitin with the experimental helix propensity scale established previously and with the statistical distribution of the residues in protein helices.

Recent advances in helix-coil theory

Peptide helices in solution form a complex mixture of all helix, all coil or, most frequently, central helices with frayed coil ends. In order to interpret experiments on helical peptides and make theoretical predictions on helices, it is therefore essential to use a helix-coil theory that takes account of this equilibrium. The original Zimm-Bragg and Lifson-Roig helix-coil theories have been greatly extended in the last 10 years to include additional interactions. These include preferences for the N-cap, N1, N2, N3 and C-cap positions, capping motifs, helix dipoles, side chain interactions and 3(10)-helix formation. These have been applied to determine energies for these preferences from experimental data and to predict the helix contents of peptides. This review discusses these newly recognised structural features of helices and how they have been included in helix-coil models.

Lifson-Roig nucleation for alpha-helices in trifluoroethanol: context has a strong effect on the helical propensity of amino acids

We have investigated the effect of substituting each of 19 common amino acids (excluding P) at the X position in the peptide acetyl-Y-VAEAK-TSXSR-VAEAK-NH(2). This very different peptide is of interest because, in previous work, we showed that nucleation in the strong alpha-helix-forming pentamers VAEAK was unable to propagate the helix structure through the sequence TSDSR, which is neither a helix-forming sequence nor a breaker, but is indifferent to helix formation. Substitution in the center of the indifferent sequence reveals an interesting measure of the helical propensity for the 19 amino acids. CD spectra were measured in various mixtures of buffer and 2,2,2-trifluoroethanol (TFE), and then analyzed for helix propensity of the amino acids using the Lifson-Roig model. However, the nucleation parameter in the Lifson-Roig model has never been measured for TFE. We have empirically found that the nucleation parameter for a solvent can be determined from the data normally used to determine only the propagation parameters. The results of the analysis of the CD show that most amino acids are excellent or good helix formers in 90% TFE, while amino acids D, W, F and G are poor helix formers for the indifferent pentamer sequence. The helix propensity of the 19 amino acids is quite different from the helix propensity measured in other peptide sequences, demonstrating the context dependence of this property. The results as a function of alcohol concentration confirm that the relative order of helical propensity of amino acids changes with solvent environment. Clearly, the prediction of alpha-helical secondary structure from protein sequence requires more than a single helical propensity for each amino acid.

Helix-coil transitions re-visited

The thermally-induced helix-coil transition in polyamino acids is a good model for determining the helix-forming propensities of amino acids but not for the two-state folding/unfolding transition in globular proteins. The equilibrium and kinetic treatments of the helix-coil transition are summarized here together with a description of applications to various types of homopolymers and copolymers. Attention is then focused on the helix-coil transition in poly-L-alanine as an example of a non-polar polyamino acid. To render such a non-polar polymer water soluble, it is necessary to introduce polar amino acids such as lysines, but care must be taken as to the location of such polar residues. If they are attached as end groups, as in a triblock copolymer, they do not perturb the helix-forming tendency of the central poly-L-alanine block significantly, but if they are introduced within the sequence of alanine residues, then the hydration properties of the lysines dominate the behavior of the resulting copolymer, thereby leading to erroneous values of the parameters characterizing the helix-forming tendency of the alanines. Neutral but polar residues, such as glutamines, also exhibit hydration-dominating properties but less so than charged lysines. Some details of the calculations for an alanine/glutamine copolymer are presented here. It is concluded that random copolymers based on a neutral water-soluble host provide reliable information about the helix-forming tendencies of amino acid residues that are introduced as guests among such neutral host residues.

Direct observation of the multistep helix formation of poly-L-glutamic acids

The helix formation dynamics of poly-L-glutamic acids (PGAs) were observed by the microsecond-resolved Fourier transform infrared (FTIR) and circular dichroism (CD) spectroscopies. The helix formation of 34-residue PGA from random coil at pH (or pD for FTIR) 8.0 was initiated by a pH jump to 4.9 using the rapid solution mixer whose mixing dead time is 50 micros. The amide I' line in the time-resolved FTIR spectra exhibited the fast (<100 micros) increase of the total helical content. The time-resolved CD spectra of the same process also showed the fast (<150 micros) formation of short helical segments (5 +/- 1 residues), which was followed by the slower (<1 ms) elongation of the short helices to longer helices (>10 residues). Similar dynamics were observed for the same pH jump of approximately 190-residue PGA, although there were additional steps that made the helix formation of approximately 190-residue PGA more complex. The observed multistep helix formation is likely caused by the strong hydrogen-bonding interactions between the protonated side chains of PGAs.

Solvation model dependency of helix-coil transition in polyalanine

Helix-coil transitions in polyalanine molecules of length 10 are studied by multi-canonical Monte Carlo simulations. The solvation effects are included by either a distance-dependent dielectric permittivity or by a term that is proportional to the solvent-accessible surface area of the peptide. We found a strong dependence of the characteristics of the helix-coil transition from the details of the solvation model.

The effect of C-terminal helix stabilization on specific DNA binding by monomeric GCN4 peptides

DNA binding by a 29-residue, monomeric, GCN4 basic region peptide, GCN4br, as well as by peptide br-C, a monomeric basic-region analogue that is helix stabilized at its C-terminal end by a Lys25. Asp29 side-chain lactam-bridged alanine-rich sequence, was studied at 25 C in an aqueous buffer containing 100 mm NaCl. Mixing of both peptides with duplex DNA containing the cAMP-responsive element (CRE) was accompanied by significant helix stabilization in the peptides, whereas mixing of the peptides with duplex DNA containing a scrambled CRE site was not. Peptide NBD-br-C was synthesized as a fluorescent probe to evaluate these peptide-DNA interactions further. Quantitative analysis of the fluorescence quenching of peptide NBD-br-C by CRE half-site DNA indicated the formation of a 1:1 complex with a dissociation constant of 1.41 +/- 0.22 microm. Competitive displacement fluorescence assays of CRE half-site binding gave dissociation constants of 0.65 +/- 0.09 microm for peptide br-C and 3.9 +/- 0.5 microM for GCN4br, which corresponds to a free energy difference of 1.1 kcal/mol that is attributed to the helix stabilization achieved in peptide br-C. This result indicates that helix initiation by the alpha-helical leucine zipper dimerization motif in native bzip proteins, such as GCN4, contributes significantly to the affinity of basic region peptides for their recognition sites on DNA. Our fluorescence assay should also prove useful for determining dissociation constants for CRE binding by other GCN4 basic region analogues under equilibrium conditions and physiological salt concentrations.

Consistent helicities from CD and template t/c data for N-templated polyalanines: progress toward resolution of the alanine helicity problem

The helicity reporting parameters t/c and [theta](222) have been measured at 2, 25, and 60 degrees C in water for the solubilized polyalanine series Ac-Hel-A(n)-(t)LInp(2)K(4)W-NH(2) of length 4 < or = n < or = 14 that bears the helix-initiating and monitoring N-cap Ac-Hel and the spaced solubilizer (t)LInp(2)K(4)W-NH(2) as a C-cap. Correlation of t/c with length shows that the helical propensity for n < or = 6 is ca. 1.0, consistent with our early reports, but that a dramatic increase in temperature dependence and helical propensity occurs for n > or = 8. A model based on hydrogen-bonding cooperativity is proposed to explain this finding, and both t/c and [theta](222) are modeled successfully by length-dependent alanine propensities at 2 degrees C of 1.03 for n = 6, 1.15, for 7 < or = n < or = 9 and 1.26 for n > or = 10. The implications of these results for the energetics of helix formation by alanine-rich peptide sequences are discussed.

Contribution of secondary structure to the heat capacity and enthalpy distribution of the unfolded state in proteins

We have recently shown that one can construct the enthalpy distribution for protein molecules from experimental knowledge of the temperature dependence of the heat capacity. For many proteins the enthalpy distribution evaluated at the midpoint of the denaturation transition (corresponding to the maximum in the heat capacity vs temperature curve) is broad and biphasic, indicating two different populations of molecules (native and unfolded) with distinctly different enthalpies. At temperatures above the denaturation point, the heat capacity for the unfolded state in many proteins is quite large and using the analysis just mentioned, we obtain a gaussian-like enthalpy distribution that is very broad. A large value of the heat capacity indicates that there are structural changes going on in the unfolded state above the transition temperature. In the present paper we investigate the origin of this large heat capacity by considering the presence of changing amounts of secondary structure (specifically, alpha-helix) in the unfolded state. For this purpose we use the empirical estimates of the Zimm-Bragg sigma and s factors for all of the native amino acids in water as determined by Scheraga and co-workers. Using myoglobin as an example, we calculate probability profiles and distribution functions for the total number of helix states in the specific-sequence molecule. Given the partition function for the specific-sequence molecule, we can then calculate a set of enthalpy moments for the molecule from which we obtain a good estimate of the enthalpy distribution in the unfolded state. This distribution turns out to be quite narrow when compared with the distribution obtained from the raw heat capacity data. We conclude that there must be other major structural changes (backbone and solvent) that are not accounted for by the inclusion of alpha-helix in the unfolded state.

Free energy distributions in proteins

Protein molecules in solution have a broad distribution of enthalpy states. A good approximation to the distribution function for enthalpy states can be calculated, using the maximum-entropy method, from the moments of the distribution that, in turn, are obtained from the experimental temperature dependence of the heat capacity. In the present paper, we show that the enthalpy probability distribution can then be formulated in terms of a free energy function that gives the free energy of the protein corresponding to a particular value of the enthalpy. By the location of the minima in this function, the free energy distribution graphically indicates the most probable values of the enthalpy for the protein. We find that the behavior of the free energy functions for proteins falls somewhere between two different cases: a two-state like function with two minima, the relative levels of the two states changing with temperature; and, a single-minimum function where the position of the minimum shifts to higher enthalpy values as the temperature is increased. We show that the temperature dependence of the free energy function can be expressed in terms of a central free energy distribution for a given, fixed temperature (which is most conveniently chosen as the temperature of the maximum in the heat capacity). The nature of this central free energy function for a given protein thus yields all of the thermodynamic behavior of that protein over the temperature range of the denaturation process.

Hydrogen-bonding classes in proteins and their contribution to the unfolding reaction

This paper proposes to assess hydrogen-bonding contributions to the protein stability, using a set of model proteins for which both X-ray structures and calorimetric unfolding data are known. Pertinent thermodynamic quantities are first estimated according to a recent model of protein energetics based on the dissolution of alkyl amides. Then it is shown that the overall free energy of hydrogen-bond formation accounts for a hydrogen-bonding propensity close to helix-forming tendencies previously found for individual amino acids. This allows us to simulate the melting curve of an alanine-rich helical 50-mer with good precision. Thereafter, hydrogen-bonding enthalpies and entropies are expressed as linear combinations of backbone-backbone, backbone-side-chain, side-chain-backbone, and side-chain-side-chain donor-acceptor contributions. On this basis, each of the four components shows a different free energy versus temperature trend. It appears that structural preference for side-chain-side-chain hydrogen bonding plays a major role in stabilizing proteins at elevated temperatures.

alpha-helix formation: discontinuous molecular dynamics on an intermediate-resolution protein model

An intermediate-resolution model of small, homogeneous peptides is introduced, and discontinuous molecular dynamics simulation is applied to study secondary structure formation. Physically, each model residue consists of a detailed three-bead backbone and a simplified single-bead side-chain. Excluded volume and hydrogen bond interactions are constructed with discontinuous (i.e., hard-sphere and square-well) potentials. Simulation results show that the backbone motion of the model is limited to realistic regions of Phi-Psi conformational space. Model polyalanine chains undergo a locally cooperative transition to form alpha-helices that are stabilized by backbone hydrogen bonding, while model polyglycine chains tend to adopt nonhelical structures. When side-chain size is increased beyond a critical diameter, steric interactions prevent formation of long alpha-helices. These trends in helicity as a function of residue type have been well documented by experimental, theoretical, and simulation studies and demonstrate the ability of the intermediate-resolution model developed in this work to accurately mimic realistic peptide behavior. The efficient algorithm used permits observation of the complete helix-coil transition within 15 min on a single-processor workstation, suggesting that simulations of very long times are possible with this model.

Electronic properties of the amino acid side chains contribute to the structural preferences in protein folding

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: alpha-helix, beta-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (esigma), and a composite of these two scales (epsilon). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for alpha-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.