Pathway Thermodynamics Highlights Kinetic Obstacles in Central Metabolism

In metabolism research, thermodynamics is usually used to determine the directionality of a reaction or the feasibility of a pathway. However, the relationship between thermodynamic potentials and fluxes is not limited to questions of directionality: thermodynamics also affects the kinetics of reactions through the flux-force relationship, which states that the logarithm of the ratio between the forward and reverse fluxes is directly proportional to the change in Gibbs energy due to a reaction (Δ_rG′). Accordingly, if an enzyme catalyzes a reaction with a Δ_rG′ of -5.7 kJ/mol then the forward flux will be roughly ten times the reverse flux. As Δ_rG′ approaches equilibrium (Δ_rG′ = 0 kJ/mol), exponentially more enzyme counterproductively catalyzes the reverse reaction, reducing the net rate at which the reaction proceeds. Thus, the enzyme level required to achieve a given flux increases dramatically near equilibrium. Here, we develop a framework for quantifying the degree to which pathways suffer these thermodynamic limitations on flux. For each pathway, we calculate a single thermodynamically-derived metric (the Max-min Driving Force, MDF), which enables objective ranking of pathways by the degree to which their flux is constrained by low thermodynamic driving force. Our framework accounts for the effect of pH, ionic strength and metabolite concentration ranges and allows us to quantify how alterations to the pathway structure affect the pathway''s thermodynamics. Applying this methodology to pathways of central metabolism sheds light on some of their features, including metabolic bypasses (e.g., fermentation pathways bypassing substrate-level phosphorylation), substrate channeling (e.g., of oxaloacetate from malate dehydrogenase to citrate synthase), and use of alternative cofactors (e.g., quinone as an electron acceptor instead of NAD). The methods presented here place another arrow in metabolic engineers'' quiver, providing a simple means of evaluating the thermodynamic and kinetic quality of different pathway chemistries that produce the same molecules.     

Kinetics and Thermodynamics of Reserpine Adsorption onto Strong Acidic Cationic Exchange Fiber

The kinetics and thermodynamics of the adsorption process of reserpine adsorbed onto the strong acidic cationic exchange fiber (SACEF) were studied by batch adsorption experiments. The adsorption capacity strongly depended on pH values, and the optimum reserpine adsorption onto the SACEF occurred at pH = 5 of reserpine solution. With the increase of temperature and initial concentration, the adsorption capacity increased. The equilibrium was attained within 20 mins. The adsorption process could be better described by the pseudo-second-order model and the Freundlich isotherm model. The calculated activation energy Ea was 4.35 kJ/mol. And the thermodynamic parameters were: 4.97<ΔH<7.44 kJ/mol, -15.29<ΔG<-11.87 kJ/mol and 41.97<ΔS<47.35 J/mol·K. The thermodynamic parameters demonstrated that the adsorption was an endothermic, spontaneous and feasible process of physisorption within the temperature range between 283 K and 323 K and the initial concentration range between 100 mg/L and 300 mg/L. All the results showed that the SACEF had a good adsorption performance for the adsorption of reserpine from alcoholic solution.     

Bioenergetic Conditions of Butyrate Metabolism by a Syntrophic, Anaerobic Bacterium in Coculture with Hydrogen-Oxidizing Methanogenic and Sulfidogenic Bacteria

The butyrate-oxidizing, proton-reducing, obligately anaerobic bacterium NSF-2 was grown in batch cocultures with either the hydrogen-oxidizing bacterium Methanospirillum hungatei PM-1 or Desulfovibrio sp. strain PS-1. Metabolism of butyrate occurred in two phases. The first phase exhibited exponential growth kinetics (phase a) and had a doubling time of 10 h. This value was independent of whether NSF-2 was cultured with a methanogen or a sulfate reducer and likely represents the maximum specific growth rate of NSF-2. This exponential growth phase was followed by a second phase with a nearly constant rate of degradation (phase b) which dominated the time course of butyrate degradation. The specific activity of H₂ uptake by the hydrogen-oxidizing bacterium controlled the bioenergetic conditions of metabolism in phase b. During this phase both the Gibbs free energy (ΔG′) and the butyrate degradation rate (v) were greater for NSF-2-Desulfovibrio sp. strain PS-1 (ΔG′ = −17.0 kJ/mol; v = 0.20 mM/h) than for NSF-2-M. hungatei PM-1 (ΔG′ = −3.8 kJ/mol, v = 0.12 mM/h). The ΔG′ value remained stable and characteristic of the two hydrogen oxidizers during phase b. The stable ΔG′ resulted from the close coupling of the rates of butyrate and H₂ oxidation. The addition of 2-bromoethanesulfonate to a NSF-2-methanogen coculture resulted in the total inhibition of butyrate degradation; the inhibition was relieved when Desulfovibrio sp. strain PS-1 was added as a new H₂ sink. When the specific activity of H₂ consumption was increased by adding higher densities of the Desulfovibrio sp. to 2-bromoethanesulfonate-inhibited NSF-2-methanogen cocultures, lower H₂ pool sizes and higher rates of butyrate degradation resulted. Thus, it is the kinetic parameters of H₂ consumption, not the type of H₂ consumer per se, that establishes the thermodynamic conditions which in turn control the rate of fatty acid degradation. The bioenergetic homeostasis we observed in phase b was a result of the kinetics of the coculture members and the feedback inhibition by hydrogen which prevents butyrate degradation rates from reaching their theoretical V_max.     

Melting studies of dangling-ended DNA hairpins: effects of end length,loop sequence and biotinylation of loop bases

The effects of 3′ single-strand dangling-ends of different lengths, sequence identity of hairpin loop, and hairpin loop biotinylation at different loop residues on DNA hairpin thermodynamic stability were investigated. Hairpins contained 16 bp stem regions and five base loops formed from the sequence, 5′-TAGTCGACGTGGTCC-N₅-GGACCACGTCGACTAG-E_n-3′. The length of the 3′ dangling-ends (E_n) was n = 13 or 22 bases. The identities of loop bases at positions 2 and 4 were varied. Biotinylation was varied at loop base positions 2, 3 or 4. Melting buffers contained 25 or 115 mM Na⁺. Average t_m values for all molecules were 73.5 and 84.0°C in 25 and 115 mM Na⁺, respectively. Average two-state parameters evaluated from van’t Hoff analysis of the melting curve shapes in 25 mM Na⁺ were ΔH_vH = 84.8 ± 15.5 kcal/mol, ΔS_vH = 244.8 ± 45.0 cal/K·mol and ΔG_vH = 11.9 ± 2.1 kcal/mol. In 115 mM Na⁺, two-state parameters were not very different at ΔH_vH = 80.42 ± 12.74 kcal/mol, ΔS_vH = 225.24 ± 35.88 cal/K·mol and ΔG_vH = 13.3 ± 2.0 kcal/mol. Differential scanning calorimetry (DSC) was performed to test the validity of the two-state assumption and evaluated van’t Hoff parameters. Thermodynamic parameters from DSC measurements (within experimental error) agreed with van’t Hoff parameters, consistent with a two-state process. Overall, dangling-end DNA hairpin stabilities are not affected by dangling-end length, loop biotinylation or sequence and vary uniformly with [Na⁺]. Consider able freedom is afforded when designing DNA hairpins as probes in nucleic acid based detection assays, such as microarrays.     

Pure-Culture Growth of Fermentative Bacteria, Facilitated by H2 Removal: Bioenergetics and H2 Production

We used an H₂-purging culture vessel to replace an H₂-consuming syntrophic partner, allowing the growth of pure cultures of Syntrophothermus lipocalidus on butyrate and Aminobacterium colombiense on alanine. By decoupling the syntrophic association, it was possible to manipulate and monitor the single organism's growth environment and determine the change in Gibbs free energy yield (ΔG) in response to changes in the concentrations of reactants and products, the purging rate, and the temperature. In each of these situations, H₂ production changed such that ΔG remained nearly constant for each organism (−11.1 ± 1.4 kJ mol butyrate⁻¹ for S. lipocalidus and −58.2 ± 1.0 kJ mol alanine⁻¹ for A. colombiense). The cellular maintenance energy, determined from the ΔG value and the hydrogen production rate at the point where the cell number was constant, was 4.6 × 10⁻¹³ kJ cell⁻¹ day⁻¹ for S. lipocalidus at 55°C and 6.2 × 10⁻¹³ kJ cell⁻¹ day⁻¹ for A. colombiense at 37°C. S. lipocalidus, in particular, seems adapted to thrive under conditions of low energy availability.     

Acetate Threshold Concentrations Suggest Varying Energy Requirements during Anaerobic Respiration by Anaeromyxobacter dehalogenans

Acetate threshold concentrations were determined under chlororespiring and Fe(III)-reducing conditions for Anaeromyxobacter dehalogenans strain 2CP-C. The acetate threshold concentrations measured were 69 ± 4, 19 ± 8, and <1 nM for chlororespiration, amorphous Fe(III) reduction, and Fe(III) citrate reduction, respectively. Residual ΔG values of −75.4 kJ/mol of electrons for chlororespiration and −41.5 kJ/mol of electrons for amorphous Fe(III) reduction were calculated at the acetate threshold concentration. By comparing threshold concentrations for different metabolisms in a single organism, this study provides insight into the metabolic use of energy under different growth conditions.     

Benzoate Fermentation by the Anaerobic Bacterium Syntrophus aciditrophicus in the Absence of Hydrogen-Using Microorganisms

The anaerobic bacterium Syntrophus aciditrophicus metabolized benzoate in pure culture in the absence of hydrogen-utilizing partners or terminal electron acceptors. The pure culture of S. aciditrophicus produced approximately 0.5 mol of cyclohexane carboxylate and 1.5 mol of acetate per mol of benzoate, while a coculture of S. aciditrophicus with the hydrogen-using methanogen Methanospirillum hungatei produced 3 mol of acetate and 0.75 mol of methane per mol of benzoate. The growth yield of the S. aciditrophicus pure culture was 6.9 g (dry weight) per mol of benzoate metabolized, whereas the growth yield of the S. aciditrophicus-M. hungatei coculture was 11.8 g (dry weight) per mol of benzoate. Cyclohexane carboxylate was metabolized by S. aciditrophicus only in a coculture with a hydrogen user and was not metabolized by S. aciditrophicus pure cultures. Cyclohex-1-ene carboxylate was incompletely degraded by S. aciditrophicus pure cultures until a free energy change (ΔG′) of −9.2 kJ/mol was reached (−4.7 kJ/mol for the hydrogen-producing reaction). Cyclohex-1-ene carboxylate, pimelate, and glutarate transiently accumulated at micromolar levels during growth of an S. aciditrophicus pure culture with benzoate. High hydrogen (10.1 kPa) and acetate (60 mM) levels inhibited benzoate metabolism by S. aciditrophicus pure cultures. These results suggest that benzoate fermentation by S. aciditrophicus in the absence of hydrogen users proceeds via a dismutation reaction in which the reducing equivalents produced during oxidation of one benzoate molecule to acetate and carbon dioxide are used to reduce another benzoate molecule to cyclohexane carboxylate, which is not metabolized further. Benzoate fermentation to acetate, CO₂, and cyclohexane carboxylate is thermodynamically favorable and can proceed at free energy values more positive than −20 kJ/mol, the postulated minimum free energy value for substrate metabolism.     

Immobilised Histidine Tagged β 2-Adrenoceptor Oriented by a Diazonium Salt Reaction and Its Application in Exploring Drug-Protein Interaction Using Ephedrine and Pseudoephedrine as Probes

A new oriented method using a diazonium salt reaction was developed for linking β ₂-adrenoceptor (β ₂-AR) on the surface of macroporous silica gel. Stationary phase containing the immobilised receptor was used to investigate the interaction between β ₂-AR and ephedrine plus pseudoephedrine by zonal elution. The isotherms of the two drugs best fit the Langmuir model. Only one type of binding site was found for ephedrine and pseudoephedrine targeting β ₂-AR. At 37 °C, the association constants during the binding were (5.94±0.05)×10³/M for ephedrine and (3.80±0.02) ×10³/M for pseudoephedrine, with the binding sites of (8.92±0.06) ×10⁻⁴ M. Thermodynamic studies showed that the binding of the two compounds to β ₂-AR was a spontaneous reaction with exothermal processes. The ΔG^θ, ΔH^θ and ΔS^θ for the interaction between ephedrine and β ₂-AR were −(22.33±0.04) kJ/mol, −(6.51±0.69) kJ/mol and 50.94±0.31 J/mol·K, respectively. For the binding of pseudoephedrine to the receptor, these values were −(21.17±0.02) kJ/mol, −(7.48±0.56) kJ/mol and 44.13±0.01 J/mol·K. Electrostatic interaction proved to be the driving force during the binding of the two drugs to β ₂-AR. The proposed immobilised method will have great potential for attaching protein to solid substrates and realizing the interactions between proteins and drugs.     

Characterization of the glutathione‐dependent reduction of the peroxiredoxin 5 homolog PfAOP from Plasmodium falciparum

Peroxiredoxins use a variety of thiols to rapidly reduce hydroperoxides and peroxynitrite. While the oxidation kinetics of peroxiredoxins have been studied in great detail, enzyme‐specific differences regarding peroxiredoxin reduction and the overall rate‐limiting step under physiological conditions often remain to be deciphered. The 1‐Cys peroxiredoxin 5 homolog PfAOP from the malaria parasite Plasmodium falciparum is an established model enzyme for glutathione/glutaredoxin‐dependent peroxiredoxins. Here, we reconstituted the catalytic cycle of PfAOP in vitro and analyzed the reaction between oxidized PfAOP and reduced glutathione (GSH) using molecular docking and stopped‐flow measurements. Molecular docking revealed that oxidized PfAOP has to adopt a locally unfolded conformation to react with GSH. Furthermore, we determined a second‐order rate constant of 6 × 10⁵ M⁻¹ s⁻¹ at 25°C and thermodynamic activation parameters ΔH ^‡, ΔS ^‡, and ΔG ^‡ of 39.8 kJ/mol, −0.8 J/mol, and 40.0 kJ/mol, respectively. The gain‐of‐function mutant PfAOP^L109M had almost identical reaction parameters. Taking into account physiological hydroperoxide and GSH concentrations, we suggest (a) that the reaction between oxidized PfAOP and GSH might be even faster than the formation of the sulfenic acid in vivo, and (b) that conformational changes are likely rate limiting for PfAOP catalysis. In summary, we characterized and quantified the reaction between GSH and the model enzyme PfAOP, thus providing detailed insights regarding the reactivity of its sulfenic acid and the versatile chemistry of peroxiredoxins.     

Removal and Recovery of Toxic Silver Ion Using Deep-Sea Bacterial Generated Biogenic Manganese Oxides

Products containing silver ion (Ag⁺) are widely used, leading to a large amount of Ag⁺-containing waste. The deep-sea manganese-oxidizing bacterium Marinobacter sp. MnI7-9 efficiently oxidizes Mn²⁺ to generate biogenic Mn oxide (BMO). The potential of BMO for recovering metal ions by adsorption has been investigated for some ions but not for Ag⁺. The main aim of this study was to develop effective methods for adsorbing and recovering Ag using BMO produced by Marinobacter sp. MnI7-9. In addition, the adsorption mechanism was determined using X-ray photoelectron spectroscopy analysis, specific surface area analysis, adsorption kinetics and thermodynamics. The results showed that BMO had a higher adsorption capacity for Ag⁺ compared to the chemical synthesized MnO₂ (CMO). The isothermal absorption curves of BMO and CMO both fit the Langmuir model well and the maximum adsorption capacities at 28°C were 8.097 mmol/g and 0.787 mmol/g, for BMO and CMO, respectively. The change in enthalpy (ΔH^θ) for BMO was 59.69 kJ/mol indicating that it acts primarily by chemical adsorption. The change in free energy (ΔG^θ) for BMO was negative, which suggests that the adsorption occurs spontaneously. Ag⁺ adsorption by BMO was driven by entropy based on the positive ΔS^θ values. The Ag⁺ adsorption kinetics by BMO fit the pseudo-second order model and the apparent activation energy of E_a is 21.72 kJ/mol. X-ray photoelectron spectroscopy analysis showed that 15.29% Ag⁺ adsorbed by BMO was transferred to Ag(0) and meant that redox reaction had happened during the adsorption. Desorption using nitric acid and Na₂S completely recovered the Ag. The results show that BMO produced by strain MnI7-9 has potential for bioremediation and reutilization of Ag⁺-containing waste.     

The First Identification of Carbohydrate Binding Modules Specific to Chitosan

Two carbohydrate binding modules (DD1 and DD2) belonging to CBM32 are located at the C terminus of a chitosanase from Paenibacillus sp. IK-5. We produced three proteins, DD1, DD2, and tandem DD1/DD2 (DD1+DD2), and characterized their binding ability. Transition temperature of thermal unfolding (T_m) of each protein was elevated by the addition of cello-, laminari-, chitin-, or chitosan-hexamer (GlcN)₆. The T_m elevation (ΔT_m) in DD1 was the highest (10.3 °C) upon the addition of (GlcN)₆ and was markedly higher than that in DD2 (1.0 °C). A synergistic effect was observed (ΔT_m = 13.6 °C), when (GlcN)₆ was added to DD1+DD2. From isothermal titration calorimetry experiments, affinities to DD1 were not clearly dependent upon chain length of (GlcN)_n; ΔG_r° values were −7.8 (n = 6), −7.6 (n = 5), −7.6 (n = 4), −7.6 (n = 3), and −7.1 (n = 2) kcal/mol, and the value was not obtained for GlcN due to the lowest affinity. DD2 bound (GlcN)_n with the lower affinities (ΔG_r° = −5.0 (n = 3) ∼ −5.2 (n = 6) kcal/mol). Isothermal titration calorimetry profiles obtained for DD1+DD2 exhibited a better fit when the two-site model was used for analysis and provided greater affinities to (GlcN)₆ for individual DD1 and DD2 sites (ΔG_r° = −8.6 and −6.4 kcal/mol, respectively). From NMR titration experiments, (GlcN)_n (n = 2∼6) were found to bind to loops extruded from the core β-sandwich of individual DD1 and DD2, and the interaction sites were similar to each other. Taken together, DD1+DD2 is specific to chitosan, and individual modules synergistically interact with at least two GlcN units, facilitating chitosan hydrolysis.     

Thermodynamic criteria for high hit rate antisense oligonucleotide design

Antisense oligonucleotides are used for therapeutic applications and in functional genomic studies. In practice, however, many of the oligonucleotides complementary to an mRNA have little or no antisense activity. Theoretical strategies to improve the ‘hit rate’ in antisense screens will reduce the cost of discovery and may lead to identification of antisense oligonucleotides with increased potency. Statistical analysis performed on data collected from more than 1000 experiments with phosphorothioate-modified oligonucleotides revealed that the oligo-probes, which form stable duplexes with RNA (ΔG^o₃₇ ≤ –30 kcal/mol) and have small self-interaction potential, are more frequently efficient than molecules that form less stable oligonucleotide–RNA hybrids or more stable self-structures. To achieve optimal statistical preference, the values for self-interaction should be (ΔG^o₃₇) ≥ –8 kcal/mol for inter-oligonucleotide pairing and (ΔG^o₃₇) ≥ –1.1 kcal/mol for intra-molecular pairing. Selection of oligonucleotides with these thermodynamic values in the analyzed experiments would have increased the ‘hit rate’ by as much as 6-fold.     

Determination of base and backbone contributions to the thermodynamics of premelting and melting transitions in B DNA

In previous papers of this series the temperature-dependent Raman spectra of poly(dA)·poly(dT) and poly(dA–dT)·poly(dA–dT) were used to characterize structurally the melting and premelting transitions in DNAs containing consecutive A·T and alternating A·T/T·A base pairs. Here, we describe procedures for obtaining thermodynamic parameters from the Raman data. The method exploits base-specific and backbone-specific Raman markers to determine separate thermodynamic contributions of A, T and deoxyribosyl-phosphate moieties to premelting and melting transitions. Key findings include the following: (i) Both poly(dA)·poly(dT) and poly(dA–dT)· poly(dA–dT) exhibit robust premelting transitions, due predominantly to backbone conformational changes. (ii) The significant van’t Hoff premelting enthalpies of poly(dA)·poly(dT) [ΔH_vH^pm = 18.0 ± 1.6 kcal·mol^–1 (kilocalories per mole cooperative unit)] and poly(dA–dT)·poly(dA–dT) (ΔH_vH^pm = 13.4 ± 2.5 kcal·mol^–1) differ by an amount (~4.6 kcal·mol^–1) estimated as the contribution from three-centered inter-base hydrogen bonding in (dA)_n·(dT)_n tracts. (iii) The overall stacking free energy of poly(dA)· poly(dT) [–6.88 kcal·mol_bp^–1 (kilocalories per mole base pair)] is greater than that of poly(dA–dT)· poly(dA–dT) (–6.31 kcal·mol_bp^–1). (iv) The difference between stacking free energies of A and T is significant in poly(dA)·poly(dT) (ΔΔG_st = 0.8 ± 0.3 kcal· mol_bp^–1), but marginal in poly(dA–dT)·poly(dA–dT) (ΔΔG_st = 0.3 ± 0.3 kcal·mol_bp^–1). (v) In poly(dA)· poly(dT), the van’t Hoff parameters for melting of A (ΔH_vH^A = 407 ± 23 kcal·mol^–1, ΔS_vH^A = 1166 ± 67 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 60.0 ± 3.2 kcal·mol^–1) are clearly distinguished from those of T (ΔH_vH^T = 185 ± 38 kcal·mol^–1, ΔS_vH^T = 516 ± 109 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 27.1 ± 5.5 kcal·mol^–1). (vi) Similar relative differences are observed in poly(dA–dT)· poly(dA–dT) (ΔH_vH^A = 333 ± 54 kcal·mol^–1, ΔS_vH^A = 961 ± 157 cal·°K^–1·mol^–1, ΔG_vH(25°C)^A = 45.0 ± 7.6 kcal· mol^–1; ΔH_vH^T = 213 ± 30 kcal·mol^–1, ΔS_vH^T = 617 ± 86 cal·°K^–1·mol^–1, ΔG_vH(25°C)^T = 29.3 ± 4.9 kcal·mol^–1). The methodology employed here distinguishes thermodynamic contributions of base stacking, base pairing and backbone conformational ordering in the molecular mechanism of double-helical B DNA formation.     

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ∼16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm_app) by 1–6°C resulting in an increase in free energy of unfolding (ΔG_unf) by 0.04–1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG_unf, (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG_unf = 1.57 kcal⁻¹) justifies the robustness of this relation.     

Kinetics of selenite uptake by mononuclear cells from peripheral human blood

The present study deals with the kinetics and thermodynamics of the uptake of⁷⁵Se-labeled SeO ₃ ^2? from incubation media to lymphocytes cultivated from eight normal individuals (14–55 years of age, two females). The uptake of SeO ₃ ^2? was evaluated on the assumption of pseudo-first-order kinetics with regard to a reacting cellular receptor pool. On the basis of the experimental observations, it was assumed that the suggested pool of receptor molecules-symbolically represented by “£H₄”—reacts with SeO ₃ ^2? in the hypothetical reaction: $$\pounds H_4 + SeO_3^{2 - } + 2H^ + \underset{{ - k_1 }}{\overset{{k_1 }}{\longleftrightarrow}}\pounds Se + 3H_2 O$$ The mean value of the change in standard free energy at 25°C was calculated to be ΔG ^o=?141.6±1.3 kJ/mol, while the corresponding mean value of the free energy of activation at 25°C was calculated to be ΔG ²⁺=?7.8±0.9 kJ/mol for the forward reaction. The calculated values of the corresponding individual changes in the respective standard enthalpies and entropies were mutually interdependent for all eight donors. ΔH ^o=?152+315ΔS ^o(kJ/mol) corresponding to the common value ΔG ^o??152 kJ/mol at 315°K. These mutual interdependencies are possibly the effect of variable conformational states (e.g., the macromolecular compactness) of the cellular receptor pools. This suggestion may furthermore be supported by the correlation traced between ΔH ^o vs the biological age in years of the donors: △H ^°?76.7?1.0 (age)kJ/mol (r = ?0.92) The calculated values of activation enthalpy ΔH ²⁺ kJ/mol and activation entropy ΔS ²⁺ (kJ/mol K) also mutually correlated linearly (r=0.998); the regression line was: △H ²⁺ = ?8.9 + 305△S²⁺ (kJ/mol) corresponding to the common value △H ²⁺ △ ?8.9 (kJ/mol) at 305°K Similarly the activation enthalpy ΔH ²⁺ vs the biological age in years correlated linearly: ΔH ²⁺=67.4?0.73(age) (kJ/mol) (r=?0.76) The range of ΔH ²⁺ studied was from 13.8 to 53.9 kJ/mol with a linearly corresponding range in ΔS ²⁺ from 73 to 205 J/mol K. The thermodynamic data reveal the selenite uptake during the hypothetical standard reaction to be exergonic and endothermic. Critical pH dependencies of the selenite uptake were explained.     

Structural and Thermodynamic Approach to Peptide Immunogenicity

In the conventional paradigm of humoral immunity, B cells recognize their cognate antigen target in its native form. However, it is well known that relatively unstable peptides bearing only partial structural resemblance to the native protein can trigger antibodies recognizing higher-order structures found in the native protein. On the basis of sound thermodynamic principles, this work reveals that stability of immunogenic proteinlike motifs is a critical parameter rationalizing the diverse humoral immune responses induced by different linear peptide epitopes. In this paradigm, peptides with a minimal amount of stability (ΔG_X<0 kcal/mol) around a proteinlike motif (X) are capable of inducing antibodies with similar affinity for both peptide and native protein, more weakly stable peptides (ΔG_X>0 kcal/mol) trigger antibodies recognizing full protein but not peptide, and unstable peptides (ΔG_X>8 kcal/mol) fail to generate antibodies against either peptide or protein. Immunization experiments involving peptides derived from the autoantigen histidyl-tRNA synthetase verify that selected peptides with varying relative stabilities predicted by molecular dynamics simulations induce antibody responses consistent with this theory. Collectively, these studies provide insight pertinent to the structural basis of immunogenicity and, at the same time, validate this form of thermodynamic and molecular modeling as an approach to probe the development/evolution of humoral immune responses.     

Influence of substituent modifications on the binding of 2-amino-1,8-naphthyridines to cytosine opposite an AP site in DNA duplexes: thermodynamic characterization

Here, we report on a significant effect of substitutions on the binding affinity of a series of 2-amino-1,8-naphthyridines, i.e., 2-amino-1,8-naphthyridine (AND), 2-amino-7-methyl-1,8-naphthyridine (AMND), 2-amino-5,7-dimethyl-1,8-naphthyridine (ADMND) and 2-amino-5,6,7-trimethyl-1,8-naphthyridine (ATMND), all of which can bind to cytosine opposite an AP site in DNA duplexes. Fluorescence titration experiments show that the binding affinity for cytosine is effectively enhanced by the introduction of methyl groups to the naphthyridine ring, and the 1:1 binding constant (10⁶ M⁻¹) follows in the order of AND (0.30) < AMND (2.7) < ADMND (6.1) < ATMND (19) in solutions containing 110 mM Na⁺ (pH 7.0, at 20°C). The thermodynamic parameters obtained by isothermal titration calorimetry experiments indicate that the introduction of methyl groups effectively reduces the loss of binding entropy, which is indeed responsible for the increase in the binding affinity. The heat capacity change (ΔC_p), as determined from temperature dependence of the binding enthalpy, is found to be significantly different between AND (−161 cal/mol K) and ATMND (−217 cal/mol K). The hydrophobic contribution appears to be a key force to explain the observed effect of substitutions on the binding affinity when the observed binding free energy (ΔG_obs) is dissected into its component terms.     

Energetic signatures of single base bulges: thermodynamic consequences and biological implications

DNA bulges are biologically consequential defects that can arise from template-primer misalignments during replication and pose challenges to the cellular DNA repair machinery. Calorimetric and spectroscopic characterizations of defect-containing duplexes reveal systematic patterns of sequence-context dependent bulge-induced destabilizations. These distinguishing energetic signatures are manifest in three coupled characteristics, namely: the magnitude of the bulge-induced duplex destabilization (ΔΔG_Bulge); the thermodynamic origins of ΔΔG_Bulge (i.e. enthalpic versus entropic); and, the cooperativity of the duplex melting transition (i.e. two-state versus non-two state). We find moderately destabilized duplexes undergo two-state dissociation and exhibit ΔΔG_Bulge values consistent with localized, nearest neighbor perturbations arising from unfavorable entropic contributions. Conversely, strongly destabilized duplexes melt in a non-two-state manner and exhibit ΔΔG_Bulge values consistent with perturbations exceeding nearest-neighbor expectations that are enthalpic in origin. Significantly, our data reveal an intriguing correlation in which the energetic impact of a single bulge base centered in one strand portends the impact of the corresponding complementary bulge base embedded in the opposite strand. We discuss potential correlations between these bulge-specific differential energetic profiles and their overall biological implications in terms of DNA recognition, repair and replication.     

Action of nicotine and analogs on acetylcholine receptors having mutations of transmitter-binding site residue αG153

A primary target for nicotine is the acetylcholine receptor channel (AChR). Some of the ability of nicotine to activate differentially AChR subtypes has been traced to a transmitter-binding site amino acid that is glycine in lower affinity and lysine in higher affinity AChRs. We studied the effects of mutations of this residue (αG153) in neuromuscular AChRs activated by nicotine and eight other agonists including nornicotine and anabasine. All of the mutations increased the unliganded gating equilibrium constant. The affinity of the resting receptor (K_d) and the net binding energy from the agonist for gating (ΔG_B) were estimated by cross-concentration fitting of single-channel currents. In all but one of the agonist/mutant combinations there was a moderate decrease in K_d and essentially no change in ΔG_B. The exceptional case was nicotine plus lysine, which showed a large, >8,000-fold decrease in K_d but no change in ΔG_B. The extraordinary specificity of this combination leads us to speculate that AChRs with a lysine at position αG153 may be exposed to a nicotine-like compound in vivo.     

Hydrogen Partial Pressures in a Thermophilic Acetate-Oxidizing Methanogenic Coculture

Hydrogen partial pressures were measured in a thermophilic coculture comprised of a eubacterial rod which oxidized acetate to H₂ and CO₂ and a hydrogenotrophic methanogen, Methanobacterium sp. strain THF. Zinder and Koch (S. H. Zinder and M. Koch, Arch. Microbiol. 138:263-272, 1984) originally predicted, on the basis of calculations of Gibbs free energies of reactions, that the H₂ partial pressure near the midpoint of growth of the coculture should be near 4 Pa (ca. 4 × 10⁻⁵ atm; ca. 0.024 μM dissolved H₂) for both organisms to be able to conserve energy for growth. H₂ partial pressures in the coculture were measured to be between 20 and 50 Pa (0.12 to 0.30 μM) during acetate utilization, approximately one order of magnitude higher than originally predicted. However, when ΔG_f (free energy of formation) values were corrected for 60°C by using the relationship ΔG_f = ΔH_f − TΔS (ΔH_f is the enthalpy or heat of formation, ΔS is the entropy value, and T is the temperature in kelvins), the predicted value was near 15 Pa, in closer agreement with the experimentally determined values. The coculture also oxidized ethanol to acetate, a more thermodynamically favorable reaction than oxidation of acetate to CO₂. During ethanol oxidation, the H₂ partial pressure reached values as high as 200 Pa. Acetate was not used until after the ethanol was consumed and the H₂ partial pressure decreased to 40 to 50 Pa. After acetate utilization, H₂ partial pressures fell to approximately 10 Pa and remained there, indicating a threshold for H₂ utilization by the methanogen. Axenic cultures of the acetate-oxidizing organism were combined with pure cultures of either Methanobacterium sp. strain THF or Methanobacterium thermoautotrophicum ΔH to form reconstituted acetate-oxidizing cocultures. The H₂ partial pressures measured in both of these reconstituted cocultures were similar to those measured in the original acetate-oxidizing rod coculture. Since M. thermoautotrophicum ΔH did not use formate as a substrate, formate is not necessarily involved in interspecies electron transfer in this coculture.