1
|
Rigobello L, Lugli F, Caporali L, Bartocci A, Fadanni J, Zerbetto F, Iommarini L, Carelli V, Ghelli AM, Musiani F. A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I. Int J Biol Macromol 2024; 273:133086. [PMID: 38871105 DOI: 10.1016/j.ijbiomac.2024.133086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 06/07/2024] [Accepted: 06/09/2024] [Indexed: 06/15/2024]
Abstract
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.
Collapse
Affiliation(s)
- Laura Rigobello
- Department of Pharmacy and Biotechnology, University of Bologna, Bologna I-40127, Italy
| | - Francesca Lugli
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna I-40126, Italy.
| | - Leonardo Caporali
- IRCCS Istituto delle Scienze Neurologiche di Bologna, Programma di Neurogenetica, Bologna I-40124, Italy
| | - Alessio Bartocci
- Department of Physics, University of Trento, Trento I-38123, Italy; INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento I-38123, Italy
| | - Jacopo Fadanni
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna I-40126, Italy
| | - Francesco Zerbetto
- Department of Chemistry "Giacomo Ciamician", University of Bologna, Bologna I-40126, Italy
| | - Luisa Iommarini
- Department of Pharmacy and Biotechnology, University of Bologna, Bologna I-40127, Italy
| | - Valerio Carelli
- IRCCS Istituto delle Scienze Neurologiche di Bologna, Programma di Neurogenetica, Bologna I-40124, Italy; Department of Biomedical and Neuromotor Sciences, University of Bologna, Bologna I-40123, Italy
| | - Anna Maria Ghelli
- Department of Pharmacy and Biotechnology, University of Bologna, Bologna I-40127, Italy; IRCCS Istituto delle Scienze Neurologiche di Bologna, Programma di Neurogenetica, Bologna I-40124, Italy
| | - Francesco Musiani
- Department of Pharmacy and Biotechnology, University of Bologna, Bologna I-40127, Italy.
| |
Collapse
|
2
|
Pirnia A, Maqdisi R, Mittal S, Sener M, Singharoy A. Perspective on Integrative Simulations of Bioenergetic Domains. J Phys Chem B 2024; 128:3302-3319. [PMID: 38562105 DOI: 10.1021/acs.jpcb.3c07335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.
Collapse
Affiliation(s)
- Adam Pirnia
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Ranel Maqdisi
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Sumit Mittal
- VIT Bhopal University, Sehore 466114, Madhya Pradesh, India
| | - Melih Sener
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| |
Collapse
|
3
|
Djurabekova A, Lasham J, Zdorevskyi O, Zickermann V, Sharma V. Long-range electron proton coupling in respiratory complex I - insights from molecular simulations of the quinone chamber and antiporter-like subunits. Biochem J 2024; 481:499-514. [PMID: 38572757 DOI: 10.1042/bcj20240009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Revised: 03/11/2024] [Accepted: 03/14/2024] [Indexed: 04/05/2024]
Abstract
Respiratory complex I is a redox-driven proton pump. Several high-resolution structures of complex I have been determined providing important information about the putative proton transfer paths and conformational transitions that may occur during catalysis. However, how redox energy is coupled to the pumping of protons remains unclear. In this article, we review biochemical, structural and molecular simulation data on complex I and discuss several coupling models, including the key unresolved mechanistic questions. Focusing both on the quinone-reductase domain as well as the proton-pumping membrane-bound domain of complex I, we discuss a molecular mechanism of proton pumping that satisfies most experimental and theoretical constraints. We suggest that protonation reactions play an important role not only in catalysis, but also in the physiologically-relevant active/deactive transition of complex I.
Collapse
Affiliation(s)
| | - Jonathan Lasham
- Department of Physics, University of Helsinki, Helsinki, Finland
| | | | - Volker Zickermann
- Institute of Biochemistry II, University Hospital, Goethe University, Frankfurt am Main, Germany
- Centre for Biomolecular Magnetic Resonance, Institute for Biophysical Chemistry, Goethe University, Frankfurt am Main, Germany
| | - Vivek Sharma
- Department of Physics, University of Helsinki, Helsinki, Finland
- HiLIFE Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| |
Collapse
|
4
|
Kim H, Saura P, Pöverlein MC, Gamiz-Hernandez AP, Kaila VRI. Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I. J Am Chem Soc 2023; 145:17075-17086. [PMID: 37490414 PMCID: PMC10416309 DOI: 10.1021/jacs.3c03086] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2023] [Indexed: 07/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Hyunho Kim
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm 10691, Sweden
| | - Patricia Saura
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm 10691, Sweden
| | | | - Ana P. Gamiz-Hernandez
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm 10691, Sweden
| | - Ville R. I. Kaila
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm 10691, Sweden
| |
Collapse
|
5
|
Melo MCR, Bernardi RC. Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike. Biophys J 2023; 122:2833-2840. [PMID: 36738105 PMCID: PMC10398237 DOI: 10.1016/j.bpj.2023.01.042] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 01/24/2023] [Accepted: 01/30/2023] [Indexed: 02/05/2023] Open
Abstract
Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.
Collapse
Affiliation(s)
- Marcelo C R Melo
- Auburn University, Department of Physics, Auburn University, Auburn, Alabama
| | - Rafael C Bernardi
- Auburn University, Department of Physics, Auburn University, Auburn, Alabama.
| |
Collapse
|
6
|
Yu J, Yan C, Dodd T, Tsai CL, Tainer JA, Tsutakawa SE, Ivanov I. Dynamic conformational switching underlies TFIIH function in transcription and DNA repair and impacts genetic diseases. Nat Commun 2023; 14:2758. [PMID: 37179334 PMCID: PMC10183003 DOI: 10.1038/s41467-023-38416-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 04/28/2023] [Indexed: 05/15/2023] Open
Abstract
Transcription factor IIH (TFIIH) is a protein assembly essential for transcription initiation and nucleotide excision repair (NER). Yet, understanding of the conformational switching underpinning these diverse TFIIH functions remains fragmentary. TFIIH mechanisms critically depend on two translocase subunits, XPB and XPD. To unravel their functions and regulation, we build cryo-EM based TFIIH models in transcription- and NER-competent states. Using simulations and graph-theoretical analysis methods, we reveal TFIIH's global motions, define TFIIH partitioning into dynamic communities and show how TFIIH reshapes itself and self-regulates depending on functional context. Our study uncovers an internal regulatory mechanism that switches XPB and XPD activities making them mutually exclusive between NER and transcription initiation. By sequentially coordinating the XPB and XPD DNA-unwinding activities, the switch ensures precise DNA incision in NER. Mapping TFIIH disease mutations onto network models reveals clustering into distinct mechanistic classes, affecting translocase functions, protein interactions and interface dynamics.
Collapse
Affiliation(s)
- Jina Yu
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Chunli Yan
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Thomas Dodd
- Department of Chemistry, Georgia State University, Atlanta, GA, USA
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA
| | - Chi-Lin Tsai
- Department of Molecular and Cellular Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - John A Tainer
- Department of Molecular and Cellular Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Susan E Tsutakawa
- Molecular Biophysics and Integrated Bioimaging, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Ivaylo Ivanov
- Department of Chemistry, Georgia State University, Atlanta, GA, USA.
- Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, GA, USA.
| |
Collapse
|
7
|
Wikström M, Djurabekova A, Sharma V. On the role of ubiquinone in the proton translocation mechanism of respiratory complex I. FEBS Lett 2023; 597:224-236. [PMID: 36180980 DOI: 10.1002/1873-3468.14506] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 09/23/2022] [Accepted: 09/23/2022] [Indexed: 01/26/2023]
Abstract
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 (QH2 ), 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.
Collapse
Affiliation(s)
- Mårten Wikström
- HiLIFE Institute of Biotechnology, University of Helsinki, Finland
| | | | - Vivek Sharma
- HiLIFE Institute of Biotechnology, University of Helsinki, Finland.,Department of Physics, University of Helsinki, Finland
| |
Collapse
|
8
|
Kampjut D, Sazanov LA. Structure of respiratory complex I – An emerging blueprint for the mechanism. Curr Opin Struct Biol 2022; 74:102350. [PMID: 35316665 PMCID: PMC7613608 DOI: 10.1016/j.sbi.2022.102350] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 01/25/2022] [Accepted: 02/08/2022] [Indexed: 11/26/2022]
Abstract
Complex I is one of the major respiratory complexes, conserved from bacteria to mammals. It oxidises NADH, reduces quinone and pumps protons across the membrane, thus playing a central role in the oxidative energy metabolism. In this review we discuss our current state of understanding the structure of complex I from various species of mammals, plants, fungi, and bacteria, as well as of several complex I-related proteins. By comparing the structural evidence from these systems in different redox states and data from mutagenesis and molecular simulations, we formulate the mechanisms of electron transfer and proton pumping and explain how they are conformationally and electrostatically coupled. Finally, we discuss the structural basis of the deactivation phenomenon in mammalian complex I.
Collapse
|
9
|
Braasch-Turi MM, Koehn JT, Kostenkova K, Van Cleave C, Ives JW, Murakami HA, Crick DC, Crans DC. Electron Transport Lipids Fold Within Membrane-Like Interfaces. Front Chem 2022; 10:827530. [PMID: 35350775 PMCID: PMC8957872 DOI: 10.3389/fchem.2022.827530] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Accepted: 01/07/2022] [Indexed: 12/23/2022] Open
Abstract
Lipoquinones, such as ubiquinones (UQ) and menaquinones (MK), function as essential lipid components of the electron transport system (ETS) by shuttling electrons and protons to facilitate the production of ATP in eukaryotes and prokaryotes. Lipoquinone function in membrane systems has been widely studied, but the exact location and conformation within membranes remains controversial. Lipoquinones, such as Coenzyme Q (UQ-10), are generally depicted simply as "Q" in life science diagrams or in extended conformations in primary literature even though specific conformations are important for function in the ETS. In this study, our goal was to determine the location, orientation, and conformation of UQ-2, a truncated analog of UQ-10, in model membrane systems and to compare our results to previously studied MK-2. Herein, we first carried out a six-step synthesis to yield UQ-2 and then demonstrated that UQ-2 adopts a folded conformation in organic solvents using 1H-1H 2D NOESY and ROESY NMR spectroscopic studies. Similarly, using 1H-1H 2D NOESY NMR spectroscopic studies, UQ-2 was found to adopt a folded, U-shaped conformation within the interface of an AOT reverse micelle model membrane system. UQ-2 was located slightly closer to the surfactant-water interface compared to the more hydrophobic MK-2. In addition, Langmuir monolayer studies determined UQ-2 resided within the monolayer water-phospholipid interface causing expansion, whereas MK-2 was more likely to be compressed out and reside within the phospholipid tails. All together these results support the model that lipoquinones fold regardless of the headgroup structure but that the polarity of the headgroup influences lipoquinone location within the membrane interface. These results have implications regarding the redox activity near the interface as quinone vs. quinol forms may facilitate locomotion of lipoquinones within the membrane. The location, orientation, and conformation of lipoquinones are critical for their function in generating cellular energy within membrane ETS, and the studies described herein shed light on the behavior of lipoquinones within membrane-like environments.
Collapse
Affiliation(s)
| | - Jordan T. Koehn
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
| | - Kateryna Kostenkova
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
| | - Cameron Van Cleave
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
| | - Jacob W. Ives
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
| | - Heide A. Murakami
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
| | - Dean C. Crick
- Cell and Molecular Biology Program, Colorado State University, Fort Collins, CO, United States
- Microbiology, Immunology, and Pathology Department, Colorado State University, Fort Collins, CO, United States
| | - Debbie C. Crans
- Chemistry Department, Colorado State University, Fort Collins, CO, United States
- Cell and Molecular Biology Program, Colorado State University, Fort Collins, CO, United States
| |
Collapse
|
10
|
Kaila VRI. Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I. Acc Chem Res 2021; 54:4462-4473. [PMID: 34894649 PMCID: PMC8697550 DOI: 10.1021/acs.accounts.1c00524] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
![]()
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.
Collapse
Affiliation(s)
- Ville R. I. Kaila
- Department of Biochemistry and Biophysics, Stockholm University, 10691 Stockholm, Sweden
| |
Collapse
|
11
|
Schimpf J, Oppermann S, Gerasimova T, Santos Seica AF, Hellwig P, Grishkovskaya I, Wohlwend D, Haselbach D, Friedrich T. Structure of the peripheral arm of a minimalistic respiratory complex I. Structure 2021; 30:80-94.e4. [PMID: 34562374 DOI: 10.1016/j.str.2021.09.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Revised: 07/09/2021] [Accepted: 09/08/2021] [Indexed: 10/20/2022]
Abstract
Respiratory complex I drives proton translocation across energy-transducing membranes by NADH oxidation coupled with (ubi)quinone reduction. In humans, its dysfunction is associated with neurodegenerative diseases. The Escherichia coli complex represents the structural minimal form of an energy-converting NADH:ubiquinone oxidoreductase. Here, we report the structure of the peripheral arm of the E. coli complex I consisting of six subunits, the FMN cofactor, and nine iron-sulfur clusters at 2.7 Å resolution obtained by cryo electron microscopy. While the cofactors are in equivalent positions as in the complex from other species, individual subunits are adapted to the absence of supernumerary proteins to guarantee structural stability. The catalytically important subunits NuoC and D are fused resulting in a specific architecture of functional importance. Striking features of the E. coli complex are scrutinized by mutagenesis and biochemical characterization of the variants. Moreover, the arrangement of the subunits sheds light on the unknown assembly of the complex.
Collapse
Affiliation(s)
- Johannes Schimpf
- Institut für Biochemie, Albert-Ludwigs-Universität Freiburg, Albertstrasse 21, 79104 Freiburg, Germany
| | - Sabrina Oppermann
- Institut für Biochemie, Albert-Ludwigs-Universität Freiburg, Albertstrasse 21, 79104 Freiburg, Germany
| | - Tatjana Gerasimova
- Institut für Biochemie, Albert-Ludwigs-Universität Freiburg, Albertstrasse 21, 79104 Freiburg, Germany; Laboratoire de Bioélectrochimie et Spectroscopie, UMR 7140, CMC, Université de Strasbourg CNRS, 4 Rue Blaise Pascal, 67081 Strasbourg, France
| | - Ana Filipa Santos Seica
- Laboratoire de Bioélectrochimie et Spectroscopie, UMR 7140, CMC, Université de Strasbourg CNRS, 4 Rue Blaise Pascal, 67081 Strasbourg, France
| | - Petra Hellwig
- Laboratoire de Bioélectrochimie et Spectroscopie, UMR 7140, CMC, Université de Strasbourg CNRS, 4 Rue Blaise Pascal, 67081 Strasbourg, France; University of Strasbourg, Institute for Advanced Studies (USIAS), 5 Allée du Général Rouvillois, 67083 Strasbourg, France
| | - Irina Grishkovskaya
- Research Institute of Molecular Pathology (IMP), Campus-Vienna-Biocenter 1, 1030 Vienna, Austria
| | - Daniel Wohlwend
- Institut für Biochemie, Albert-Ludwigs-Universität Freiburg, Albertstrasse 21, 79104 Freiburg, Germany
| | - David Haselbach
- Research Institute of Molecular Pathology (IMP), Campus-Vienna-Biocenter 1, 1030 Vienna, Austria
| | - Thorsten Friedrich
- Institut für Biochemie, Albert-Ludwigs-Universität Freiburg, Albertstrasse 21, 79104 Freiburg, Germany.
| |
Collapse
|
12
|
Van Cleave C, Koehn JT, Pereira CS, Haase AA, Peters BJ, Croslow SW, McLaughlin KG, Werst KR, Goach AL, Crick DC, Arantes GM, Crans DC. Interactions of Truncated Menaquinones in Lipid Monolayers and Bilayers. Int J Mol Sci 2021; 22:9755. [PMID: 34575937 PMCID: PMC8470443 DOI: 10.3390/ijms22189755] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Revised: 09/02/2021] [Accepted: 09/05/2021] [Indexed: 11/28/2022] Open
Abstract
Menaquinones (MK) are hydrophobic molecules that consist of a naphthoquinone headgroup and a repeating isoprenyl side chain and are cofactors used in bacterial electron transport systems to generate cellular energy. We have previously demonstrated that the folded conformation of truncated MK homologues, MK-1 and MK-2, in both solution and reverse micelle microemulsions depended on environment. There is little information on how MKs associate with phospholipids in a model membrane system and how MKs affect phospholipid organization. In this manuscript, we used a combination of Langmuir monolayer studies and molecular dynamics (MD) simulations to probe these questions on truncated MK homologues, MK-1 through MK-4 within a model membrane. We observed that truncated MKs reside farther away from the interfacial water than ubiquinones are are located closer to the phospholipid tails. We also observed that phospholipid packing does not change at physiological pressure in the presence of truncated MKs, though a difference in phospholipid packing has been observed in the presence of ubiquinones. We found through MD simulations that for truncated MKs, the folded conformation varied, but MKs location and association with the bilayer remained unchanged at physiological conditions regardless of side chain length. Combined, this manuscript provides fundamental information, both experimental and computational, on the location, association, and conformation of truncated MK homologues in model membrane environments relevant to bacterial energy production.
Collapse
Affiliation(s)
- Cameron Van Cleave
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
| | - Jordan T. Koehn
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
| | - Caroline Simões Pereira
- Department of Biochemistry, Institutio de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, São Paulo 05508-900, SP, Brazil; (C.S.P.); (G.M.A.)
| | - Allison A. Haase
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
| | - Benjamin J. Peters
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
| | - Seth W. Croslow
- Department of Chemistry, Monmouth College, Monmouth, IL 61462, USA; (S.W.C.); (K.G.M.); (A.L.G.)
| | - Kyle G. McLaughlin
- Department of Chemistry, Monmouth College, Monmouth, IL 61462, USA; (S.W.C.); (K.G.M.); (A.L.G.)
| | - Katarina R. Werst
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
| | - Audra L. Goach
- Department of Chemistry, Monmouth College, Monmouth, IL 61462, USA; (S.W.C.); (K.G.M.); (A.L.G.)
| | - Dean C. Crick
- Cell and Molecular Biology Program, Colorado State University, Fort Collins, CO 80523, USA;
- Department of Microbiology, Immunology and Pathology, Colorado State University, Fort Collins, CO 80523, USA
| | - Guilherme Menegon Arantes
- Department of Biochemistry, Institutio de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, São Paulo 05508-900, SP, Brazil; (C.S.P.); (G.M.A.)
| | - Debbie C. Crans
- Department of Chemistry, Colorado State University, Fort Collins, CO 80523, USA; (C.V.C.); (J.T.K.); (A.A.H.); (B.J.P.); (K.R.W.)
- Cell and Molecular Biology Program, Colorado State University, Fort Collins, CO 80523, USA;
| |
Collapse
|
13
|
Poor Person's pH Simulation of Membrane Proteins. Methods Mol Biol 2021. [PMID: 34302678 DOI: 10.1007/978-1-0716-1468-6_12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/01/2023]
Abstract
pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and Vo domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations.
Collapse
|
14
|
Abstract
The electron transport chain of mitochondria is initiated by the respiratory complex I that converts chemical energy into a proton motive force to power synthesis of adenosine triphosphate. On a chemical level, complex I catalyzes elementary electron and proton transfer processes that couple across large molecular distances of >300 Å. However, under low oxygen concentrations, the respiratory chain operates in reverse mode and produces harmful reactive oxygen species. To avoid cell damage, the mitochondrial complex I transitions into a deactive state that inhibits turnover by molecular principles that remain elusive. By combining large-scale molecular simulations with cryo-electron microscopy data, we show here that complex I deactivation blocks the communication between proton pumping and redox modules by conformational and hydration changes. 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.
Collapse
|
15
|
Kaur D, Khaniya U, Zhang Y, Gunner MR. Protein Motifs for Proton Transfers That Build the Transmembrane Proton Gradient. Front Chem 2021; 9:660954. [PMID: 34211960 PMCID: PMC8239185 DOI: 10.3389/fchem.2021.660954] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 05/31/2021] [Indexed: 11/13/2022] Open
Abstract
Biological membranes are barriers to polar molecules, so membrane embedded proteins control the transfers between cellular compartments. Protein controlled transport moves substrates and activates cellular signaling cascades. In addition, the electrochemical gradient across mitochondrial, bacterial and chloroplast membranes, is a key source of stored cellular energy. This is generated by electron, proton and ion transfers through proteins. The gradient is used to fuel ATP synthesis and to drive active transport. Here the mechanisms by which protons move into the buried active sites of Photosystem II (PSII), bacterial RCs (bRCs) and through the proton pumps, Bacteriorhodopsin (bR), Complex I and Cytochrome c oxidase (CcO), are reviewed. These proteins all use water filled proton transfer paths. The proton pumps, that move protons uphill from low to high concentration compartments, also utilize Proton Loading Sites (PLS), that transiently load and unload protons and gates, which block backflow of protons. PLS and gates should be synchronized so PLS proton affinity is high when the gate opens to the side with few protons and low when the path is open to the high concentration side. Proton transfer paths in the proteins we describe have different design features. Linear paths are seen with a unique entry and exit and a relatively straight path between them. Alternatively, paths can be complex with a tangle of possible routes. Likewise, PLS can be a single residue that changes protonation state or a cluster of residues with multiple charge and tautomer states.
Collapse
Affiliation(s)
- Divya Kaur
- Department of Chemistry, The Graduate Center, City University of New York, New York, NY, United States.,Department of Physics, City College of New York, New York, NY, United States
| | - Umesh Khaniya
- Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| | - Yingying Zhang
- Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| | - M R Gunner
- Department of Chemistry, The Graduate Center, City University of New York, New York, NY, United States.,Department of Physics, City College of New York, New York, NY, United States.,Department of Physics, The Graduate Center, City University of New York, New York, NY, United States
| |
Collapse
|
16
|
A conserved arginine residue is critical for stabilizing the N2 FeS cluster in mitochondrial complex I. J Biol Chem 2021; 296:100474. [PMID: 33640456 PMCID: PMC8042128 DOI: 10.1016/j.jbc.2021.100474] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 02/18/2021] [Accepted: 02/24/2021] [Indexed: 11/22/2022] Open
Abstract
Respiratory complex I (NADH:ubiquinone oxidoreductase), the first enzyme of the electron-transport chain, captures the free energy released by NADH oxidation and ubiquinone reduction to translocate protons across an energy-transducing membrane and drive ATP synthesis during oxidative phosphorylation. The cofactor that transfers the electrons directly to ubiquinone is an iron–sulfur cluster (N2) located in the NDUFS2/NUCM subunit. A nearby arginine residue (R121), which forms part of the second coordination sphere of the N2 cluster, is known to be posttranslationally dimethylated but its functional and structural significance are not known. Here, we show that mutations of this arginine residue (R121M/K) abolish the quinone-reductase activity, concomitant with disappearance of the N2 signature from the electron paramagnetic resonance (EPR) spectrum. Analysis of the cryo-EM structure of NDUFS2-R121M complex I at 3.7 Å resolution identified the absence of the cubane N2 cluster as the cause of the dysfunction, within an otherwise intact enzyme. The mutation further induced localized disorder in nearby elements of the quinone-binding site, consistent with the close connections between the cluster and substrate-binding regions. Our results demonstrate that R121 is required for the formation and/or stability of the N2 cluster and highlight the importance of structural analyses for mechanistic interpretation of biochemical and spectroscopic data on complex I variants.
Collapse
|
17
|
Pamplona R, Jové M, Mota-Martorell N, Barja G. Is the NDUFV2 subunit of the hydrophilic complex I domain a key determinant of animal longevity? FEBS J 2021; 288:6652-6673. [PMID: 33455045 DOI: 10.1111/febs.15714] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Revised: 12/02/2020] [Accepted: 01/14/2021] [Indexed: 12/18/2022]
Abstract
Complex I, a component of the electron transport chain, plays a central functional role in cell bioenergetics and the biology of free radicals. The structural and functional N module of complex I is one of the main sites of the generation of free radicals. The NDUFV2 subunit/N1a cluster is a component of this module. Furthermore, the rate of free radical production is linked to animal longevity. In this review, we explore the hypothesis that NDUFV2 is the only conserved core subunit designed with a regulatory function to ensure correct electron transfer and free radical production, that low gene expression and protein abundance of the NDUFV2 subunit is an evolutionary adaptation needed to achieve a longevity phenotype, and that these features are determinants of the lower free radical generation at the mitochondrial level and a slower rate of aging of long-lived animals.
Collapse
Affiliation(s)
- Reinald Pamplona
- Department of Experimental Medicine, University of Lleida-Lleida Biomedical Research Institute (UdL-IRBLleida), Lleida, Spain
| | - Mariona Jové
- Department of Experimental Medicine, University of Lleida-Lleida Biomedical Research Institute (UdL-IRBLleida), Lleida, Spain
| | - Natalia Mota-Martorell
- Department of Experimental Medicine, University of Lleida-Lleida Biomedical Research Institute (UdL-IRBLleida), Lleida, Spain
| | - Gustavo Barja
- Department of Genetics, Physiology and Microbiology, Complutense University of Madrid, Madrid, Spain
| |
Collapse
|
18
|
Röpke M, Saura P, Riepl D, Pöverlein MC, Kaila VRI. Functional Water Wires Catalyze Long-Range Proton Pumping in the Mammalian Respiratory Complex I. J Am Chem Soc 2020; 142:21758-21766. [PMID: 33325238 PMCID: PMC7785131 DOI: 10.1021/jacs.0c09209] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The respiratory complex I is a gigantic
(1 MDa) redox-driven proton
pump that reduces the ubiquinone pool and generates proton motive
force to power ATP synthesis in mitochondria. Despite resolved molecular
structures and biochemical characterization of the enzyme from multiple
organisms, its long-range (∼300 Å) proton-coupled electron
transfer (PCET) mechanism remains unsolved. We employ here microsecond
molecular dynamics simulations to probe the dynamics of the mammalian
complex I in combination with hybrid quantum/classical (QM/MM) free
energy calculations to explore how proton pumping reactions are triggered
within its 200 Å wide membrane domain. Our simulations predict
extensive hydration dynamics of the antiporter-like subunits in complex
I that enable lateral proton transfer reactions on a microsecond time
scale. We further show how the coupling between conserved ion pairs
and charged residues modulate the proton transfer dynamics, and how
transmembrane helices and gating residues control the hydration process.
Our findings suggest that the mammalian complex I pumps protons by
tightly linked conformational and electrostatic coupling principles.
Collapse
Affiliation(s)
- Michael Röpke
- Center for Integrated Protein Science Munich at the Department of Chemistry, Technical University of Munich, Lichtenbergstrasse 4, D85748 Garching, Germany
| | - Patricia Saura
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Daniel Riepl
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Maximilian C Pöverlein
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Ville R I Kaila
- Department of Biochemistry and Biophysics, Stockholm University, SE-106 91 Stockholm, Sweden.,Center for Integrated Protein Science Munich at the Department of Chemistry, Technical University of Munich, Lichtenbergstrasse 4, D85748 Garching, Germany
| |
Collapse
|
19
|
Haapanen O, Reidelbach M, Sharma V. Coupling of quinone dynamics to proton pumping in respiratory complex I. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2020; 1861:148287. [PMID: 32777307 DOI: 10.1016/j.bbabio.2020.148287] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 07/09/2020] [Accepted: 08/03/2020] [Indexed: 12/12/2022]
Abstract
Respiratory complex I (NADH:quinone oxidoreductase) plays a central role in generating the proton electrochemical gradient in mitochondrial and bacterial membranes, which is needed to generate ATP. Several high-resolution structures of complex I have been determined, revealing its intricate architecture and complementing the biochemical and biophysical studies. However, the molecular mechanism of long-range coupling between ubiquinone (Q) reduction and proton pumping is not known. Computer simulations have been applied to decipher the dynamics of Q molecule in the ~30 Å long Q tunnel. In this short report, we discuss the binding and dynamics of Q at computationally predicted Q binding sites, many of which are supported by structural data on complex I. We suggest that the binding of Q at these sites is coupled to proton pumping by means of conformational rearrangements in the conserved loops of core subunits.
Collapse
Affiliation(s)
- Outi Haapanen
- Department of Physics, University of Helsinki, Finland
| | | | - Vivek Sharma
- Department of Physics, University of Helsinki, Finland; HiLIFE Institute of Biotechnology, University of Helsinki, Finland.
| |
Collapse
|
20
|
Li F, Shrivastava IH, Hanlon P, Dagda RK, Gasanoff ES. Molecular Mechanism by which Cobra Venom Cardiotoxins Interact with the Outer Mitochondrial Membrane. Toxins (Basel) 2020; 12:E425. [PMID: 32605112 PMCID: PMC7404710 DOI: 10.3390/toxins12070425] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2020] [Revised: 06/23/2020] [Accepted: 06/26/2020] [Indexed: 12/29/2022] Open
Abstract
Cardiotoxin CTII from Najaoxiana cobra venom translocates to the intermembrane space (IMS) of mitochondria to disrupt the structure and function of the inner mitochondrial membrane. At low concentrations, CTII facilitates ATP-synthase activity, presumably via the formation of non-bilayer, immobilized phospholipids that are critical in modulating ATP-synthase activity. In this study, we investigated the effects of another cardiotoxin CTI from Najaoxiana cobra venom on the structure of mitochondrial membranes and on mitochondrial-derived ATP synthesis. By employing robust biophysical methods including 31P-NMR and 1H-NMR spectroscopy, we analyzed the effects of CTI and CTII on phospholipid packing and dynamics in model phosphatidylcholine (PC) membranes enriched with 2.5 and 5.0 mol% of cardiolipin (CL), a phospholipid composition that mimics that in the outer mitochondrial membrane (OMM). These experiments revealed that CTII converted a higher percentage of bilayer phospholipids to a non-bilayer and immobilized state and both cardiotoxins utilized CL and PC molecules to form non-bilayer structures. Furthermore, in order to gain further understanding on how cardiotoxins bind to mitochondrial membranes, we employed molecular dynamics (MD) and molecular docking simulations to investigate the molecular mechanisms by which CTII and CTI interactively bind with an in silico phospholipid membrane that models the composition similar to the OMM. In brief, MD studies suggest that CTII utilized the N-terminal region to embed the phospholipid bilayer more avidly in a horizontal orientation with respect to the lipid bilayer and thereby penetrate at a faster rate compared with CTI. Molecular dynamics along with the Autodock studies identified critical amino acid residues on the molecular surfaces of CTII and CTI that facilitated the long-range and short-range interactions of cardiotoxins with CL and PC. Based on our compiled data and our published findings, we provide a conceptual model that explains a molecular mechanism by which snake venom cardiotoxins, including CTI and CTII, interact with mitochondrial membranes to alter the mitochondrial membrane structure to either upregulate ATP-synthase activity or disrupt mitochondrial function.
Collapse
Affiliation(s)
- Feng Li
- STEM Program, Science Department, Chaoyang KaiWen Academy, Yard 46, 3rd Baoquan Street, Chaoyang District, Beijing 100018, China; (F.L.); (P.H.); (E.S.G.)
| | - Indira H. Shrivastava
- Department of Computational and System Biology, University of Pittsburgh, Pittsburgh, PA 15260, USA;
- Department of Environmental and Occupational Health, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | - Paul Hanlon
- STEM Program, Science Department, Chaoyang KaiWen Academy, Yard 46, 3rd Baoquan Street, Chaoyang District, Beijing 100018, China; (F.L.); (P.H.); (E.S.G.)
| | - Ruben K. Dagda
- Reno School of Medicine, Department of Pharmacology, University of Nevada, Reno, NV 89557, USA
| | - Edward S. Gasanoff
- STEM Program, Science Department, Chaoyang KaiWen Academy, Yard 46, 3rd Baoquan Street, Chaoyang District, Beijing 100018, China; (F.L.); (P.H.); (E.S.G.)
| |
Collapse
|
21
|
Khaniya U, Gupta C, Cai X, Mao J, Kaur D, Zhang Y, Singharoy A, Gunner MR. Hydrogen bond network analysis reveals the pathway for the proton transfer in the E-channel of T. thermophilus Complex I. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2020; 1861:148240. [PMID: 32531220 DOI: 10.1016/j.bbabio.2020.148240] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2019] [Revised: 05/19/2020] [Accepted: 06/03/2020] [Indexed: 10/24/2022]
Abstract
Complex I, NADH-ubiquinone oxidoreductase, is the first enzyme in the mitochondrial and bacterial aerobic respiratory chain. It pumps four protons through four transiently open pathways from the high pH, negative, N-side of the membrane to the positive, P-side driven by the exergonic transfer of electrons from NADH to a quinone. Three protons transfer through subunits descended from antiporters, while the fourth, E-channel is unique. The path through the E-channel is determined by a network analysis of hydrogen bonded pathways obtained by Monte Carlo sampling of protonation states, polar hydrogen orientation and water occupancy. Input coordinates are derived from molecular dynamics trajectories comparing oxidized, reduced (dihydro) and no menaquinone-8 (MQ). A complex proton transfer path from the N- to the P-side is found consisting of six clusters of highly connected hydrogen-bonded residues. The network connectivity depends on the presence of quinone and its redox state, supporting a role for this cofactor in coupling electron and proton transfers. The N-side is more organized with MQ-bound complex I facilitating proton entry, while the P-side is more connected in the apo-protein, facilitating proton exit. Subunit Nqo8 forms the core of the E channel; Nqo4 provides the N-side entry, Nqo7 and then Nqo10 join the pathway in the middle, while Nqo11 contributes to the P-side exit.
Collapse
Affiliation(s)
- Umesh Khaniya
- Department of Physics, City College of New York, New York 10031, USA; Department of Physics, The Graduate Center, City University of New York, New York 10016, USA
| | - Chitrak Gupta
- School of Molecular Sciences, Arizona State University, Tempe, AZ, USA; Biodesign Institute, Arizona State University, Tempe, AZ, USA
| | - Xiuhong Cai
- Department of Physics, City College of New York, New York 10031, USA; Department of Physics, The Graduate Center, City University of New York, New York 10016, USA
| | - Junjun Mao
- Department of Physics, City College of New York, New York 10031, USA
| | - Divya Kaur
- Department of Physics, City College of New York, New York 10031, USA; Department of Chemistry, The Graduate Center, City University of New York, New York 10016, USA
| | - Yingying Zhang
- Department of Physics, City College of New York, New York 10031, USA; Department of Physics, The Graduate Center, City University of New York, New York 10016, USA
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, AZ, USA; Biodesign Institute, Arizona State University, Tempe, AZ, USA
| | - M R Gunner
- Department of Physics, City College of New York, New York 10031, USA; Department of Physics, The Graduate Center, City University of New York, New York 10016, USA; Department of Chemistry, The Graduate Center, City University of New York, New York 10016, USA.
| |
Collapse
|