Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand.
Ontology highlight
ABSTRACT:
SUBMITTER: Jiang H
PROVIDER: S-EPMC9305431 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Nitrogenase is the only enzyme in nature that can fix N2 from the air. The active cofactor of the leading form of this enzyme contains seven irons and one molybdenum connected by sulfide bridges. In several recent experimental studies, it has been suggested that the cofactor is very flexible, and might lose one of its sulfides during catalysis. In this study, the possible loss of a sulfide has been investigated by model calculations. In previous studies, we have shown that there should be four activation steps before catalysis starts, and this study is based on that finding. It was found here that, after the four reductions in the activation steps, a sulfide will become very loosely bound and can be released in a quite exergonic step with a low barrier. The binding of N2 has no part in that release. In our previous studies, we suggested that the central carbide should be protonated three times after the four activation steps. With the new finding, there will instead be a loss of a sulfide, as the barrier for the loss is much lower than the ones for protonating the carbide. Still, it is suggested here that the carbide will be protonated anyway, but only with one proton, in the E3 to E4 step. A very complicated transition state for H2 formation involving a large structural change was obtained. The combined step, with a loss of H2 and binding of N2 , is calculated to be endergonic by +2.3 kcal mol-1 ; this is in excellent agreement with experiments in which an easily reversible step has been found.
| S-EPMC9303661 | biostudies-literature
Project description:BackgroundEstimating the rate constants of a biochemical reaction system with known stoichiometry from noisy time series measurements of molecular concentrations is an important step for building predictive models of cellular function. Inference techniques currently available in the literature may produce rate constant values that defy necessary constraints imposed by the fundamental laws of thermodynamics. As a result, these techniques may lead to biochemical reaction systems whose concentration dynamics could not possibly occur in nature. Therefore, development of a thermodynamically consistent approach for estimating the rate constants of a biochemical reaction system is highly desirable.ResultsWe introduce a Bayesian analysis approach for computing thermodynamically consistent estimates of the rate constants of a closed biochemical reaction system with known stoichiometry given experimental data. Our method employs an appropriately designed prior probability density function that effectively integrates fundamental biophysical and thermodynamic knowledge into the inference problem. Moreover, it takes into account experimental strategies for collecting informative observations of molecular concentrations through perturbations. The proposed method employs a maximization-expectation-maximization algorithm that provides thermodynamically feasible estimates of the rate constant values and computes appropriate measures of estimation accuracy. We demonstrate various aspects of the proposed method on synthetic data obtained by simulating a subset of a well-known model of the EGF/ERK signaling pathway, and examine its robustness under conditions that violate key assumptions. Software, coded in MATLAB®, which implements all Bayesian analysis techniques discussed in this paper, is available free of charge at http://www.cis.jhu.edu/~goutsias/CSS%20lab/software.html.ConclusionsOur approach provides an attractive statistical methodology for estimating thermodynamically feasible values for the rate constants of a biochemical reaction system from noisy time series observations of molecular concentrations obtained through perturbations. The proposed technique is theoretically sound and computationally feasible, but restricted to quantitative data obtained from closed biochemical reaction systems. This necessitates development of similar techniques for estimating the rate constants of open biochemical reaction systems, which are more realistic models of cellular function.
| S-EPMC3248051 | biostudies-literature
Project description:This paper describes a comparative study of the decomposition of two nitrosyl iron complexes (NICs) with penicillamine thiolic ligands [Fe2(SC5H11NO2)2(NO)4]SO4 ·5H2O (I) and glutathione- (GSH-) ligands [Fe2(SC10H17N3O6)2(NO)4]SO4 ·2H2O (II), which spontaneously evolve to NO in aqueous medium. NO formation was measured by a sensor electrode and by spectrophotometric methods by measuring the formation of a hemoglobin- (Hb-) NO complex. The NO evolution reaction rate from (I)??k 1 = (4.6 ± 0.1)·10(-3)?s(-1) and the elimination rate constant of the penicillamine ligand k 2 = (1.8 ± 0.2)·10(-3)?s(-1) at 25°C in 0.05?M phosphate buffer, ?pH 7.0, was calculated using kinetic modeling based on the experimental data. Both reactions are reversible. Spectrophotometry and mass-spectrometry methods have firmly shown that the penicillamine ligand is exchanged for GS(-) during decomposition of 1.5·10(-4)?M (I) in the presence of 10(-3)?M GSH, with 76% yield in 24?h. As has been established, such behaviour is caused by the resistance of (II) to decomposition due to the higher affinity of iron to GSH in the complex. The discovered reaction may impede S-glutathionylation of the essential enzyme systems in the presence of (I) and is important for metabolism of NIC, connected with its antitumor activity.
| S-EPMC3984828 | biostudies-literature
Project description:Reaction centers (RCs) are integral membrane proteins that undergo a series of electron transfer reactions during the process of photosynthesis. In the Q(A) site of RCs from Rhodobacter sphaeroides, ubiquinone-10 is reduced, by a single electron transfer, to its semiquinone. The neutral quinone and anionic semiquinone have similar affinities, which is required for correct in situ reaction thermodynamics. A previous study showed that despite similar affinities, anionic quinones associate and dissociate from the Q(A) site at rates ≈10(4) times slower than neutral quinones indicating that anionic quinones encounter larger binding barriers (Madeo, J.; Gunner, M. R. Modeling binding kinetics at the Q(A) site in bacterial reaction centers. Biochemistry 2005, 44, 10994-11004). The present study investigates these barriers computationally, using steered molecular dynamics (SMD) to model the unbinding of neutral ground state ubiquinone (UQ) and its reduced anionic semiquinone (SQ(-)) from the Q(A) site. In agreement with experiment, the SMD unbinding barrier for SQ(-) is larger than for UQ. Multi Conformational Continuum Electrostatics (MCCE), used here to calculate the binding energy, shows that SQ(-) and UQ have comparable affinities. In the Q(A) site, there are stronger binding interactions for SQ(-) compared to UQ, especially electrostatic attraction to a bound non-heme Fe(2+). These interactions compensate for the higher SQ(-) desolvation penalty, allowing both redox states to have similar affinities. These additional interactions also increase the dissociation barrier for SQ(-) relative to UQ. Thus, the slower SQ(-) dissociation rate is a direct physical consequence of the additional binding interactions required to achieve a Q(A) site affinity similar to that of UQ. By a similar mechanism, the slower association rate is caused by stronger interactions between SQ(-) and the polar solvent. Thus, stronger interactions for both the unbound and bound states of charged and highly polar ligands can slow their binding kinetics without a conformational gate. Implications of this for other systems are discussed.
| S-EPMC3202297 | biostudies-literature
Project description:Nitrogenases catalyze the reduction of dinitrogen (N2) gas to ammonium at a complex heterometallic cofactor. This most commonly occurs at the FeMo cofactor (FeMoco), a [Mo-7Fe-9S-C] cluster whose exact reactivity and substrate-binding mode remain unknown. Alternative nitrogenases replace molybdenum with either vanadium or iron and differ in reactivity, most prominently in the ability of vanadium nitrogenase to reduce CO to hydrocarbons. Here we report the 1.35-Å structure of vanadium nitrogenase from Azotobacter vinelandii. The 240-kDa protein contains an additional α-helical subunit that is not present in molybdenum nitrogenase. The FeV cofactor (FeVco) is a [V-7Fe-8S-C] cluster with a homocitrate ligand to vanadium. Unexpectedly, it lacks one sulfide ion compared to FeMoco, which is replaced by a bridging ligand, likely a μ-1,3-carbonate. The anion fits into a pocket within the protein that is obstructed in molybdenum nitrogenase, and its different chemical character helps to rationalize the altered chemical properties of this unique N2- and CO-fixing enzyme.
| S-EPMC5563456 | biostudies-literature
Project description:A reaction based fluorescence turn-on strategy for hydrogen sulfide (H(2)S) was developed. This strategy was based on a H(2)S-specific Michael addition-cyclization sequence. Other biological thiols such as cysteine and glutathione did not pursue the reaction and therefore did not turn on the fluorescence/consume the substrates. The probes showed good selectivity and sensitivity for hydrogen sulfide.
| S-EPMC3336739 | biostudies-literature
Project description:Lithium alloys are synthesized by reactions between lithium metal and group 14 elements, such as carbon, silicon, germanium, and tin. The nitrogenation and denitrogenation properties are investigated by thermal and structural analyses. All alloys dissociate the nitrogen triple bond of gaseous molecules to form atomic state as nitrides below 500 °C, which is lower than those required for conventional thermochemical and catalytic processes on nitride syntheses. For all alloys except for germanium, it is indicated that nanosized lithium nitride is formed as the product. The denitrogenation (nitrogen desorption) reaction by lithium nitride and metals, which is an ideal opposite reaction of nitrogenation, occurs by heating up to 600 °C to form lithium alloys. Among them, the lithium-tin alloy is a potential material to control the dissociation and recombination of nitrogen below 500 °C by the reversible reaction with the largest amount of utilizable lithium in the alloy phase. The nitrogenation and denitrogenation reactions of the lithium alloys at lower temperature are realized by the high reactivity with nitrogen and mobility of lithium. The above reactions based on lithium alloys are adapted to the ammonia synthesis. As a result, ammonia can be synthesized below 500 °C under 0.5 MPa of pressure. Therefore, the reaction using lithium alloys is recognized as a pseudocatalyst for the ammonia synthesis.
| S-EPMC6640966 | biostudies-literature
Project description:Nitrogenase is the only enzyme that can convert N2 to NH3 . Crystallographic structures have indicated that one of the sulfide ligands of the active-site FeMo cluster, S2B, can be replaced by an inhibitor, like CO and OH- , and it has been suggested that it may be displaced also during the normal reaction. We have investigated possible proton transfer pathways within the FeMo cluster during the conversion of N2 H2 to two molecules of NH3 , assuming that the protons enter the cluster at the S3B, S4B or S5A sulfide ions and are then transferred to the substrate. We use combined quantum mechanical and molecular mechanical (QM/MM) calculations with the TPSS and B3LYP functionals. The calculations indicate that the barriers for these reactions are reasonable if the S2B ligand remains bound to the cluster, but they become prohibitively high if S2B has dissociated. This suggests that it is unlikely that S2B reversibly dissociates during the normal reaction cycle.
| S-EPMC9804283 | biostudies-literature
Project description:The pepsin folding mechanism involves a prosegment (PS) domain that catalyzes folding, which is then removed, resulting in a kinetically trapped native state. Although native pepsin (Np) is kinetically stable, it is irreversibly denatured due to a large folding barrier, and in the absence of the PS it folds to a more thermodynamically stable denatured state, termed refolded pepsin (Rp). This system serves as a model to understand the nature of kinetic barriers and folding transitions between compact states. Quasielastic neutron scattering (QENS) was used to characterize and compare the flexibility of Np, as a kinetically trapped state, with that of Rp, as a thermodynamically stable fold. Additionally, the dynamics of Np were compared with those of a partially unfolded form and a thermally stabilized, inhibitor-bound form. QENS revealed length-scale-dependent differences between Np and Rp on a picosecond timescale and indicated greater flexibility in Np, leading to the conclusion that kinetic stabilization likely does not correspond to reduced internal dynamics. Furthermore, large differences were observed upon inhibition, indicating that QENS of proteins in solution may prove useful for examining the role of conformational entropy changes in ligand binding.
| S-EPMC3183798 | biostudies-other
Project description:Coarse-graining of fully atomistic molecular dynamics simulations is a long-standing goal in order to allow the description of processes occurring on biologically relevant timescales. For example, the prediction of pathways, rates and rate-limiting steps in protein-ligand unbinding is crucial for modern drug discovery. To achieve the enhanced sampling, we perform dissipation-corrected targeted molecular dynamics simulations, which yield free energy and friction profiles of molecular processes under consideration. Subsequently, we use these fields to perform temperature-boosted Langevin simulations which account for the desired kinetics occurring on multisecond timescales and beyond. Adopting the dissociation of solvated sodium chloride, trypsin-benzamidine and Hsp90-inhibitor protein-ligand complexes as test problems, we reproduce rates from molecular dynamics simulation and experiments within a factor of 2-20, and dissociation constants within a factor of 1-4. Analysis of friction profiles reveals that binding and unbinding dynamics are mediated by changes of the surrounding hydration shells in all investigated systems.
| S-EPMC7286908 | biostudies-literature