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 433 records · Page 24Linked to original sources
PubMed1996

Structure-based thermodynamic scale of alpha-helix propensities in amino acids.

A structural parameterization of the folding energetics has been used to predict the effect of single amino acid mutations at exposed locations in alpha-helices. The results have been used to derive a structure-based thermodynamic scale of alpha-helix propensities for amino acids. The structure-based thermodynamic analysis was performed for four different systems for which structural and experimental thermodynamic data are available: T4 lysozyme [Blaber et al (1994) J. Mol. Biol.235, 600-624], barnase [Horovitz et al. (1992) J.Mol.Biol.227,560-568], a synthetic leucine zipper [O'Neil & Degrado (1990) Science 250, 646-651], and a synthetic peptide [Lyu et al. (1990) Science 250, 669-673]. These studies have permitted the optimization of the set of solvent-accessible surface areas (ASA) for all amino acids in the unfolded state. It is shown that a single set of structure/thermodynamic parameters accounts well for all the experimental data sets of helix propensities. For T4 lysozyme, the average value of the absolute difference between predicted and experimental delta G values is 0.09 kcal/mol, for barnase 0.14 kcal/mol, for the synthetic coiled-coil 0.11 kcal/mol, and for the synthetic peptide 0.08 kcal/mol. In addition, this approach predicts well the overall stability of the proteins and rationalizes the differences in alpha-helix propensities between amino acids. The excellent agreement observed between predicted and experimental delta G values for all amino acids validates the use of this structural parameterization in free energy calculations for folding or binding.

Amino Acids↗
PubMed1996

Structural and thermodynamic characterization of the interaction of the SH3 domain from Fyn with the proline-rich binding site on the p85 subunit of PI3-kinase.

The interaction of the Fyn SH3 domain with the p85 subunit of PI3-kinase is investigated using structural detail and thermodynamic data. The solution structure complex of the SH3 domain with a proline-rich peptide mimic of the binding site on the p85 subunit is described. This indicates that the peptide binds as a poly(L-proline) type II helix. Circular dichroism spectroscopic studies reveal that in the unbound state the peptide exhibits no structure. Thermodynamic data for the binding of this peptide to the SH3 domain suggest that the weak binding (approximately 31 microM) of this interaction is, in part, due to the entropically unfavorable effect of helix formation (delta S0 = -78 J.mol-1.K-1). Binding of the SH3 domain to the intact p85 subunit (minus its own SH3 domain) is tighter, and the entropic and enthalpic contributions are very different from those given by the peptide interaction (delta S0 = +252 J.mol-1.K-1; delta H0 = +44 kJ.mol-1). From these dramatically different thermodynamic measurements we are able to conclude that the interaction of the proline-rich peptide does not effectively mimic the interaction of the intact p85 subunit with the SH3 domain and suggest that other interactions could be important.

Amino Acid Sequence↗
PubMed1997

Thermodynamic stability of annexin V E17G: equilibrium parameters from an irreversible unfolding reaction.

Conformational stability of the membrane-binding protein annexin V E17G has been determined by high-sensitivity differential scanning microcalorimetry (DSC) measurements and by isothermal, guanidinium hydrochloride (GdnHCl)-induced unfolding studies. Wild-type annexin V and the E17G mutant protein studied here are structurally almost identical. Therefore, it can be expected that the present results will not deviate significantly from the stability data of the wild-type molecule. Thermal unfolding is irreversible, while GdnHCl unfolding shows a high degree of reversibility. We were able to demonstrate that characteristic features of annexin V E17G unfolding permit us to extract from the kinetically controlled heat capacity curves thermodynamic equilibrium parameters at the high heating rates. The thermodynamic quantities obtained from the DSC studies in phosphate buffer at pH 7.0 are as follows: t1/2 = 54.7 degrees C (heating rate of 2.34 K min-1), delta H0 = 690 kJ mol-1, and delta Cp = 10.3 kJ mol-1 K-1 which correspondends to a value of delta G0D (20 degrees C) of 53.4 kJ mol-1. When compared on a per gram basis, these thermodynamic parameters classify annexin V E17G as a marginally stable protein. This conclusion is consistent with structural and functional features of the protein that require conformational adaptability for hinge-bending motions and pore formation on interaction with membranes. We observed a large difference between the change in the Gibbs energy value derived from the heat capacity studies and that determined from the GdnHCl unfolding curve. The difference appears to stem from a specific interaction of the protein with the denaturant that results in both a low half-denaturation concentration C1/2 of 1.74 M and a small slope (6.0 kJ L mol-2) of the delta Gapp versus [GdnHCl] plot. The extraordinary interaction of annexin V with GdnHCl is also manifested in the enormous depression of the transition temperature delta t1/2 (= 18 degrees C) when the GdnHCl concentration is increased from 0 to 1 M. "Regular" proteins experience an average decrease in the transition temperature of 8 +/- 2 degrees C per 1 M change in the concentration of GdnHCl.

Annexin A5↗
PubMed1997

Reexamination of the hyper thermodynamic stability of asymmetric cyanomet valency hybrid hemoglobin, (alpha+CN-beta+CN-)(alpha beta): no preferentially populating asymmetric hybrid at equilibrium.

It has been reported that hybridization of the equimolar mixture of cyanomethemoglobin and deoxyhemoglobin through dimer exchange reaction results in establishment of an approximately binomial (1:2:1) equilibrium distribution of these parental hemoglobins and their hybrid molecule, (alpha+CN-beta+CN-)(alpha beta), within several days under anaerobic conditions at pH 7.4, 21.5 degrees C, leading to a hyper (i.e., about 170 times) thermodynamic stability of (alpha+CN-beta+CN-)(alpha beta) relative to the stability of the other diliganded species at pH 7.4, 21.5 degrees C [Daugherty, M. A., Shea, M. A., Johnson, J. A., LiCata, V. J., Turner, G. J., & Ackers, G. K. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 1110-1114]. To examine whether the published "binomiality" for this deoxy-cyanomet hybrid system really reflects the thermodynamic stability of (alpha+CN-beta+CN-)(alpha beta), we used a binomial (1:2:1) equilibrium distribution of the equimolar mixture of cyanomethemoglobin and fully oxygenated hemoglobin as a starting condition, and then this system was rapidly deoxygenated. We found that the relative population of the hybrid was reduced to 8.6% of the total upon deoxygenation. It was also found that the hybridization experiment under anaerobic conditions was not allowed to continue for a long time due to a valency exchange reaction between deoxy and cyanomet derivatives. For instance, a 48 h incubation resulted in the oxidation of 44% of Fe2+ to Fe3+ hemes in the original deoxyhemoglobin and the reduction of 42% of Fe3+ to Fe2+ hemes in the original cyanomethemoglobin. These results suggest that a real distribution of the deoxy-cyanomet hybrid system at equilibrium is fairly far from 1:2:1 (binomial distribution), and the thermodynamic stability of (alpha+CN-beta+CN-)(alpha beta) is less than one-tenth of the hyperstability previously reported. In addition, most of the previous results on deoxy-cyanomet valency hybrids placed under long anaerobic conditions should be subject to reexamination due to possible valency exchange reactions.

Adult↗
PubMed1997

Investigation of the structural basis for thermodynamic stabilities of tandem GU wobble pairs: NMR structures of (rGGAGUUCC)2 and (rGGAUGUCC)2.

The symmetric, tandem GU mismatch motifs, and , which only differ in the mismatch order, have an average difference in thermodynamic stability of 2 kcal/mol at 37 degrees C. Thermodynamic studies of duplexes containing these motifs indicate the effect is largely localized to the mismatches and adjacent base pairs. The three-dimensional structures of two representative duplexes, (rGGAGUUCC)2 and (rGGAUGUCC)2, were determined by two-dimensional NMR and a simulated annealing protocol. Local deviations are similar to other intrahelical GU mismatches with little effect on backbone torsion angles and a slight overtwisting between the base pair 5' of the G of the mismatch and the mismatch itself. Comparisons of the resulting stacking patterns along with electrostatic potential maps suggest that interactions between highly negative electrostatic regions between base pairs may play a role in the observed thermodynamic differences.

Base Composition↗
PubMed1997

Thermodynamic and structural analysis of phosphotyrosine polypeptide binding to Grb2-SH2.

A thermodynamic analysis using isothermal titration calorimetry (ITC) has been performed to examine the binding interaction between the SH2 (Src homology 2) domain of growth factor receptor binding protein 2 (Grb2-SH2) and one of its phosphotyrosine (pY) polypeptide ligands. Interaction of the Shc-derived phosphotyrosine hexapeptide Ac-SpYVNVQ-NH2 with Grb2-SH2 was both enthalpically and entropically favorable (DeltaH = -7.55 kcal mol-1, -TDeltaS = -1.46 kcal mol-1 , DeltaG = -9.01 kcal mol-1, T = 20 degrees C). ITC experiments using five alanine-substituted peptides were performed to examine the role of each side chain in binding. The results were consistent with homology models of the Grb2-SH2-Shc hexapeptide complex which identified several possible hydrogen bonds between Grb2-SH2 and the phosphotyrosine and conserved asparagine(+2) side chains of the Shc hexapeptide. These studies also demonstrated that the hydrophobic valine(+1) side chain contributes significantly to the favorable entropic component of binding. The thermodynamic and structural data are consistent with a Grb2-SH2 recognition motif of pY-hydrophobic-N-X (where X is any amino acid residue). The measured heat capacity of binding (DeltaCp = -146 cal mol-1 K-1) was very similar to computed values using semiempirical estimates (DeltaCp = -106 to -193 cal mol-1 K-1) derived from apolar and polar accessible surface area values calculated from several homology models of the Grb2-SH2-Shc hexapeptide complex. The homology model which most closely reproduced the measured DeltaCp value is also the model which had the lowest RMS deviation from the subsequently determined crystal structure. Calculations based on the thermodynamic data and these semiempirical estimates indicated that the binding event involves burial of nearly comparable apolar (677 A2) and polar (609 A2) surface areas.

Adaptor Proteins, Signal Transducing↗
PubMed1998

The impact of an exocyclic cytosine adduct on DNA duplex properties: significant thermodynamic consequences despite modest lesion-induced structural alterations.

The exocyclic base adduct 3,N4-deoxyethenocytosine (epsilonC) is a common DNA lesion that can arise from carcinogen exposure and/or as a biproduct of cellular processes. We have examined the thermal and thermodynamic impact of this lesion on DNA duplex properties, as well as the structural alterations imparted by the lesion. For these studies, we used calorimetric and spectroscopic techniques to investigate a family of 13-mer DNA duplexes of the form (5'CGCATGNGTACGC3')x(3'GCGTACNCATGCG5'), where the central NxN base pair represents the four standard Watson-Crick base pairs (corresponding to four control duplexes), and where either one of the N bases has been replaced by epsilonC, yielding eight test duplexes. Studies on these 12 duplexes permit us to assess the impact of the epsilonC lesion as a function of sequence context. Our spectroscopic and calorimetric data allow us to reach the following conclusions: (i) The epsilonC lesion imparts a large penalty on duplex stability, with sequence context only modestly modulating the extent of this lesion-induced destabilization. This result contrasts with our recent studies of duplexes with abasic sites, where sequence context was found to be the predominant determinant of thermodynamic damage. (ii) For the epsilonC-containing duplexes, sequence context effects are most often observed in the enthalpic contribution to lesion-induced duplex destabilization. However, due to compensating entropies, the free energy changes associated with this lesion-induced duplex destablization are nearly independent of sequence context. (iii) Despite significant lesion-induced changes in duplex energetics, our spectroscopic probes detect only modest lesion-induced changes in duplex structure. In fact, the overall duplex maintains a global B-form conformation, in agreement with NMR structural data. We discuss possible interpretations of the apparent disparity between the severe thermodynamic and relatively mild structural impacts of the epsilonC lesion on duplex properties. We also note and discuss the implications of empirical correlations between biophysical and biological properties of lesion-containing duplexes.

Base Sequence↗
PubMed1998

Monitoring the sizes of denatured ensembles of staphylococcal nuclease proteins: implications regarding m values, intermediates, and thermodynamics.

Fluorescence and size-exclusion chromatography (SEC) are used to monitor urea denaturation of wild-type staphylococcal nuclease (SN) as well as the m+ and m- mutants A69T and V66W, respectively. It is found that the SEC partition coefficient, 1/Kd, is directly proportional to the Stokes radii of proteins. From the Stokes radii, the denatured ensembles of the three proteins are found to be highly compact in the limit of low urea concentration and expand significantly with increasing urea concentration. The m values from fluorescence-detected denaturation of the SN proteins are generally considered to reflect the relative sizes of denatured ensembles. However, the rank order of m values of the SN proteins studied do not correspond to the rank order of denatured ensemble sizes detected by 1/Kd, suggesting that m values reflect more than just surface area increases on denaturation. SEC provides two complementary ways to demonstrate the existence of intermediates in urea denaturation and illustrates that V66W undergoes a three-state transition. Fluorescence-detected urea denaturations of A69T and wt SN do not correspond with 1/Kd-detected denaturation profiles, a result that would ordinarily mean that the transitions are non-two-state. However, this interpretation fails to recognize the rapidly changing size and thermodynamic character of the denatured ensembles of these proteins both within and outside of the transition zone. The implications of the changing sizes and thermodynamic character of the denatured ensembles for SN proteins are manifold, requiring a reconsideration of the thermodynamics of proteins whose denatured ensembles behave as those of SN proteins.

Alanine↗
PubMed1999

Comparing the thermodynamic stabilities of a related thermophilic and mesophilic enzyme.

Several models have been proposed to explain the high temperatures required to denature enzymes from thermophilic organisms; some involve greater maximum thermodynamic stability for the thermophile, and others do not. To test these models, we reversibly melted two analogous protein domains in a two-state manner. E2cd is the isolated catalytic domain of cellulase E2 from the thermophile Thermomonospora fusca. CenAP30 is the analogous domain of the cellulase CenA from the mesophile Cellulomonas fimi. When reversibly denatured in a common buffer, the thermophilic enzyme E2cd had a temperature of melting (Tm) of 72.2 degrees C, a van't Hoff enthalpy of unfolding (DeltaHVH) of 190 kcal/mol, and an entropy of unfolding (DeltaSu) of 0.55 kcal/(mol*K); the mesophilic enzyme CenAP30 had a Tm of 56.4 degrees C, a DeltaHVH of 107 kcal/mol, and a DeltaSu of 0. 32 kcal/(mol*K). The higher DeltaHVH and DeltaSu values for E2cd suggest that its free energy of unfolding (DeltaGu) has a steeper dependence on temperature at the Tm than CenAP30. This result supports models that predict a greater maximum thermodynamic stability for thermophilic enzymes than for their mesophilic counterparts. This was further explored by urea denaturation. Under reducing conditions at 30 degrees C, E2cd had a concentration of melting (Cm) of 5.2 M and a DeltaGu of 11.2 kcal/mol; CenAP30 had a Cm of 2.6 M and a DeltaGu of 4.3 kcal/mol. Under nonreducing conditions, the Cm and DeltaGu of CenAP30 were increased to 4.5 M and 10.8 kcal/mol at 30 degrees C; the Cm for E2cd was increased to at least 7.4 M at 32 degrees C. We were unable to determine a DeltaGu value for E2cd under nonreducing conditions due to problems with reversibility. These data suggest that E2cd attains its greater thermal stability (DeltaTm = 15.8 degrees C) through a greater thermodynamic stability (DeltaDeltaGu = 6.9 kcal/mol) compared to its mesophilic analogue CenAP30.

Actinomycetales↗
PubMed1999

Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A.A, C.C, G.G, and T.T mismatches.

Thermodynamic measurements are reported for 51 DNA duplexes with A.A, C.C, G.G, and T.T single mismatches in all possible Watson-Crick contexts. These measurements were used to test the applicability of the nearest-neighbor model and to calculate the 16 unique nearest-neighbor parameters for the 4 single like with like base mismatches next to a Watson-Crick pair. The observed trend in stabilities of mismatches at 37 degrees C is G.G > T.T approximately A.A > C.C. The observed stability trend for the closing Watson-Crick pair on the 5' side of the mismatch is G.C >/= C.G >/= A.T >/= T.A. The mismatch contribution to duplex stability ranges from -2.22 kcal/mol for GGC.GGC to +2.66 kcal/mol for ACT.ACT. The mismatch nearest-neighbor parameters predict the measured thermodynamics with average deviations of DeltaG degrees 37 = 3.3%, DeltaH degrees = 7. 4%, DeltaS degrees = 8.1%, and TM = 1.1 degrees C. The imino proton region of 1-D NMR spectra shows that G.G and T.T mismatches form hydrogen-bonded structures that vary depending on the Watson-Crick context. The data reported here combined with our previous work provide for the first time a complete set of thermodynamic parameters for molecular recognition of DNA by DNA with or without single internal mismatches. The results are useful for primer design and understanding the mechanism of triplet repeat diseases.

Base Pair Mismatch↗
PubMed1999

Thermodynamic instability of human lambda 6 light chains: correlation with fibrillogenicity.

Certain types of human light chains have the propensity to deposit pathologically as amyloid fibrils as evidenced by the preferential association of monoclonal lambda 6 proteins with AL amyloidosis. However, the molecular features that render such proteins amyloidogenic have not been elucidated. Based upon the demonstrated relationship between the thermodynamic stability of light chains and their propensity to aggregate in vitro, we have initiated studies where the thermodynamic properties and fibrillogenic potential of two recombinant (r) V lambda 6 molecules were compared. The first protein was generated from cDNA cloned from marrow-derived plasma cells from a patient (Wil) who had AL amyloidosis and renal amyloid deposits; the second was from a patient (Jto) with multiple myeloma in whom the lambda 6 protein was deposited not as amyloid but in the form of renal tubular casts. The thermodynamic stabilities of rV lambda 6Wil and -Jto were determined from chaotropic and thermal denaturation studies. Based upon the Delta GH2O, Delta H, Delta G25 degrees C, Tm, and Cm values, the rV lambda 6Wil was less stable than its nonamyloidogenic counterpart, rV lambda 6Jto. Measurement of fibril formation using a novel in vitro fibril forming assay demonstrated that although both rV lambda 6 proteins formed fibrils in vitro, Wil had a shorter lag time and exhibited faster kinetics under physiologic conditions. Comparative amino acid sequence analyses of these two components and other lambda 6 amyloid-associated light chains revealed that the Jto protein had certain primary structural features that we posit contributed to its increased stability and thus rendered this protein nonamyloidogenic. Our studies provide the first evidence that stabilizing interactions within the V L domain can influence the kinetics of light chain fibrillogenicity.

Amino Acid Sequence↗
PubMed2000

Thermodynamic analysis of cavity creating mutations in an engineered leucine zipper and energetics of glycerol-induced coiled coil stabilization.

Protein stability in vitro can be influenced either by introduction of mutations or by changes in the chemical composition of the solvent. Recently, we have characterized the thermodynamic stability and the rate of folding of the engineered dimeric leucine zipper A(2), which has a strengthened hydrophobic core [Dürr, E., Jelesarov, I., and Bosshard, H. R. (1999) Biochemistry 38, 870-880]. Here we report on the energetic consequences of a cavity introduced by Leu/Ala substitution at the tightly packed dimeric interface and how addition of 30% glycerol affects the folding thermodynamics of A(2) and the cavity mutants. Folding could be described by a two-state transition from two unfolded monomers to a coiled coil dimer. Removal of six methylene groups by Leu/Ala substitutions destabilized the dimeric coiled coil by 25 kJ mol(-1) at pH 3.5 and 25 degrees C in aqueous buffer. Destabilization was purely entropic at around room temperature and became increasingly enthalpic at elevated temperatures. Mutations were accompanied by a decrease of the unfolding heat capacity by 0.5 kJ K(-1) mol(-1). Addition of 30% glycerol increased the free energy of folding of A(2) and the cavity mutants by 5-10 kJ mol(-1) and lowered the unfolding heat capacity by 25% for A(2) and by 50% for the Leu/Ala mutants. The origin of the stabilizing effect of glycerol varied with temperature. Stabilization of the parent leucine zipper A(2) was enthalpic with an unfavorable entropic component between 0 and 100 degrees C. In the case of cavity mutants, glycerol induced enthalpic stabilization below 50 degrees C and entropic stabilization above 50 degrees C. The effect of glycerol could not be accounted for solely by the enthalpy and entropy of transfer or protein surface from water to glycerol/water mixture. We propose that in the presence of glycerol the folded coiled coil dimer is better packed and displays less intramolecular fluctuations, leading to enhanced enthalpic interactions and to an increase of the entropy of folding. This work demonstrates that mutational and solvent effects on protein stability can be thermodynamically complex and that it may not be sufficient to only analyze changes of enthalpy and entropy at the unfolding temperature (T(m)) to understand the mechanisms of protein stabilization.

Acids↗
PubMed2006

Thermodynamic analysis of the impact of the surfactant-protein interactions on the molecular parameters and surface behavior of food proteins.

This paper reports on the thermodynamics of the interactions between surfactants (anionic, CITREM, SSL; nonionic, PGE; zwitterionic, phospholipids) and food proteins (sodium caseinate, legumin) depending on the chemical structure and molecular state (individual molecules, micelles) of the surfactants and the molecular parameters (conformation, molar mass, charge) of the proteins under changes of pH in the range from 7.2 to 5.0 and temperature from 293 to 323 K. The marked effect of the protein-surfactant interactions on the molecular parameters (the weight-average molar mass, the gyration and hydrodynamic radii) and the thermodynamic affinity of the proteins for an aqueous medium were determined by a combination of static and dynamic laser light scattering. Thermodynamically justified schematic sketches of the molecular mechanisms of the complex formation between like-charged proteins and surfactants have been proposed. In response to the complex formation between the proteins and the surfactants, the more stable and fine foams have been detected generally.

Anions↗
PubMedDate not supplied

Kinetic and thermodynamic investigation on ascorbate oxidase activity and stability of a Cucurbita maxima extract.

The kinetic and thermodynamic properties of ascorbate oxidase (AO) activity and stability of a Cucurbita maxima extract were investigated. Activity tests performed at 25 degrees C using initial ascorbic acid concentration in the range 50-750 M allowed estimating the Michaelis constant for this substrate (Km = 126 microM) and the maximum initial rate of ascorbic acid oxidation (A0,max = 1.57 mM min-1). The main thermodynamic parameters of the enzyme reaction (DeltaH* = 10.3 kJ mol-1; DeltaG* = 87.2 kJ mol-1; DeltaS* = -258 J mol-1 K-1) were estimated through activity tests performed at 25-48 C. Within such a temperature range, no decrease in the initial reaction rate was detected. The long-term thermostability of the raw extract was then investigated by means of residual activity tests carried out at 10-70 degrees C, which allowed estimating the thermodynamic parameters of the irreversible enzyme inactivation as well (DeltaH*D = 51.7 kJ mol-1; DeltaG*D = 103 kJ mol-1; S*D = -160 J mol-1 K-1). Taking into account the specific rate of AO inactivation determined at different temperatures, we also estimated the enzyme half-life (1047 min at 10 degrees C and 21.2 min at 70 degrees C) and predicted the integral activity of a continuous system using this enzyme preparation. This work should be considered as a preliminary attempt to characterize the AO activity of a C. maxima extract before its concentration by liquid-liquid extraction techniques.

Ascorbate Oxidase↗
PubMedDate not supplied

Prediction of ligand-receptor binding thermodynamics by free energy force field (FEFF) 3D-QSAR analysis: application to a set of peptidometic renin inhibitors.

A methodology is presented and applied in which the accurate estimation of ligand-receptor binding thermodynamics is achieved by formulating the calculation as a QSAR problem. When the receptor geometry is known, the free energy force field (FEFF) ligand-receptor binding energy terms can be calculated and used as independent variables in constructing FEFF 3D-QSARs. The FEFF 3D-QSAR analysis of a series of transition state inhibitors of renin was carried out. From a statistical analysis of the free energy contributions to the binding process, FEFF 3D-QSARs were constructed that reveal the change in solvation free energy upon binding and the intramolecular vacuum internal energy of the ligand in the unbound state are the most significant FEFF terms in determining the binding free energy, delta G. Other terms, such as ligand stretching, bending, and torsion energy changes, the intermolecular van der Waals interaction energy, and change in ligand conformational entropy upon binding, are also found to make significant contributions in some FEFF 3D-QSAR delta G models and in delta H and delta S binding models. Overall, a relatively small number of the thermodynamic contributions to the ligand-receptor binding process dominates the thermodynamics of binding in a given model.

Amino Acid Sequence↗
PubMed2006

A stochastic regression approach to analyzing thermodynamic uncertainty in chemical speciation modeling.

Chemical speciation modeling is a vital tool for assessing the bioavailability of inorganic species, yet significant uncertainties in thermodynamic parameters and model form limit its potential for decision-making. In this paper we present a novel method for the quantification of thermodynamic parameter uncertainty and ionic strength correction model uncertainty using Bayesian Markov Chain Monte Carlo (MCMC) estimation methods. These methods allow for the inclusion of correlation modeling, which has not been present in previous work. The MCMC simulations are used to model a natural river water to determine the uncertainty in the calculated environmental speciation of ethylenediamenetetraacetate, a chelating agent that has attracted considerable environmental interest. The results indicate that incorporating correlation among related thermodynamic parameters into the uncertainty model is necessary to correctly quantify the overall system uncertainty. This result indicates the superiority of MCMC estimation methods overtraditional Monte Carlo methods when available data are used to estimate parameter uncertainty in systems with closely related model parameters.

Bayes Theorem↗
PubMed2001

Stereochemical, structural, and thermodynamic origins of stability differences between stereoisomeric benzo[a]pyrene diol epoxide deoxyadenosine adducts in a DNA mutational hot spot sequence.

Benzo[a]pyrene (BP), a prototype polycyclic aromatic hydrocarbon (PAH), can be metabolically activated to the enantiomeric benzo[a]pyrene diol epoxides (BPDEs), (+)-(7R,8S,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene and the (-)-(7S,8R,9R,10S) enantiomer. These can react with adenine residues in DNA, to produce the stereoisomeric 10S (+)- and 10R (-)-trans-anti-[BP]-N(6)-dA adducts. High-resolution NMR solution studies indicate that in DNA duplexes the 10R (-) adduct is intercalated on the 5'-side of the modified adenine, while the 10S (+) adduct is disordered, exhibits multiple adduct conformations, and is positioned on the 3'-side of the modified adenine. Duplexes containing the 10S (+) adduct positioned at A within codon 61 of the human N-ras sequence CAA are thermodynamically less stable and more easily excised by human DNA repair enzymes than those containing the 10R (-) adduct. However, the molecular origins of these differences are not understood and represent a fascinating opportunity for elucidating structure-function relationships. We have carried out a computational investigation to uncover the structural and thermodynamic origins of these effects in the 11-mer duplex sequence d(CGGACAAGAAG).d(CTTCTTGTCCG) by performing a 2-ns molecular dynamics simulation using NMR solution structures as the basis for the starting models. Then, we applied the MM-PBSA (molecular mechanics Poisson-Boltzmann surface area) method to compute free energy differences between the stereoisomeric adducts. The 10R (-) isomer is more stable by approximately 13 kcal/mol, of which approximately 10 kcal/mol is enthalpic, which agrees quite well with their observed differences in thermodynamic stability. The lower stability of the 10S (+) adduct is due to diminished stacking by the BP moiety in the intercalation pocket, more helix unwinding, and a diminished quality of Watson-Crick base pairing. The latter stems from conformational heterogeneity involving a syn-anti equilibrium of the glycosidic bond in the modified adenine residue. The lower stability and conformational heterogeneity of the 10S (+) adduct may play a role in its enhanced susceptibility to nucleotide excision repair.

7,8-Dihydro-7,8-dihydroxybenzo(a)pyrene 9,10-oxide↗
PubMed2003

Studies into the thermodynamic origin of negative cooperativity in ion-pairing molecular recognition.

The association of synthetic receptors to target guests often proceeds through the cooperative action of multiple binding forces. An investigation into the thermodynamic origin of cooperativity in ion-pairing host-guest binding in water is described. The binding affinities of 1,2,3,4-butanetetracarboxylate, tricarballate, glutarate, and acetate to a C(3)(v) symmetric metallo-host (1) are characterized in terms of the binding constants (K(a)) and the thermodynamic parameters deltaG degrees, deltaH degrees, and deltaS degrees, as determined by isothermal titration calorimetry (ITC). These values are used to determine the individual contributions of the binding interaction to the overall binding. Several ways to view the combination of the individual binding events that make up the whole are analyzed, all of which lead to the conclusion of negative cooperativity. Combined, the data were used to evaluate the thermodynamic origin of negative cooperativity for this series of guests, revealing that entropy is the largest contributing factor. An interpretation of this result focuses upon differences in the number of water molecules displaced upon binding.

Acetic Acid↗