RESUMO
Botulinum neurotoxins (BoNTs) are the most potent toxins known and are used to treat an increasing number of medical disorders. All BoNTs are naturally co-expressed with a protective partner protein (NTNH) with which they form a 300 kDa complex, to resist acidic and proteolytic attack from the digestive tract. We have previously identified a new botulinum neurotoxin serotype, BoNT/X, that has unique and therapeutically attractive properties. We present the cryo-EM structure of the BoNT/X-NTNH/X complex and the crystal structure of the isolated NTNH protein. Unexpectedly, the BoNT/X complex is stable and protease-resistant at both neutral and acidic pH and disassembles only in alkaline conditions. Using the stabilizing effect of NTNH, we isolated BoNT/X and showed that it has very low potency both in vitro and in vivo. Given the high catalytic activity and translocation efficacy of BoNT/X, low activity of the full toxin is likely due to the receptor-binding domain, which presents very weak ganglioside binding and exposed hydrophobic surfaces.
RESUMO
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.
Assuntos
Hidrogênio , Hidrogenase , Simulação de Dinâmica Molecular , Pyrococcus furiosus , Hidrogenase/química , Hidrogenase/metabolismo , Hidrogênio/química , Hidrogênio/metabolismo , Pyrococcus furiosus/enzimologia , Prótons , Teoria da Densidade Funcional , Domínio Catalítico , OxirreduçãoRESUMO
Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III2IV2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria.
Assuntos
Microscopia Crioeletrônica , Mycobacterium smegmatis , Mycobacterium smegmatis/metabolismo , Mycobacterium smegmatis/enzimologia , Transporte de Elétrons , Oxirredução , Proteínas de Bactérias/metabolismo , Proteínas de Bactérias/química , Prótons , Complexo III da Cadeia de Transporte de Elétrons/metabolismo , Complexo III da Cadeia de Transporte de Elétrons/química , Oxigênio/metabolismo , Complexo IV da Cadeia de Transporte de Elétrons/metabolismo , Complexo IV da Cadeia de Transporte de Elétrons/química , Domínio Catalítico , Modelos MolecularesRESUMO
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.
Assuntos
Simulação de Dinâmica Molecular , Prótons , Termodinâmica , Água , Água/química , Teoria Quântica , Teoria da Densidade FuncionalRESUMO
Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.
Assuntos
Proteínas de Bactérias , Entomoplasmataceae , Ribonucleotídeo Redutases , Transporte de Elétrons , Prótons , Ribonucleotídeo Redutases/química , Cristalografia por Raios X/métodos , Entomoplasmataceae/enzimologia , Domínio Catalítico , Proteínas de Bactérias/químicaRESUMO
Directional ion transport across biological membranes plays a central role in many cellular processes. Elucidating the molecular determinants for vectorial ion transport is key to understanding the functional mechanism of membrane-bound ion pumps. The extensive investigation of the light-driven proton pump bacteriorhodopsin from Halobacterium salinarum(HsBR) enabled a detailed description of outward proton transport. Although the structure of inward-directed proton pumping rhodopsins is very similar to HsBR, little is known about their protonation pathway, and hence, the molecular reasons for the vectoriality of proton translocation remain unclear. Here, we employ a combined experimental and theoretical approach to tracking protonation steps in the light-driven inward proton pump xenorhodopsin from Nanosalina sp. (NsXeR). Time-resolved infrared spectroscopy reveals the transient deprotonation of D220 concomitantly with deprotonation of the retinal Schiff base. Our molecular dynamics simulations support a proton release pathway from the retinal Schiff base via a hydrogen-bonded water wire leading to D220 that could provide a putative gating point for the proton release and with allosteric interactions to the retinal Schiff base. Our findings support the key role of D220 in mediating proton release to the cytoplasmic side and provide evidence that this residue is not the primary proton acceptor of the proton transiently released by the retinal Schiff base.
RESUMO
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.
Assuntos
Complexo I de Transporte de Elétrons , Prótons , Complexo I de Transporte de Elétrons/metabolismo , Oxirredução , Transporte de Elétrons , Quinonas , Bombas de Próton/metabolismo , CatáliseRESUMO
Recent breakthroughs in neural network-based structure prediction methods, such as AlphaFold2 and RoseTTAFold, have dramatically improved the quality of computational protein structure prediction. These models also provide statistical confidence scores that can estimate uncertainties in the predicted structures, but it remains unclear to what extent these scores are related to the intrinsic conformational dynamics of proteins. Here, we compare AlphaFold2 prediction scores with explicit large-scale molecular dynamics simulations of 28 one- and two-domain proteins with varying degrees of flexibility. We demonstrate a strong correlation between the statistical prediction scores and the explicit motion derived from extensive atomistic molecular dynamics simulations and further derive an elastic network model based on the statistical scores of AlphFold2 (AF-ENM), which we benchmark in combination with coarse-grained molecular dynamics simulations. We show that our AF-ENM method reproduces the global protein dynamics with improved accuracy, providing a powerful way to derive effective molecular dynamics using neural network-based structure prediction models.
Assuntos
Simulação de Dinâmica Molecular , Proteínas , Proteínas/química , Redes Neurais de Computação , Conformação Molecular , Movimento (Física) , Conformação ProteicaRESUMO
Electron bifurcation is a fundamental energy coupling mechanism widespread in microorganisms that thrive under anoxic conditions. These organisms employ hydrogen to reduce CO2, but the molecular mechanisms have remained enigmatic. The key enzyme responsible for powering these thermodynamically challenging reactions is the electron-bifurcating [FeFe]-hydrogenase HydABC that reduces low-potential ferredoxins (Fd) by oxidizing hydrogen gas (H2). By combining single-particle cryo-electron microscopy (cryoEM) under catalytic turnover conditions with site-directed mutagenesis experiments, functional studies, infrared spectroscopy, and molecular simulations, we show that HydABC from the acetogenic bacteria Acetobacterium woodii and Thermoanaerobacter kivui employ a single flavin mononucleotide (FMN) cofactor to establish electron transfer pathways to the NAD(P)+ and Fd reduction sites by a mechanism that is fundamentally different from classical flavin-based electron bifurcation enzymes. By modulation of the NAD(P)+ binding affinity via reduction of a nearby iron-sulfur cluster, HydABC switches between the exergonic NAD(P)+ reduction and endergonic Fd reduction modes. Our combined findings suggest that the conformational dynamics establish a redox-driven kinetic gate that prevents the backflow of the electrons from the Fd reduction branch toward the FMN site, providing a basis for understanding general mechanistic principles of electron-bifurcating hydrogenases.
Assuntos
Elétrons , Hidrogenase , Hidrogenase/química , NAD/metabolismo , Microscopia Crioeletrônica , Ferredoxinas/química , Oxirredução , Hidrogênio/química , Transporte de ElétronsRESUMO
Zn2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn2+ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn2+ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-ß-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn2+ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn2+ site geometries that could guide the further development of biomolecular Zn2+ models.
Assuntos
Metaloproteínas , Zinco , Zinco/química , Ligantes , Benchmarking , Sítios de Ligação , Metaloproteínas/químicaRESUMO
Redox reactions are central to biochemistry and are both controlled by and induce protein structural changes. Here, we describe structural rearrangements and crosstalk within the Bacillus cereus ribonucleotide reductase R2b-NrdI complex, a di-metal carboxylate-flavoprotein system, as part of the mechanism generating the essential catalytic free radical of the enzyme. Femtosecond crystallography at an X-ray free electron laser was utilized to obtain structures at room temperature in defined redox states without suffering photoreduction. Together with density functional theory calculations, we show that the flavin is under steric strain in the R2b-NrdI protein complex, likely tuning its redox properties to promote superoxide generation. Moreover, a binding site in close vicinity to the expected flavin O2 interaction site is observed to be controlled by the redox state of the flavin and linked to the channel proposed to funnel the produced superoxide species from NrdI to the di-manganese site in protein R2b. These specific features are coupled to further structural changes around the R2b-NrdI interaction surface. The mechanistic implications for the control of reactive oxygen species and radical generation in protein R2b are discussed.
Assuntos
Ribonucleotídeo Redutases , Cristalografia por Raios X , Flavinas/metabolismo , Oxirredução , Ribonucleotídeo Redutases/química , SuperóxidosRESUMO
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 O2 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 a3 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.
Assuntos
Complexo IV da Cadeia de Transporte de Elétrons , Prótons , Aerobiose , Transporte de Elétrons , Complexo IV da Cadeia de Transporte de Elétrons/química , Oxigênio/metabolismo , Água/químicaRESUMO
Dioxygenases catalyze stereoselective oxygen atom transfer in metabolic pathways of biological, industrial, and pharmaceutical importance, but their precise chemical principles remain controversial. The α-ketoglutarate (αKG)-dependent dioxygenase AsqJ synthesizes biomedically active quinolone alkaloids via desaturation and subsequent epoxidation of a carbon-carbon bond in the cyclopeptin substrate. Here, we combine high-resolution X-ray crystallography with enzyme engineering, quantum-classical (QM/MM) simulations, and biochemical assays to describe a peroxidic intermediate that bridges the substrate and active site metal ion in AsqJ. Homolytic cleavage of this moiety during substrate epoxidation generates an activated high-valent ferryl (FeIV = O) species that mediates the next catalytic cycle, possibly without the consumption of the metabolically valuable αKG cosubstrate. Our combined findings provide an important understanding of chemical bond activation principles in complex enzymatic reaction networks and molecular mechanisms of dioxygenases.
Assuntos
Dioxigenases , Carbono , Catálise , Domínio Catalítico , Dioxigenases/química , Ácidos Cetoglutáricos/metabolismo , Oxigênio/químicaRESUMO
The mitochondrial electron transport chain maintains the proton motive force that powers adenosine triphosphate (ATP) synthesis. The energy for this process comes from oxidation of reduced nicotinamide adenine dinucleotide (NADH) and succinate, with the electrons from this oxidation passed via intermediate carriers to oxygen. Complex IV (CIV), the terminal oxidase, transfers electrons from the intermediate electron carrier cytochrome c to oxygen, contributing to the proton motive force in the process. Within CIV, protons move through the K and D pathways during turnover. The former is responsible for transferring two protons to the enzyme's catalytic site upon its reduction, where they eventually combine with oxygen and electrons to form water. CIV is the main site for respiratory regulation, and although previous studies showed that steroid binding can regulate CIV activity, little is known about how this regulation occurs. Here, we characterize the interaction between CIV and steroids using a combination of kinetic experiments, structure determination, and molecular simulations. We show that molecules with a sterol moiety, such as glyco-diosgenin and cholesteryl hemisuccinate, reversibly inhibit CIV. Flash photolysis experiments probing the rapid equilibration of electrons within CIV demonstrate that binding of these molecules inhibits proton uptake through the K pathway. Single particle cryogenic electron microscopy (cryo-EM) of CIV with glyco-diosgenin reveals a previously undescribed steroid binding site adjacent to the K pathway, and molecular simulations suggest that the steroid binding modulates the conformational dynamics of key residues and proton transfer kinetics within this pathway. The binding pose of the sterol group sheds light on possible structural gating mechanisms in the CIV catalytic cycle.
Assuntos
Diosgenina , Complexo IV da Cadeia de Transporte de Elétrons , Esteroides , Animais , Sítios de Ligação , Domínio Catalítico/efeitos dos fármacos , Bovinos , Diosgenina/farmacologia , Transporte de Elétrons , Complexo IV da Cadeia de Transporte de Elétrons/antagonistas & inibidores , Complexo IV da Cadeia de Transporte de Elétrons/química , Oxirredução , Oxigênio/metabolismo , Conformação Proteica , Prótons , Esteroides/química , Esteroides/farmacologia , EsteróisRESUMO
The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a ß-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.
Assuntos
Proteínas de Choque Térmico HSP90 , Chaperonas Moleculares , Proteínas de Choque Térmico HSP90/metabolismo , Modelos Moleculares , Chaperonas Moleculares/metabolismo , Simulação de Dinâmica Molecular , Conformação Proteica , Dobramento de ProteínaRESUMO
Photosystem II (PSII) catalyzes light-driven water oxidization, releasing O2 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 Mn4O5Ca 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 Tyrz (Tyr161) lowers the reaction barrier for the water-mediated proton transfer from a Ca2+-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 Mn4O5Ca cluster lowers the proton transfer barrier through the water ligand sphere of the Mn4O5Ca cluster to Asp61 via a similar ion-pair dissociation process, while the resulting Mn-bound oxo/oxyl species leads to O2 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.
Assuntos
Complexo de Proteína do Fotossistema II , Prótons , Elétrons , Oxirredução , Oxigênio/química , Complexo de Proteína do Fotossistema II/química , Água/químicaRESUMO
Photosystem II (PSII), the water/plastoquinone photo-oxidoreductase, plays a key energy input role in the biosphere. [Formula: see text], the reduced semiquinone form of the nonexchangeable quinone, is often considered capable of a side reaction with O2, forming superoxide, but this reaction has not yet been demonstrated experimentally. Here, using chlorophyll fluorescence in plant PSII membranes, we show that O2 does oxidize [Formula: see text] at physiological O2 concentrations with a t1/2 of 10 s. Superoxide is formed stoichiometrically, and the reaction kinetics are controlled by the accessibility of O2 to a binding site near [Formula: see text], with an apparent dissociation constant of 70 ± 20 µM. Unexpectedly, [Formula: see text] could only reduce O2 when bicarbonate was absent from its binding site on the nonheme iron (Fe2+) and the addition of bicarbonate or formate blocked the O2-dependant decay of [Formula: see text] These results, together with molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations, indicate that electron transfer from [Formula: see text] to O2 occurs when the O2 is bound to the empty bicarbonate site on Fe2+ A protective role for bicarbonate in PSII was recently reported, involving long-lived [Formula: see text] triggering bicarbonate dissociation from Fe2+ [Brinkert et al, Proc. Natl. Acad. Sci. U.S.A. 113, 12144-12149 (2016)]. The present findings extend this mechanism by showing that bicarbonate release allows O2 to bind to Fe2+ and to oxidize [Formula: see text] This could be beneficial by oxidizing [Formula: see text] and by producing superoxide, a chemical signal for the overreduced state of the electron transfer chain.
Assuntos
Bicarbonatos/metabolismo , Oxigênio/metabolismo , Complexo de Proteína do Fotossistema II/metabolismo , Chlamydomonas reinhardtii/metabolismo , Clorofila/metabolismo , Transporte de Elétrons/fisiologia , Formiatos/metabolismo , Oxirredução , Quinonas/metabolismo , Spinacia oleracea/metabolismoRESUMO
Biological energy conversion is catalyzed by membrane-bound proteins that transduce chemical or light energy into energy forms that power endergonic processes in the cell. At a molecular level, these catalytic processes involve elementary electron-, proton-, charge-, and energy-transfer reactions that take place in the intricate molecular machineries of cell respiration and photosynthesis. Recent developments in structural biology, particularly cryo-electron microscopy (cryoEM), have resolved the molecular architecture of several energy transducing proteins, but detailed mechanistic principles of their charge transfer reactions still remain poorly understood and a major challenge for modern biochemical research. To this end, multiscale molecular simulations provide a powerful approach to probe mechanistic principles on a broad range of time scales (femtoseconds to milliseconds) and spatial resolutions (101-106 atoms), although technical challenges also require balancing between the computational accuracy, cost, and approximations introduced within the model. Here we discuss how the combination of atomistic (aMD) and hybrid quantum/classical molecular dynamics (QM/MM MD) simulations with free energy (FE) sampling methods can be used to probe mechanistic principles of enzymes responsible for biological energy conversion. We present mechanistic explorations of long-range proton-coupled electron transfer (PCET) dynamics in the highly intricate respiratory chain enzyme Complex I, which functions as a redox-driven proton pump in bacterial and mitochondrial respiratory chains by catalyzing a 300 Å fully reversible PCET process. This process is initiated by a hydride (H-) transfer between NADH and FMN, followed by long-range (>100 Å) electron transfer along a wire of 8 FeS centers leading to a quinone biding site. The reduction of the quinone to quinol initiates dissociation of the latter to a second membrane-bound binding site, and triggers proton pumping across the membrane domain of complex I, in subunits up to 200 Å away from the active site. Our simulations across different size and time scales suggest that transient charge transfer reactions lead to changes in the internal hydration state of key regions, local electric fields, and the conformation of conserved ion pairs, which in turn modulate the dynamics of functional steps along the reaction cycle. Similar functional principles, which operate on much shorter length scales, are also found in some unrelated proteins, suggesting that enzymes may employ conserved principles in the catalysis of biological energy transduction processes.
Assuntos
Complexo I de Transporte de Elétrons , Prótons , Respiração Celular , Microscopia Crioeletrônica , Transporte de Elétrons , Complexo I de Transporte de Elétrons/metabolismo , ElétronsRESUMO
The membrane-bound hydrogenase (Mbh) is a redox-driven Na+/H+ transporter that employs the energy from hydrogen gas (H2) production to catalyze proton pumping and Na+/H+ exchange across cytoplasmic membranes of archaea. Despite a recently resolved structure of this ancient energy-transducing enzyme [Yu et al. Cell 2018, 173, 1636-1649], the molecular principles of its redox-driven ion-transport mechanism remain puzzling and of major interest for understanding bioenergetic principles of early cells. Here we use atomistic molecular dynamics (MD) simulations in combination with data clustering methods and quantum chemical calculations to probe principles underlying proton reduction as well as proton and sodium transport in Mbh from the hyperthermophilic archaeon Pyrococcus furiosus. We identify putative Na+ binding sites and proton pathways leading across the membrane and to the NiFe-active center as well as conformational changes that regulate ion uptake. We suggest that Na+ binding and protonation changes at a putative ion-binding site couple to proton transfer across the antiporter-like MbhH subunit by modulating the conformational state of a conserved ion pair at the subunit interface. Our findings illustrate conserved coupling principles within the complex I superfamily and provide functional insight into archaeal energy transduction mechanisms.
Assuntos
Proteínas Arqueais/química , Hidrogenase/química , Trocadores de Sódio-Hidrogênio/química , Proteínas Arqueais/metabolismo , Catálise , Domínio Catalítico , Hidrogenase/metabolismo , Transporte de Íons , Simulação de Dinâmica Molecular , Ligação Proteica , Prótons , Pyrococcus furiosus/enzimologia , Sódio/química , Sódio/metabolismo , Trocadores de Sódio-Hidrogênio/metabolismo , Água/química , Água/metabolismoRESUMO
Cellular respiration is powered by membrane-bound redox enzymes that convert chemical energy into an electrochemical proton gradient and drive the energy metabolism. By combining large-scale classical and quantum mechanical simulations with cryo-electron microscopy data, we resolve here molecular details of conformational changes linked to proton pumping in the mammalian complex I. Our data suggest that complex I deactivation blocks water-mediated proton transfer between a membrane-bound quinone site and proton-pumping modules, decoupling the energy-transduction machinery. We identify a putative gating region at the interface between membrane domain subunits ND1 and ND3/ND4L/ND6 that modulates the proton transfer by conformational changes in transmembrane helices and bulky residues. The region is perturbed by mutations linked to human mitochondrial disorders and is suggested to also undergo conformational changes during catalysis of simpler complex I variants that lack the "active"-to-"deactive" transition. Our findings suggest that conformational changes in transmembrane helices modulate the proton transfer dynamics by wetting/dewetting transitions and provide important functional insight into the mammalian respiratory complex I.