Molecular thermodynamic model to predict the alpha-helical secondary structure of polypeptide chains in solution.

The native state of a protein molecule in aqueous solutions represents one of the lowest states of Gibbs energy [Anfinsen, C.B. (1973) Science 181, 223-230]. Much progress has been made about the rules of protein folding [King, J. (1989) Chem. Eng. News 67, 32-54] and the dominant forces in protein folding [Dill, K.A. (1990) Biochemistry 29, 7133-7155]. However, the quantitative contributions of different Gibbs energy terms to protein stability remains a controversial issue [Moult, J., & Unger, R. (1991) Biochemistry 30, 3816-3824]. A molecular thermodynamic model has been proposed for the Gibbs energy of folding a residue in aqueous homopolypeptides from a random-coiled state to either the alpha-helix state or the beta-sheet state [Chen, C.-C., Zhu, Y., King, J.A., & Evans, L.B. (1992) Biopolymers 32, 1375-1392]. In this work, we present a generalization of the molecular thermodynamic model for the Gibbs energy of folding natural and synthetic heteropolypeptides from random-coiled conformations into alpha-helical conformations. The generalized model incorporates the intrinsic folding potential due to residue-solvent interactions, the cooperative folding effect due to residue-residue interactions, and the location and length of alpha-helices. The utility of the model was demonstrated by examining the stability of alpha-helical conformations of a number of natural polypeptides including C-peptide (residues 1-13) and S-peptide (residues 1-20) of RNase A (bovine pancreatic ribonuclease A), the P alpha fragment in BPTI (bovine pancreatic trypsin inhibitor), and synthetic polypeptides (the copolymers of different amino acid residues) including alanine-based peptides (16 or 17 residues long) in water. The computed Gibbs energies correspond well with the experimental data on helicity. The results also accounted for the effects of amino acid substitution and temperature on the stability of alpha-helical conformations of the test polypeptides.     

Free energy perturbation study on a Trp-binding mutant (Ser88-->Cys) of the trp-repressor.

The Ser88-->Cys mutant of the trp-repressor showed a lower affinity for the corepressor than the wild-type repressor [delta delta G = 1.7 +/- 0.3 kcal/mol, Chou and Matthews (1989) J. Biol. Chem., 264, 18314-18319]. A molecular dynamics/free energy cycle perturbation study was performed to understand the origin of the decreased affinity. A value (delta delta G = 1.58 +/- 0.28 kcal/mol) comparable with the experimental value was obtained by the simulation. Free energy component analysis revealed that destabilization of the van der Waals interaction between Ser88 and Trp109 (corepressor) mainly contributed to the decreased affinity of the mutant. The rotational transition of the hydroxyl (sulfhydryl) group of Ser88 (Cys88) during the simulations affected the contributions of Arg84 and water to the free energy change in the aporepressor and those of Arg84 and Trp109 to that in the holorepressor. However, the contributions from different residues compensated each other, and the total free energy changes were almost invariable in the various simulations.     

Differential scanning calorimetric study of the thermal unfolding of mutant forms of phage T4 lysozyme.

In two recent papers, we reported the effects of several point mutations on the thermodynamics of the thermal unfolding of the lysozyme of phage T4 as determined by differential scanning calorimetry. The mutants studied were R96H [Kitamura, S., & Sturtevant, J.M. (1989) Biochemistry 28, 3788-3792] and T157 replaced by A, E, I, L, N, R, and V [Connelly, P., Ghosaini, L., Hu, C.-Q., Kitamura, S., Tanaka, A., & Sturtevant, J.M. (1991) Biochemistry 30, 1887-1891]. Here we report the results of a similar study of the single mutations A82P, A93P, and G113A and the double mutation C54T:C97A. The three single mutants all show small apparent stabilization at pH 2.5 and 46.2 degrees C (the denaturational temperature of the wild-type protein), amounting to -0.5 +/- 0.4 kcal mol-1 in free energy, whereas the double mutant shows a weak apparent destabilization, +0.8 +/- 0.4 kcal mol-1. As in all our previous studies of mutant proteins, the enthalpy changes produced by these mutations are in general of much larger magnitude than the corresponding free energy changes and frequently of opposite sign.     

Strength and co-operativity of contributions of surface salt bridges to protein stability

Many of the interactions that stabilize proteins are co-operative and cannot be reduced to a sum of pairwise interactions. Such interactions may be analysed by protein engineering methods using multiple thermodynamic cycles comprising wild-type protein and all combinations of mutants in the interacting residues. There is a triad of charged residues on the surface of barnase, comprising residues Asp8, Asp12 and Arg110, that interact by forming two exposed salt bridges. The three residues have been mutated to alanine to give all the single, double and triple mutants. The free energies of unfolding of wild-type and the seven mutant proteins have been determined and the results analysed to give the contributions of the residues in the two salt bridges to protein stability. It is possible to isolate the energies of forming the salt bridges relative to the solvation of the separated ions by water. In the intact triad, the apparent contribution to the stabilization energy of the protein of the salt bridge between Asp12 and Arg110 is -1.25 kcal mol-1, whereas that of the salt bridge between Asp8 with Arg110 is -0.98 kcal mol-1. The strengths of the two salt bridges are coupled: the energy of each is reduced by 0.77 kcal mol-1 when the other is absent. The salt-linked triad, relative to alanine residues at the same positions, does not contribute to the stability of the protein since the favourable interactions of the salt bridges are more than offset by other electrostatic and non-electrostatic energy terms. Salt-linked triads occur in other proteins, for example, haemoglobin, where the energy of only the salt-bridge term is important and so the coupling of salt bridges could be of general importance to the stability and function of proteins.     

Sequence context modulates the stability of a GxxxG-mediated transmembrane helix-helix dimer

To quantify the relationship between sequence and transmembrane dimer stability, a systematic mutagenesis and thermodynamic study of the protein-protein interaction residues in the glycophorin A transmembrane helix-helix dimer was carried out. The results demonstrate that the glycophorin A transmembrane sequence dimerizes when its GxxxG motif is abolished by mutation to large aliphatic residues, suggesting that the sequence encodes an intrinsic propensity to self-associate independent of a GxxxG motif. In the presence of an intact GxxxG motif, the glycophorin A dimer stability can be modulated over a span of -0.5 kcal mol(-1) to +3.2 kcal mol(-1) by mutating the surrounding sequence context. Thus, these flanking residues play an active role in determining the transmembrane dimer stability. To assess the structural consequences of the thermodynamic effects of mutations, molecular models of mutant transmembrane domains were constructed, and a structure-based parameterization of the free energy change due to mutation was carried out. The changes in association free energy for glycophorin A mutants can be explained primarily by changes in packing interactions at the protein-protein interface. The energy cost of removing favorable van der Waals interactions was found to be 0.039 kcal mol(-1) per A2 of favorable occluded surface area. The value corresponds well with estimates for mutations in bacteriorhodopsin as well as for those mutations in the interiors of soluble proteins that create packing defects.     

Role of arginine 129 in heparin binding and activation of antithrombin

The contribution of Arg(129) of the serpin, antithrombin, to the mechanism of allosteric activation of the protein by heparin was determined from the effect of mutating this residue to either His or Gln. R129H and R129Q antithrombins bound pentasaccharide and full-length heparins containing the antithrombin recognition sequence with similar large reductions in affinity ranging from 400- to 2500-fold relative to the control serpin, corresponding to a loss of 28-35% of the binding free energy. The salt dependence of pentasaccharide binding showed that the binding defect of the mutant serpin resulted from the loss of approximately 2 ionic interactions, suggesting that Arg(129) binds the pentasaccharide cooperatively with other residues. Rapid kinetic studies showed that the mutation minimally affected the initial low affinity binding of heparin to antithrombin, but greatly affected the subsequent conformational activation of the serpin leading to high affinity heparin binding, although not enough to disfavor activation. Consistent with these findings, the mutant antithrombin was normally activated by heparin for accelerated inhibition of factor Xa and thrombin. These results support an important role for Arg(129) in an induced-fit mechanism of heparin activation of antithrombin wherein conformational activation of the serpin positions Arg(129) and other residues for cooperative interactions with the heparin pentasaccharide so as to lock the serpin in the activated state.     

Understanding the thermodynamic stability of an RNA hairpin and its mutant.

The thermodynamic stability of RNA hairpin loops has been a subject of considerable interest in the recent past (Wimberly et al., 1991). There have been experimental reports indicating that the hairpins with a C(UUCG)G loop sequence are thermodynamically very stable (Wimberly et al., 1991). We used the solution structure of GGAC(UUCG)GUCC (Cheong et al., 1990; Varani et al., 1991) as the starting conformation in our attempt to understand its thermodynamic stability. We carried out molecular dynamics/free energy simulations to understand the basis for the destabilization of the C(UUCG)G loop by mutating cytosine (C7)-->uracil. Because of the limited length of simulation and the presence of kinetic barriers (solvent intervention) to the uracil-->cytosine mutation, all of our computed free energy differences are based on multiple forward simulations. Based on these calculations we find that the cytosine-->uracil mutation in the loop destabilizes it by approximately 1.5kcal/mol relative to that of the reference state, an A-form RNA but with cytosine (C7) looped out. This is the same sign and magnitude as that observed in the thermodynamic studies carried out by Varani et al.(1991). We have carried out free energy component analysis to understand the effect of mutating the cytosine residue to uracil on the thermodynamic stability of the C(UUCG)G hairpin loops. Our calculations show that the most significant contribution to the stability is from the phosphate group linking U5 and U6, which favors the cytosine residue over uracil by about 6.0 kcal/mol. The residues U5, U6, and G8 in the loop region also contribute significantly to the stability. The contributions from the salt and solvent compensate each other, indicating the dynamic nature of interactions of the environment with the nucleic acid system and the coupling between these two components.     

Energetics of sequence-specific protein-DNA association: computational analysis of integrase Tn916 binding to its target DNA

The N-terminal domain of the bacterial integrase Tn916 specifically recognizes the 11 bp DNA target site by positioning the face of a three-stranded beta-sheet into the major groove. Binding is linked to structural adaptation. We have characterized INT-DBD binding to DNA in detail by calorimetry [Milev, S., Gorfe, A., Karshikoff, A., Clubb, R. T., Bosshard, H. R., and Jelesarov, I. (2003) Biochemistry 42, 3481-3491]. Our thermodynamic analysis has indicated that the major driving force of association is the hydrophobic effect while polar interactions contribute less. To gain more comprehensive information about the binding process, we performed a computational analysis of the binding free energy and report here the results. A hybrid molecular mechanics/continuum approach was followed. The total binding free energy is predicted with reasonable accuracy. The calculations confirm that nonpolar effects stabilize the protein-DNA complex while electrostatics opposes binding. Structural changes optimizing surface complementarity are costly in terms of energy. The energetic consequences from the replacement of nine DNA-contacting residues by alanine were investigated. The calculations correctly predict the binding affinity decrease of eight mutations and the destabilizing effect of one wild-type residue. Bulky side chains stabilize the wild-type complex through packing interactions and favorable nonpolar dehydration, but the net nonpolar energy changes do not correlate with the relative affinity loss upon mutation. Discrete protein-DNA electrostatic interactions may be net stabilizing or net destabilizing depending on the local environment. In contrast to nonpolar energy changes, the magnitude of the electrostatic free energy ranks the mutations according to the experimentally measured DeltaDeltaG. Free energy decomposition analysis from a structural perspective leads to detailed information about the thermodynamic strategy used by INT-DBD for sequence-specific DNA binding.     

Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Gly --> Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.     

The effect of surface charge balance on thermodynamic stability and kinetics of refolding of firefly luciferase

Thermodynamic stability and refolding kinetics of firefly luciferase and three representative mutants with depletion of negative charge on a flexible loop via substitution of Glu by Arg (ER mutant) or Lys (EK mutant) as well as insertion of another Arg in ER mutants (ERR mutant) was investigated. According to thermodynamic studies, structural stability of ERR and ER mutants are enhanced compared to WT protein, whereas, these mutants become prone to aggregation at higher temperatures. Accordingly, it was concluded that enhanced structural stability of mutants depends on more compactness of folded state, whereas aggregation at higher temperatures in mutants is due to weakening of intermolecular repulsive electrostatic interactions and increase of intermolecular hydrophobic interactions. Kinetic results indicate that early events of protein folding are accelerated in mutants.     

Enhanced thermodynamic stabilities of yeast iso-1-cytochromes c with amino acid replacements at positions 52 and 102.

We have determined the structures and thermodynamic stabilities of the wild type Asn-52 and unusually thermostable mutant Ile-52 yeast iso-1-cytochromes c (Das, G., Hickey, D. R. McLendon, D., McLendon, G., and Sherman, F. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 496-499). Although both structures were similar, Water-166, buried within the wild type protein, is excluded from the Ile-52 mutant, which substantially reorganizes the local hydrogen bonding. Wild type Cys-102 was replaced with alanine or serine to eliminate dimerization in vitro. The Cys-102 (wild type), Ala-102, and Ser-102 proteins were equally stable, whereas the chemically modified Cys-102-SCH3 was less stable. The order of stability observed with replacements at positions 52 and 102 was as follows: Ile-52 Ala-102 greater than Ala-52 Ala-102 greater than Asn-52 Ala-102 ("normal") greater than Gly-52 Ala-102. No significant stabilization was attributed to potential energy interactions expressed as helix-forming propensities of replacements at position 52. A high correlation between differences in free energy changes and transfer free energies suggests hydrophobic interactions are the main factor for enhancing stability in the Ile-52 mutant. Additional possible contributions to the thermostability of the Ile-52 variant are energetic effects due to packing and hydrogen bonding changes surrounding position 52.     

Protein-protein recognition: an experimental and computational study of the R89K mutation in Raf and its effect on Ras binding

Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras-Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [gamma-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9+/-1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein-protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 A away, namely residues 37-41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2-3 specific water molecules. This analysis will provide insights into the large experimental database of Ras-Raf mutations.     

Pairwise electrostatic interactions between alpha-neurotoxins and gamma, delta, and epsilon subunits of the nicotinic acetylcholine receptor

alpha-Neurotoxins bind with high affinity to alpha-gamma and alpha-delta subunit interfaces of the nicotinic acetylcholine receptor. Since this high affinity complex likely involves a van der Waals surface area of approximately 1200 A(2) and 25-35 residues on the receptor surface, analysis of side chains should delineate major interactions and the orientation of bound alpha-neurotoxin. Three distinct regions on the gamma subunit, defined by Trp(55), Leu(119), Asp(174), and Glu(176), contribute to alpha-toxin affinity. Of six charge reversal mutations on the three loops of Naja mossambica mossambica alpha-toxin, Lys(27) --> Glu, Arg(33) --> Glu, and Arg(36) --> Glu in loop II reduce binding energy substantially, while mutations in loops I and III have little effect. Paired residues were analyzed by thermodynamic mutant cycles to delineate electrostatic linkages between the six alpha-toxin charge reversal mutations and three key residues on the gamma subunit. Large coupling energies were found between Arg(33) at the tip of loop II and gammaLeu(119) (-5.7 kcal/mol) and between Lys(27) and gammaGlu(176) (-5.9 kcal/mol). gammaTrp(55) couples strongly to both Arg(33) and Lys(27), whereas gammaAsp(174) couples minimally to charged alpha-toxin residues. Arg(36), despite strong energetic contributions, does not partner with any gamma subunit residues, perhaps indicating its proximity to the alpha subunit. By analyzing cationic, neutral and anionic residues in the mutant cycles, interactions at gamma176 and gamma119 can be distinguished from those at gamma55.     

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.     

Protein stabilization by salt bridges: concepts, experimental approaches and clarification of some misunderstandings

Salt bridges in proteins are bonds between oppositely charged residues that are sufficiently close to each other to experience electrostatic attraction. They contribute to protein structure and to the specificity of interaction of proteins with other biomolecules, but in doing so they need not necessarily increase a protein's free energy of unfolding. The net electrostatic free energy of a salt bridge can be partitioned into three components: charge-charge interactions, interactions of charges with permanent dipoles, and desolvation of charges. Energetically favorable Coulombic charge-charge interaction is opposed by often unfavorable desolvation of interacting charges. As a consequence, salt bridges may destabilize the structure of the folded protein. There are two ways to estimate the free energy contribution of salt bridges by experiment: the pK(a) approach and the mutation approach. In the pK(a) approach, the contribution of charges to the free energy of unfolding of a protein is obtained from the change of pK(a) of ionizable groups caused by altered electrostatic interactions upon folding of the protein. The pK(a) approach provides the relative free energy gained or lost when ionizable groups are being charged. In the mutation approach, the coupling free energy between interacting charges is obtained from a double mutant cycle. The coupling free energy is an indirect and approximate measure of the free energy of charge-charge interaction. Neither the pK(a) approach nor the mutation approach can provide the net free energy of a salt bridge. Currently, this is obtained only by computational methods which, however, are often prone to large uncertainties due to simplifying assumptions and insufficient structural information on which calculations are based. This state of affairs makes the precise thermodynamic quantification of salt bridge energies very difficult. This review is focused on concepts and on the assessment of experimental methods and does not cover the vast literature.     

Electrostatic stabilization of a thermophilic cold shock protein.

The cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus differs from its mesophilic homolog Bs-CspB from Bacillus subtilis by 15.8 kJ mol(-1) in the Gibbs free energy of denaturation (DeltaG(D)). The two proteins vary in sequence at 12 positions but only two of them, Arg3 and Leu66 of Bc-Csp, which replace Glu3 and Glu66 of Bs-CspB, are responsible for the additional stability of Bc-Csp. These two positions are near the ends of the protein chain, but close to each other in the three-dimensional structure. The Glu3Arg exchange alone changed the stability by more than 11 kJ mol(-1). Here, we elucidated the molecular origins of the stability difference between the two proteins by a mutational analysis. Electrostatic contributions to stability were characterized by measuring the thermodynamic stabilities of many variants as a function of salt concentration. Double and triple mutant analyses indicate that the stabilization by the Glu3Arg exchange originates from three sources. Improved hydrophobic interactions of the aliphatic moiety of Arg3 contribute about 4 kJ mol(-1). Another 4 kJ mol(-1) is gained from the relief of a pairwise electrostatic repulsion between Glu3 and Glu66, as in the mesophilic protein, and 3 kJ mol(-1) originate from a general electrostatic stabilization by the positive charge of Arg3, which is not caused by a pairwise interaction. Mutations of all potential partners for an ion pair within a radius of 10 A around Arg3 had only marginal effects on stability. The Glu3-->Arg3 charge reversal thus optimizes ionic interactions at the protein surface by both local and global effects. However, it cannot convert the coulombic repulsion with another Glu residue into a corresponding attraction. Avoidance of unfavorable coulombic repulsions is probably a much simpler route to thermostability than the creation of stabilizing surface ion pairs, which can form only at the expense of conformational entropy.     

Basic interdomain boundary residues in calmodulin decrease calcium affinity of sites I and II by stabilizing helix-helix interactions

Calmodulin is an EF-hand calcium-binding protein (148 a.a.) essential in intracellular signal transduction. Its homologous N- and C-terminal domains are separated by a linker that appears disordered in NMR studies. In a study of an N-domain fragment of Paramecium CaM (PCaM1-75), the addition of linker residues 76 to 80 (MKEQD) raised the Tm by 9 degrees C and lowered calcium binding by 0.54 kcal/mol (Sorensen et al., [Biochemistry 2002;41:15-20]), showing that these tether residues affect energetics as well as being a barrier to diffusion. To determine the individual contributions of residues 74 through 80 (RKMKEQD) to stability and calcium affinity, we compared a nested series of 7 fragments (PCaM1-74 to PCaM1-80). For the first 4, PCaM1-74 through PCaM1-77, single amino acid additions at the C-terminus corresponded to stepwise increases in thermostability and decreases in calcium affinity with a net change of 13.5 degrees C in Tm and 0.55 kcal/mol in free energy. The thermodynamic properties of fragments PCaM1-77 through PCaM1-80 were nearly identical. We concluded that the 3 basic residues in the sequence from 74 to 77 (RKMK) are critical to the increased stability and decreased calcium affinity of the longer N-domain fragments. Comparisons of NMR (HSQC) spectra of 15N-PCaM1-74 and 15N-PCaM1-80 and analysis of high-resolution structural models suggest these residues are latched to amino acids in helix A of CaM. The addition of residues E78, Q79, and D80 had a minimal effect on sites I and II, but they may contribute to the mechanism of energetic communication between the domains.     

Human platelet factor 4 subunit association/dissociation thermodynamics and kinetics

Platelet factor 4 (PF4) monomers (7800 daltons) form dimers and tetramers in varying molar ratios under certain solution conditions [Mayo, K. H., & Chen, M. J. (1989) Biochemistry 28, 9469]. The presence of a simplified aromatic region (one Tyr and two His) and resolved monomer, dimer, and tetramer Y60 3,5 ring proton resonances makes study of PF4 aggregate association/dissociation thermodynamics and kinetics possible. PF4 protein subunit association/dissociation equilibrium thermodynamic parameters have been derived by 1H NMR (500MHz) resonance line-fitting analysis of steady-state Y60 3,5 ring proton resonance monomer-dimer-tetramer populations as a function of temperature from 10 to 40 degrees C. Below 10 degrees C and above 40 degrees C, resonance broadening and overlap severely impaired analysis. Enthalpic and entropic contributions to dimer association Gibb's free energy [-5.1 kcal/mol (30 degrees C)] are +2.5 +/- 1 kcal/mol and +26 +/- 7 eu, respectively, and for tetramer association Gibb's free energy [-5.7 kcal/mol (30 degrees C)], they are -7.5 +/- 1 kcal/mol and -7 +/- 3 eu, respectively. These thermodynamic parameters are consistent with low dielectric medium electrostatic/hydrophobic interactions governing dimer formation and hydrogen bonding governing tetramer formation. Association/dissociation kinetic parameters, i.e., steady-state jump rates, have been derived from exchange-induced line-width increases and from 1H NMR (500 MHz) saturation-transfer and spin-lattice (Tl) relaxation experiments. From dissociation jump rates and equilibrium constants, association rate constants were estimated. For dimer and tetramer equilibria at 30 degrees C, unimolecular dissociation rate constants are 35 +/- 10 s-1 for dimer dissociation and 6 +/- 2 s-1 for tetramer dissociation.(ABSTRACT TRUNCATED AT 250 WORDS)     

Mutations in alpha-helical solvent-exposed sites of eglin c have long-range effects: evidence from molecular dynamics simulations

Eglin c is a small protease inhibitor whose structural and thermodynamic properties have been well studied. Previous thermodynamic measurements on mutants at solvent-accessible positions in the protein's helix have shown the unexpected result that the data could be best fit by the inclusion of residue- and position-specific parameters to the model. To explore the origins of this surprising result, long molecular dynamics simulations in explicit solvent have been performed. These simulations indicate specific long-range interactions between the solvent-exposed residues in the eglin c alpha-helix and binding loop, an unexpected observation for such a small protein. The residues involved in the interaction are on opposite sides of the protein, about 25 A apart. Simulations of alanine substitutions at the solvent-exposed helix positions, arginine 22, glutamic acid 23, threonine 26, and leucine 27, show both small and large perturbations of eglin c dynamics. Two mutations exhibit large impacts on the long-range helix-loop interactions. Previous stability measurements (Yi et al., Biochemistry 2003;42:7594-7603) had indicated that an alanine substitution at position 27 was less stabilizing than at other solvent-exposed positions in the helix. The L27A mutation effects observed in these simulations suggest that the position-dependent loss of stability measured in wet bench experiments is derived from changes in dynamics that involve long-range interactions; thus, these simulations support the hypothesis that solvent-exposed positions in helices are not always equivalent.