Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I

Electronic Polarization Leads to a Drier Dewetted State for Hydrophobic Gating in the Big Potassium Channel

In the hydrophobic gating mechanism proposed for some ion channels, ion permeation is not blocked by the physical dimension of the channel pore but by its dewetted state which constitutes the energetic bottleneck. A major source of uncertainty in the mechanism is that the dewetted state was not observed in experiments and only probed in simulations using nonpolarizable force fields, which do not accurately represent the properties of confined water. Here we analyze hydration of the central cavity in the pore-gate domain of the Big Potassium channel using molecular dynamics and grand canonical Monte Carlo simulations with enhanced sampling techniques. Including polarization leads to a much drier dewetted state and a higher barrier for the transition to the wet state, suggesting more effective hydrophobic gating. The simulations also identify two backbone carbonyls at the bottom of the selectivity filter as good candidates for characterizing the dewetted state using infrared spectroscopies.

QM/MM Free Energy Calculations of Long-Range Biological Protonation Dynamics by Adaptive and Focused Sampling

Water-mediated proton transfer reactions are central for catalytic processes in a wide range of biochemical systems, ranging from biological energy conversion to chemical transformations in the metabolism. Yet, the accurate computational treatment of such complex biochemical reactions is highly challenging and requires the application of multiscale methods, in particular hybrid quantum/classical (QM/MM) approaches combined with free energy simulations. Here, we combine the unique exploration power of new advanced sampling methods with density functional theory (DFT)-based QM/MM free energy methods for multiscale simulations of long-range protonation dynamics in biological systems. In this regard, we show that combining multiple walkers/well-tempered metadynamics with an extended system adaptive biasing force method (MWE) provides a powerful approach for exploration of water-mediated proton transfer reactions in complex biochemical systems. We compare and combine the MWE method also with QM/MM umbrella sampling and explore the sampling of the free energy landscape with both geometric (linear combination of proton transfer distances) and physical (center of excess charge) reaction coordinates and show how these affect the convergence of the potential of mean force (PMF) and the activation free energy. We find that the QM/MM-MWE method can efficiently explore both direct and water-mediated proton transfer pathways together with forward and reverse hole transfer mechanisms in the highly complex proton channel of respiratory Complex I, while the QM/MM-US approach shows a systematic convergence of selected long-range proton transfer pathways. In this regard, we show that the PMF along multiple proton transfer pathways is recovered by combining the strengths of both approaches in a QM/MM-MWE/focused US (FUS) scheme and reveals new mechanistic insight into the proton transfer principles of Complex I. Our findings provide a promising basis for the quantitative multiscale simulations of long-range proton transfer reactions in biological systems.

Mechanistic Principles of Hydrogen Evolution in the Membrane-Bound Hydrogenase

The membrane-bound hydrogenase (Mbh) from Pyrococcus furiosus is an archaeal member of the Complex I superfamily. It catalyzes the reduction of protons to H2 gas powered by a [NiFe] active site and transduces the free energy into proton pumping and Na+/H+ exchange across the membrane. Despite recent structural advances, the mechanistic principles of H2 catalysis and ion transport in Mbh remain elusive. Here, we probe how the redox chemistry drives the reduction of the proton to H2 and how the catalysis couples to conformational dynamics in the membrane domain of Mbh. By combining large-scale quantum chemical density functional theory (DFT) and correlated ab initio wave function methods with atomistic molecular dynamics simulations, we show that the proton transfer reactions required for the catalysis are gated by electric field effects that direct the protons by water-mediated reactions from Glu21L toward the [NiFe] site, or alternatively along the nearby His75L pathway that also becomes energetically feasible in certain reaction steps. These local proton-coupled electron transfer (PCET) reactions induce conformational changes around the active site that provide a key coupling element via conserved loop structures to the ion transport activity. We find that H2 forms in a heterolytic proton reduction step, with spin crossovers tuning the energetics along key reaction steps. On a general level, our work showcases the role of electric fields in enzyme catalysis and how these effects are employed by the [NiFe] active site of Mbh to drive PCET reactions and ion transport.

Photochromic radical states in 3D covalent organic frameworks with zyg topology for enhanced photocatalysis

Covalent-organic frameworks (COFs) with photoinduced donor-acceptor (D-A) radical pairs show enhanced photocatalytic activity in principle. However, achieving long-lived charge separation in COFs proves challenging due to the rapid charge recombination. Here, we develop a novel strategy by combining [6 + 4] nodes to construct zyg-type 3D COFs, first reported in COF chemistry. This structure type exhibits a fused Olympic-rings-like shape, which provides a platform for stabilizing the photoinduced D-A radical pairs. The zyg-type COFs containing catalytically active moieties such as triphenylamine and phenothiazine (PTZ) show superior photocatalytic production rates of hydrogen peroxide (H2O2). Significantly, the photochromic radical states of these COFs show up to 400% enhancement in photocatalytic activity compared to the parent states, achieving a remarkable H2O2 synthesis rate of 3324 μmol g-1 h-1, which makes the PTZ-COF one of the best crystalline porous photocatalysts in H2O2 production. This work will shed light on the synthesis of efficient 3D COF photocatalysts built on topologies that can facilitate photogenerating D-A radical pairs for enhanced photocatalysis.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Mechanism of mitochondrial oxidative phosphorylation disorder in male infertility

ABSTRACT

Male infertility has become a global concern, accounting for 20-70% of infertility. Dysfunctional spermatogenesis is the most common cause of male infertility; thus, treating abnormal spermatogenesis may improve male infertility and has attracted the attention of the medical community. Mitochondria are essential organelles that maintain cell homeostasis and normal physiological functions in various ways, such as mitochondrial oxidative phosphorylation (OXPHOS). Mitochondrial OXPHOS transmits electrons through the respiratory chain, synthesizes adenosine triphosphate (ATP), and produces reactive oxygen species (ROS). These mechanisms are vital for spermatogenesis, especially to maintain the normal function of testicular Sertoli cells and germ cells. The disruption of mitochondrial OXPHOS caused by external factors can result in inadequate cellular energy supply, oxidative stress, apoptosis, or ferroptosis, all inhibiting spermatogenesis and damaging the male reproductive system, leading to male infertility. This article summarizes the latest pathological mechanism of mitochondrial OXPHOS disorder in testicular Sertoli cells and germ cells, which disrupts spermatogenesis and results in male infertility. In addition, we also briefly outline the current treatment of spermatogenic malfunction caused by mitochondrial OXPHOS disorders. However, relevant treatments have not been fully elucidated. Therefore, targeting mitochondrial OXPHOS disorders in Sertoli cells and germ cells is a research direction worthy of attention. We believe this review will provide new and more accurate ideas for treating male infertility.

Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields

Prediction of the hydration levels of protein cavities and active sites is important to both mechanistic analysis and ligand design. Due to the unique microscopic environment of these buried water molecules, a polarizable model is expected to be crucial for an accurate treatment of protein internal hydration in simulations. Here we adapt a nonequilibrium candidate Monte Carlo approach for conducting grand canonical Monte Carlo simulations with the Drude polarizable force field. The GPU implementation enables the efficient sampling of internal cavity hydration levels in biomolecular systems. We also develop an enhanced sampling approach referred to as B-walking, which satisfies detailed balance and readily combines with grand canonical integration to efficiently calculate quantitative binding free energies of water to protein cavities. Applications of these developments are illustrated in a solvent box and the polar ligand binding site in trypsin. Our simulation results show that including electronic polarization leads to a modest but clear improvement in the description of water position and occupancy compared to the crystal structure. The B-walking approach enhances the range of water sampling in different chemical potential windows and thus improves the accuracy of water binding free energy calculations.

Enhanced Sampling of Buried Charges in Free Energy Calculations Using Replica Exchange with Charge Tempering

Buried ionizable groups in proteins often play important structural and functional roles. However, it is generally challenging to study the detailed molecular mechanisms solely based on experimental measurements. Free energy calculations using atomistic simulations, on the other hand, complement experimental studies and can provide high temporal and spatial resolution information that can lead to mechanistic insights. Nevertheless, it is also well recognized that sufficient sampling of such atomistic simulations can be challenging, considering that structural changes related to the buried charges may be very slow. In the present study, we describe a simple but effective enhanced sampling technique called replica exchange with charge tempering (REChgT) with a novel free energy method, multisite λ dynamics (MSλD), to study two systems containing buried charges, pKa prediction of a small molecule, orotate, in complex with the dihydroorotate dehydrogenase, and relative stability of a Glu-Lys pair buried in the hydrophobic core of two variants of Staphylococcal nuclease. Compared to the original MSλD simulations, the usage of REChgT dramatically increases sampling in both conformational and alchemical spaces, which directly translates into a significant reduction of wall time to converge the free energy calculations. This study highlights the importance of sufficient sampling toward developing improved free energy methods.

Photo-oxidation of methanol in complexes with pyrido[2,3-b]pyrazine: a nonadiabatic molecular dynamics study

Excited-state Proton-Coupled Electron Transfer (PCET) constitutes a key step in the photo-oxidation of small, electron-rich systems possessing acidic hydrogen atoms, such as water or alcohols, which can play a vital role in green hydrogen production. In this contribution, we employ ab initio quantum-chemical methods and on-the-fly nonadiabatic molecular dynamics simulations to study the mechanism and the photodynamics of PCET in 1 : 1 pyrido[2,3-b]pyrazine complexes with methanol. We find the process to be ultrafast and efficient when the intramolecular hydrogen bond is formed with one of the β-positioned nitrogen atoms. The complex exhibiting a hydrogen bond with an isolated nitrogen site, on the contrary, shows much lower reactivity. We explain this effect with the stabilization of the reactive ππ* charge-transfer electronic state in the former case.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.

Mutation at the entrance of the quinone cavity severely disrupts quinone binding in respiratory complex I

In all resolved structures of complex I, there exists a tunnel-like Q-chamber for ubiquinone binding and reduction. The entrance to the Q-chamber in ND1 subunit forms a narrow bottleneck, which is rather tight and requires thermal conformational changes for ubiquinone to get in and out of the binding chamber. The substitution of alanine with threonine at the bottleneck (AlaThr MUT), associated with 3460/ND1 mtDNA mutation in human complex I, is implicated in Leber's Hereditary Optic Neuropathy (LHON). Here, we show the AlaThr MUT further narrows the Q-chamber entrance cross-section area by almost 30%, increasing the activation free energy barrier of quinone passage by approximately 5 kJ mol-1. This severely disrupts quinone binding and reduction as quinone passage through the bottleneck is slowed down almost tenfold. Our estimate of the increase in free energy barrier is entirely due to the bottleneck narrowing, leading to a reduction of the transition state entropy between WT and MUT, and thus more difficult quinone passage. Additionally, we investigate details of possible water exchange between the Q-chamber and membrane. We find water exchange is dynamic in WT but may be severely slowed in MUT. We propose that LHON symptoms caused by 3460/ND1 mtDNA mutation are due to slowed quinone binding. This leads to an increased production of reactive oxidative species due to upstream electron backup at the FMN site of complex I, thus resulting in a mt bioenergetic defect.

Mitochondrial Cristae Morphology Reflecting Metabolism, Superoxide Formation, Redox Homeostasis, and Pathology

Significance: Mitochondrial (mt) reticulum network in the cell possesses amazing ultramorphology of parallel lamellar cristae, formed by the invaginated inner mitochondrial membrane. Its non-invaginated part, the inner boundary membrane (IBM) forms a cylindrical sandwich with the outer mitochondrial membrane (OMM). Crista membranes (CMs) meet IBM at crista junctions (CJs) of mt cristae organizing system (MICOS) complexes connected to OMM sorting and assembly machinery (SAM). Cristae dimensions, shape, and CJs have characteristic patterns for different metabolic regimes, physiological and pathological situations. Recent Advances: Cristae-shaping proteins were characterized, namely rows of ATP-synthase dimers forming the crista lamella edges, MICOS subunits, optic atrophy 1 (OPA1) isoforms and mitochondrial genome maintenance 1 (MGM1) filaments, prohibitins, and others. Detailed cristae ultramorphology changes were imaged by focused-ion beam/scanning electron microscopy. Dynamics of crista lamellae and mobile CJs were demonstrated by nanoscopy in living cells. With tBID-induced apoptosis a single entirely fused cristae reticulum was observed in a mitochondrial spheroid. Critical Issues: The mobility and composition of MICOS, OPA1, and ATP-synthase dimeric rows regulated by post-translational modifications might be exclusively responsible for cristae morphology changes, but ion fluxes across CM and resulting osmotic forces might be also involved. Inevitably, cristae ultramorphology should reflect also mitochondrial redox homeostasis, but details are unknown. Disordered cristae typically reflect higher superoxide formation. Future Directions: To link redox homeostasis to cristae ultramorphology and define markers, recent progress will help in uncovering mechanisms involved in proton-coupled electron transfer via the respiratory chain and in regulation of cristae architecture, leading to structural determination of superoxide formation sites and cristae ultramorphology changes in diseases. Antioxid. Redox Signal. 39, 635-683.

Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I

Complex I is a redox-driven proton pump that drives electron transport chains and powers oxidative phosphorylation across all domains of life. Yet, despite recently resolved structures from multiple organisms, it still remains unclear how the redox reactions in Complex I trigger proton pumping up to 200 Å away from the active site. Here, we show that the proton-coupled electron transfer reactions during quinone reduction drive long-range conformational changes of conserved loops and trans-membrane (TM) helices in the membrane domain of Complex I from Yarrowia lipolytica. We find that the conformational switching triggers a π → α transition in a TM helix (TM3ND6) and establishes a proton pathway between the quinone chamber and the antiporter-like subunits, responsible for proton pumping. Our large-scale (>20 μs) atomistic molecular dynamics (MD) simulations in combination with quantum/classical (QM/MM) free energy calculations show that the helix transition controls the barrier for proton transfer reactions by wetting transitions and electrostatic effects. The conformational switching is enabled by re-arrangements of ion pairs that propagate from the quinone binding site to the membrane domain via an extended network of conserved residues. We find that these redox-driven changes create a conserved coupling network within the Complex I superfamily, with point mutations leading to drastic activity changes and mitochondrial disorders. On a general level, our findings illustrate how catalysis controls large-scale protein conformational changes and enables ion transport across biological membranes.

A self-sustaining serpentinization mega-engine feeds the fougerite nanoengines implicated in the emergence of guided metabolism

The demonstration by Ivan Barnes et al. that the serpentinization of fresh Alpine-type ultramafic rocks results in the exhalation of hot alkaline fluids is foundational to the submarine alkaline vent theory (AVT) for life's emergence to its 'improbable' thermodynamic state. In AVT, such alkaline fluids ≤ 150°C, bearing H2 > CH4 > HS--generated and driven convectively by a serpentinizing exothermic mega-engine operating in the ultramafic crust-exhale into the iron-rich, CO2> > > NO3--bearing Hadean ocean to result in hydrothermal precipitate mounds comprising macromolecular ferroferric-carbonate oxyhydroxide and minor sulfide. As the nanocrystalline minerals fougerite/green rust and mackinawite (FeS), they compose the spontaneously precipitated inorganic membranes that keep the highly contrasting solutions apart, thereby maintaining redox and pH disequilibria. They do so in the form of fine chimneys and chemical gardens. The same disequilibria drive the reduction of CO2 to HCOO- or CO, and the oxidation of CH4 to a methyl group-the two products reacting to form acetate in a sequence antedating the 'energy-producing' acetyl coenzyme-A pathway. Fougerite is a 2D-layered mineral in which the hydrous interlayers themselves harbor 2D solutions, in effect constricted to ~ 1D by preferentially directed electron hopping/tunneling, and proton Gröthuss 'bucket-brigading' when subject to charge. As a redox-driven nanoengine or peristaltic pump, fougerite forces the ordered reduction of nitrate to ammonium, the amination of pyruvate and oxalate to alanine and glycine, and their condensation to short peptides. In turn, these peptides have the flexibility to sequester the founding inorganic iron oxyhydroxide, sulfide, and pyrophosphate clusters, to produce metal- and phosphate-dosed organic films and cells. As the feed to the hydrothermal mound fails, the only equivalent sustenance on offer to the first autotrophs is the still mildly serpentinizing upper crust beneath. While the conditions here are very much less bountiful, they do offer the similar feed and disequilibria the survivors are accustomed to. Sometime during this transition, a replicating non-ribosomal guidance system is discovered to provide the rules to take on the incrementally changing surroundings. The details of how these replicating apparatuses emerged are the hard problem, but by doing so the progenote archaea and bacteria could begin to colonize what would become the deep biosphere. Indeed, that the anaerobic nitrate-respiring methanotrophic archaea and the deep-branching Acetothermia presently comprise a portion of that microbiome occupying serpentinizing rocks offers circumstantial support for this notion. However, the inescapable, if jarring conclusion is drawn that, absent fougerite/green rust, there would be no structured channelway to life.

H2O2 selectively damages the binuclear iron-sulfur cluster N1b of respiratory complex I

NADH:ubiquinone oxidoreductase, respiratory complex I, plays a major role in cellular energy metabolism by coupling electron transfer with proton translocation. Electron transfer is catalyzed by a flavin mononucleotide and a series of iron-sulfur (Fe/S) clusters. As a by-product of the reaction, the reduced flavin generates reactive oxygen species (ROS). It was suggested that the ROS generated by the respiratory chain in general could damage the Fe/S clusters of the complex. Here, we show that the binuclear Fe/S cluster N1b is specifically damaged by H₂O₂, however, only at high concentrations. But under the same conditions, the activity of the complex is hardly affected, since N1b can be easily bypassed during electron transfer.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

The decrease in diurnal oxygen production in Elodea under the influence of high geomagnetic variability: the role of light, temperature and atmospheric pressure

Epidemiological studies have indicated adverse effects of geomagnetic disturbance on human health, including increased mortality. There is evidence from plant and animal studies that help to elucidate this interaction. This study tests the hypothesis that geomagnetic disturbance affects living systems, by modifying the metabolic process of photosynthesis, in the natural environment.Continuous 24-h measurements of dissolved oxygen in flasks containing Holtfreiter's solution and strands of healthy Elodea were recorded from May 1996, until September 1998, in an electromagnetically quiet, purpose built, garden shed environment, without mains electricity. Sensormeter recordings of oxygen, light, temperature and air pressure were uploaded weekly to a PC. The hourly total geomagnetic field measurements were obtained from the nearest observatory.Significant decrease in oxygen (diurnal volume of oxygen divided by plant mass and diurnal light), (O/WL), was found on days of high geomagnetic field variability throughout 11 recorded months of the year 1997. This result was independent of temperature and atmospheric pressure. No significant decrease in O/WL during high geomagnetic variability was found for the 7 months recorded in 1996. The 1996 and 1997 data both showed a significant decrease in the diurnal time lag between peak light and peak oxygen for diurnal high geomagnetic variability compared with low geomagnetic variability. Cross correlation analysis for 1997 and 1998 data showed a decrease in positive correlation of oxygen with light in high geomagnetic variability, compared with low geomagnetic variability, and increased positive correlation with the geomagnetic field instead. These experiments support a hypothesis of high geomagnetic field variability as a weak zeitgeber, and a metabolic depressant for photosynthetic oxygen production in plants.

On the role of ubiquinone in the proton translocation mechanism of respiratory complex I

Complex I converts oxidoreduction energy into a proton electrochemical gradient across the inner mitochondrial or bacterial cell membrane. This gradient is the primary source of energy for aerobic synthesis of ATP. Oxidation of reduced nicotinamide adenine dinucleotide (NADH) by ubiquinone (Q) yields NAD⁺ and ubiquinol (QH₂ ), which is tightly coupled to translocation of four protons from the negatively to the positively charged side of the membrane. Electrons from NADH oxidation reach the iron-sulfur centre N2 positioned near the bottom of a tunnel that extends circa 30 Å from the membrane domain into the hydrophilic domain of the complex. The tunnel is occupied by ubiquinone, which can take a distal position near the N2 centre or proximal positions closer to the membrane. Here, we review important structural, kinetic and thermodynamic properties of ubiquinone that define its role in complex I function. We suggest that this function exceeds that of a mere substrate or electron acceptor and propose that ubiquinone may be the redox element of complex I coupling electron transfer to proton translocation.

Single protonation of the reduced quinone in respiratory complex I drives four-proton pumping

Complex I is a key proton-pumping enzyme in bacterial and mitochondrial respiratory electron transport chains. Using quantum chemistry and electrostatic calculations, we have examined the pKa of the reduced quinone QH-/QH2 in the catalytic cavity of complex I. We find that pKa (QH-/QH2) is very high, above 20. This means that the energy of a single protonation reaction of the doubly reduced quinone (i.e. the reduced semiquinone QH-) is sufficient to drive four protons across the membrane with a potential of 180 mV. Based on these calculations, we propose a possible scheme of redox-linked proton pumping by complex I. The model explains how the energy of the protonation reaction can be divided equally among four pumping units of the pump, and how a single proton can drive translocation of four additional protons in multiple pumping blocks.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Knockdown of NDUFC1 inhibits cell proliferation, migration, and invasion of hepatocellular carcinoma

Background

NADH: ubiquinone oxidoreductase subunit C1(NDUFC1) encodes a subunit of the Complex I, which may support the structural stability of Complex I and assist in its biogenesis. The expression and functional roles of NDUFC1 in hepatocellular carcinoma (HCC) remain unknown.

Result

We knocked down the expression of NDUFC1 in HCC cell lines to explore the effects of NDUFC1 downregulation on HCC in vitro. MTT assay determined that downregulation of NDUFC1 significantly inhibited cell proliferation. Flow cytometry with (propidium iodide) PI staining indicated silencing of NDUFC1 arrested cell cycle of BEL-7404 cells at G2 phase and SK-HEP-1 cells at S/G2 phase. Annexin V-PI double staining and flow cytometric analysis showed that the downregulation of NDUFC1 significantly increased the population of apoptotic cells. Wound-healing assay and transwell assay indicated that the downregulation of NDUFC1 suppressed the migration and invasion of HCC cells. According to the detection of complex1 activity, we found that the activity of NDUFC1 silenced group decreased, whereas the content of ROS increased. Furthermore, combined with bioinformatics analysis of senescence-related genes, we found that the silence of NDUFC1 in HCC could induce senescence and inhibit autophagy. In addition, NDUFC1 could correlate positively with cancer-related pathways, among which the p53 pathways and the PI3K/Akt/mTOR pathways. Finally, NDUFC1 is high expression in HCC specimens. High NDUFC1 expression was associated with poor prognosis and was an independent risk factor for reduced overall survival (OS).

Conclusions

Our study indicated, for the first time, that NDUFC1 is an independent risk factor for the poor prognosis of HCC patients. NDUFC1 may promote tumor progression by inhibiting mitochondrial Complex I and up-regulating ROS through multiple cancer-related and senescence-related pathways of HCC, including p53 pathways and PI3K/Akt/mTOR pathways. We suppose that NDUFC1 might be a potential target for the mitochondrial metabolism therapy of HCC.

Molecular Principles of Redox-Coupled Protonation Dynamics in Photosystem II

Photosystem II (PSII) catalyzes light-driven water oxidization, releasing O₂ into the atmosphere and transferring the electrons for the synthesis of biomass. However, despite decades of structural and functional studies, the water oxidation mechanism of PSII has remained puzzling and a major challenge for modern chemical research. Here, we show that PSII catalyzes redox-triggered proton transfer between its oxygen-evolving Mn₄O₅Ca cluster and a nearby cluster of conserved buried ion-pairs, which are connected to the bulk solvent via a proton pathway. By using multi-scale quantum and classical simulations, we find that oxidation of a redox-active Tyr_z (Tyr161) lowers the reaction barrier for the water-mediated proton transfer from a Ca²⁺-bound water molecule (W3) to Asp61 via conformational changes in a nearby ion-pair (Asp61/Lys317). Deprotonation of this W3 substrate water triggers its migration toward Mn1 to a position identified in recent X-ray free-electron laser (XFEL) experiments [Ibrahim et al. Proc. Natl. Acad. Sci. USA 2020, 117, 12,624-12,635]. Further oxidation of the Mn₄O₅Ca cluster lowers the proton transfer barrier through the water ligand sphere of the Mn₄O₅Ca cluster to Asp61 via a similar ion-pair dissociation process, while the resulting Mn-bound oxo/oxyl species leads to O₂ formation by a radical coupling mechanism. The proposed redox-coupled protonation mechanism shows a striking resemblance to functional motifs in other enzymes involved in biological energy conversion, with an interplay between hydration changes, ion-pair dynamics, and electric fields that modulate the catalytic barriers.

Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants

Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.

The clinically relevant triple mutation in the mtND1 gene inactivates Escherichia coli complex I

NADH:ubiquinone oxidoreductase (respiratory complex I) plays a major role in cellular energy metabolism. Complex I deficiencies are the most common cause of mitochondrial dysfunction. Patients suffering from a variety of neurodegenerative diseases carry numerous mutations in the mitochondrially encoded subunits of the complex. The biochemical consequences of these mutations are largely unknown because these genes are difficult to access experimentally. Here, we use Escherichia coli as a model system to characterize the effect of a 7 bp inversion in mtND1 (m.3902-3908inv7) that results in a triple mutation. The triple mutant grew poorly but contained a normal amount of the stably assembled variant. The variant showed no enzymatic activity, which might contribute to the deleterious effect of the mutation in humans.