Prediction of the thermodynamics of protein unfolding: the helix-coil transition of poly(L-alanine)

Water-mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water-mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins-which is a well-accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free-energy calculation for a thermodynamic cycle of a protein folding using a liquid-state density functional theory. Here, the solvation-free energy for a leucine zipper formation was examined in the coiled-coil protein GCN4-p1, a typical model for hydrophobic interactions, which demonstrated that water-mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.

Quantifying the Native Energetics Stabilizing Bacteriorhodopsin by Single-Molecule Force Spectroscopy

We quantified the equilibrium (un)folding free energy ΔG_{0} of an eight-amino-acid region starting from the fully folded state of the model membrane-protein bacteriorhodopsin using single-molecule force spectroscopy. Analysis of equilibrium and nonequilibrium data yielded consistent, high-precision determinations of ΔG_{0} via multiple techniques (force-dependent kinetics, Crooks fluctuation theorem, and inverse Boltzmann analysis). We also deduced the full 1D projection of the free-energy landscape in this region. Importantly, ΔG_{0} was determined in bacteriorhodopsin's native bilayer, an advance over traditional results obtained by chemical denaturation in nonphysiological detergent micelles.

Antimicrobial peptides: Promising alternatives in the post feeding antibiotic era

Antimicrobial peptides (AMPs), critical components of the innate immune system, are widely distributed throughout the animal and plant kingdoms. They can protect against a broad array of infection-causing agents, such as bacteria, fungi, parasites, viruses, and tumor cells, and also exhibit immunomodulatory activity. AMPs exert antimicrobial activities primarily through mechanisms involving membrane disruption, so they have a lower likelihood of inducing drug resistance. Extensive studies on the structure-activity relationship have revealed that net charge, hydrophobicity, and amphipathicity are the most important physicochemical and structural determinants endowing AMPs with antimicrobial potency and cell selectivity. This review summarizes the recent advances in AMPs development with respect to characteristics, structure-activity relationships, functions, antimicrobial mechanisms, expression regulation, and applications in food, medicine, and animals.

More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.

Extended solvent-contact model for protein solvation: test cases for dipeptides

Solvation effects are critically important in the structural stabilization and functional optimization of proteins. Here, we propose a new solvation free energy function for proteins, and test its applicability in predicting the solvation free energies of dipeptides. The present solvation model involves the improvement of the previous solvent-contact model assuming that the molecular solvation free energy could be given by the sum over the individual atomic contributions. In addition to the existing solvent-contact term, the modified solvation free energy function includes the self-solvation term that reflects the effects of intramolecular interactions in the solute molecule on solute-solvent interactions. Four kinds of atomic parameters should be determined in this solvation model: atomic fragmental volume, maximum atomic occupancy, atomic solvation, and atomic self-solvation parameters. All of these parameters for 16 atom types are optimized with a standard genetic algorithm in such a way to minimize the difference between the solvation free energies of dipeptides obtained from high-level quantum chemical calculations and those predicted by the solvation free energy function. The solvation free energies of dipeptides estimated from the new solvation model are in good agreement with the quantum chemical results. Therefore, the optimized solvation free energy function is expected to be useful for examining the structural and energetic features of proteins in aqueous solution.

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).

Thermal stability of hydrophobic helical oligomers: a lattice simulation study in explicit water

We investigate the thermal stability of helical hydrophobic oligomers using a three-dimensional, water-explicit lattice model and the Wang-Landau Monte Carlo method. The degree of oligomer helicity is controlled by the parameter ε(mm) < 0, which mimics monomer-monomer hydrogen bond interactions leading to the formation of helical turns in atomistic proteins. We vary |ε(mm)| between 0 and 4.5 kcal/mol and therefore investigate systems ranging from flexible homopolymers (i.e., those with no secondary structure) to helical oligomers that are stable over a broad range of temperatures. We find that systems with |ε(mm)| ≤ 2.0 kcal/mol exhibit a broad thermal unfolding transition at high temperature, leading to an ensemble of random coils. In contrast, the structure of conformations involved in a second, low-temperature, transition is strongly dependent on |ε(mm)|. Weakly helical oligomers are observed when |ε(mm)| ≤ 1.0 kcal/mol and exhibit a low-temperature, cold-unfolding-like transition to an ensemble of strongly water-penetrated globular conformations. For higher |ε(mm)| (1.7 kcal/mol ≤ |ε(mm)| ≤ 2.0 kcal/mol), cold unfolding is suppressed, and the low-temperature conformational transition becomes a "crystallization", in which a "molten" helix is transformed into a defect-free helix. The molten helix preserves ≥50% of the helical contacts observed in the "crystal" at a lower temperature. When |ε(mm)| = 4.5 kcal/mol, we find that conformational transitions are largely suppressed within the range of temperatures investigated.

Mutational analysis of m-values as a strategy to identify cold-resistant substructures of the protein ensemble

Characterizing the native ensemble of protein is an important yet difficult objective of structural biology. The structural dynamics of protein macromolecules play key roles in biological function, but the short lifetimes and low population of near-native states of the protein ensemble limit their ability to be studied directly. In part to address such issues, it was shown recently that the cooperative substructures that populate a protein ensemble could be ascertained by NMR methods performed at very cold temperatures. What is presented here is an argument that these same substructures can also be determined by denaturant-induced unfolding studies performed on protein at room temperature. Data supporting this argument are given for Staphylococcal nuclease, chymotrypsin inhibitor 2, and ubiquitin. The observation of an agreement between the thermodynamics of the protein ensemble simulated under very cold temperatures to the apparent sensitivity of the ensemble to chemical denaturants at room temperature also suggests that the overall structural-thermodynamic character of an ensemble is surprisingly robust and preserved even in the presence of strong denaturing conditions.

Anti-group A streptococcal vaccine epitope: structure, stability, and its ability to interact with HLA class II molecules

Streptococcus pyogenes infections remain a health problem in several countries due to poststreptococcal sequelae. We developed a vaccine epitope (StreptInCor) composed of 55 amino acids residues of the C-terminal portion of the M protein that encompasses both T and B cell protective epitopes. The nuclear magnetic resonance (NMR) structure of the StreptInCor peptide showed that the structure was composed of two microdomains linked by an 18-residue α-helix. A chemical stability study of the StreptInCor folding/unfolding process using far-UV circular dichroism showed that the structure was chemically stable with respect to pH and the concentration of urea. The T cell epitope is located in the first microdomain and encompasses 11 out of the 18 α-helix residues, whereas the B cell epitope is in the second microdomain and showed no α-helical structure. The prediction of StreptInCor epitope binding to different HLA class II molecules was evaluated based on an analysis of the 55 residues and the theoretical possibilities for the processed peptides to fit into the P1, P4, P6, and P9 pockets in the groove of several HLA class II molecules. We observed 7 potential sites along the amino acid sequence of StreptInCor that were capable of recognizing HLA class II molecules (DRB1*, DRB3*, DRB4*, and DRB5*). StreptInCor-overlapping peptides induced cellular and humoral immune responses of individuals bearing different HLA class II molecules and could be considered as a universal vaccine epitope.

Hydration and conformational mechanics of single, end-tethered elastin-like polypeptides

We investigated the effect of temperature, ionic strength, solvent polarity, and type of guest residue on the force-extension behavior of single, end-tethered elastin-like polypeptides (ELPs), using single molecule force spectroscopy (SMFS). ELPs are stimulus-responsive polypeptides that contain repeats of the five amino acids Val-Pro-Gly-Xaa-Gly (VPGXG), where Xaa is a guest residue that can be any amino acid with the exception of proline. We fitted the force-extension data with a freely jointed chain (FJC) model which allowed us to resolve small differences in the effective Kuhn segment length distributions that largely arise from differences in the hydrophobic hydration behavior of ELP. Our results agree qualitatively with predictions from recent molecular dynamics simulations and demonstrate that hydrophobic hydration modulates the molecular elasticity for ELPs. Furthermore, our results show that SMFS, when combined with our approach for data analysis, can be used to study the subtleties of polypeptide-water interactions and thus provides a basis for the study of hydrophobic hydration in intrinsically unstructured biomacromolecules.

Conformational sampling of peptides in cellular environments

Biological systems provide a complex environment that can be understood in terms of its dielectric properties. High concentrations of macromolecules and cosolvents effectively reduce the dielectric constant of cellular environments, thereby affecting the conformational sampling of biomolecules. To examine this effect in more detail, the conformational preference of alanine dipeptide, poly-alanine, and melittin in different dielectric environments is studied with computer simulations based on recently developed generalized Born methodology. Results from these simulations suggest that extended conformations are favored over alpha-helical conformations at the dipeptide level at and below dielectric constants of 5-10. Furthermore, lower-dielectric environments begin to significantly stabilize helical structures in poly-alanine at epsilon = 20. In the more complex peptide melittin, different dielectric environments shift the equilibrium between two main conformations: a nearly fully extended helix that is most stable in low dielectrics and a compact, V-shaped conformation consisting of two helices that is preferred in higher dielectric environments. An additional conformation is only found to be significantly populated at intermediate dielectric constants. Good agreement with previous studies of different peptides in specific, less-polar solvent environments, suggest that helix stabilization and shifts in conformational preferences in such environments are primarily due to a reduced dielectric environment rather than specific molecular details. The findings presented here make predictions of how peptide sampling may be altered in dense cellular environments with reduced dielectric response.

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).

Solvation Model Based on Order Parameters and a Fast Sampling Method for the Calculation of the Solvation Free Energies of Peptides

An analytical solvation model is proposed as a function of an order parameter, which represents the local arrangement of water molecules in the first solvation shell of peptide atoms. The model is combined with a fast sampling method, rotational isomeric state Monte Carlo, to sample efficiently the torsional degrees of freedom on a peptide backbone. This order parameter solvation model is shown to reproduce without ad hoc fitting parameters the solvation free energies of single amino acids and tripeptides with slightly better accuracy than the generalized Born model but with several orders of magnitude improvement in efficiency. This method is a potential candidate for efficiently and accurately tackling some important issues in biophysical chemistry that are related to solvation, for example, protein folding, ligand binding, etc. Our results also present fundamental new insights into solvation. Specifically, the local water geometry, represented in this work by a properly defined order parameter, carries the majority, if not all, of the energetic information of solvation, including solute-solvent interactions and solvent reorganization in the presence of the solute.

Quantitative analysis of helix-coil transition of block copolypeptide, Glu12-Ala12, by combined use of CD and NMR spectroscopy

To investigate helix-coil transition mechanisms, conformations of Glu12-Ala12, EA, in aqueous solution have been studied in detail over the pH range from 2 to 8 and the temperature range from 20 to 60 degrees C using CD and NMR spectroscopy. The 750-MHz NMR spectra displayed excellent dispersion of the backbone amide proton signals, and permitted essentially complete sequence-specific resonance assignments. These assignments, together with short- and medium-range nuclear Overhauser effect (NOE) constraints and coupling constants, enable us to analyze conformational characteristics of all the residues in the EA peptide individually. A combined use of CD and NMR techniques reveals that the EA peptide assumes a stable alpha-helix from Glu12 to Ala19 in 0.1 M NaCl solution at 20 degrees C above pH 7. The alpha-helix is getting longer as decreasing pH. Below pH 4, the peptide assumes the longest alpha-helix from Glu3 to Ala23. The important observation of the present study is that the helix-coil transition occurs stepwise, residue by residue, from both the N- and C-termini of the alpha-helix. No conformational equilibrium between the helical and random-coil states is detected for the residues in the central region of the alpha-helix. Quantitative analysis of temperature-induced helix-to-coil transitions at various pHs provides a pH-independent residual enthalpy change delta H(r) = 0.95 kcal res(-1). Similar values have been reported for a 50-residue alanine-rich peptide (1.2 kcal res(-1)), poly-L-glutamate (1.1 kcal res(-1)), poly-L-lysine (1.1 kcal res(-1)), and poly-L-alanine (0.86 kcal res(-1)). Those investigations, along with our present result, suggest that delta H(r) is mainly determined by the transformation of the backbone associated with the disruption of the intramolecular hydrogen bond. These results should increase our understanding of the helix-coil transition.

Understanding the alpha-helix to coil transition in polypeptides using network rigidity: predicting heat and cold denaturation in mixed solvent conditions

Thermodynamic stability in polypeptides is described using a novel Distance Constraint Model (DCM). Here, microscopic interactions are represented as constraints. A topological arrangement of constraints define a mechanical framework. Each constraint in the framework is associated with an enthalpic and entropic contribution. All accessible topological arrangements of distance constraints form an ensemble of mechanical frameworks, each representing a microstate of the polypeptide. A partition function is calculated exactly using a transfer matrix approach, where in many respects the DCM is similar to the Lifson-Roig model. The crucial difference is that the effect of network rigidity is explicitly calculated for each mechanical framework in the ensemble. Network rigidity is a mechanical interaction that provides a mechanism for long-range molecular cooperativity and enables a proper treatment of the nonadditivity of a microscopic free energy decomposition. Accounting for (1) helix <--> coil conformation changes along the backbone similar to the Lifson-Roig model, (2) i to i + 4 hydrogen-bond formation <--> breaking similar to the Zimm-Bragg model, and (3) structured <--> unstructured solvent interaction (hydration effects), a six-parameter DCM describes normal and inverted helix-coil transitions in polypeptides. Under suitable mixed solvent conditions heat and cold denaturation is predicted. Model parameters are fitted to experimental data showing different degrees of cold denaturation in monomeric polypeptides in aqueous hexafluoroisopropanol (HFIP) solution at various HFIP concentrations. By assuming a linear HFIP concentration dependence (up to 6% by mole fraction) on model parameters, all essential experimentally observed features are captured.

Direct access to the cooperative substructure of proteins and the protein ensemble via cold denaturation

The modern view of protein thermodynamics predicts that proteins undergo cold-induced unfolding. Unfortunately, the properties of proteins and water conspire to prevent the detailed observation of this fundamental process. Here we use protein encapsulation to allow cold denaturation of the protein ubiquitin to be monitored by high-resolution NMR at temperatures approaching -35 degrees C. The cold-induced unfolding of ubiquitin is found to be highly noncooperative, in distinct contrast to the thermal melting of this and other proteins. These results demonstrate the potential of cold denaturation as a means to dissect the cooperative substructures of proteins and to provide a rigorous framework for testing statistical thermodynamic treatments of protein stability, dynamics and function.

Exploring Flory's isolated-pair hypothesis: statistical mechanics of helix-coil transitions in polyalanine and the C-peptide from RNase A

To evaluate Flory's isolated-pair hypothesis in the context of helical peptides, we explore equilibrium conformations of alpha-helix-forming polypeptides as a function of temperature by means of replica exchange molecular dynamics in conjunction with the CHARMM/GB implicit solvent force field and the weighted histogram analysis method. From these simulations, Zimm-Bragg parameters, s and sigma, of Ac-Alan-NMe are computed as a function of temperature. The values obtained for s(T) and sigma(T) remain unchanged along the length of the polypeptide except for very short chains and yield results consistent with measurements based on short helix-forming peptides but suggest larger s values than anticipated from polymer-based measurements. From direct estimates of the density of states for Ac-Alan-NMe (n = 3-20) and peptide constructs based on the C peptide from RNase A, the conformational entropy is calculated versus temperature. The calculated S(T) shows a clear proportionality to the chain length over a wide range of temperature. This is observed in polypeptides with both significantly branched and simple methyl (alanine) side chains. These results provide evidence for the validity of Flory's isolated pair hypothesis, at least in the context of helical peptides and helix-to-coil transitions in these peptides.

Thermodynamics of the coil-alpha-helix transition of amphipathic peptides in a membrane environment: the role of vesicle curvature

The binding of peptides or proteins to a bilayer membrane is often coupled with a random coil-->alpha-helix transition. Knowledge of the energetics of this membrane-induced folding event is essential for the understanding of the mechanism of membrane activity. In a recent study [Wieprecht et al., J. Mol. Biol. 294 (1999) 785-794], we have developed an approach which allows an analysis of the energetics of membrane-induced folding. We have systematically varied the helix content of the amphipathic peptide magainin-2-amide by synthesizing analogs where two adjacent amino acid residues were substituted by their corresponding D-enantiomers and have measured their binding to small unilamellar vesicles (SUVs). Correlation of the binding parameters with the helicities allowed the evaluation of the thermodynamic parameters of helix formation. Since SUVs (30 nm in diameter) are characterized by a non-ideal lipid packing due to their high membrane curvature, we have now extended our studies to large unilamellar vesicles (LUVs) (100 nm in diameter) with a lipid packing close to planar membranes. While the free energy of binding was similar for SUVs and LUVs, the binding enthalpies and entropies were distinctly different for the two membrane systems. The thermodynamic parameters of the coil-helix transition were nevertheless not affected by the vesicle size. Helix formation at the membrane surface of LUVs (SUVs) was characterized by an enthalpy change of -0.8 (-0.7) kcal/mol per residue, an entropy change of-2.3 (-1.9) cal/mol K per residue, and a free energy change of -0.12 (-0.14) kcal/mol per residue. Helix formation accounted for approximately 50% of the free energy of binding underlining its major role as a driving force for membrane-binding.

The enthalpy of the alanine peptide helix measured by isothermal titration calorimetry using metal-binding to induce helix formation

The goal of this study is to use the model system described earlier to make direct measurements of the enthalpy of helix formation at different temperatures. For this we studied model alanine peptides in which helix formation can be triggered by metal (La(3+)) binding. The heat of La(3+) interaction with the peptides at different temperatures is measured by isothermal titration calorimetry. Circular dichroism spectroscopy is used to follow helix formation. Peptides of increasing length (12-, 16-, and 19-aa residues) that contain a La(3+)-binding loop followed by helices of increasing length, are used to separate the heat of metal binding from the enthalpy of helix formation. We demonstrate that (i) the enthalpy of helix formation is -0.9 +/- 0.1 kcal/mol; (ii) the enthalpy of helix formation is independent of the peptide length; (iii) the enthalpy of helix formation does not depend significantly on temperature in the range from 5 to 45 degrees C, suggesting that the heat capacity change on helix formation is very small. Thus, the use of metal binding to induce helix formation has an enormous potential for measuring various thermodynamic properties of alpha-helices.

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.

Role of helix formation for the retention of peptides in reversed-phase high-performance liquid chromatography

In order to get insight into the role of helix formation for retention in reversed-phase HPLC, we have studied the isocratic retention behavior of amphipathic and non-amphipathic potentially helical model peptides. Plots of the logarithmic capacity factor in absence of organic solvent (ln k0) versus l/T were used to derive the enthalpy, deltaH0, the free energy, deltaG0, the entropy of interaction, deltaS0, and the heat capacity change, deltaCp. Retention of all peptides was accompanied by negative deltaCp revealing that hydrophobic interactions play a large role independent of peptide sequence and secondary structure. deltaH0 was negative for the amphipathic analogs and was attributed mainly to helix formation of these peptides upon interaction with the stationary phase. In contrast, deltaH0 was considerably less exothermic or even endothermic for the non-amphipathic analogs. The differences in helix formation between the individual analogs were quantified on the basis of thermodynamic data of helix formation previously derived for peptides in a hydrophobic environment. Correlation of the helicity with the free energy of stationary phase interaction revealed that helix formation accounts for approximately 40-70% of deltaG0, and is hence in addition to the hydrophobic effect a major driving force of retention.

Thermodynamics of the alpha-helix-coil transition of amphipathic peptides in a membrane environment: implications for the peptide-membrane binding equilibrium

Amphipathic alpha-helices are the membrane binding motif in many proteins. The corresponding peptides are often random coil in solution but are folded into an alpha-helix upon interaction with the membrane. The energetics of this ubiquitous folding process are still a matter of conjecture. Here, we present a new method to quantitatively analyze the thermodynamics of peptide folding at the membrane interface. We have systematically varied the helix content of a given amphipathic peptide when bound to the membrane and have correlated the thermodynamic binding parameters determined by isothermal titration calorimetry with the alpha-helix content obtained by circular dichroism spectroscopy. The peptides investigated were the antibiotic magainin 2 amide and three analogs in which two adjacent amino acid residues were substituted by their d-enantiomers. The thermodynamic parameters controlling the alpha-helix formation were found to be linearly related to the helicity of the membrane-bound peptides. Helix formation at the membrane surface is characterized by an enthalpy change of DeltaH(helix) approximately -0.7 kcal/mol per residue, an entropy change of DeltaS(helix) approximately -1.9 cal/molK residue and a free energy change of DeltaG(helix)=-0.14 kcal/mol residue. Helix formation is a strong driving force of peptide insertion into the membrane and accounts for about 50 % of the free energy of binding. An increase in temperature entails an unfolding of the membrane-bound helix. The temperature dependence can be described with the Zimm-Bragg theory and the enthalpy of unfolding agrees with that deduced from isothermal titration calorimetry.

Monte Carlo simulations of beta-hairpin folding at constant temperature

Monte Carlo simulations were applied to beta-hairpin folding of a valine-based peptide. Two valine residues in the middle of the peptide were substituted with glycine, to serve as turn residues. Unlike lattice model simulations, structure prediction methods, and unfolding simulations, our simulations used an atom-based model, constant temperature (274 K), and non-beta-hairpin initial conformations. Based on the concept of solvent reference, the effective energy function simplified the solvent calculation and overcame the multiple minima problem. Driven by the hydrophobic interaction, the peptide first folded into a compact U-shaped conformation with a central turn, in analogy to the initial collapse with simultaneous nucleation in protein folding. The peptide units in the U-shaped conformation then reoriented, gradually forming hydrogen bonds in the beta-hairpin pattern from the beta-turn to the ends of the strands. With the same energy function, an alanine-based peptide folded into helix-dominated structures. The basic structure types (alpha-helix or beta-hairpin) that formed during the simulations depended upon the amino acid sequence. Compared with helix, beta-hairpin folding is driven mainly by the hydrophobic interaction. Hydrogen bonding is necessary to maintain the ordered secondary structure.

Thermal stability of Escherichia coli ribonuclease HI and its active site mutants in the presence and absence of the Mg2+ ion. Proposal of a novel catalytic role for Glu48

Escherichia coli ribonuclease HI, which requires divalent cations (Mg2+ or Mn2+) for activity, was thermostabilized by 2.6-3.0 kcal/mol in the presence of the Mg2+, Mn2+, or Ca2+ ion, probably because the negative charge repulsion around the active site was canceled upon the binding of these metal ions. The dissociation constants were determined to be 0.71 mM for Mg2+, 0.035 mM for Mn2+, and 0.16 mM for Ca2+. Likewise, various active site mutants at Asp10, Glu48, Asp70, or Asp134 were thermostabilized by 0.4-3.0 kcal/mol in the presence of the Mg2+ ion, suggesting that this ion binds to these mutant proteins as well. The dissociation constants of Mg2+ were determined to be 9.8 mM for D10N, 1.1 mM for E48Q, 18.8 mM for D70N, and 1.8 mM for D134N. Thus, the mutation of Asp10 or Asp70 to Asn considerably impairs the Mg2+ binding, whereas the mutation of Glu48 to Gln or Asp134 to Asn does not. Comparison of the thermal stability of the mutant proteins with that of the wild-type protein in the absence of the Mg2+ ion suggests that the negative charge repulsion between Asp10 and Asp70 is responsible for the binding of the metal cofactor. Glu48 may be required to anchor a water molecule, which functions as a general acid.

Effects of pH and Ca2+ on heterodimer and heterotetramer formation by chromogranin A and chromogranin B

The two major proteins of the secretory vesicles of neuroendocrine cells, chromogranin A (CGA) and chromogranin B (CGB), have been shown to undergo pH- and Ca2+-dependent conformational changes and aggregation and have been suggested to play essential roles during secretory vesicle biogenesis in the trans-Golgi network. CGA has been shown to exist primarily in a tetrameric state at pH 5.5 and primarily in a dimeric state at pH 7.5, and CGB has been shown to exist in a monomeric state at both pH 5.5 and pH 7.5. Using purified CGA and CGB, it recently has been shown that CGA interacts with CGB at pH 5.5 (Yoo, S. H.(1996) J. Biol. Chem. 271, 1558-1565). In expanding this investigation, we have studied the temperature dependence of the pH-dependent interaction of CGA and CGB by analytical ultracentrifugation and found that two molecules of CGA bound to two molecules of CGB at pH 5.5 with DeltaG0 values of -43.6 kcal/mol in the absence of Ca2+ at 37 degrees C and -40.3 kcal/mol in the presence of 0.1 mM Ca2+. However, one molecule of CGA bound to one molecule of CGB at pH 7.5 with DeltaG0 values of -13.6 kcal/mol in the absence of Ca2+ at 37 degrees C. The magnitude of DeltaG0 values increased with increasing temperatures at both pH values. However, the values for enthalpy and entropy changes decreased with increasing temperatures in both pH levels, suggesting formation of more ordered structures. In the absence of Ca2+ at pH 5. 5, the heterotetramerization reaction at 37 degrees C was entropically driven, whereas in the presence of Ca2+ (0.1 mM) the heterotetramerization was virtually an enthalpic reaction. On the other hand, the heterodimer formation in the absence of Ca2+ at pH 7. 5 showed large negative enthalpy and entropy changes at 37 degrees C, indicating an enthalpic interaction compensated by entropic changes. In view of the interaction of tetrameric CGA with tetrameric inositol 1,4,5-trisphosphate (IP3) receptor and the existence of heterotetrameric IP3 receptor in the cell, the heterotetramer formation by CGA and CGB not only raises the possibility of interaction between the heterotetrameric chromogranin and heterotetrameric IP3 receptor but also appears to reflect their important roles in the cell.

Folding simulations of alanine-based peptides with lysine residues

The folding of short alanine-based peptides with different numbers of lysine residues is simulated at constant temperature (274 K) using the rigid-element Monte Carlo method. The solvent-referenced potential has prevented the multiple-minima problem in helix folding. From various initial structures, the peptides with three lysine residues fold into helix-dominated conformations with the calculated average helicity in the range of 60-80%. The peptide with six lysine residues shows only 8-14% helicity. These results agree well with experimental observations. The intramolecular electrostatic interaction of the charged lysine side chains and their electrostatic hydration destabilize the helical conformations of the peptide with six lysine residues, whereas these effects on the peptides with three lysine residues are small. The simulations provide insight into the helix-folding mechanism, including the beta-bend intermediate in helix initiation, the (i, i + 3) hydrogen bonds, the asymmetrical helix propagation, and the asymmetrical helicities in the N- and C-terminal regions. These findings are consistent with previous studies.

Thermodynamic study of the pH-dependent interaction of chromogranin A with an intraluminal loop peptide of the inositol 1,4,5-trisphosphate receptor

The secretory vesicles of adrenal chromaffin cells have previously been identified as a major inositol 1,4,5-trisphosphate (IP3)-sensitive Ca2+ store, and their Ca2+ store role has been attributed to the presence of chromogranin A, a high capacity, low affinity Ca2+ binding protein. Chromogranin A has since been shown to exist primarily in a dimeric state at pH 7.5 and primarily in a tetrameric state at the intravesicular pH of 5.5 and has also been shown to interact with the membrane proteins of secretory vesicles at pH 5.5, including a 260-kDa protein reactive to IP3 receptor antibody [Yoo, S. H. (1994) J. Biol. Chem. 269, 12001-12006]. In a recent study, chromogranin A was shown to interact with one of the intraluminal loop regions of the IP3 receptor at pH 5.5 but not at pH 7.5 [Yoo, S. H., & Lewis, M. S. (1994) FEBS Lett. 341, 28-32]. To gain further insight, we have studied the temperature dependence of the pH-dependent interaction of chromogranin A with the intraluminal peptide of the the IP3 receptor by analytical ultracentrifugation, using multiwavelength scan analysis, and found that four molecules of the intraluminal domain peptide of the IP3 receptor bound to each chromogranin A tetramer with delta Go values ranging from -23.6 to -27.6 kcal mol-1 in the absence and presence of 35 mN Ca2+.(ABSTRACT TRUNCATED AT 250 WORDS)

Estimation of changes in side chain configurational entropy in binding and folding: general methods and application to helix formation

Theoretical estimations of changes in side chain configurational entropy are essential for understanding the different contributions to the overall thermodynamic behavior of important biological processes like folding and binding. The configurational entropy of any given side chain in any particular protein can be evaluated from the complete energy profile of the side chain. Calculations of the energy profiles can be performed using the side chain single bond dihedrals as the only independent variables as long as the structures at each value of the dihedrals are allowed to relax through small changes in the valence bond angles. The probabilities of different side chain conformers obtained from these energy profiles are very similar to the conformer populations obtained by analysis of side chain preferences in the proteins of the Protein Data Bank. Also, side chain conformational entropies obtained from the energy profiles agree extremely well with those obtained from the Protein Data Bank conformer populations. Changes in side chain configurational entropy in binding and folding can be computed as differences in conformational entropy because, in most cases, the frequency of the rotational oscillation around the energy minimum of any given conformer does not appear to change significantly in the reactions. Changes of side chain conformational entropy calculated in this way were compared with experimental values. The only available experimental data--the effect of side chain substitution on the stability of alpha-helices--were used for this comparison. The experimental values were corrected to subtract the solvent contributions. This comparison yields an excellent agreement between calculated and experimental values, validating not only the theoretical estimates but also the separability of the entropic contributions into configurational terms and solvation related terms.

Thermal unfolding of tetrameric melittin: comparison with the molten globule state of cytochrome c

Whereas melittin at micromolar concentrations is unfolded under conditions of low salt at neutral pH, it transforms to a tetrameric alpha-helical structure under several conditions, such as high peptide concentration, high anion concentration, or alkaline pH. The anion- and pH-dependent stabilization of the tetrameric structure is similar to that of the molten globule state of several acid-denatured proteins, suggesting that tetrameric melittin might be a state similar to the molten globule state. To test this possibility, we studied the thermal unfolding of tetrameric melittin using far-UV CD and differential scanning calorimetry. The latter technique revealed a broad but distinct heat absorption peak. The heat absorption curves were consistent with the unfolding transition observed by CD and were explainable by a 2-state transition mechanism between the tetrameric alpha-helical state and the monomeric unfolded state. From the peptide or salt-concentration dependence of unfolding, the heat capacity change upon unfolding was estimated to be 5 kJ (mol of tetramer)-1 K-1 at 30 degrees C and decreased with increasing temperature. The observed change in heat capacity was much smaller than that predicted from the crystallographic structure (9.2 kJ (mol of tetramer)-1 K-1), suggesting that the hydrophobic residues of tetrameric melittin in solution are exposed in comparison with the crystallographic structure. However, the results also indicate that the structure is more ordered than that of a typical molten globule state. We consider that the conformation is intermediate between the molten globule state and the native state of globular proteins.

Elucidating the folding problem of helical peptides using empirical parameters

Using an empirical analysis of experimental data we have estimated a set of energy contributions which accounts for the stability of isolated alpha-helices. With this database and an algorithm based on statistical mechanics, we describe the average helical behaviour in solution of 323 peptides and the helicity per residue of those peptides analyzed by nuclear magnetic resonance. Moreover the algorithm successfully detects the alpha-helical tendency, in solution, of a peptide corresponding to a beta-strand of ubiquitin.

Helix folding simulations with various initial conformations

Using a solvent-referenced energy calculation, a 16-residue peptide with alanine side chains folded into predominantly alpha-helical conformations during constant temperature (274 K) simulations. From different initial conformations, helical conformations were reached and the multiple energy minima did not become a serious problem. Under the same conditions, the simulation did not indiscriminately fold a sequence such as polyglycine into stable helices. Interesting observations from the simulations relate to the folding mechanism. The electrostatic interactions between the successive amides favored extended conformations (or beta strands) and caused energy barriers to helix folding. beta-bends were observed as intermediates during helix nucleation. The helix propagation toward the C-terminus seemed faster than that toward the N-terminus. In helical conformations, hydrogen bond oscillation between the i,i+ 4 and the i,i+3 patterns was observed. The i,i+3 hydrogen bonds occurred more frequently during helix propagation and deformation near both ends of the helical segment.

Helix-forming tendencies of nonpolar amino acids predicted by Monte Carlo simulated annealing

Monte Carlo simulated annealing is applied to the study of the alpha-helix-forming tendencies of seven nonpolar amino acids, Ala, Leu, Met, Phe, Ile, Val, and Gly. Homooligomers of 10 amino acids are used and the helix tendency is calculated by folding alpha-helicies from completely random initial conformations. The results of the simulation imply that Met, Ala, and Leu are helix formers and that Val, Ile, and Gly are helix breakers, while Phe comes in between the two groups. The differences between helix formers and breakers turned out to be large in agreement with the recent experiments with short peptides. It is argued from the energy distributions of the obtained conformations that the helix tendency is small for the helix breakers because of steric hindrance of side chains. Homoglycine is shown to favor a random coil conformation. The beta-strand tendencies of the same homooligomers are also considered, and they are shown to agree with the frequencies of amino acids in beta-sheet from the protein data base.

De novo design and creation of a stable artificial protein

Protein de novo design has been performed, as an exercise of the inverse folding problem. A beta/alpha-barrel protein was designed and synthesized using the Escherichia coli expression system for the structural characterization. A tertiary model with a two-fold symmetry was built, based upon the geometrical parameters extracted from X-ray crystal structures of several beta/alpha-barrel proteins. Amino acid frequencies at each position on the alpha- and beta-structures were investigated, and an amino acid sequence with 201 residues was designed. The associated gene was chemically synthesized and the fusion protein with human growth hormone was expressed in Escherichia coli. The purified protein after being cleaved and refolded was found to be stable and globular with the large amount of secondary structures. However, it has similar characteristics to the molten globules of natural proteins, with loose packing of side-chains. The approach for the tight packing is discussed.

Thermodynamic characterization of an artificially designed amphiphilic alpha-helical peptide containing periodic prolines: observations of high thermal stability and cold denaturation

To investigate the structural stability of proteins, we analyzed the thermodynamics of an artificially designed 30-residue peptide. The designed peptide, NH2-EELLPLAEALAPLLEALLPLAEALAPLLKK-COOH (PERI COIL-1), with prolines at i + 7 positions, forms a pentameric alpha-helical structure in aqueous solution. The thermal denaturation curves of the CD at 222 nm (pH 7.5) show an unusual cold denaturation occurring well above 0 degrees C and no thermal denaturation is observable under 90 degrees C. This conformational change is reversible and depends on peptide concentration. A 2-state model between the monomeric denatured state (5D) and the pentameric helical state (H5) was sufficient to analyze 5 thermal denaturation curves of PERI COIL-1 with concentrations between 23 and 286 microM. The analysis was carried out by a nonlinear least-squares method using 3 fitting parameters: the midpoint temperature, Tm, the enthalpy change, delta H(Tm), and the heat capacity change, delta Cp. The association number (n = 5) was determined by sedimentation equilibrium and was not used as a fitting parameter. The heat capacity change suggests that the hydrophobic residues are buried in the helical state and exposed in the denatured one, as it occurs normally for natural globular proteins. On the other hand, the enthalpy and the entropy changes have values close to those found for coiled-coils and are quite distinct from typical values reported for natural globular proteins. In particular, the enthalpy change extrapolated at 110 degrees C is about 3 kJ/mol per amino acid residue, i.e., half of the value found for globular proteins.(ABSTRACT TRUNCATED AT 250 WORDS)