1
|
Wang P, Demaray J, Moroz S, Stuchebrukhov AA. Searching for proton transfer channels in respiratory complex I. Biophys J 2024:S0006-3495(24)00518-6. [PMID: 39095988 DOI: 10.1016/j.bpj.2024.07.041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 06/19/2024] [Accepted: 07/30/2024] [Indexed: 08/04/2024] Open
Abstract
We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.
Collapse
Affiliation(s)
- Panyue Wang
- Department of Chemistry, University of California at Davis, Davis, California
| | - Jackson Demaray
- Department of Chemistry, University of California at Davis, Davis, California
| | - Stanislav Moroz
- Department of Chemistry, University of California at Davis, Davis, California
| | | |
Collapse
|
2
|
Emperador S, Habbane M, López-Gallardo E, Del Rio A, Llobet L, Mateo J, Sanz-López AM, Fernández-García MJ, Sánchez-Tocino H, Benbunan-Ferreiro S, Calabuig-Goena M, Narvaez-Palazón C, Fernández-Vega B, González-Iglesias H, Urreizti R, Artuch R, Pacheu-Grau D, Bayona-Bafaluy P, Montoya J, Ruiz-Pesini E. Identification and characterization of a new pathologic mutation in a large Leber hereditary optic neuropathy pedigree. Orphanet J Rare Dis 2024; 19:148. [PMID: 38582886 PMCID: PMC10999093 DOI: 10.1186/s13023-024-03165-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2023] [Accepted: 03/30/2024] [Indexed: 04/08/2024] Open
Abstract
BACKGROUND Most patients suffering from Leber hereditary optic neuropathy carry one of the three classic pathologic mutations, but not all individuals with these genetic alterations develop the disease. There are different risk factors that modify the penetrance of these mutations. The remaining patients carry one of a set of very rare genetic variants and, it appears that, some of the risk factors that modify the penetrance of the classical pathologic mutations may also affect the phenotype of these other rare mutations. RESULTS We describe a large family including 95 maternally related individuals, showing 30 patients with Leber hereditary optic neuropathy. The mutation responsible for the phenotype is a novel transition, m.3734A > G, in the mitochondrial gene encoding the ND1 subunit of respiratory complex I. Molecular-genetic, biochemical and cellular studies corroborate the pathogenicity of this genetic change. CONCLUSIONS With the study of this family, we confirm that, also for this very rare mutation, sex and age are important factors modifying penetrance. Moreover, this pedigree offers an excellent opportunity to search for other genetic or environmental factors that additionally contribute to modify penetrance.
Collapse
Affiliation(s)
- Sonia Emperador
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Mouna Habbane
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Laboratoire Biologie Et Santé, Faculté Des Sciences Ben M'Sick, Hassan II University of Casablanca, 20670, Casablanca, Morocco
| | - Ester López-Gallardo
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Alejandro Del Rio
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
| | - Laura Llobet
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain
- Certest Biotec, 50840-San Mateo de Gállego, Zaragoza, Spain
| | - Javier Mateo
- Servicio de Oftalmología, Hospital Clínico Universitario Lozano Blesa, 50009, Zaragoza, Spain
| | - Ana María Sanz-López
- Servicio de Oftalmología, Hospital Universitario de Toledo, 45004, Toledo, Spain
| | | | | | - Sol Benbunan-Ferreiro
- Servicio de Oftalmología. Hospital Universitario Río Hortega, 47012, Valladolid, Spain
| | - María Calabuig-Goena
- Servicio de Oftalmología. Hospital Universitario Río Hortega, 47012, Valladolid, Spain
| | | | | | - Hector González-Iglesias
- Instituto de Productos Lácteos de Asturias, Consejo Superior de Investigaciones Científicas (IPLA-CSIC), 33300-Villaviciosa, Asturias, Spain
| | - Roser Urreizti
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
- Departament de Bioquímica Clínica, Institut de Recerca Sant Joan de Déu, 08950, Barcelona, Spain
| | - Rafael Artuch
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
- Departament de Bioquímica Clínica, Institut de Recerca Sant Joan de Déu, 08950, Barcelona, Spain
| | - David Pacheu-Grau
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Pilar Bayona-Bafaluy
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain
- Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, 50018, Zaragoza, Spain
| | - Julio Montoya
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain.
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain.
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain.
| | - Eduardo Ruiz-Pesini
- Departamento de Bioquímica, Biología Molecular y Celular, Universidad de Zaragoza, 50009- and 50013, Zaragoza, Spain.
- Instituto de Investigación Sanitaria (IIS) de Aragón, 50009, Zaragoza, Spain.
- Centro de Investigaciones Biomédicas en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, 28029, Madrid, Spain.
| |
Collapse
|
3
|
Lasham J, Djurabekova A, Zickermann V, Vonck J, Sharma V. Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures. J Phys Chem B 2024; 128:2304-2316. [PMID: 38430110 DOI: 10.1021/acs.jpcb.3c07421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2024]
Abstract
Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.
Collapse
Affiliation(s)
- Jonathan Lasham
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
| | - Amina Djurabekova
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
| | - Volker Zickermann
- Institute of Biochemistry II, University Hospital, Goethe University, 60590 Frankfurt am Main, Germany
- Centre for Biomolecular Magnetic Resonance, Institute for Biophysical Chemistry, Goethe University, 60438 Frankfurt am Main, Germany
| | - Janet Vonck
- Department of Structural Biology, Max Planck Institute of Biophysics, 60438 Frankfurt am Main, Germany
| | - Vivek Sharma
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland
- HiLIFE Institute of Biotechnology, University of Helsinki, 00014 Helsinki, Finland
| |
Collapse
|
4
|
Pathak B, Kesari S, Patwari GN. Enticing a Proton using Single Ammonia Molecule as Bait. J Phys Chem B 2024; 128:1022-1028. [PMID: 38240575 DOI: 10.1021/acs.jpcb.3c06761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
In microhydrated acid-solvent clusters, deprotonation of an acid is assisted by a critical number of solvent molecules and a solvent electric field. Born-Oppenheimer molecular dynamics simulations reveal that trifluoroacetic acid undergoes spontaneous proton transfer in water clusters, with the critical number being five. Acetic acid and phenol, on the other hand, do not dissociate even in the presence of a large number of water molecules (in excess of 40). The addition of a single ammonia molecule to the water cluster, which interacts directly with the protic group, lowers the critical number of solvent water molecules required for proton transfer to three and seven in the case of acetic acid and phenol, respectively. The population of the undissociated and the proton-transferred structures get dispersed to form separate islands on the electric field versus the O-H distance representation with the cusp representing the critical values. The critical electric fields for the spontaneous proton transfer are around 254, 237, and 318 MV cm-1 for trifluoroacetic acid, acetic acid, and phenol, respectively. In the case of phenol, the free energy profiles suggest that proton transfer to the ammonia moiety embedded in water promotes proton transfer efficiently due to the higher basicity of ammonia and enhanced hydrogen bonding network of solvent water, vis-à-vis phenol-ammonia clusters.
Collapse
Affiliation(s)
- Bijaya Pathak
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Shaivi Kesari
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - G Naresh Patwari
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
5
|
Grba DN, Chung I, Bridges HR, Agip ANA, Hirst J. Investigation of hydrated channels and proton pathways in a high-resolution cryo-EM structure of mammalian complex I. SCIENCE ADVANCES 2023; 9:eadi1359. [PMID: 37531432 PMCID: PMC10396290 DOI: 10.1126/sciadv.adi1359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 07/03/2023] [Indexed: 08/04/2023]
Abstract
Respiratory complex I, a key enzyme in mammalian metabolism, captures the energy released by reduction of ubiquinone by NADH to drive protons across the inner mitochondrial membrane, generating the proton-motive force for ATP synthesis. Despite remarkable advances in structural knowledge of this complicated membrane-bound enzyme, its mechanism of catalysis remains controversial. In particular, how ubiquinone reduction is coupled to proton pumping and the pathways and mechanisms of proton translocation are contested. We present a 2.4-Å resolution cryo-EM structure of complex I from mouse heart mitochondria in the closed, active (ready-to-go) resting state, with 2945 water molecules modeled. By analyzing the networks of charged and polar residues and water molecules present, we evaluate candidate pathways for proton transfer through the enzyme, for the chemical protons for ubiquinone reduction, and for the protons transported across the membrane. Last, we compare our data to the predictions of extant mechanistic models, and identify key questions to answer in future work to test them.
Collapse
|
6
|
Wei RJ, Khaniya U, Mao J, Liu J, Batista VS, Gunner MR. Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers. PHOTOSYNTHESIS RESEARCH 2023; 156:101-112. [PMID: 36307598 DOI: 10.1007/s11120-022-00973-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Accepted: 10/03/2022] [Indexed: 06/16/2023]
Abstract
Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the QB site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to QB. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with QB•- are compared with earlier results obtained with QB [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.
Collapse
Affiliation(s)
- Rongmei Judy Wei
- Ph.D. Program in Chemistry, The Graduate Center, City University of New York, New York, NY, 10016, USA
- Department of Physics, City College of New York, New York, NY, 10031, USA
| | - Umesh Khaniya
- Department of Physics, City College of New York, New York, NY, 10031, USA
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - Junjun Mao
- Department of Physics, City College of New York, New York, NY, 10031, USA
| | - Jinchan Liu
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT, 06520, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, CT, 06520, USA
| | - M R Gunner
- Ph.D. Program in Chemistry, The Graduate Center, City University of New York, New York, NY, 10016, USA.
- Department of Physics, City College of New York, New York, NY, 10031, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
7
|
Tunnel dynamics of quinone derivatives and its coupling to protein conformational rearrangements in respiratory complex I. BIOCHIMICA ET BIOPHYSICA ACTA. BIOENERGETICS 2023; 1864:148951. [PMID: 36509126 DOI: 10.1016/j.bbabio.2022.148951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Revised: 11/30/2022] [Accepted: 12/02/2022] [Indexed: 12/13/2022]
Abstract
Respiratory complex I in mitochondria and bacteria catalyzes the transfer of electrons from NADH to quinone (Q). The free energy available from the reaction is used to pump protons and to establish a membrane proton electrochemical gradient, which drives ATP synthesis. Even though several high-resolution structures of complex I have been resolved, how Q reduction is linked with proton pumping, remains unknown. Here, microsecond long molecular dynamics (MD) simulations were performed on Yarrowia lipolytica complex I structures where Q molecules have been resolved in the ~30 Å long Q tunnel. MD simulations of several different redox/protonation states of Q reveal the coupling between the Q dynamics and the restructuring of conserved loops and ion pairs. Oxidized quinone stabilizes towards the N2 FeS cluster, a binding mode not previously described in Yarrowia lipolytica complex I structures. On the other hand, reduced (and protonated) species tend to diffuse towards the Q binding sites closer to the tunnel entrance. Mechanistic and physiological relevance of these results are discussed.
Collapse
|
8
|
Alkhaldi HA, Vik SB. Subunits E-F-G of E. coli Complex I can form an active complex when expressed alone, but in time-delayed assembly co-expression of B-CD-E-F-G is optimal. BIOCHIMICA ET BIOPHYSICA ACTA. BIOENERGETICS 2022; 1863:148593. [PMID: 35850264 PMCID: PMC9783743 DOI: 10.1016/j.bbabio.2022.148593] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/06/2022] [Revised: 06/15/2022] [Accepted: 07/11/2022] [Indexed: 12/27/2022]
Abstract
Respiratory Complex I from E. coli is a proto-type of the mitochondrial enzyme, consisting of a 6-subunit peripheral arm (B-CD-E-F-G-I) and a 7-subunit membrane arm. When subunits E-F-G (N-module), were expressed alone they formed an active complex as determined by co-immunoprecipitation and native gel electrophoresis. When co-expressed with subunits B and CD, only a complex of E-F-G was found. When these five subunits were co-expressed with subunit I and two membrane subunits, A and H, a complex of B-CD-E-F-G-I was membrane-bound, constituting the N- and Q-modules. Assembly of Complex I was also followed by splitting the genes between two plasmids, in three different groupings, and expressing them simultaneously, or with time-delay of expression from one plasmid. When the B-CD-E-F-G genes were co-expressed after a time-delay, assembly was over 90 % of that when the whole operon was expressed together. In summary, E-F-G was the only soluble subcomplex detected in these studies, but assembly was not optimal when these subunits were expressed either first or last. Co-expression of subunits B and CD with E-F-G provided a higher level of assembly, indicating that integrated assembly of N- and Q-modules provides a more efficient pathway.
Collapse
Affiliation(s)
- Hind A Alkhaldi
- Department of Biological Sciences, Southern Methodist University, Dallas, TX 75275-0376, USA
| | - Steven B Vik
- Department of Biological Sciences, Southern Methodist University, Dallas, TX 75275-0376, USA.
| |
Collapse
|
9
|
Kubo S, Takada S. Rotational Mechanism of FO Motor in the F-Type ATP Synthase Driven by the Proton Motive Force. Front Microbiol 2022; 13:872565. [PMID: 35783438 PMCID: PMC9243769 DOI: 10.3389/fmicb.2022.872565] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 04/20/2022] [Indexed: 11/13/2022] Open
Abstract
In FOF1 ATP synthase, driven by the proton motive force across the membrane, the FO motor rotates the central rotor and induces conformational changes in the F1 motor, resulting in ATP synthesis. Recently, many near-atomic resolution structural models have been obtained using cryo-electron microscopy. Despite high resolution, however, static information alone cannot elucidate how and where the protons pass through the FO and how proton passage is coupled to FO rotation. Here, we review theoretical and computational studies based on FO structure models. All-atom molecular dynamics (MD) simulations elucidated changes in the protonation/deprotonation of glutamate—the protein-carrier residue—during rotation and revealed the protonation states that form the “water wire” required for long-range proton hopping. Coarse-grained MD simulations unveiled a free energy surface based on the protonation state and rotational angle of the rotor. Hybrid Monte Carlo and MD simulations showed how proton transfer is coupled to rotation.
Collapse
Affiliation(s)
- Shintaroh Kubo
- Department of Anatomy and Cell Biology, McGill University, Montreal, QC, Canada
- *Correspondence: Shintaroh Kubo
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- Shoji Takada
| |
Collapse
|
10
|
Wei RJ, Zhang Y, Mao J, Kaur D, Khaniya U, Gunner MR. Comparison of proton transfer paths to the Q A and Q B sites of the Rb. sphaeroides photosynthetic reaction centers. PHOTOSYNTHESIS RESEARCH 2022; 152:153-165. [PMID: 35344134 DOI: 10.1007/s11120-022-00906-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2021] [Accepted: 02/18/2022] [Indexed: 06/14/2023]
Abstract
The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the QA site cycles between quinone, QA, and anionic semiquinone, QA·-, being reduced once and never binding protons. In the QB site, ubiquinone is reduced twice by QA·-, binds two protons and is released into the membrane as the quinol, QH2. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. QA is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to QB. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as AspM17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake.
Collapse
Affiliation(s)
- Rongmei Judy Wei
- Ph.D. Program in Chemistry, The Graduate Center, City University of New York, New York, NY, 10016, USA
- Department of Physics, City College of New York, New York, NY, 10031, USA
| | - Yingying Zhang
- Department of Physics, City College of New York, New York, NY, 10031, USA
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - Junjun Mao
- Department of Physics, City College of New York, New York, NY, 10031, USA
| | - Divya Kaur
- Department of Chemistry, Brock University, 500 Glenridge Avenue, St. Catharines, ON, L2S 3A1, Canada
| | - Umesh Khaniya
- Department of Physics, City College of New York, New York, NY, 10031, USA
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - M R Gunner
- Ph.D. Program in Chemistry, The Graduate Center, City University of New York, New York, NY, 10016, USA.
- Department of Physics, City College of New York, New York, NY, 10031, USA.
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
11
|
Khaniya U, Mao J, Wei RJ, Gunner MR. Characterizing Protein Protonation Microstates Using Monte Carlo Sampling. J Phys Chem B 2022; 126:2476-2485. [PMID: 35344367 PMCID: PMC8997239 DOI: 10.1021/acs.jpcb.2c00139] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Proteins are polyelectrolytes with acidic and basic amino acids Asp, Glu, Arg, Lys, and His, making up ≈25% of the residues. The protonation state of residues, cofactors, and ligands defines a "protonation microstate". In an ensemble of proteins some residues will be ionized and others neutral, leading to a mixture of protonation microstates rather than in a single one as is often assumed. The microstate distribution changes with pH. The protein environment also modifies residue proton affinity so microstate distributions change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ΔG, ΔH, and ΔS between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. At pH 7 the lysozyme charge ranges from 6 to 10, with 24 accepted protonation microstates, while RCs have ≈50,000. A weighted Pearson correlation analysis shows coupling between residue protonation states in RCs and how they change when the quinone in the QB site is reduced. Protonation microstates can be used to define input MD parameters and provide insight into the motion of protons coupled to reactions.
Collapse
Affiliation(s)
- Umesh Khaniya
- Department of Physics, City College of New York, New York, New York 10031, United States.,Department of Physics, The Graduate Center, City University of New York, New York, New York 10016, United States
| | - Junjun Mao
- Department of Physics, City College of New York, New York, New York 10031, United States
| | - Rongmei Judy Wei
- Department of Physics, City College of New York, New York, New York 10031, United States.,Department of Chemistry, The Graduate Center, City University of New York, New York, New York 10016, United States
| | - M R Gunner
- Department of Physics, City College of New York, New York, New York 10031, United States.,Department of Physics, The Graduate Center, City University of New York, New York, New York 10016, United States.,Department of Chemistry, The Graduate Center, City University of New York, New York, New York 10016, United States
| |
Collapse
|
12
|
Al-Attar S, Rendon J, Sidore M, Duneau JP, Seduk F, Biaso F, Grimaldi S, Guigliarelli B, Magalon A. Gating of Substrate Access and Long-Range Proton Transfer in Escherichia coli Nitrate Reductase A: The Essential Role of a Remote Glutamate Residue. ACS Catal 2021. [DOI: 10.1021/acscatal.1c03988] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Sinan Al-Attar
- Laboratoire de Chimie Bactérienne (UMR7283), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Julia Rendon
- Laboratoire de Bioénergétique et Ingénierie des Protéines (UMR7281), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Marlon Sidore
- Laboratoire d’Ingénierie des Systèmes Macromoléculaires (UMR7255), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Jean-Pierre Duneau
- Laboratoire d’Ingénierie des Systèmes Macromoléculaires (UMR7255), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Farida Seduk
- Laboratoire de Chimie Bactérienne (UMR7283), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Frédéric Biaso
- Laboratoire de Bioénergétique et Ingénierie des Protéines (UMR7281), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Stéphane Grimaldi
- Laboratoire de Bioénergétique et Ingénierie des Protéines (UMR7281), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Bruno Guigliarelli
- Laboratoire de Bioénergétique et Ingénierie des Protéines (UMR7281), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| | - Axel Magalon
- Laboratoire de Chimie Bactérienne (UMR7283), IMM, IM2B, Aix Marseille Université, CNRS, 13402 Marseille, France
| |
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
|
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
|
15
|
Zhang F, Vik SB. Analysis of the assembly pathway for membrane subunits of Complex I reveals that subunit L (ND5) can assemble last in E. coli. BBA ADVANCES 2021; 1. [DOI: 10.1016/j.bbadva.2021.100027] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
|
16
|
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
|
17
|
Gutiérrez-Fernández J, Kaszuba K, Minhas GS, Baradaran R, Tambalo M, Gallagher DT, Sazanov LA. Key role of quinone in the mechanism of respiratory complex I. Nat Commun 2020; 11:4135. [PMID: 32811817 PMCID: PMC7434922 DOI: 10.1038/s41467-020-17957-0] [Citation(s) in RCA: 62] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 07/28/2020] [Indexed: 01/12/2023] Open
Abstract
Complex I is the first and the largest enzyme of respiratory chains in bacteria and mitochondria. The mechanism which couples spatially separated transfer of electrons to proton translocation in complex I is not known. Here we report five crystal structures of T. thermophilus enzyme in complex with NADH or quinone-like compounds. We also determined cryo-EM structures of major and minor native states of the complex, differing in the position of the peripheral arm. Crystal structures show that binding of quinone-like compounds (but not of NADH) leads to a related global conformational change, accompanied by local re-arrangements propagating from the quinone site to the nearest proton channel. Normal mode and molecular dynamics analyses indicate that these are likely to represent the first steps in the proton translocation mechanism. Our results suggest that quinone binding and chemistry play a key role in the coupling mechanism of complex I. Complex I (NADH:ubiquinone oxidoreductase) is the first enzyme of the respiratory chain in bacteria and mitochondria. Here, the authors present cryo-EM and crystal structures of T. thermophilus complex I in different conformational states and further analyse them by Normal Mode Analysis and molecular dynamics simulations and conclude that quinone redox reactions are important for the coupling mechanism of complex I.
Collapse
Affiliation(s)
| | - Karol Kaszuba
- Institute of Science and Technology Austria, Am Campus 1, A-3400, Klosterneuburg, Austria
| | - Gurdeep S Minhas
- Medical Research Council Mitochondrial Biology Unit, Keith Peters Building, Hills rd, Cambridge, CB2 0XY, UK.,Sosei Heptares, Steinmetz Building, Granta Park, Cambridge, CB21 6DG, UK
| | - Rozbeh Baradaran
- Medical Research Council Mitochondrial Biology Unit, Keith Peters Building, Hills rd, Cambridge, CB2 0XY, UK.,Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, 17165, Solna, Sweden
| | - Margherita Tambalo
- Institute of Science and Technology Austria, Am Campus 1, A-3400, Klosterneuburg, Austria
| | - David T Gallagher
- Medical Research Council Mitochondrial Biology Unit, Keith Peters Building, Hills rd, Cambridge, CB2 0XY, UK
| | - Leonid A Sazanov
- Institute of Science and Technology Austria, Am Campus 1, A-3400, Klosterneuburg, Austria.
| |
Collapse
|
18
|
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
|
19
|
Steiner J, Sazanov L. Structure and mechanism of the Mrp complex, an ancient cation/proton antiporter. eLife 2020; 9:59407. [PMID: 32735215 PMCID: PMC7419157 DOI: 10.7554/elife.59407] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 07/30/2020] [Indexed: 11/19/2022] Open
Abstract
Multiple resistance and pH adaptation (Mrp) antiporters are multi-subunit Na+ (or K+)/H+ exchangers representing an ancestor of many essential redox-driven proton pumps, such as respiratory complex I. The mechanism of coupling between ion or electron transfer and proton translocation in this large protein family is unknown. Here, we present the structure of the Mrp complex from Anoxybacillus flavithermus solved by cryo-EM at 3.0 Å resolution. It is a dimer of seven-subunit protomers with 50 trans-membrane helices each. Surface charge distribution within each monomer is remarkably asymmetric, revealing probable proton and sodium translocation pathways. On the basis of the structure we propose a mechanism where the coupling between sodium and proton translocation is facilitated by a series of electrostatic interactions between a cation and key charged residues. This mechanism is likely to be applicable to the entire family of redox proton pumps, where electron transfer to substrates replaces cation movements.
Collapse
Affiliation(s)
- Julia Steiner
- Institute of Science and Technology Austria, Klosterneuburg, Austria
| | - Leonid Sazanov
- Institute of Science and Technology Austria, Klosterneuburg, Austria
| |
Collapse
|