Calorimetric study of the order-disorder conformational transition in succinoglycan

Thermally induced order-disorder conformational transition in succinoglycan was studied using the method of high-sensitivity differential scanning microcalorimetry within the range of polysaccharide concentrations from 0.1 to 3.5 mg mL⁻¹ at NaCl concentrations 0, 0.01, and 0.1M. The positions and shapes of the excess heat capacity curves depended substantially on both the NaCl and polysaccharide concentrations. At low polysaccharide concentrations in salt-free solution the experimental curves were closely approximated by the two-state model suggesting the transition mechanism to be of the single helix-coil type. With increasing polysaccharide and/or NaCl concentration, the experimental curves changed significantly in symmetry, which indicated a changing transition mechanism. At high polysaccharide concentrations or in the presence of the salt, the order-disorder transition of succinoglycan was shown to include two stages: the cooperative dissociation of the helix dimer and subsequent two-state melting of the helix monomer. The dependence of thermodynamic parameters for the dissociation and melting of helix structures in succinoglycan on NaCl and polysaccharide concentrations was obtained by fitting the experimental excess heat capacity curves. The cooperativity parameter σ for the single helix-coil transition as well as the average length of the helix segment of succinoglycan were calculated. Some features of succinoglycan ordering in solution are discussed. © 1996 John Wiley & Sons, Inc.     

The nature of stacking interations in polynucleotides. Molecular states in Oligo- and polyribocytidylic acids by relaxation analysis.

The dynamics of the helix-coil transition of single-stranded poly(C) (polyribocytidylate) and CpC (cytidyly(3'-5')cytosine) was investigated by an improved cable temperature-jump technique. The single-strand relaxation was characterized by following the ultraviolet (uv) absorbance changes at 248 and 280 nm. Poly(C) and CpC showed single relaxation processes with amplitudes corresponding to those expected from equilibrium melting curves. The relaxation time contants in the range of 25-100 ns were independent of the nucleotide concentration, but strongly dependent upon temperature. Using thermodynanic parameters obtained from circular dichroism (CD) and uv absorbance melting curves, the following rate constants k (at 20 degrees C, 1.05 M ionic strength, pH 7) and activation enthalpies EA were calculated for poly (C): helix formation kR = 1.11 X 10(-7) s-1 (EAR = 2.6 kcal); helix dissociation kD = 2.1 X 10(6) s-1 (EAD = 11.9 kcal). The rate constants obtained for CpC were higher by a factor of about 2 in kR and 12 in kD, whereas the activation enthalpies closely corresponded to those found for the polymer. In addition to the single-stranded helix-coil relaxation, poly(C) and CpC exhibit a relaxation process with a time constant below 25 ns and maximum amplitudes at wavelengths lambda greater than or equal to 285 nm. The same process is found in cytidine and is attributed to hydration equilibria. The hydration reaction can be considered to be in equilibrium during the entire time range of the helix-coil transition and thus the data obtained for the helix-coil transition can be described by a simple two-state model. The rate parameters indicate the existence of relatively high energy barriers in the helix-coil transition and provide strong evidence evidence against an oscillating dimer model. If there is an ensemble of substates for one of the states (as may be expected for the coil form), the energy difference between the populated substates is small compared with the energy difference between the major conformational states.     

Predicting helical segments in proteins by a helix-coil transition theory with parameters derived from a structural database of proteins

A novel helix-coil transition theory has been developed. This new theory contains more types of interactions than similar theories developed earlier. The parameters of the models were obtained from a database of 351 nonhomologous proteins. No manual adjustment of the parameters was performed. The interaction parameters obtained in this manner were found to be physically meaningful, consistent with current understanding of helix stabilizing/destabilizing interactions. Novel insights into helix stabilizing/destabilizing interactions have also emerged from this analysis. The theory developed here worked well in sorting out helical residues from amino acid sequences. If the theory was forced to make prediction on every residue of a given amino acid sequence, its performance was the best among ten other secondary structural prediction algorithms in distinguishing helical residues from nonhelical ones. The theory worked even better if one only required it to make prediction on residues that were “predictable” (identifiable by the theory); >90% predictive reliability could be achieved. The helical residues or segments identified by the helix-coil transition theory can be used as secondary structural contraints to speed up the prediction of the three-dimensional structure of a protein by reducing the dimension of a computational protein folding problem. Possible further improvements of this helix-coil transition theory are also discussed. Proteins 28:344–359, 1997. © 1997 Wiley-Liss, Inc.     

Evaluating Molecular Mechanical Potentials for Helical Peptides and Proteins

Multiple variants of the AMBER all-atom force field were quantitatively evaluated with respect to their ability to accurately characterize helix-coil equilibria in explicit solvent simulations. Using a global distributed computing network, absolute conformational convergence was achieved for large ensembles of the capped A₂₁ and F_s helical peptides. Further assessment of these AMBER variants was conducted via simulations of a flexible 164-residue five-helix-bundle protein, apolipophorin-III, on the 100 ns timescale. Of the contemporary potentials that had not been assessed previously, the AMBER-99SB force field showed significant helix-destabilizing tendencies, with beta bridge formation occurring in helical peptides, and unfolding of apolipophorin-III occurring on the tens of nanoseconds timescale. The AMBER-03 force field, while showing adequate helical propensities for both peptides and stabilizing apolipophorin-III, (i) predicts an unexpected decrease in helicity with ALA→ARG⁺ substitution, (ii) lacks experimentally observed 3₁₀ helical content, and (iii) deviates strongly from average apolipophorin-III NMR structural properties. As is observed for AMBER-99SB, AMBER-03 significantly overweighs the contribution of extended and polyproline backbone configurations to the conformational equilibrium. In contrast, the AMBER-99φ force field, which was previously shown to best reproduce experimental measurements of the helix-coil transition in model helical peptides, adequately stabilizes apolipophorin-III and yields both an average gyration radius and polar solvent exposed surface area that are in excellent agreement with the NMR ensemble.     

The elastic rod model for DNA and its application to the tertiary structure of DNA minicircles in mononucleosomes.

Explicit solutions to the equations of equilibrium in the theory of the elastic rod model for DNA are employed to develop a procedure for finding the configuration that minimizes the elastic energy of a minicircle in a mononucleosome with specified values of the minicircle size N in base pairs, the extent w of wrapping of DNA about the histone core particle, the helical repeat h(0)b of the bound DNA, and the linking number Lk of the minicircle. The procedure permits a determination of the set Y(N, w, h(0)b) of integral values of Lk for which the minimum energy configuration does not involve self-contact, and graphs of writhe versus w are presented for such values of Lk. For the range of N of interest here, 330 < N < 370, the set Y(N, w, h(0)b) is of primary importance: when Lk is not in Y(N, w, h(0)b), the configurations compatible with Lk have elastic energies high enough to preclude the occurrence of an observable concentration of topoisomer Lk in an equilibrium distribution of topoisomers. Equilibrium distributions of Lk, calculated by setting differences in the free energy of the extranucleosomal loop equal to differences in equilibrium elastic energy, are found to be very close to Gaussian when computed under the assumption that w is fixed, but far from Gaussian when it is assumed that w fluctuates between two values. The theoretical results given suggest a method by which one may calculate DNA-histone binding energies from measured equilibrium distributions of Lk.     

Long α helices projecting from the membrane as the dimer interface in the voltage-gated H+ channel

The voltage-gated H⁺ channel (Hv) is a H⁺-permeable voltage-sensor domain (VSD) protein that consists of four transmembrane segments (S1–S4). Hv assembles as a dimeric channel and two transmembrane channel domains function cooperatively, which is mediated by the coiled-coil assembly domain in the cytoplasmic C terminus. However, the structural basis of the interdomain interactions remains unknown. Here, we provide a picture of the dimer configuration based on the analyses of interactions among two VSDs and a coiled-coil domain. Systematic mutations of the linker region between S4 of VSD and the coiled-coil showed that the channel gating was altered in the helical periodicity with the linker length, suggesting that two domains are linked by helices. Cross-linking analyses revealed that the two S4 helices were situated closely in the dimeric channel. The interaction interface between the two S4 and the assembly interface of the coiled-coil domain were aligned in the same direction based on the phase angle calculation along α helices. Collectively, we propose that continuous helices stretching from the transmembrane to the cytoplasmic region in the dimeric interface regulate the channel activation in the Hv dimer.     

Monomer concentrations and dimerization constants in crystallizing lysozyme solutions by dialysis kinetics.

Dialysis kinetics measurements have been made to study the effect of ionic strength on the dimerization of lysozyme in acidic solutions that lead to the growth of tetragonal lysozyme crystals. Using glutaraldehyde cross-linked dimers of lysozyme, we have determined that both monomers and dimers can escape from 25,000 molecular weight cutoff dialysis membranes with velocity constants of 5.1 x 10(-7) and 1.0 x 10(-7) s(-1) for the monomer and dimer species, respectively. The flux from 25K MWCO membranes has been measured for lysozyme in pH 4.0 buffered solutions of 1, 3, 4, 5, and 7% NaCl over a wide range of protein concentrations. Assuming that dimerization is the first step in crystallization, a simple monomer to dimer equilibrium was used to model the flux rates. Dimerization constants calculated at low protein concentrations were 265, 750, 1212, and 7879 M(-1) for 3, 4, 5, and 7% NaCl, respectively. These values indicate that dimerization increases with the ionic strength of the solution suggesting that aggregation is moderated by electrostatic interactions. At high protein concentrations and high supersaturation, the dimerization model does not describe the data well. However, the Li model that uses a pathway of monomer <-> dimer <-> tetramer <-> octamer <-> 16-mer fits the measured flux data remarkably well suggesting the presence of higher order aggregates in crystallizing solutions.     

Contractility and conformation

Contractility in fibers can arise from changes of macromolecular conformation caused by changes in some thermodynamic variable such as temperature, pH, or solvent composition. Illustrations are given of contractile processes in fibers and of changes in macromolecular conformation in dilute solution. These may involve order-disorder transitions, e.g. of the type exhibited by the helix-coil transition. A statistical mechanical treatment of the helix-coil transition involves the assignment of statistical weights to various states and the proper counting of these statistical weights in the formation and evaluation of the partition function; the thermodynamic properties of the system are derivable from the partition function. The counting procedure involved in the consideration of the α-helix and random coil is described. In addition, the factors affecting the relative stabilities of various helical conformations are discussed. These considerations of macromolecular conformation provide a basis for discussing contractile mechanisms in which changes of conformation are involved.     

A thermodynamic scale for leucine zipper stability and dimerization specificity: e and g interhelical interactions.

The leucine zipper is a dimeric coiled-coil protein structure composed of two amphipathic alpha-helices with the hydrophobic surfaces interacting to create the dimer interface. This structure has been found to mediate the dimerization of two abundant classes of DNA binding proteins: the bZIP and bHLH-Zip proteins. Several workers have reported that amino acids in the e and g positions of the coiled coil can modulate dimerization stability and specificity. Using the bZIP protein VBP as a host molecule, we report a thermodynamic scale (delta delta G) for 27 interhelical interactions in 35 proteins between amino acids in the g and the following e positions (g<==>e') of a leucine zipper coiled coil. We have examined the four commonly occurring amino acids in the e and g positions of bZIP proteins, lysine (K), arginine (R), glutamine (Q), glutamic acid (E), as well as the only other remaining charged amino acid aspartic acid (D), and finally alanine (A) as a reference amino acid. These results indicate that E<==>R is the most stable interhelical pair, being 0.35 kcal/mol more stable than E<==>K. A thermodynamic cycle analysis shows that the E<==>R pair is 1.33 kcal/mol more stable than A<==>A with -1.14 kcal/mol of coupling energy (delta delta Gint) coming from the interaction of E with R. The E<==>K coupling energy is only -0.14 kcal/mol. E interacts with more specificity than Q. The R<==>R pair is less stable than the K<==>K by 0.24 kcal/mol. R interacts with more specificity than K. Q forms more stable pairs with the basic amino acids K and R rather than with E. Changing amino acids in the e position to A creates bZIP proteins that form tetramers.     

Detection and characterization of the intermediate on the folding pathway of human alpha-lactalbumin

To discuss the relation between the folding mechanism and the chemical structure of proteins, the reversible unfolding reactions of human alpha-lactalbumin by acidification and by guanidine hydrochloride at 25 degrees C are studied by means of circular dichroism, difference spectra and pH-jump measurements and are compared with those for bovine alpha-lactalbumin. As shown previously for bovine alpha-lactalbumin, the folding process at neutral pH is not explained by a simple two-state mechanism but involves an intermediate form that has the same amount of helical structures as the native form. The transition between the intermediate and the fully denatured states is too rapid to be measured and corresponds to the helix-coil transition of the backbone. One of the differences of human alpha-lactalbumin from the bovine protein is the remarkable stability of the intermediate at neutral pH, which can be explained by differences in the primary chemical structure. Another difference is the existence at acid pH of an additional helical form, which is more helical than the native form. The transition from this to the intermediate or to the fully denatured one also is shown to resemble the helix-coil transition. The following folding scheme of human alpha-lactalbumin is proposed: formula: (see text). Here N is the native form, and the intermediate is a macroscopic state distributed around the state A3 at neutral pH, while the distribution in the acid and fully denautured states shifts toward Am and A-n, respectively.     

Effect of sequence-specific interactions on the stability of helical conformations in polypeptides

A modification of the Zimm–Bragg two-state model for the helix–coil transition in polypeptides, which considers the effect of charge–dipole, charge–charge, and other specific interactions on helix stability, is presented. The new model introduces a series of adjustable parameters whose values are estimated by fitting to recent spectroscopic results on medium-sized peptides. This formalism, based on traditional two-state helix–coil transition models, provides a framework in which data on the helix contents of peptides of specific sequence can be rationalized by a statistical mechanical theory.     

Parameters of helix-coil transition theory for alanine-based peptides of varying chain lengths in water.

Thermal unfolding curves have been measured for a series of short alanine-based peptides that contain repeating sequences and varying chain lengths. Standard helix-coil theory successfully fits the observed transition curves, even for these short peptides. The results provide values for sigma, the helix nucleation constant, delta H0, the enthalpy change on helix formation, and for s (0 degree C), the average helix propagation parameter at 0 degree C. The enthalpy change agrees with the value determined calorimetrically. The success of helix-coil theory in describing the unfolding transitions of short peptides in water indicates that helical propensities, or s values, can be determined from substitution experiments in short alanine-based peptides.     

Spectral Analysis of Naturally Occurring Methylxanthines (Theophylline,Theobromine and Caffeine) Binding with DNA

Nucleic acids exist in a dynamic equilibrium with a number of molecules that constantly interact with them and regulate the cellular activities. The inherent nature of the structure and conformational integrity of these macromolecules can lead to altered biological activity through proper targeting of nucleic acids binding ligands or drug molecules. We studied the interaction of naturally occurring methylxanthines such as theophylline, theobromine and caffeine with DNA, using UV absorption and Fourier transform infrared (FTIR) spectroscopic methods, and especially monitored their binding affinity in the presence of Mg²⁺ and during helix-coil transitions of DNA by temperature (T_m) or pH melting profiles. The study indicates that all these molecules effectively bind to DNA in a dose dependent manner. The overall binding constants of DNA-theophylline = 3.5×10³ M⁻¹, DNA-theobromine = 1.1×10³ M⁻¹, and DNA-Caffeine = 3.8×10³ M⁻¹. On the other hand T_m/pH melting profiles showed 24–35% of enhanced binding activity of methylxanthines during helix-coil transitions of DNA rather than to its native double helical structure. The FTIR analysis divulged that theophylline, theobromine and caffeine interact with all the base pairs of DNA (A-T; G-C) and phosphate group through hydrogen bond (H-bond) interaction. In the presence of Mg²⁺, methylxanthines altered the structure of DNA from B to A-family. However, the B-family structure of DNA remained unaltered in DNA-methylxanthines complexes or in the absence of Mg²⁺. The spectral analyses indicated the order of binding affinity as “caffeine≥theophylline>theobromine” to the native double helical DNA, and “theophylline≥theobromine>caffeine to the denatured form of DNA and in the presence of divalent metal ions.     

Chemical denaturation: potential impact of undetected intermediates in the free energy of unfolding and m-values obtained from a two-state assumption.

The chemical unfolding transition of a protein was simulated, including the presence of an intermediate (I) in equilibrium with the native (N) and unfolded (U) states. The calculations included free energies of unfolding, DeltaGuw, in the range of 1.4 kcal/mol to 10 kcal/mol and three different global m-values. The simulations included a broad range of equilibrium constants for the N left arrow over right arrow I process. The dependence of the N <--> I equilibrium on the concentration of denaturant was also included in the simulations. Apparent DeltaGuw and m-values were obtained from the simulated unfolding transitions by fitting the data to a two-state unfolding process. The potential errors were calculated for two typical experimental situations: 1) the unfolding is monitored by a physical property that does not distinguish between native and intermediate states (case I), and 2) the physical property does not distinguish between intermediate and unfolded states (case II). The results obtained indicated that in the presence of an intermediate, and in both experimental situations, the free energy of unfolding and the m-values could be largely underestimated. The errors in DeltaGuw and m-values do not depend on the m-values that characterize the global N <--> U transition. They are dependent on the equilibrium constant for the N <--> I transition and its characteristic m1-value. The extent of the underestimation increases for higher energies of unfolding. Including no random error in the simulations, it was estimated that the underestimation in DeltaGuw could range between 25% and 35% for unfolding transitions of 3-10 kcal/mol (case I). In case II, the underestimation in DeltaGuw could be even larger than in case I. In the same energy range, a 50% error in the m-value could also take place. The fact that most of the mutant proteins are characterized by both a lower m-value and a lower stability than the wild-type protein suggests that in some cases the results could have been underestimated due to the application of the two-state assumption.     

Binding of alkaline cations to the double-helical form of gramicidin.

Gramicidin is a polypeptide antibiotic that forms monovalent cation-specific channels in membrane environments. In organic solvents and in lipids containing unsaturated fatty acid chains, it forms a double-helical "pore" structure, in which two monomers are intertwined. This form of gramicidin can bind two cations inside its lumen, and the crystal structures of both an ion complex and an ion-free form have been determined. In this study, we have used circular dichroism (CD) spectroscopy to examine the binding mechanism and the binding constants (K1 and K2) of cations to gramicidin in the double helical form in methanol solution. The dramatic change in optical rotation in the far-ultraviolet CD spectrum of gramicidin provides a useful tool for monitoring the binding. The binding mechanism appears to involve a large conformation change associated with the binding of ions to the first of the two sites. The calculated values for the K1 binding constants for alkaline cations are considerably smaller than the K2 binding constants. The order of binding affinity for alkaline cations is similar to that for the helical dimer "channel" form of gramicidin, i.e., Cs+ approximately Rb+ > > K+ > Li+, but in comparison to the helical dimer form, the binding to double-helical dimers is dominated by a cation size-dependent conformational change in the gramicidin structure.     

A calorimetric characterization of the salt dependence of the stability of the GCN4 leucine zipper.

The effects of different salts (LiCl, NaCl, ChoCl, KF, KCl, and KBr) on the structural stability of a 33-residue peptide corresponding to the leucine zipper region of GCN4 have been studied by high-sensitivity differential scanning calorimetry. These experiments have allowed an estimation of the salt dependence of the thermodynamic parameters that define the stability of the coiled coil. Independent of the nature of the salt, a destabilization of the coiled coil is always observed upon increasing salt concentration up to a maximum of approximately 0.5 M, depending on the specific cation or anion. At higher salt concentrations, this effect is reversed and a stabilization of the leucine zipper is observed. The effect of salt concentration is primarily entropic, judging from the lack of a significant salt dependence of the transition enthalpy. The salt dependence of the stability of the peptide is complex, suggesting the presence of specific salt effects at high salt concentrations in addition to the nonspecific electrostatic effects that are prevalent at lower salt concentrations. The data is consistent with the existence of specific interactions between anions and peptide with an affinity that follows a reverse size order (F- > Cl- > Br-). Under all conditions studied, the coiled coil undergoes reversible thermal unfolding that can be well represented by a reaction of the form N2<==>2U, indicating that the unfolding is a two-state process in which the helices are only stable when they are in the coiled coil conformation.     

Pathway and Mechanism of pH Dependent Human Hemoglobin Tetramer-Dimer-Monomer Dissociations

Hemoglobin dissociation is of great interest in protein process and clinical medicine as well as in artificial blood research. However, the pathway and mechanisms of pH-dependent human Hb dissociation are not clear, whether Hb would really dissociate into monomers is still a question. Therefore, we have conducted a multi-technique investigation on the structure and function of human Hb versus pH. Here we demonstrate that tetramer hemoglobin can easily dissociate into dimer in abnormal pH and the tetramer → dimer dissociation is reversible if pH returns to normal physiological value. When the environmental pH becomes more acidic (<6.5) or alkaline (>8.0), Hb can further dissociate from dimer to monomer. The proportion of monomers increases while the fraction of dimers decreases as pH declines from 6.2 to 5.4. The dimer → monomer dissociation is accompanied with series changes of protein structure thus it is an irreversible process. The structural changes in the dissociated Hbs result in some loss of their functions. Both the Hb dimer and monomer cannot adequately carry and release oxygen to the tissues in circulation. These findings provide a comprehensive understanding on the pH-dependent protein transitions of human Hb, give guideline to explain complex protein processes and the means to control protein dissociation or re-association reaction. They are also of practical value in clinical medicine, blood preservation and blood substitute development.     

How p53 binds DNA as a tetramer.

The p53 tumor suppressor protein is a tetramer that binds sequence-specifically to a DNA consensus sequence consisting of two consecutive half-sites, with each half-site being formed by two head-to-head quarter-sites (--><-- --><--). Each p53 subunit binds to one quarter-site, resulting in all four DNA quarter-sites being occupied by one p53 tetramer. The tetramerization domain forms a symmetric dimer of dimers, and two contrasting models have the two DNA-binding domains of each dimer bound to either consecutive or alternating quarter-sites. We show here that the two monomers within a dimer bind to a half-site (two consecutive quarter-sites), but not to separated (alternating) quarter-sites. Tetramers bind similarly, with the two dimers within each tetramer binding to pairs of half-sites. Although one dimer within the tetramer is sufficient for binding to one half-site in DNA, concurrent interaction of the second dimer with a second half-site in DNA drastically enhances binding affinity (at least 50-fold). This cooperative dimer-dimer interaction occurs independently of tetramerization and is a primary mechanism responsible for the stabilization of p53 DNA binding. Based on these findings, we present a model of p53 binding to the consensus sequence, with the tetramer binding DNA as a pair of clamps.     

A General Framework for Thermodynamically Consistent Parameterization and Efficient Sampling of Enzymatic Reactions

Kinetic models provide the means to understand and predict the dynamic behaviour of enzymes upon different perturbations. Despite their obvious advantages, classical parameterizations require large amounts of data to fit their parameters. Particularly, enzymes displaying complex reaction and regulatory (allosteric) mechanisms require a great number of parameters and are therefore often represented by approximate formulae, thereby facilitating the fitting but ignoring many real kinetic behaviours. Here, we show that full exploration of the plausible kinetic space for any enzyme can be achieved using sampling strategies provided a thermodynamically feasible parameterization is used. To this end, we developed a General Reaction Assembly and Sampling Platform (GRASP) capable of consistently parameterizing and sampling accurate kinetic models using minimal reference data. The former integrates the generalized MWC model and the elementary reaction formalism. By formulating the appropriate thermodynamic constraints, our framework enables parameterization of any oligomeric enzyme kinetics without sacrificing complexity or using simplifying assumptions. This thermodynamically safe parameterization relies on the definition of a reference state upon which feasible parameter sets can be efficiently sampled. Uniform sampling of the kinetics space enabled dissecting enzyme catalysis and revealing the impact of thermodynamics on reaction kinetics. Our analysis distinguished three reaction elasticity regions for common biochemical reactions: a steep linear region (0> ΔG_r >-2 kJ/mol), a transition region (-2> ΔG_r >-20 kJ/mol) and a constant elasticity region (ΔG_r <-20 kJ/mol). We also applied this framework to model more complex kinetic behaviours such as the monomeric cooperativity of the mammalian glucokinase and the ultrasensitive response of the phosphoenolpyruvate carboxylase of Escherichia coli. In both cases, our approach described appropriately not only the kinetic behaviour of these enzymes, but it also provided insights about the particular features underpinning the observed kinetics. Overall, this framework will enable systematic parameterization and sampling of enzymatic reactions.