Predicting Thermodynamic Properties of PBXTHs with New Quantum Topological Indexes

Novel group quantitative structure-property relationship (QSPR) models on the thermodynamic properties of PBXTHs were presented, by the multiple linear regression (MLR) analysis method. Four thermodynamic properties were studied: the entropy (S^θ), the standard enthalpy of formation (Δ_fH^θ), the standard Gibbs energy of formation (Δ_fG^θ), and the relative standard Gibbs energy of formation (Δ_RG^θ). The results by the formula indicate that the calculated and predicted data in this study are in good agreement with those in literature and the deviation is within the experimental errors. To validate the estimation reliability for internal samples and the predictive ability for other samples, leave-one-out (LOO) cross validation (CV) and external validation were performed, and the results show that the models are satisfactory.     

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.     

Characterization of the kinetic and thermodynamic landscape of RNA folding using a novel application of isothermal titration calorimetry

A novel isothermal titration calorimetry (ITC) method was applied to investigate RNA helical packing driven by the GAAA tetraloop–receptor interaction in magnesium and potassium solutions. Both the kinetics and thermodynamics were obtained in individual ITC experiments, and analysis of the kinetic data over a range of temperatures provided Arrhenius activation energies (ΔH^‡) and Eyring transition state entropies (ΔS^‡). The resulting rich dataset reveals strongly contrasting kinetic and thermodynamic profiles for this RNA folding system when stabilized by potassium versus magnesium. In potassium, association is highly exothermic (ΔH_25°C = −41.6 ± 1.2 kcal/mol in 150 mM KCl) and the transition state is enthalpically barrierless (ΔH^‡ = −0.6 ± 0.5). These parameters are sigificantly positively shifted in magnesium (ΔH_25°C = −20.5 ± 2.1 kcal/mol, ΔH^‡ = 7.3 ± 2.2 kcal/mol in 0.5 mM MgCl₂). Mixed salt solutions approximating physiological conditions exhibit an intermediate thermodynamic character. The cation-dependent thermodynamic landscape may reflect either a salt-dependent unbound receptor conformation, or alternatively and more generally, it may reflect a small per-cation enthalpic penalty associated with folding-coupled magnesium uptake.     

Thermodynamic analysis of F1‐ATPase rotary catalysis using high‐speed imaging

F₁-ATPase (F₁) is a rotary motor protein fueled by ATP hydrolysis. Although the mechanism for coupling rotation and catalysis has been well studied, the molecular details of individual reaction steps remain elusive. In this study, we performed high-speed imaging of F₁ rotation at various temperatures using the total internal reflection dark-field (TIRDF) illumination system, which allows resolution of the F₁ catalytic reaction into elementary reaction steps with a high temporal resolution of 72 µs. At a high concentration of ATP, F₁ rotation comprised distinct 80° and 40° substeps. The 80° substep, which exhibited significant temperature dependence, is triggered by the temperature-sensitive reaction, whereas the 40° substep is triggered by ATP hydrolysis and the release of inorganic phosphate (P_i). Then, we conducted Arrhenius analysis of the reaction rates to obtain the thermodynamic parameters for individual reaction steps, that is, ATP binding, ATP hydrolysis, P_i release, and TS reaction. Although all reaction steps exhibited similar activation free energy values, ΔG^‡ = 53–56 kJ mol⁻¹, the contributions of the enthalpy (ΔH^‡), and entropy (ΔS^‡) terms were significantly different; the reaction steps that induce tight subunit packing, for example, ATP binding and TS reaction, showed high positive values of both ΔH^‡ and ΔS^‡. The results may reflect modulation of the excluded volume as a function of subunit packing tightness at individual reaction steps, leading to a gain or loss in water entropy.     

Determination of the secondary structure of group II bulge loops using the fluorescent probe 2-aminopurine

Eleven RNA hairpins containing 2-aminopurine (2-AP) in either base-paired or single nucleotide bulge loop positions were optically melted in 1 M NaCl; and, the thermodynamic parameters ΔH°, ΔS°, ΔG°₃₇, and T_M for each hairpin were determined. Substitution of 2-AP for an A (adenosine) at a bulge position (where either the 2-AP or A is the bulge) in the stem of a hairpin, does not affect the stability of the hairpin. For group II bulge loops such as AA/U, where there is ambiguity as to which of the A residues is paired with the U, hairpins with 2-AP substituted for either the 5′ or 3′ position in the hairpin stem have similar stability. Fluorescent melts were performed to monitor the environment of the 2-AP. When the 2-AP was located distal to the hairpin loop on either the 5′ or 3′ side of the hairpin stem, the change in fluorescent intensity upon heating was indicative of an unpaired nucleotide. A database of phylogenetically determined RNA secondary structures was examined to explore the presence of naturally occurring bulge loops embedded within a hairpin stem. The distribution of bulge loops is discussed and related to the stability of hairpin structures.     

Specificity of RSG-1.2 peptide binding to RRE-IIB RNA element of HIV-1 over Rev peptide is mainly enthalpic in origin

Rev is an essential HIV-1 regulatory protein which binds to the Rev responsive element (RRE) present within the env gene of HIV-1 RNA genome. This binding facilitates the transport of the RNA to the cytoplasm, which in turn triggers the switch between viral latency and active viral replication. Essential components of this complex have been localized to a minimal arginine rich Rev peptide and stem IIB region of RRE. A synthetic peptide known as RSG-1.2 binds with high binding affinity and specificity to the RRE-IIB than the Rev peptide, however the thermodynamic basis of this specificity has not yet been addressed. The present study aims to probe the thermodynamic origin of this specificity of RSG-1.2 over Rev Peptide for RRE-IIB. The temperature dependent melting studies show that RSG-1.2 binding stabilizes the RRE structure significantly (ΔT _m = 4.3°C), in contrast to Rev binding. Interestingly the thermodynamic signatures of the binding have also been found to be different for both the peptides. At pH 7.5, RSG-1.2 binds RRE-IIB with a K_a = 16.2±0.6×10⁷ M⁻¹ where enthalpic change ΔH = −13.9±0.1 kcal/mol is the main driving force with limited unfavorable contribution from entropic change TΔS = −2.8±0.1 kcal/mol. A large part of ΔH may be due to specific stacking between U72 and Arg15. In contrast binding of Rev (K_a = 3.1±0.4×10⁷ M⁻¹) is driven mainly by entropy (ΔH = 0 kcal/mol and TΔS = 10.2±0.2 kcal/mol) which arises from major conformational changes in the RNA upon binding.     

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.     

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.     

Group contribution method for thermodynamic analysis of complex metabolic networks

A new, to our knowledge, group contribution method based on the group contribution method of Mavrovouniotis is introduced for estimating the standard Gibbs free energy of formation (Δ_fG′°) and reaction (Δ_rG′°) in biochemical systems. Gibbs free energy contribution values were estimated for 74 distinct molecular substructures and 11 interaction factors using multiple linear regression against a training set of 645 reactions and 224 compounds. The standard error for the fitted values was 1.90 kcal/mol. Cross-validation analysis was utilized to determine the accuracy of the methodology in estimating Δ_rG′° and Δ_fG′° for reactions and compounds not included in the training set, and based on the results of the cross-validation, the standard error involved in these estimations is 2.22 kcal/mol. This group contribution method is demonstrated to be capable of estimating Δ_rG′° and Δ_fG′° for the majority of the biochemical compounds and reactions found in the iJR904 and iAF1260 genome-scale metabolic models of Escherichia coli and in the Kyoto Encyclopedia of Genes and Genomes and University of Minnesota Biocatalysis and Biodegradation Database. A web-based implementation of this new group contribution method is available free at http://sparta.chem-eng.northwestern.edu/cgi-bin/GCM/WebGCM.cgi.     

Engineered Bi-Histidine Metal Chelation Sites Map the Structure of the Mechanical Unfolding Transition State of an Elastomeric Protein Domain GB1

Determining the structure of the transition state is critical for elucidating the mechanism behind how proteins fold and unfold. Due to its high free energy, however, the transition state generally cannot be trapped and studied directly using traditional structural biology methods. Thus, characterizing the structure of the transition state that occurs as proteins fold and unfold remains a major challenge. Here, we report a novel (to our knowledge) method that uses engineered bi-histidine (bi-His) metal-binding sites to directly map the structure of the mechanical unfolding transition state of proteins. This method is adapted from the traditional ψ-value analysis, which uses engineered bi-His metal chelation sites to probe chemical (un)folding transition-state structure. The ϕ^M2+_U-value is defined as ΔΔG_‡-N/ΔΔG_U-N, which is the energetic effects of metal chelation by the bi-His site on the unfolding energy barrier (ΔG_‡-N) relative to its thermodynamic stability (ΔG_U-N) and can be used to obtain information about the transition state in the mutational site. As a proof of principle, we used the small protein GB1 as a model system and set out to map its mechanical unfolding transition-state structure. Using single-molecule atomic force microscopy and spectrofluorimetry, we directly quantified the effect of divalent metal ion binding on the mechanical unfolding free energy and thermodynamic stability of GB1, which allowed us to quantify ϕ^M2+_U-values for different sites in GB1. Our results enabled us to map the structure of the mechanical unfolding transition state of GB1. Within GB1’s mechanical unfolding transition state, the interface between force-bearing β-strands 1 and 4 is largely disrupted, and the first β-hairpin is partially disordered while the second β-hairpin and the α-helix remain structured. Our results demonstrate the unique application of ψ-value analysis in elucidating the structure of the transition state that occurs during the mechanical unfolding process, offering a potentially powerful new method for investigating the design of novel elastomeric proteins.     

Distal Histidine Stabilizes Bound O2 and Acts as a Gate for Ligand Entry in Both Subunits of Adult Human Hemoglobin

The role of the distal histidine in regulating ligand binding to adult human hemoglobin (HbA) was re-examined systematically by preparing His(E7) to Gly, Ala, Leu, Gln, Phe, and Trp mutants of both Hb subunits. Rate constants for O₂, CO, and NO binding were measured using rapid mixing and laser photolysis experiments designed to minimize autoxidation of the unstable apolar E7 mutants. Replacing His(E7) with Gly, Ala, Leu, or Phe causes 20–500-fold increases in the rates of O₂ dissociation from either Hb subunit, demonstrating unambiguously that the native His(E7) imidazole side chain forms a strong hydrogen bond with bound O₂ in both the α and β chains (ΔG_{His(E7)H-bond} ≈ −8 kJ/mol). As the size of the E7 amino acid is increased from Gly to Phe, decreases in k_O2′, k_NO′, and calculated bimolecular rates of CO entry (k_entry′) are observed. Replacing His(E7) with Trp causes further decreases in k_O2′, k_NO′, and k_entry′ to 1–2 μm⁻¹ s⁻¹ in β subunits, whereas ligand rebinding to αTrp(E7) subunits after photolysis is markedly biphasic, with fast k_O2′, k_CO′, and k_NO′ values ≈150 μm⁻¹ s⁻¹ and slow rate constants ≈0.1 to 1 μm⁻¹ s⁻¹. Rapid bimolecular rebinding to an open α subunit conformation occurs immediately after photolysis of the αTrp(E7) mutant at high ligand concentrations. However, at equilibrium the closed αTrp(E7) side chain inhibits the rate of ligand binding >200-fold. These data suggest strongly that the E7 side chain functions as a gate for ligand entry in both HbA subunits.     

Thermodynamic and kinetic basis for recognition and repair of 8-oxoguanine in DNA by human 8-oxoguanine-DNA glycosylase

We have used a stepwise increase in ligand complexity approach to estimate the relative contributions of the nucleotide units of DNA containing 7,8-dihydro-8-oxoguanine (oxoG) to its total affinity for human 8-oxoguanine DNA glycosylase (OGG1) and construct thermodynamic models of the enzyme interaction with cognate and non-cognate DNA. Non-specific OGG1 interactions with 10–13 nt pairs within its DNA-binding cleft provides approximately 5 orders of magnitude of its affinity for DNA (ΔG° approximately −6.7 kcal/mol). The relative contribution of the oxoG unit of DNA (ΔG° approximately −3.3 kcal/mol) together with other specific interactions (ΔG° approximately −0.7 kcal/mol) provide approximately 3 orders of magnitude of the affinity. Formation of the Michaelis complex of OGG1 with the cognate DNA cannot account for the major part of the enzyme specificity, which lies in the k_cat term instead; the rate increases by 6–7 orders of magnitude for cognate DNA as compared with non-cognate one. The k_cat values for substrates of different sequences correlate with the DNA twist, while the K_M values correlate with ΔG° of the DNA fragments surrounding the lesion (position from −6 to +6). The functions for predicting the K_M and k_cat values for different sequences containing oxoG were found.     

Mechanism of Inhibition by C-terminal α-Helices of the ? Subunit of Escherichia coli FoF1-ATP Synthase

The ϵ subunit of bacterial F_oF₁-ATP synthase (F_oF₁), a rotary motor protein, is known to inhibit the ATP hydrolysis reaction of this enzyme. The inhibitory effect is modulated by the conformation of the C-terminal α-helices of ϵ, and the “extended” but not “hairpin-folded” state is responsible for inhibition. Although the inhibition of ATP hydrolysis by the C-terminal domain of ϵ has been extensively studied, the effect on ATP synthesis is not fully understood. In this study, we generated an Escherichia coli F_oF₁ (EF_oF₁) mutant in which the ϵ subunit lacked the C-terminal domain (F_oF₁^ϵΔC), and ATP synthesis driven by acid-base transition (ΔpH) and the K⁺-valinomycin diffusion potential (ΔΨ) was compared in detail with that of the wild-type enzyme (F_oF₁^ϵWT). The turnover numbers (k_cat) of F_oF₁^ϵWT were severalfold lower than those of F_oF₁^ϵΔC. F_oF₁^ϵWT showed higher Michaelis constants (K_m). The dependence of the activities of F_oF₁^ϵWT and F_oF₁^ϵΔC on various combinations of ΔpH and ΔΨ was similar, suggesting that the rate-limiting step in ATP synthesis was unaltered by the C-terminal domain of ϵ. Solubilized F_oF₁^ϵWT also showed lower k_cat and higher K_m values for ATP hydrolysis than the corresponding values of F_oF₁^ϵΔC. These results suggest that the C-terminal domain of the ϵ subunit of EF_oF₁ slows multiple elementary steps in both the ATP synthesis/hydrolysis reactions by restricting the rotation of the γ subunit.     

Molecular-Level Thermodynamic Switch Controls Chemical Equilibrium in Sequence-Specific Hydrophobic Interaction of 35 Dipeptide Pairs

Applying the Planck-Benzinger methodology, the sequence-specific hydrophobic interactions of 35 dipeptide pairs were examined over a temperature range of 273–333 K, based on data reported by Nemethy and Scheraga in 1962. The hydrophobic interaction in these sequence-specific dipeptide pairs is highly similar in its thermodynamic behavior to that of other biological systems. The results imply that the negative Gibbs free energy change minimum at a well-defined stable temperature, T_s, where the bound unavailable energy, TΔS^o = 0, has its origin in the sequence-specific hydrophobic interactions, are highly dependent on details of molecular structure. Each case confirms the existence of a thermodynamic molecular switch wherein a change of sign in ΔCp^o(T)_reaction (change in specific heat capacity of reaction at constant pressure) leads to true negative minimum in the Gibbs free energy change of reaction, ΔG^o(T)_reaction, and hence a maximum in the related equilibrium constant, K_eq. Indeed, all interacting biological systems examined to date by Chun using the Planck-Benzinger methodology have shown such a thermodynamic switch at the molecular level, suggesting its existence may be universal.     

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

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

M?ssbauer spectroscopy of the nitrogenase proteins from Klebsiella pneumoniae. Structural assignments and mechanistic conclusions

The Mo–Fe protein and the Fe protein which together constitute the nitrogenase of Klebsiella pneumoniae were prepared from bacteria grown in ⁵⁷Fe-enriched medium. The Mössbauer spectrum of the Mo–Fe protein, as isolated in the presence of Na₂S₂O₄, showed that the protein contained three iron species, called M4, M5 and M6. The area of the spectrum associated with species M4, with δ=0.65mm/s and ΔE=3.05mm/s at 4.2°K, corresponded to two iron atoms/molecule of protein and it is interpreted as being due to a high-spin ferrous, spin-coupled pair of iron atoms. The iron atoms of species M4 may be involved in the quaternary structure of the protein. Species M5, with δ=0.61mm/s and ΔE=0.83mm/s at 77°K, corresponded to eight iron atoms/molecule of protein and is interpreted as being due to Fe₄S₄ or Fe₂S₂ low-spin ferrous iron clusters. Species M6, with δ=0.37mm/s and ΔE=0.71mm/s at 77°K, also corresponded to eight iron atoms/molecule of protein and, at 4.2°K, became a broad shallow absorption, characteristic of magnetic hyperfine interaction. Oxidation of the Mo–Fe protein with the redox dye Lauth''s Violet did not affect the activity of the protein but changed species M4, M5 and M6 into the species M1 (δ=0.37mm/s, ΔE=0.75mm/s at 77°K, broad magnetic component at 4.2°K) and M2 (δ=0.35mm/s, ΔE=0.9mm/s at 4.2°K). In the presence of the Fe protein, Na₂S₂O₄, ATP and Mg²⁺, the M6 component of the Mo–Fe protein was replaced by species M7 with δ=0.46mm/s, ΔE=1.04mm/s at 4.2°K. The change in Mössbauer parameters associated with the M6 → M7 transformation was very similar to the change observed on reduction of the high-potential Fe protein from Chromatium vinosum. In contrast, Na₂S₂O₄-reduced Fe protein contained only one type of iron cluster (F4). Species F4 had δ=0.50mm/s, ΔE=0.9mm/s at 195°K, and at 4.2°K broadened in a manner characteristic of a magnetic hyperfine interaction, associated with half-integral spin, equally distributed over all four atoms of the Fe protein. The Mössbauer spectra of the Mo–Fe and the Fe protein under argon were unaffected by the reducible substrates N₂ and C₂H₂ and the inhibitor CO in the presence of ATP, Mg²⁺ and Na₂S₂O₄. A number of Mössbauer spectral species associated with inactivated Mo–Fe and Fe proteins are described and discussed.     

dGPredictor: Automated fragmentation method for metabolic reaction free energy prediction and de novo pathway design

Group contribution (GC) methods are conventionally used in thermodynamics analysis of metabolic pathways to estimate the standard Gibbs energy change (Δ_rG′^o) of enzymatic reactions from limited experimental measurements. However, these methods are limited by their dependence on manually curated groups and inability to capture stereochemical information, leading to low reaction coverage. Herein, we introduce an automated molecular fingerprint-based thermodynamic analysis tool called dGPredictor that enables the consideration of stereochemistry within metabolite structures and thus increases reaction coverage. dGPredictor has comparable prediction accuracy compared to existing GC methods and can capture Gibbs energy changes for isomerase and transferase reactions, which exhibit no overall group changes. We also demonstrate dGPredictor’s ability to predict the Gibbs energy change for novel reactions and seamless integration within de novo metabolic pathway design tools such as novoStoic for safeguarding against the inclusion of reaction steps with infeasible directionalities. To facilitate easy access to dGPredictor, we developed a graphical user interface to predict the standard Gibbs energy change for reactions at various pH and ionic strengths. The tool allows customized user input of known metabolites as KEGG IDs and novel metabolites as InChI strings (https://github.com/maranasgroup/dGPredictor).     

Corynebacterium glutamicum tailored for efficient isobutanol production

We recently engineered Corynebacterium glutamicum for aerobic production of 2-ketoisovalerate by inactivation of the pyruvate dehydrogenase complex, pyruvate:quinone oxidoreductase, transaminase B, and additional overexpression of the ilvBNCD genes, encoding acetohydroxyacid synthase, acetohydroxyacid isomeroreductase, and dihydroxyacid dehydratase. Based on this strain, we engineered C. glutamicum for the production of isobutanol from glucose under oxygen deprivation conditions by inactivation of l-lactate and malate dehydrogenases, implementation of ketoacid decarboxylase from Lactococcus lactis, alcohol dehydrogenase 2 (ADH2) from Saccharomyces cerevisiae, and expression of the pntAB transhydrogenase genes from Escherichia coli. The resulting strain produced isobutanol with a substrate-specific yield (Y_P/S) of 0.60 ± 0.02 mol per mol of glucose. Interestingly, a chromosomally encoded alcohol dehydrogenase rather than the plasmid-encoded ADH2 from S. cerevisiae was involved in isobutanol formation with C. glutamicum, and overexpression of the corresponding adhA gene increased the Y_P/S to 0.77 ± 0.01 mol of isobutanol per mol of glucose. Inactivation of the malic enzyme significantly reduced the Y_P/S, indicating that the metabolic cycle consisting of pyruvate and/or phosphoenolpyruvate carboxylase, malate dehydrogenase, and malic enzyme is responsible for the conversion of NADH+H⁺ to NADPH+H⁺. In fed-batch fermentations with an aerobic growth phase and an oxygen-depleted production phase, the most promising strain, C. glutamicum ΔaceE Δpqo ΔilvE ΔldhA Δmdh(pJC4ilvBNCD-pntAB)(pBB1kivd-adhA), produced about 175 mM isobutanol, with a volumetric productivity of 4.4 mM h⁻¹, and showed an overall Y_P/S of about 0.48 mol per mol of glucose in the production phase.