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 703 records · Page 39Linked to original sources

Thermodynamic investigation of hirudin binding to the slow and fast forms of thrombin: evidence for folding transitions in the inhibitor and protease coupled to binding.

Temperature dependent studies of the interaction of the clotting enzyme thrombin with the potent natural inhibitor hirudin reveal a large negative heat capacity change of -1.7(+/- 0.2) kcal/mol per K associated with the formation of the thrombin-hirudin complex, independent of the allosteric state of the enzyme. Binding of N-terminal fragments of hirudin (hir1-49 and hir1-43) is characterized by heat capacity changes of -1.2(+/- 0.1) and -0.9(+/- 0.1) kcal/mol per K, respectively. The magnitude of these heat capacity changes is unprecedented for protease-inhibitor interactions. A thermodynamic analysis based on observed heat capacity and entropy changes predicts that binding is accompanied by substantial coupled folding transitions in both hirudin and thrombin. In the absence of a structure of free thrombin, analysis of differences in the predicted number of residues which fold upon binding hirudin and its fragments leads to the following structural model: three surface loops in thrombin (W60d, W148 and fibrinogen binding loops) are disordered in the free state and fold upon formation of the thrombin-hirudin complex. Molecular dynamics simulations, run over a time scale of 5 ps, are consistent with the hypothesis of large scale coupled folding transitions in both hirudin and thrombin upon formation of the complex. Comparison of the thermodynamics for the interaction of hirudin with the slow and fast forms of thrombin allows dissection of the coupling free energy for allosteric switching. The coupling free energy for the slow-->fast transition increases linearly, in absolute value, with temperature. The coupling enthalpy and entropy terms for hirudin were found to be delta Hoc = 12(+/- 1) kcal/mol and delta Soc = 47(+/- 4) cal/mol per K. Preferential interaction with the fast form is therefore due to the balance of two opposite forces, both quite large in magnitude. The contribution of enthalpic effects opposes the slow-->fast transition and stabilizes binding to the slow form. The contribution of entropic effects favors the slow-->fast transition and stabilizes binding to the fast form. In the physiological temperature range the entropic effects prevail and result in preferential binding of hirudin to the fast form. The region of thrombin recognizing the N-terminal domain of hirudin contains most of the residues that are energetically linked to the slow-->fast transition. This region is part of the "allosteric core" of thrombin and includes the W60d loop, shaping the specificity site S2, and the Na+ binding loop connecting the last two beta-strands of the B chain.

Allosteric Regulation↗

Thermodynamic stability effects of single peptide bond hydrolysis in protein inhibitors of serine proteinases.

Kunitz type soybean trypsin inhibitor (STI) and basic pancreatic trypsin inhibitor (BPTI) were used as model proteins to measure the thermodynamic consequences of single peptide bond hydrolysis. For each inhibitor the reactive site cleaved form was prepared and, additionally, each inhibitor was selectively cleaved at a Met residue. Selective cleavage generally led to reducing of the T(den) value from 7 K up to 75 K. For STI cleaved at Met84 a slight stabilization (increase of T(den) by 1.0 K) was observed. In terms of deltaG(den), the difference between the most extreme cleavage effect was 11.44 kcal/mole, much larger than resulting from the theoretical effects of crosslinks. It was found that hydrolysis of a single peptide bond affects not only entropy, but also enthalpy and heat capacity parameters. Moreover, the sign of change is opposite for two inhibitors: deltaH(den) and deltaS(den) increase for both cleaved forms of STI, while they decrease for two nicked forms of BPTI. To understand the stability effects, a thermodynamic cycle analysis was applied based on comparison of stabilities of intact and cleaved protein with peptide bond hydrolysis equilibria in native and denatured states. The cycle revealed a good agreement of the theoretical effect of crosslink elimination with a free energy difference for hydrolysis of a single peptide bond in a denatured protein. It appears that hydrolysis constants for single peptide bonds in a native protein span over at least 20 orders of magnitude. They are very low for peptide bonds placed in alpha-helices and very high if cleavage reaction leads to a formation of a new secondary structure element.

Aprotinin↗

Thermodynamics of the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site.

What are the thermodynamic consequences of the stepwise conversion of a highly specific (consensus) protein-DNA interface to one that is nonspecific? How do the magnitudes of key favorable contributions to complex stability (burial of hydrophobic surfaces and reduction of DNA phosphate charge density) change as the DNA sequence of the specific site is detuned? To address these questions we investigated the binding of lac repressor (LacI) to a series of 40 bp fragments carrying symmetric (consensus) and variant operator sequences over a range of temperatures and salt concentrations. Variant DNA sites contained symmetrical single and double base-pair substitutions at positions 4 and/or 5 [sequence: see text] in each 10 bp half site of the symmetric lac operator (Osym). Non-specific interactions were examined using a 40 bp non-operator DNA fragment. Disruption of the consensus interface by a single symmetrical substitution reduces the observed equilibrium association constant (K(obs)) for Osym by three to four orders of magnitude; double symmetrical substitutions approach the six orders in magnitude difference between specific and non-specific binding to a 40 bp fragment. At these adjacent positions in the consensus site, the free energy effects of multiple substitutions are non-additive: the first reduces /deltaG(obs)o/ by 3 to 5 kcal mol(-1), approximately halfway to the non-specific level, whereas the second is less deleterious, reducing /deltaG(obs)o/ by less than 3 kcal mol(-1). Variant-specific dependences of K(obs) on temperature and salt concentration characterize these LacI-operator interactions. In general, binding constants and standard free energies of binding both exhibit characteristic extrema near 290 K. As a consequence, both the enthalpic and entropic contributions to stability of Osym and variant complexes change from positive (i.e. entropy driven) at lower temperatures to negative (i.e. enthalpy driven) at higher temperatures, indicating that the heat capacity change upon binding, deltaC(obs)o, is large and negative. In general, /deltaC(obs)o/ decreases as the specificity and stability of the variant complex decreases. Stabilities of complexes of LacI with Osym and all variant operators are strongly [salt]-dependent. Binding constants for the variant complexes exhibit a power-dependence on [salt] that is larger in magnitude (i.e. more negative) than for Osym, but no obvious trend relates changes in contributions from the polyelectrolyte effect and the observed reductions in stability (delta deltaG(obs)o). These variant-specific thermodynamic signatures provide novel insights into the consequences of converting a consensus interface to a less specific one; such insights are not obtained from comparisons at the level of delta deltaG(obs)o. We propose that this variant-specific behavior arises from a strong effect of operator sequence on the extent of induced conformational changes in the protein (and possibly also in the DNA site) which accompany binding.

Bacterial Proteins↗

Thermodynamics of denaturation of mutants of barnase with disulfide crosslinks.

We have measured the effects of disulfide crosslinks on the thermodynamics of denaturation of three mutants of barnase that contain cystine and the corresponding single and double cysteine mutants. At first sight, the data are consistent with the hypothesis that disulfide crosslinks stabilise proteins through entropic destabilisation of the denatured state, but the decreases in the entropy of denaturation are larger than predicted and are accompanied by decreases in the enthalpy of denaturation. These effects are not a unique feature of the disulfide crosslink and are observed in a range of non-crosslinked mutants of barnase as part of a general enthalpy-entropy compensation phenomenon. Similarly, effects on the heat capacity change for denaturation (delta C(p)d), determined from the slope of the enthalpy of denaturation versus temperature, are not confined to mutants with disulfide crosslinks. The value of delta C(p)d is lower in four stabilised mutants than in wild-type barnase, irrespective of the presence of a disulfide crosslink, while the delta C(p)d remains unchanged in a destabilised mutant containing a disulfide. The variation in delta C(p)d may result from an inherent temperature-dependence of delta C(p)d, since it is measured for each mutant over a different temperature range. The thermodynamics of denaturation of the disulfide mutant with a crosslink between positions 70 and 92 change anomalously with pH but in a similar way to that of the D93N mutant of barnase, which lacks the D93-R69 salt-bridge present in the wild-type. This finding confirms initial observations in the X-ray structure of this disulfide mutant that the salt-bridge has been disrupted by the introduced crosslink.

Bacterial Proteins↗

Human alpha-thrombin inhibition by the highly selective compounds N-ethoxycarbonyl-D-Phe-Pro-alpha-azaLys p-nitrophenyl ester and N-carbobenzoxy-Pro-alpha-azaLys p-nitrophenyl ester: a kinetic, thermodynamic and X-ray crystallographic study.

Kinetics, thermodynamics and structural aspects of human alpha-thrombin (thrombin) inhibition by newly synthesized low molecular weight derivatives of alpha-azalysine have been investigated. The thrombin catalyzed hydrolysis of N-ethoxycarbonyl-D-Phe-Pro-alpha-azaLys p-nitrophenyl ester (Eoc-D-Phe-Pro-azaLys-ONp) and N-carbobenzoxy-Pro-alpha-azaLys p-nitrophenyl ester (Cbz-Pro-azaLys-ONp) was investigated at pH 6.2 and 21.0 degrees C, and analyzed in parallel with that of N-alpha-(N,N-dimethylcarbamoyl)-alpha-azalysine p-nitrophenyl ester (Dmc-azaLys-ONp). Decarboxylation following the enzymatic hydrolysis of these p-nitrophenyl esters gave the corresponding 1-peptidyl-2(4-aminobutyl) hydrazines (peptidyl-Abh) showing properties of thrombin competitive inhibitors. Therefore, thermodynamics for the reversible binding of D-Phe-Pro-Abh, Cbz-Pro-Abh and Dmc-Abh to thrombin was examined. These results are consistent with the minimum four-step catalytic mechanism for product inhibition of serine proteinases. Eoc-D-Phe-Pro-azaLys-ONp and Eoc-D-Phe-Pro-Abh display a sub-micromolar affinity for thrombin together with a high selectivity versus homologous plasmatic and pancreatic serine proteinases acting on cationic substrates. The three-dimensional structures of the reversible non-covalent thrombin:Eoc-D-Phe-Pro-Abh and thrombin:Cbz-Pro-Abh complexes have been determined by X-ray crystallography at 2.0 A resolution (R-factor = 0.169 and 0.179, respectively), and analyzed in parallel with that of the thrombin:Dmc-azaLys acyl-enzyme adduct. Both Eoc-D-Phe-Pro-Abh and Cbz-Pro-Abh competitive inhibitors are accommodated in the thrombin active center, spanning the region between the aryl binding site and the S1 primary specificity subsite.

Animals↗

Thermodynamics and kinetics of unfolding of the thermostable trimeric adenylate kinase from the archaeon Sulfolobus acidocaldarius.

The thermal stability of adenylate kinase from the thermoacidophilic archaeon Sulfolobus acidocaldarius was characterized comprehensively using denaturant-induced unfolding, differential scanning calorimetry, circular dichroism spectroscopy, and enzymological inactivation studies. The thermally induced unfolding of the protein is irreversible due to aggregation, whereas the unfolding induced by guanidinium chloride is reversible. The protein is known to be a homotrimer in its native state and we established that it unfolds upon dissociation in the case of denaturant unfolding. We measured the thermodynamic stability of the protein in a temperature range from 5 to 70 degrees C using denaturant unfolding. The protein has a maximum of stability (intrinsic free energy) of 31 kcal/mol-trimer (130 kJ/mol-trimer) at 32 degrees C (based on the linear extrapolation model). The heat capacity change upon unfolding DeltaCp and the m-value were considered to be constant in this temperature range and calculated to be 2.86 kcal/mol-trimer (11.9 kJ/mol-trimer) and 5.67 kcal/mol-trimer M (23.7 kJ/mol-trimer M), respectively. The influence of trimerization on thermodynamic stability was investigated. The several interrelated aspects of thermal stability such as unfolding kinetics, the temperature-dependence of the free energy, and the concentration and temperature-dependencies of the fraction of denatured protein are described quantitatively. The properties of the Gibbs-Helmholtz function of the adenylate kinase from S. acidocaldarius, in particular, and of oligomeric proteins, in general terms, are discussed and compared with the properties of the analogous function for monomeric proteins. Moreover, we discuss methodological aspects: we obtained the analytical expression of the denaturant-unfolding isotherm for homotrimeric proteins; we include a formula Appendix containing the derivations of the expressions used.

Adenylate Kinase↗

Thermodynamics of the unfolding of the cold-shock protein from Thermotoga maritima.

Proteins from (hyper-)thermophiles are known to exhibit high intrinsic stabilities. Commonly, their thermodynamic characterization is impeded by irreversible side reactions of the thermal analysis or calorimetrical problems. Small single-domain proteins are suitable candidates to overcome these obstacles. Here, the thermodynamics of the thermal denaturation of the recombinant cold-shock protein (Csp) from the hyperthermophilic bacterium Thermotoga maritima (Tm) was studied by differential scanning calorimetry. The unfolding transition can be described over a broad pH range (3.5-8.5) by a reversible two-state process. Maximum stability (DeltaG (25 degrees C)=6.5 kcal/mol) was observed at pH 5-6 where Tm Csp unfolds with a melting temperature at 95 degrees C. The heat capacity difference between the native and the denatured states is 1.1(+/-0.1) kcal/(mol K). At pH 7, thermal denaturation occurs at 82 degrees C. The corresponding free energy profile has its maximum at 30 degrees C with DeltaGN-->U=4.8(+/-0.5) kcal/mol. At the optimal growth temperature of T. maritima (80 degrees C), Tm Csp in the absence of ligands is only marginally stable, with a free energy of stabilization not far beyond the thermal energy. With the known stabilizing effect of nucleic acids in mind, this suggests a highly dynamical interaction of Tm Csp with its target molecules.

Bacterial Proteins↗

Folding of a model three-helix bundle protein: a thermodynamic and kinetic analysis.

The kinetics and thermodynamics of an off-lattice model for a three-helix bundle protein are investigated as a function of a bias gap parameter that determines the energy difference between native and non-native contacts. A simple dihedral potential is used to introduce the tendency to form right-handed helices. For each value of the bias parameter, 100 trajectories of up to one microsecond are performed. Such statistically valid sampling of the kinetics is made possible by the use of the discrete molecular dynamics method with square-well interactions. This permits much faster simulations for off-lattice models than do continuous potentials. It is found that major folding pathways can be defined, although ensembles with considerable structural variation are involved. The large gap models generally fold faster than those with a smaller gap. For the large gap models, the kinetic intermediates are non-obligatory, while both obligatory and non-obligatory intermediates are present for small gap models. Certain large gap intermediates have a two-helix microdomain with one helix extended outward (as in domain-swapped dimers); the small gap intermediates have more diverse structures. The importance of studying the kinetic, as well as the thermodynamics, of folding for an understanding of the mechanism is discussed and the relation between kinetic and equilibrium intermediates is examined. It is found that the behavior of this model system has aspects that encompass both the "new" view and the "old" view of protein folding.

Algorithms↗

Anticooperative ligand binding properties of recombinant ferric Vitreoscilla homodimeric hemoglobin: a thermodynamic, kinetic and X-ray crystallographic study.

Thermodynamics and kinetics for cyanide, azide, thiocyanate and imidazole binding to recombinant ferric Vitreoscilla sp. homodimeric hemoglobin (Vitreoscilla Hb) have been determined at pH 6.4 and 7.0, and 20.0 degrees C, in solution and in the crystalline state. Moreover, the three-dimensional structures of the diligated thiocyanate and imidazole derivatives of recombinant ferric Vitreoscilla Hb have been determined by X-ray crystallography at 1.8 A (Rfactor=19.9%) and 2.1 A (Rfactor=23.8%) resolution, respectively. Ferric Vitreoscilla Hb displays an anticooperative ligand binding behaviour in solution. This very unusual feature can only be accounted for by assuming ligand-linked conformational changes in the monoligated species, which lead to the observed 300-fold decrease in the affinity of cyanide, azide, thiocyanate and imidazole for the monoligated ferric Vitreoscilla Hb with respect to that of the fully unligated homodimer. In the crystalline state, thermodynamics for azide and imidazole binding to ferric Vitreoscilla Hb may be described as a simple process with an overall ligand affinity for the homodimer corresponding to that for diligation in solution. These data suggest that the ligand-free homodimer, observed in the crystalline state, is constrained in a low affinity conformation whose ligand binding properties closely resemble those of the monoligated species in solution. From the kinetic viewpoint, anticooperativity is reflected by the 300-fold decrease of the second-order rate constant for cyanide and imidazole binding to the monoligated ferric Vitreoscilla Hb with respect to that for ligand association to the ligand-free homodimer in solution. On the other hand, values of the first-order rate constant for cyanide and imidazole dissociation from the diligated and monoligated derivatives of ferric Vitreoscilla Hb in solution are closely similar. As a whole, ligand binding and structural properties of ferric Vitreoscilla Hb appear to be unique among all Hbs investigated to date.

Crystallography, X-Ray↗

DNA sequence context modulates the impact of a cisplatin 1,2-d(GpG) intrastrand cross-link on the conformational and thermodynamic properties of duplex DNA.

The anticancer activity of cisplatin derives from its ability to bind and cross-link DNA, with the major adduct being the 1,2-d(GpG) intrastrand cross-link. Here, the consequences of this adduct on the conformation, thermal stability, and energetics of duplex DNA are assessed, and the modulation of these parameters by the sequence context of the adduct is evaluated. The properties of a family of 15-mer DNA duplexes containing a single 1,2-d(GpG) cis-¿Pt(NH(3))(2)¿(2+) intrastrand cross-link are probed in different sequence contexts where the flanking base-pairs are systematically varied from T.A to C.G to A.T. By using a combination of spectroscopic and calorimetric techniques, the structural, thermal, and thermodynamic properties of each duplex, both with and without the cross-link, are characterized. Circular dichroism spectroscopic data reveal that the cross-link alters the structure of the host duplex in a manner consistent with a shift from a B-like to an A-like conformation. Thermal denaturation data reveal that the cross-link induces substantial thermal and thermodynamic destabilization of the host duplex. Significantly, the magnitudes of these cross-link-induced effects on duplex structure, thermal stability, and energetics are influenced by the bases that flank the adduct. The presence of flanking A.T base-pairs, relative to T.A or C.G base-pairs, enhances the extent of cross-link-induced alteration to an A-like conformation and dampens the extent of cross-link-induced duplex destabilization. These results are discussed in terms of available structural data, and in terms of the selective recognition of cisplatin-DNA adducts by HMG-domain proteins.

Base Pairing↗

A thermodynamic and structural analysis of DNA minor-groove complex formation.

As part of an effort to develop a better understanding of the structural and thermodynamic principles of DNA minor groove recognition, we have investigated complexes of three diphenylfuran dications with the d(CGCGAATTCGCG)(2) duplex. The parent compound, furamidine (DB75), has two amidine substituents while DB244 has cyclopentyl amidine substituents and DB226 has 3-pentyl amidines. The structure for the DB244-DNA complex is reported here and is compared to the structure of the DB75 complex. Crystals were not obtained with DB226 but information from the DB75 and DB244 structures as well as previous NMR results on DB226 indicate that all three compounds bind in the minor groove at the AATT site of the duplex. DB244 and DB75 penetrate to the floor of the groove and form hydrogen bonds with T8 on one strand and T20 on the opposite strand while DB226 forms a complex with fewer interactions. Binding studies by surface plasmon resonance (SPR) yield -delta G degrees values in the order DB244>DB75>DB226 that are relatively constant with temperature. The equilibrium binding constants for DB244 are 10-20 times greater than that for DB226. Isothermal titration calorimetric (ITC) experiments indicate that, in contrast to delta G degrees, delta H degrees varies considerably with temperature to yield large negative delta Cp degrees values. The thermodynamic results, analyzed in terms of structures of the DNA complexes, provide an explanation of why DB244 binds more strongly to DNA than DB75, while DB266 binds more weakly. All three compounds have a major contribution to binding from hydrophobic interactions but the hydrophobic term is most favorable for DB244. DB244 also has strong contributions from molecular interactions in its DNA complex and all of these factors combine to give it the largest-delta G degrees for binding. Although the factors that influence the energetics of minor groove interactions are varied and complex, results from the literature coupled with those on the furan derivatives indicate that there are some common characteristics for minor groove recognition by unfused heterocyclic cations that can be used in molecular design.

Base Sequence↗

Molecular-level thermodynamic and kinetic parameters for the self-assembly of apoferritin molecules into crystals.

The self-assembly of apoferritin molecules into crystals is a suitable model for protein crystallization and aggregation; these processes underlie several biological and biomedical phenomena, as well as for protein and virus self-assembly. We use the atomic force microscope in situ, during the crystallization of apoferritin to visualize and quantify at the molecular level the processes responsible for crystal growth. To evaluate the governing thermodynamic parameters, we image the configuration of the incorporation sites, "kinks", on the surface of a growing crystal. We show that the kinks are due to thermal fluctuations of the molecules at the crystal-solution interface. This allows evaluation of the free energy of the intermolecular bond phi=3.0 k(B)T=7.3 kJ/mol. The crystallization free energy, extracted from the protein solubility, is -42 kJ/mol. Published determinations of the second virial coefficient and the protein solubility between 0 and 40 degrees C revealed that the enthalpy of crystallization is close to zero. Analyses based on these three values suggest that the main component in the crystallization driving force is the entropy gain of the water molecules bound to the protein molecules in solution and released upon crystallization. Furthermore, monitoring the incorporation of individual molecules in to the kinks, we determine the characteristic frequency of attachment of individual molecules at one set of conditions. This allows a correlation between the mesoscopic kinetic coefficient for growth and the molecular-level thermodynamic and kinetic parameters determined here. We found that step growth velocity, scaled by the molecular size, equals the product of the kink density and attachment frequency, i.e. the latter pair are the molecular-level parameters for self-assembly of the molecules into crystals.

Apoferritins↗

Thermodynamics of trimer-of-hairpins formation by the SIV gp41 envelope protein.

The gp41 envelope protein mediates the entry of primate immunodeficiency viruses into target cells by promoting the fusion of viral and cellular membranes. The structure of the gp41 ectodomain core represents a trimer of identical helical hairpins in which a central trimeric coiled-coil made up of three amino-terminal helices is wrapped in an outer layer of three antiparallel carboxyl-terminal helices. Triggering formation of this fusion-active gp41 conformation appears to cause close membrane apposition and thus overcome the activation energy barrier for lipid bilayer fusion. We present a detailed description of the folding thermodynamics of the simian immunodeficiency virus (SIV) gp41 core by using a recombinant trimeric N34(L6)C28 model. Differential scanning calorimetry and spectroscopic experiments on denaturant-induced and thermal unfolding indicate that the free energy of association of three 68 residue N34(L6)C28 peptides to a trimer-of-hairpins is 76 kJ mol(-1) at pH 7.0 and 25 degrees C in physiological buffer. The associated enthalpy change, Delta H(unf), is 177 kJ mol(-1), while the entropy of unfolding, Delta S(unf), is 0.32 kJ K(-1) mol(-1). The temperature of maximal stability is close to 20 degrees C. The unfolding heat capacity increment is approximately 9 kJ K(-1) mol(-1) (approximately 45 J K(-1) mol residue(-1)), which is lower than expected for unfolding of the trimer to an extended and fully hydrated polypeptide chain. Replacement by isoleucine of the polar residues Thr582 or Thr586 buried in the N-terminal trimeric coiled-coil interface leads to very strong stabilization of the trimer-of-hairpins, 30-35 kJ mol(-1). Single-point mutations in the central coiled-coil thus strongly stabilize the gp41 core structure. These thermodynamic characteristics may be important for the refolding of the gp41 envelope protein into its fusion-active conformation during membrane fusion.

Amino Acid Sequence↗

Identification of thermodynamically relevant interactions between EF-Tu and backbone elements of tRNA.

A set of 45 different tRNAs, each containing a single deoxynucleotide substitution covering the upper half of the molecule was used in conjunction with a high-throughput ribonuclease protection assay to investigate the thermodynamic role of 2' hydroxyl groups in stabilizing a complex with elongation factor Tu (EF-Tu) from Thermus thermophilus. Five distinct 2' hydroxyl groups were identified where substitution with a proton resulted in an approximately tenfold decrease in the binding affinity. The same five 2' hydroxyl groups reduced the affinity of the interaction with the nearly identical Thermus aquaticus EF-Tu. Four of these 2' hydroxyl groups were observed to form hydrogen bonds in a co-crystal structure of tRNA(Phe) and T. aquaticus EF-Tu, while the fifth 2' hydroxyl group can be associated with an intramolecular hydrogen bond in the tRNA. However, four additional hydrogen bonds to 2' hydroxyl groups observed in the crystal structure show no thermodynamic effect upon disruption. Some of these discrepancies may be reconciled based on the unbound structures of the protein and RNA.

Alanine↗

The osmophobic effect: natural selection of a thermodynamic force in protein folding.

Intracellular organic osmolytes are present in certain organisms adapted to harsh environments and these osmolytes protect intracellular macromolecules against the denaturing environmental stress. In natural selection of organic osmolytes as protein stabilizers, it appears that the osmolyte property selected for is the unfavorable interaction between the osmolyte and the peptide backbone, a solvophobic thermodynamic force that we call the osmophobic effect. Because the peptide backbone is highly exposed to osmolyte in the denatured state, the osmophobic effect preferentially raises the free energy of the denatured state, shifting the equilibrium in favor of the native state. By focusing the solvophobic force on the denatured state, the native state is left free to function relatively unfettered by the presence of osmolyte. The osmophobic effect is a newly uncovered thermodynamic force in nature that complements the well-recognized hydrophobic interactions, hydrogen bonding, electrostatic and dispersion forces that drive protein folding. In organisms whose survival depends on the intracellular presence of osmolytes that can counteract denaturing stresses, the osmophobic effect is as fundamental to protein folding as these well-recognized forces.

Amino Acid Sequence↗

Comparative thermodynamics for monomer and dimer sequence-dependent binding of a heterocyclic dication in the DNA minor groove.

Phenylamidine cationic groups linked by a furan ring (furamidine) and related symmetric diamidine compounds bind as monomers in the minor groove of AT sequences of DNA. DB293, an unsymmetric derivative with one of the phenyl rings of furamidine replaced with a benzimidazole, can bind to AT sequences as a monomer but binds more strongly to GC-containing minor-groove DNA sites as a stacked dimer. The dimer-binding mode has high affinity, is highly cooperative and sequence selective. In order to develop a better understanding of the correlation between structural and thermodynamic aspects of DNA molecular recognition, DB293 was used as a model to compare the binding of minor-groove agents with AT and mixed sequence DNA sites. Isothermal titration calorimetry and surface plasmon resonance results clearly show that the binding of DB293 and other related compounds into the minor groove of AT sequences is largely entropy-driven while the binding of DB293 as a dimer into the minor groove of GC-containing sequences is largely enthalpy-driven. At 25 degrees C, for example, the AT binding has DeltaG degrees, DeltaH degrees and TDeltaS degrees values of -9.6, -3.6 and 6.0 kcal/mol while the values for dimer binding to a GC-containing site are -9.0, -10.9 and -1.9 kcal/mol (per mol of bound compound), respectively. These results show that the thermodynamic components for binding of compounds of this type to DNA are very dependent on the structure, solvation and sequence of the DNA binding site.

Base Sequence↗

The incorporation of gel pressure into the irreversible thermodynamic equation of fluid flow in order to explain biological tissue swelling.

After the Donnan osmotic potential derived theory of corneal stromal swelling was suggested (Hodson, 1971) it later was proposed, and verified experimentally, (Elliott et al., 1980) that the theory could be applied even in situations where there were no bounding cellular membranes attached to the corneal stroma. This was a paradox. The central problem with swelling in the absence of membranes (i.e. with a gel boundary whose reflexion coefficient is essentially zero) is that the driving force is not readily understandable as osmotic (but see Oster & Peskin, 1992) and yet the magnitude of the swelling fitted quantitatively with Donnan theory. We believe that we have resolved this paradox with a new and expanded irreversible thermodynamic relationship for fluid flows into and across biological tissues. In the course of the derivation, fresh concepts arise: for example the relationship of salt disparity is described which forbids diffusible salt-generated chemical and osmotic potentials to be simultaneously at equilibrium in the presence of ionised macromolecules and this relationship is developed to generate a new intrinsic thermodynamic property which is termed gel pressure and which drives fluid flow. Gel pressure provides a theoretical basis for biological tissue swelling. Microscopically, gel pressure is identified as the electrostatic potential developed by the mutual repulsion of the fixed matrix charges.

Animals↗

Non-equilibrium thermodynamics of molecular evolution.

The evolution of the information complexity of a large database of protein sequences is investigated. The information entropy for protein sequences is determined from their algorithmic complexity and is found to change with evolutionary time at a constant rate. The information content of changed residues is always lower than the content of conserved residues. This indicates that sequences are becoming less random throughout evolution. It also shows that the system is being driven toward minimal complexity production. The change in information content per amino acid substitution is virtually identical for all the protein sequences studied. These results are interpreted with a statistical mechanical theory that ties sequence information to the thermodynamics of protein structure. Sequence evolution is viewed as a means to drive the system to minimum thermodynamic entropy production in a stable, non-equilibrium state. This theory provides a physical framework for understanding molecular evolution and incorporates features of both the neutralist and selectionist models.

Amino Acid Sequence↗