Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “THERMODYNAMICS”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10Linked to original sources
PubMed1999

Kinetic and thermodynamic aspects of lipid translocation in biological membranes.

A theoretical analysis of the lipid translocation in cellular bilayer membranes is presented. We focus on an integrative model of active and passive transport processes determining the asymmetrical distribution of the major lipid components between the monolayers. The active translocation of the aminophospholipids phosphatidylserine and phosphatidylethanolamine is mathematically described by kinetic equations resulting from a realistic ATP-dependent transport mechanism. Concerning the passive transport of the aminophospholipids as well as of phosphatidylcholine, sphingomyelin, and cholesterol, two different approaches are used. The first treatment makes use of thermodynamic flux-force relationships. Relevant forces are transversal concentration differences of the lipids as well as differences in the mechanical states of the monolayers due to lateral compressions. Both forces, originating primarily from the operation of an aminophospholipid translocase, are expressed as functions of the lipid compositions of the two monolayers. In the case of mechanical forces, lipid-specific parameters such as different molecular surface areas and compression force constants are taken into account. Using invariance principles, it is shown how the phenomenological coefficients depend on the total lipid amounts. In a second approach, passive transport is analyzed in terms of kinetic mechanisms of carrier-mediated translocation, where mechanical effects are incorporated into the translocation rate constants. The thermodynamic as well as the kinetic approach are applied to simulate the time-dependent redistribution of the lipid components in human red blood cells. In the thermodynamic model the steady-state asymmetrical lipid distribution of erythrocyte membranes is simulated well under certain parameter restrictions: 1) the time scales of uncoupled passive transbilayer movement must be different among the lipid species; 2) positive cross-couplings of the passive lipid fluxes are needed, which, however, may be chosen lipid-unspecifically. A comparison of the thermodynamic and the kinetic approaches reveals that antiport mechanisms for passive lipid movements may be excluded. Simulations with kinetic symport mechanisms are in qualitative agreement with experimental data but show discrepancies in the asymmetrical distribution for sphingomyelin.

Biological Transport↗
PubMed1998

On the use of deuterated phospholipids for infrared spectroscopic studies of monomolecular films: a thermodynamic analysis of single and binary component phospholipid monolayers.

A thermodynamic study of monolayer mixing behavior was performed using binary lipid mixtures of DPPC + DOPG (7:1 mol:mol) as well as a mixture containing the phosphocholine acyl chain perdeuterated analog of DPPC, DPPC-d62 (i.e. 7:1 mol:mol DPPC-d62 + DOPG). An analysis was performed on these isotherms that calculated the Helmholtz excess free energy of mixing to determine potential thermodynamic differences in the fully protiated DPPC:DOPG monolayer versus the identical monolayer substituted with acyl chain perdeuterated DPPC (i.e. DPPC-d62). The main conclusions of this study are that: (1) the temperature-dependent thermodynamic phase behavior of DPPC-d62 differs substantially from that of DPPC, and these differences appear magnified when these molecules are studied as monomolecular films at the A/W interface as opposed to bulk phase systems; and (2) the calculated excess free energy of mixing values, deltaFxs(A), for the deuterium-containing 7:1 DPPC-d62:DOPG binary monolayer film show a type of near-ideal mixing behavior above 20 degrees C that is consistent with the 'squeezing-out' of the DOPG component in the binary mixture containing the DPPC-d62 component. This type of behavior does not occur in the graph of the excess free energy of mixing for the fully protiated 7:1 DPPC:DOPG binary monolayer film. The use of acyl chain perdeuterated DPPC (i.e. DPPC-d62) in binary monomolecular mixtures with other phospholipids produces a real and measurable difference on the thermodynamic properties of the monolayer when compared to the case of the fully protiated monomolecular film. In particular, the presence of DPPC-d62 in a 7:1 mol:mol DPPC-d62:DOPG binary monomolecular film may over-state the fluid nature of the monolayer at any given temperature, thereby leading to an over-estimate of the amount of material potentially 'squeezed out' of the monomolecular film. These results have implications for a commonly used IR spectroscopic method that relies on the incorporation of deuterium-labeled phospholipids into monolayers for vibrational spectroscopic analysis of individual components in a multicomponent monomolecular film.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed2003

Solvation effects and driving forces for protein thermodynamic and kinetic cooperativity: how adequate is native-centric topological modeling?

What energetic and solvation effects underlie the remarkable two-state thermodynamics and folding/unfolding kinetics of small single-domain proteins? To address this question, we investigate the folding and unfolding of a hierarchy of continuum Langevin dynamics models of chymotrypsin inhibitor 2. We find that residue-based additive Gō-like contact energies, although native-centric, are by themselves insufficient for protein-like calorimetric two-state cooperativity. Further native biases by local conformational preferences are necessary for protein-like thermodynamics. Kinetically, however, even models with both contact and local native-centric energies do not produce simple two-state chevron plots. Thus a model protein's thermodynamic cooperativity is not sufficient for simple two-state kinetics. The models tested appear to have increasing internal friction with increasing native stability, leading to chevron rollovers that typify kinetics that are commonly referred to as non-two-state. The free energy profiles of these models are found to be sensitive to the choice of native contacts and the presumed spatial ranges of the contact interactions. Motivated by explicit-water considerations, we explore recent treatments of solvent granularity that incorporate desolvation free energy barriers into effective implicit-solvent intraprotein interactions. This additional feature reduces both folding and unfolding rates vis-à-vis that of the corresponding models without desolvation barriers, but the kinetics remain non-two-state. Taken together, our observations suggest that interaction mechanisms more intricate than simple Gō-like constructs and pairwise additive solvation-like contributions are needed to rationalize some of the most basic generic protein properties. Therefore, as experimental constraints on protein chain models, requiring a consistent account of protein-like thermodynamic and kinetic cooperativity can be more stringent and productive for some applications than simply requiring a model heteropolymer to fold to a target structure.

Calorimetry↗
PubMed2003

Evidence of a thermal unfolding dimeric intermediate for the Escherichia coli histone-like HU proteins: thermodynamics and structure.

The Escherichia coli histone-like HU protein pool is composed of three dimeric forms: two homodimers, EcHUalpha(2) and EcHUbeta(2), and a heterodimer, EcHUalphabeta. The relative abundance of these dimeric forms varies during cell growth and in response to environmental changes, suggesting that each dimer plays different physiological roles. Here, differential scanning calorimetry and circular dichroism (CD) were used to study the thermal stability of the three E.coli HU dimers and show that each of them has its own thermodynamic signature. Unlike the other HU proteins studied so far, which melt through a single step (N(2)<-->2D), this present thermodynamic study shows that the three E.coli dimers melt according to a two-step mechanism (N(2)<-->I(2)<-->2D). The native dimer, N(2), melts partially into a dimeric intermediate, I(2), which in turn yields the unfolded monomers, D. In addition, the crystal structure of the EcHUalpha(2) dimer has been solved. Comparative thermodynamic and structural analysis between EcHUalpha(2) and the HU homodimer from Bacillus stearothermophilus suggests that the E.coli dimer is constituted by two subdomains of different energetic properties. The CD study indicates that the intermediate, I(2), corresponds to an HU dimer having partly lost its alpha-helices. The partially unfolded dimer I(2) is unable to complex with high-affinity, single-stranded break-containing DNA. These structural, thermodynamic and functional results suggest that the N(2)<-->I(2) equilibrium plays a central role in the physiology of E.coli HU. The I(2) molecular species seems to be the EcHUbeta(2) preferential conformation, possibly related to its role in the E.coli cold-shock adaptation. Besides, I(2) might be required in E.coli for the HU chain exchange, which allows the heterodimer formation from homodimers.

Amino Acid Sequence↗
PubMed2000

Analysis of protein self-association under conditions of the thermodynamic non-ideality.

Analysis of protein-protein interactions in highly concentrated solutions requires a consideration of the non-ideality in such solutions which is expressed by the virial coefficients. Different equations are presented to estimate effects of the thermodynamic non-ideality on the macromolecular interaction of self-associating proteins in sedimentation equilibrium experiments. Usually the influence of thermodynamic non-ideal behavior are described by concentration power series. The convergence of such power series is limited at high solute concentration. When expressing the thermodynamic non-ideality by an activity power series this disadvantage can be minimized. The developed centrifuge equations are the basis for a global analysis to estimate equilibrium constants and the corresponding thermodynamic activities of the reactants. Based on fit analysis of synthetic concentration profiles it was established that marked deviations from the expected association constants are observed for proteins with strong association forces between solute molecules. Considerable differences were also observed in weakly interacting systems. This was due to the excluded volume of the protein which is similar in magnitude to the binding constant. For interactions with moderate affinities values extremely close to the true binding values were obtained, as confirmed by experimental results with concanavalin A.

Concanavalin A↗
PubMed2002

Thermodynamics of reactions catalyzed by PABA synthase.

Microcalorimetry and high-performance liquid chromatography (HPLC) have been used to conduct a thermodynamic investigation of reactions catalyzed by PABA synthase, the enzyme located at the first step in the shikimic acid metabolic pathway leading from chorismate to 4-aminobenzoate (PABA). The overall biochemical reaction catalyzed by the PabB and PabC components of PABA synthase is: chorismate(aq)+ammonia(aq)=4-aminobenzoate(aq)+pyruvate(aq)+H(2)O(l). This reaction can be divided into two partial reactions involving the intermediate 4-amino-4-deoxychorismate (ADC): chorismate(aq)+ammonia(aq)=ADC(aq)+H(2)O(l) and ADC(aq)=4-aminobenzoate(aq)+pyruvate(aq). Microcalorimetric measurements were performed on all three of these reactions at a temperature of 298.15 K and pH values in the range 8.72-8.77. Equilibrium measurements were performed on the first partial (ADC synthase) reaction at T=298.15 K and at pH=8.78. The saturation molality of 4-aminobenzoate(cr) in water is (0.00382+/-0.0004) mol kg(-1) at T=298.15 K. The results of the equilibrium and calorimetric measurements were analyzed in terms of a chemical equilibrium model that accounts for the multiplicity of ionic states of the reactants and products. These calculations gave thermodynamic quantities at the temperature 298.15 K and an ionic strength of zero for chemical reference reactions involving specific ionic forms. For the reaction: chorismate(2-)(aq)+NH(4)(+)(aq)=ADC(-)(aq)+H(2)O(l), K=(10.8+/-4.2) and Delta(r)H(m)(o)=-(35+/-15) kJ mol(-1). For the reaction: ADC(-)(aq)=4-aminobenzoate(-)(aq)+pyruvate(-)(aq)+H(+)(aq), Delta(r)H(m)(o)=-(139+/-23) kJ mol(-1). For the reaction: chorismate(2-)(aq)+NH(4)(+)(aq)=4-aminobenzoate(-)(aq)+pyruvate(-)(aq)+H(2)O(l)+H(+)(aq), Delta(r)H(m)(o)=-(174+/-6) kJ mol(-1). Thermodynamic cycle calculations were used to calculate thermodynamic quantities for three additional reactions that utilize L-glutamine rather than ammonia and that are pertinent to this branch point of the shikimic acid pathway. The quantities obtained in this study permit the calculation of the position of equilibrium of these reactions as a function of temperature, pH, and ionic strength. Values of the apparent equilibrium constants and the standard transformed Gibbs energy changes Delta(r)G'(m)(o) under approximately physiological conditions are given.

Calorimetry↗
PubMed2003

Thermodynamic characterization of the interaction of human cyclophilin 18 with cyclosporin A.

Isothermal titration calorimetry (ITC) was used to investigate thermodynamic parameters of the cyclosporin A (CsA)-cyclophilin 18 (hCyp18) association reaction. We have calculated the thermodynamic parameters (enthalpy, entropy, heat capacity, and free energy of binding) of the CsA/hCyp18 complexation. All but two methods described in the literature underestimate the affinity to hCyp18 of CsA. We found that the association constant (1.1.10(8) M(-1) at 10 degrees C) of CsA to hCyp18 is in close agreement with the reciprocal of the reported inhibitory constant of the peptidylprolyl cis/trans isomerase activity of hCyp18. Interpretation of the thermodynamic parameters in buffered solution of water, 30% glycerol and D(2)O leads to the conclusion that the highly specific binding of CsA to hCyp18 is mainly mediated through hydrogen bonding and to a lesser degree through hydrophobic interaction. Furthermore, the pH dependence of the association constant was determined and analyzed according to a single proton linkage model, resulting in a pK(a) value of 5.7 in free hCyp18 and below 4.5 in the CsA complexed form. Titration experiments using different single component buffers possessing different heats of ionization allowed us to estimate that statistically half a proton is transferred upon CsA binding from the binding interface of hCyp18 to the buffer at pH 5.5. No proton transfer was detected at pH 7.5. The thermodynamic results are discussed in relation to the published X-ray and NMR structure of the free and CsA complexed hCyp18.

Algorithms↗
PubMed2000

Thermodynamics of reactions catalyzed by anthranilate synthase.

Microcalorimetry and high performance liquid chromatography have been used to conduct a thermodynamic investigation of reactions catalyzed by anthranilate synthase, the enzyme located at the first step in the biosynthetic pathway leading from chorismate to tryptophan. One of the overall biochemical reactions catalyzed by anthranilate synthase is: chorismate(aq) + ammonia(aq) = anthranilate(aq) + pyruvate(aq) + H2O(l). This reaction can be divided into two partial reactions involving the intermediate 2-amino-4-deoxyisochorismate (ADIC): chorismate(aq) + ammonia(aq) = ADIC(aq) + H2O(l) and ADIC(aq) = anthranilate(aq) + pyruvate(aq). The native anthranilate synthase and a mutant form of it that is deficient in ADIC lyase activity but has ADIC synthase activity were used to study the overall ammonia-dependent reaction and the first of the above two partial reactions, respectively. Microcalorimetric measurements were performed on the overall reaction at a temperature of 298.15 K and pH 7.79. Equilibrium measurements were performed on the first partial (ADIC synthase) reaction at temperatures ranging from 288.15 to 302.65 K, and at pH values from 7.76 to 8.08. The results of the equilibrium and calorimetric measurements were analyzed in terms of a chemical equilibrium model that accounts for the multiplicity of ionic states of the reactants and products. These calculations gave thermodynamic quantities at the temperature 298.15 K and an ionic strength of zero for chemical reference reactions involving specific ionic forms. For the reaction: chorismate2-(aq) + NH4+(aq) = anthranilate-(aq) + pyruvate-(aq) + H+(aq) + H2O(l), delta rHmo = -(116.3 +/- 5.4) kJ mol-1. For the reaction: chorismate2-(aq) + NH4+(aq) = ADIC-(aq) + H2O(l), K = (20.3 +/- 4.5) and delta rHmo = (7.5 +/- 0.6) kJ mol-1. Thermodynamic cycle calculations were used to calculate thermodynamic quantities for three additional reactions that are pertinent to this branch point of the chorismate pathway. The quantities obtained in this study permit the calculation of the position of equilibrium of these reactions as a function of temperature, pH, and ionic strength. Values of the apparent equilibrium constants and the standard transformed Gibbs energy changes delta rG'mo under approximately physiological conditions are given.

Anthranilate Synthase↗
PubMed1993

Sensitive fluorescence-based thermodynamic and kinetic measurements of DNA hybridization in solution.

Kinetic and thermodynamic constants associated with DNA hybridization were determined in solution using fluorescence measurements and complementary fluorophore-labeled oligomers. One oligomer was labeled with a 5'-terminal fluorescein, and the other was labeled with a 3'-terminal rhodamine. The juxtaposition of the two labels in double-stranded complexes results in a strong quenching of the fluorescein emission, thereby providing the means for distinguishing single-stranded DNA from double-stranded DNA. Since measurements were based on fluorescence, DNA denaturation and association could be monitored routinely at strand concentrations 100-1000-fold lower than permitted by absorbance hypochromicity measurements. To determine if fluorescence quenching mirrored base pair formation, temperature profiles of DNA association and dissociation were constructed from both absorbance hypochromicity and fluorescence quenching measurements at a number of different DNA concentrations. Analyses of these profiles using the "all-or-none" model of hybridization provided thermodynamic data which were statistically indistinguishable between the two measurement methods, thus validating the use of fluorescence quenching in thermodynamic studies of oligomers. The effects of fluorophore attachment on the thermodynamic properties of the DNA strands were investigated by analyzing the melting curves of different combinations of unlabeled and labeled complementary oligomers. The presence of both labels was found to stabilize the double-stranded DNA by about -1.5 kcal in delta G degrees 298, primarily due to the fluorescein label. Association and dissociation rate constants were determined by fluorescence measurements at different temperatures, and linear Arrhenius plots were obtained. The fluorescence measurements provided a unique "label dilution" method for measuring dissociation rate constants of oligomers based upon the dynamic association and dissociation of complementary DNA strands at constant temperature. Association rate measurements were simplified since relatively low concentrations of complementary oligomers could be mixed, thereby reducing hybridization rates and eliminating the need for rapid mixing and measurement techniques.

Base Sequence↗
PubMed1993

Thermodynamics of single-stranded RNA and DNA interactions with oligolysines containing tryptophan. Effects of base composition.

We have examined the thermodynamics of binding of a series of oligolysines (net charge z = +2 to +10) containing one, two, or three tryptophans to several single-stranded (ss) homo-polynucleotides [poly(A), poly(C), poly(I), poly(dU), poly(dT)] and duplex (ds) DNA in order to investigate the effects of peptide charge, tryptophan content, and polynucleotide base and sugar type. Equilibrium association constants, Kobs, were measured as a function of monovalent salt concentration (KCH3CO2) and temperature by monitoring the quenching of the peptide tryptophan fluorescence upon interaction with the polynucleotides, from which the dependence of delta G(o)obs, delta H(o)obs, and delta S(o)obs on [KCH3CO2] was obtained. As observed previously with poly(U) [Mascotti, D.P., & Lohman, T.M. (1992) Biochemistry 31, 8932], the dependence of delta G(o)obs on [K+] for peptide binding to each polynucleotide is entirely entropic in origin (i.e., delta H(o)obs is independent of [K+]), consistent with the conclusion that Kobs increases with decreasing salt concentration due to the favorable increase in entropy resulting from the displacement of bound cations (K+) from the nucleic acid upon formation of the complex. For each ss polynucleotide, we find that significantly less than one potassium ion is released thermodynamically per net positive peptide charge, as determined from the value of delta log Kobs/delta log [K+]. Interestingly, (-delta log Kobs/delta log [K+])/z decreases with increasing peptide charge for poly(A), poly(C), and poly(dT), contrary to the behavior observed with poly(U) and ds-DNA, which may reflect a significant release of bound water upon formation of peptide complexes with these ss homo-polynucleotides or an increased binding of K+ to the ss polynucleotide with increasing [K+]. Alternatively, there may be conformational differences between the bound states of oligolysines of low charge, relative to oligolysines of higher charge. However, in all cases, peptides with z < +4 display different thermodynamics of binding than peptides with z > +4. The presence of tryptophan (Trp) within these peptides does not influence the salt dependence of Kobs for binding to poly(A), poly(C), or poly(dT). However, the Trp content of the peptide does contribute significantly to the thermodynamics of these interactions: Trp interactions result in a favorable contribution to delta H(o)obs, but an unfavorable contribution to delta S(o)obs, with little effect on delta G(o)obs due to entropy-enthalpy compensations. Oligolysines containing Trp also display a small, but significant, dependence of Kobs on base composition, with Kobs decreasing in the order poly(I) >> poly(dT) approximately poly(U) approximately poly(A) >> poly(C).

Amino Acid Sequence↗
PubMed1992

Use of liquid hydrocarbon and amide transfer data to estimate contributions to thermodynamic functions of protein folding from the removal of nonpolar and polar surface from water.

This extension of the liquid hydrocarbon model seeks to quantify the thermodynamic contributions to protein stability from the removal of nonpolar and polar surface from water. Thermodynamic data for the transfer of hydrocarbons and organic amides from water to the pure liquid phase are analyzed to obtain contributions to the thermodynamics of folding from the reduction in water-accessible surface area. Although the removal of nonpolar surface makes the dominant contribution to the standard heat capacity change of folding (delta C0fold), here we show that inclusion of the contribution from removal of polar surface allows a quantitative prediction of delta C0fold within the uncertainty of the calorimetrically determined value. Moreover, analysis of the contribution of polar surface area to the enthalpy of transfer of liquid amides provides a means of estimating the contributions from changes in nonpolar and polar surface area as well as other factors to the enthalpy of folding (delta H0fold). In addition to estimates of delta H0fold, this extension of the liquid hydrocarbon model provides a thermodynamic explanation for the observation [Privalov, P. L., & Khechinashvili, N. N. (1974) J. Mol. Biol. 86, 665-684] that the specific enthalpy of folding (cal g-1) of a number of globular proteins converges to a common value at approximately 383 K. Because amounts of nonpolar and polar surface area buried by these proteins upon folding are found to be linear functions of molar mass, estimates of both delta C0fold and delta H0fold may be obtained given only the molar mass of the protein of interest.(ABSTRACT TRUNCATED AT 250 WORDS)

Amides↗
PubMed2001

Thermodynamics of metal ion binding. 1. Metal ion binding by wild-type carbonic anhydrase.

Understanding the energetic consequences of molecular structure in aqueous solution is a prerequisite to the rational design of synthetic motifs with predictable properties. Such properties include ligand binding and the collapse of polymer chains into discrete three-dimensional structures. Despite advances in macromolecular structure determination, correlations of structure with high-resolution thermodynamic data remain limited. Here we compare thermodynamic parameters for the binding of Zn(II), Cu(II), and Co(II) to human carbonic anhydrase II. These calorimetrically determined values are interpreted in terms of high-resolution X-ray crystallographic data. While both zinc and cobalt are bound with a 1:1 stoichiometry, CAII binds two copper ions. Considering only the high-affinity site, there is a diminution in the enthalpy of binding through the series Co(II) --> Zn(II) --> Cu(II) that mirrors the enthalpy of hydration; this observation reinforces the notion that the thermodynamics of solute association with water is at least as important as the thermodynamics of solute-solute interaction and that these effects must be considered when interpreting association in aqueous solution. Additionally, DeltaC(p) data suggest that zinc binding to CAII proceeds with a greater contribution from desolvation than does binding of either copper or cobalt, suggesting Nature optimizes binding by optimizing desolvation.

Apoenzymes↗
PubMed2001

Thermodynamics of metal ion binding. 2. Metal ion binding by carbonic anhydrase variants.

The ability to construct molecular motifs with predictable properties in aqueous solution requires an extensive knowledge of the relationships between structure and energetics. The design of metal binding motifs is currently an area of intense interest in the bioorganic community. To date synthetic motifs designed to bind metal ions lack the remarkable affinities observed in biological systems. To better understand the structural basis of metal ion affinity, we report here the thermodynamics of binding of divalent zinc ions to wild-type and mutant carbonic anhydrases and the interpretation of these parameters in terms of structure. Mutations were made both to the direct His ligand at position 94 and to indirect, or second-shell, ligands Gln-92, Glu-117, and Thr-199. The thermodynamics of ligand binding by several mutant proteins is complicated by the development of a second zinc binding site on mutation; such effects must be considered carefully in the interpretation of thermodynamic data. In all instances modification of the protein produces a complex series of changes in both the enthalpy and entropy of ligand binding. In most cases these effects are most readily rationalized in terms of ligand and protein desolvation, rather than in terms of changes in the direct interactions of ligand and protein. Alteration of second-shell ligands, thought to function primarily by orienting the direct ligands, produces profoundly different effects on the enthalpy of binding, depending on the nature of the residue. These results suggest a range of activities for these ligands, contributing both enthalpic and entropic effects to the overall thermodynamics of binding. Together, our results demonstrate the importance of understanding relationships between structure and hydration in the construction of novel ligands and biological polymers.

Asparagine↗
PubMed1991

A thermodynamic analysis of RNA transcript elongation and termination in Escherichia coli.

In the first part of this paper we present a thermodynamic analysis of the elongation phase of transcription in Escherichia coli. The stability of the elongation complex is described by a "free energy of formation" function (delta G zero f) that is a sum of terms for forming (i) a locally denatured 17-base-pair DNA "bubble"; (ii) a constant-length hybrid between the 3'-terminal 12-nucleotide residues of the RNA transcript and the corresponding region of the DNA template strand; and (iii) a set of binding interactions between the polymerase and certain DNA and RNA residues within and near the "transcription bubble". The transcriptional elongation complex is very stable at most positions along a natural DNA template and moves in a highly processive fashion. At these positions, the delta G zero f function provides a quantitative measure of the stability of the elongation complex. Besides allowing for the polymerization of the RNA transcript, the elongation complex also serves to define the context within which transcript termination occurs. In the second part of the paper the thermodynamic analysis is extended to discriminate between template positions at which the elongation complex is stable and positions at which it is rendered relatively unstable by the presence of a string of rU residues at the 3'-terminus of the RNA together with the formation of a specific RNA hairpin just upstream of this point. Most factor-independent (intrinsic) termination events are thermodynamically disallowed at the former positions and are thermodynamically allowed at the latter positions. The extended form of the analysis closely predicts the exact sites of termination at a number of intrinsic terminators (and attenuators) in the E. coli genome. It also correctly predicts bidirectional function for a number of bidirectional terminators. In some cases it may identify terminators that are similar to the intrinsic type but that require additional protein factors, unusual polymerase-nucleic acid interactions, or rate-limiting conformational changes in order to function. Finally, it successfully locates intrinsic terminators within a number of E. coli operons and discriminates between these terminators and the surrounding DNA sequence.

Algorithms↗
PubMed2002

Flavin thermodynamics explain the oxygen insensitivity of enteric nitroreductases.

Bacterial nitroreductases are NAD(P)H-dependent flavoenzymes which catalyze the oxygen-insensitive reduction of nitroaromatics, quinones, and riboflavin derivatives. Despite their broad substrate specificity, their reactivity is very specific for two-electron, not one-electron, chemistry. We now describe the thermodynamic properties of the flavin mononucleotide cofactor of Enterobacter cloacae nitroreductase (NR), determined under a variety of solution conditions. The two-electron redox midpoint potential of NR is -190 mV at pH 7.0, and both the pH dependence of the midpoint potential and the optical spectrum of the reduced enzyme indicate that the transition is from neutral oxidized flavin to anionic flavin hydroquinone. The one-electron-reduced semiquinone states of both the free enzyme and an NR-substrate analogue complex are strongly suppressed based on optical spectroscopy and electron paramagnetic resonance measurements. This can explain the oxygen insensitivity of NR and its homologues, as it makes the execution of one-electron chemistry thermodynamically unfavorable. Therefore, we have established a chemical basis for the recent finding that a nitroreductase is a member of the soxRS oxidative defense regulon in Escherichia coli [Liochev, S. I., Hausladen, A., Fridovich, I. (1999) Proc. Natl. Acad. Sci. U.S.A. 96 (7), 3537-3539]. We also report binding affinities for the FMN cofactor in all three oxidation states either determined fluorometrically or calculated using thermodynamic cycles. Thus, we provide a detailed picture of the thermodynamics underlying the unusual activity of NR.

Apoenzymes↗
PubMed2004

Thermodynamic parameters based on a nearest-neighbor model for DNA sequences with a single-bulge loop.

All 64 possible thermodynamic parameters for a single-bulge loop in the middle of a sequence were derived from optical melting studies. The relative stability of a single bulge depended on both the type of bulged base and its flanking base pairs. The contribution of the single bulge to helix stability ranged from 3.69 kcal/mol for a TAT bulge to -1.05 kcal/mol for an ACC bulge. Thermodynamics for 10 sequences with a GTG bulge were determined to test the applicability of the nearest-neighbor model to a single-bulge loop. Thermodynamic parameters for the GTG bulge and Watson-Crick base pairs predict, DeltaH degrees, DeltaS degrees, and T(M)(50 microM) values with average deviations of 3.0%, 4.3%, 4.7%, and 0.9 degrees C, respectively. The prediction accuracy was within the limits of what can be expected for a nearest-neighbor model. This certified that the thermodynamics for single-bulge loops can be estimated adequately using a nearest-neighbor model.

DNA↗
PubMed2004

Positive selection dictates the choice between kinetic and thermodynamic protein folding and stability in subtilases.

Subtilisin E (SbtE) is a member of the ubiquitous superfamily of serine proteases called subtilases and serves as a model for understanding propeptide-mediated protein folding mechanisms. Unlike most proteins that adopt thermodynamically stable conformations, the native state of SbtE is trapped into a kinetically stable conformation. While kinetic stability offers distinct functional advantages to the native state, the constraints that dictate the selection between kinetic and thermodynamic folding and stability remain unknown. Using highly conserved subtilases, we demonstrate that adaptive evolution of sequence dictates selection of folding pathways. Intracellular and extracellular serine proteases (ISPs and ESPs, respectively) constitute two subfamilies within the family of subtilases that have highly conserved sequences, structures, and catalytic activities. Our studies on the folding pathways of subtilisin E (SbtE), an ESP, and its homologue intracellular serine protease 1 (ISP1), an ISP, show that although topology, contact order, and hydrophobicity that drive protein folding reactions are conserved, ISP1 and SbtE fold through significantly different pathways and kinetics. While SbtE absolutely requires the propeptide to fold into a kinetically trapped conformer, ISP1 folds to a thermodynamically stable state more than 1 million times faster and independent of a propeptide. Furthermore, kinetics establish that ISP1 and SbtE fold through different intermediate states. An evolutionary analysis of folding constraints in subtilases suggests that observed differences in folding pathways may be mediated through positive selection of specific residues that map mostly onto the protein surface. Together, our results demonstrate that closely related subtilases can fold through distinct pathways and mechanisms, and suggest that fine sequence details can dictate the choice between kinetic and thermodynamic folding and stability.

Amino Acid Sequence↗
PubMed2004

The molecular code for hemoglobin allostery revealed by linking the thermodynamics and kinetics of quaternary structural change. 1. Microstate linear free energy relations.

A novel model linking the thermodynamics and kinetics of hemoglobin's allosteric (R --> T) and ligand binding reactions is applied to photolysis data for human HbCO. To describe hemoglobin's kinetics at the microscopic level of structural transitions and ligand-binding events for individual [ij]-ligation microstates ((ij)R --> (ij)T, (ij)R + CO --> ((i)(+1))(k)R, and (ij)T + CO --> ((i)(+1))(k)T), the model calculates activation energies, (ij)DeltaG(++), from previously measured cooperative free energies of the equilibrium microstates (Huang, Y., and Ackers, G. K. (1996) Biochemistry 35, 704-718) by using linear free energy relations ((ij)DeltaG(++) - (01)DeltaG(++) = alpha[(ij)DeltaG - (01)DeltaG], where the parameter alpha, describing the variation of activation energy with reaction energy perturbation, can depend on the natures of both the reaction and the perturbation). The alpha value measured here for the allosteric dynamics, 0.21 +/- 0.03, corresponds closely to values observed previously, strongly suggesting that the thermodynamic microstate energies directly underlie the allosteric kinetics (as opposed to the alpha((ij)DeltaG(RT)) serving merely as arbitrary fitting parameters). Besides systematizing the study of hemoglobin kinetics, the utility of the microstate linear free energy model lies in the ability to test microscopic aspects of allosteric dynamics such as the "symmetry rule" for quaternary change deduced previously from thermodynamic evidence (Ackers, G. K., et al. (1992) Science 255, 54-63). Reflecting a remarkably detailed correspondence between thermodynamics and kinetics, we find that a kinetic model that includes the large free energy splitting between doubly ligated T microstates implied by the symmetry rule fits the data significantly better than one that does not.

Allosteric Regulation↗