1
|
Ohnuki J, Okazaki KI. Integration of AlphaFold with Molecular Dynamics for Efficient Conformational Sampling of Transporter Protein NarK. J Phys Chem B 2024. [PMID: 39066727 DOI: 10.1021/acs.jpcb.4c02726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Transporter proteins carry their substrate across the cell membrane by changing their conformation. Thus, conformational dynamics are crucial for transport function. However, clarifying the complete transport cycle is challenging even with the current structural biology approach. Molecular dynamics (MD) simulation is a computational approach that can provide the time-resolved conformational dynamics of transporter proteins in atomic details but suffers from a high computational cost. Here, we integrate state-of-the-art protein structure prediction AI, AlphaFold2 (AF2), with MD simulation to reduce the computational cost. Focusing on the transporter protein NarK, we first show that AF2 sampled broad conformations of NarK, including the inward-open, occluded, and outward-open states. We also applied the coevolution-informed mutation in AF2, identifying state-shifting mutations. Then, we show that MD simulations from AF2-generated outward-open conformation, which is experimentally unresolved, captured the essence of the conformational state. We also found that MD simulations from AF2-generated intermediates showed transient dynamics like a transition state connecting two conformational states. This study paves the way for efficient conformational sampling of transporter proteins.
Collapse
Affiliation(s)
- Jun Ohnuki
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8585, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi 444-8585, Japan
| | - Kei-Ichi Okazaki
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi 444-8585, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
2
|
Pasquadibisceglie A, Quadrotta V, Polticelli F. In Silico Analysis of the Structural Dynamics and Substrate Recognition Determinants of the Human Mitochondrial Carnitine/Acylcarnitine SLC25A20 Transporter. Int J Mol Sci 2023; 24:ijms24043946. [PMID: 36835358 PMCID: PMC9961348 DOI: 10.3390/ijms24043946] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 02/02/2023] [Accepted: 02/09/2023] [Indexed: 02/18/2023] Open
Abstract
The Carnitine-Acylcarnitine Carrier is a member of the mitochondrial Solute Carrier Family 25 (SLC25), known as SLC25A20, involved in the electroneutral exchange of acylcarnitine and carnitine across the inner mitochondrial membrane. It acts as a master regulator of fatty acids β-oxidation and is known to be involved in neonatal pathologies and cancer. The transport mechanism, also known as "alternating access", involves a conformational transition in which the binding site is accessible from one side of the membrane or the other. In this study, through a combination of state-of-the-art modelling techniques, molecular dynamics, and molecular docking, the structural dynamics of SLC25A20 and the early substrates recognition step have been analyzed. The results obtained demonstrated a significant asymmetry in the conformational changes leading to the transition from the c- to the m-state, confirming previous observations on other homologous transporters. Moreover, analysis of the MD simulations' trajectories of the apo-protein in the two conformational states allowed for a better understanding of the role of SLC25A20 Asp231His and Ala281Val pathogenic mutations, which are at the basis of Carnitine-Acylcarnitine Translocase Deficiency. Finally, molecular docking coupled to molecular dynamics simulations lend support to the multi-step substrates recognition and translocation mechanism already hypothesized for the ADP/ATP carrier.
Collapse
Affiliation(s)
| | | | - Fabio Polticelli
- Department of Sciences, University of Roma Tre, 00146 Rome, Italy
- National Institute of Nuclear Physics, Roma Tre Section, 00146 Rome, Italy
- Correspondence:
| |
Collapse
|
3
|
Giangregorio N, Pierri CL, Tonazzi A, Incampo G, Tragni V, De Grassi A, Indiveri C. Proline/Glycine residues of the PG-levels guide conformational changes along the transport cycle in the mitochondrial carnitine/acylcarnitine carrier (SLC25A20). Int J Biol Macromol 2022; 221:1453-1465. [PMID: 36122779 DOI: 10.1016/j.ijbiomac.2022.09.135] [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: 07/19/2022] [Revised: 09/05/2022] [Accepted: 09/15/2022] [Indexed: 11/19/2022]
Abstract
Mitochondrial carnitine/acylcarnitine carrier (CAC) is a member of the mitochondrial carrier (MC) family and imports acylcarnitine into the mitochondrial matrix in exchange for carnitine, playing a pivotal role in carnitine shuttle, crucial for fatty acid oxidation. The crystallized structure of CAC has not been solved yet, however, the availability of several in vitro/in silico studies, also based on the crystallized structures of the ADP/ATP carrier in the cytosolic-conformation and in the matrix-conformation, has made possible to confirm the hypothesis of the single-binding centered-gated pore mechanism for all the members of the MC family. In addition, our recent bioinformatics analyses allowed quantifying in silico the importance of protein residues of MC substrate binding region, of those involved in the formation of the matrix and cytosolic gates, and of those belonging to the Pro/Gly (PG) levels, proposed to be crucial for the tilting/kinking/bending of the six MC transmembrane helices, funneling the substrate translocation pathway. Here we present a combined in silico/in vitro analysis employed for investigating the role played by a group of 6 proline residues and 6 glycine residues, highly conserved in CAC, belonging to MC PG-levels. Residues of the PG-levels surround the similarly located MC common substrate binding region, and were proposed to lead conformational changes and substrate translocation, following substrate binding. For our analysis, we employed 3D molecular modeling approaches, alanine scanning site-directed mutagenesis and in vitro transport assays. Our analysis reveals that P130 (H3), G268 (H6) and G220 (H5), mutated in alanine, affect severely CAC transport activity (mutant catalytic efficiency lower than 5 % compared to the wild type CAC), most likely due to their major role in triggering CAC conformational changes, following carnitine binding. Notably, P30A (H1) and G121A (H3) CAC mutants, increase the carnitine uptake up to 217 % and 112 %, respectively, compared to the wild type CAC.
Collapse
Affiliation(s)
- Nicola Giangregorio
- CNR Institute of Biomembranes, Bioenergetics and Molecular Biotechnologies (IBIOM), Via Amendola 122/O, 70126 Bari, Italy.
| | - Ciro Leonardo Pierri
- Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Italy, Via E. Orabona, 4, 70126 Bari, Italy.
| | - Annamaria Tonazzi
- CNR Institute of Biomembranes, Bioenergetics and Molecular Biotechnologies (IBIOM), Via Amendola 122/O, 70126 Bari, Italy
| | - Giovanna Incampo
- Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Italy, Via E. Orabona, 4, 70126 Bari, Italy
| | - Vincenzo Tragni
- CNR Institute of Biomembranes, Bioenergetics and Molecular Biotechnologies (IBIOM), Via Amendola 122/O, 70126 Bari, Italy; Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Italy, Via E. Orabona, 4, 70126 Bari, Italy
| | - Anna De Grassi
- Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, Italy, Via E. Orabona, 4, 70126 Bari, Italy
| | - Cesare Indiveri
- CNR Institute of Biomembranes, Bioenergetics and Molecular Biotechnologies (IBIOM), Via Amendola 122/O, 70126 Bari, Italy; Department DiBEST (Biologia, Ecologia, Scienze della Terra) Unit of Biochemistry and Molecular Biotechnology, University of Calabria, Via Bucci 4C, 87036 Arcavacata di Rende, Italy
| |
Collapse
|
4
|
Mavridou V, King MS, Tavoulari S, Ruprecht JJ, Palmer SM, Kunji ERS. Substrate binding in the mitochondrial ADP/ATP carrier is a step-wise process guiding the structural changes in the transport cycle. Nat Commun 2022; 13:3585. [PMID: 35739110 PMCID: PMC9226169 DOI: 10.1038/s41467-022-31366-5] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Accepted: 06/14/2022] [Indexed: 02/02/2023] Open
Abstract
Mitochondrial ADP/ATP carriers import ADP into the mitochondrial matrix and export ATP to the cytosol to fuel cellular processes. Structures of the inhibited cytoplasmic- and matrix-open states have confirmed an alternating access transport mechanism, but the molecular details of substrate binding remain unresolved. Here, we evaluate the role of the solvent-exposed residues of the translocation pathway in the process of substrate binding. We identify the main binding site, comprising three positively charged and a set of aliphatic and aromatic residues, which bind ADP and ATP in both states. Additionally, there are two pairs of asparagine/arginine residues on opposite sides of this site that are involved in substrate binding in a state-dependent manner. Thus, the substrates are directed through a series of binding poses, inducing the conformational changes of the carrier that lead to their translocation. The properties of this site explain the electrogenic and reversible nature of adenine nucleotide transport.
Collapse
Affiliation(s)
- Vasiliki Mavridou
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| | - Martin S. King
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| | - Sotiria Tavoulari
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| | - Jonathan J. Ruprecht
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| | - Shane M. Palmer
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| | - Edmund R. S. Kunji
- grid.5335.00000000121885934Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Keith Peters Building, Hills Road, Cambridge, CB2 0XY UK
| |
Collapse
|
5
|
Principal Component Analysis and Related Methods for Investigating the Dynamics of Biological Macromolecules. J 2022. [DOI: 10.3390/j5020021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Principal component analysis (PCA) is used to reduce the dimensionalities of high-dimensional datasets in a variety of research areas. For example, biological macromolecules, such as proteins, exhibit many degrees of freedom, allowing them to adopt intricate structures and exhibit complex functions by undergoing large conformational changes. Therefore, molecular simulations of and experiments on proteins generate a large number of structure variations in high-dimensional space. PCA and many PCA-related methods have been developed to extract key features from such structural data, and these approaches have been widely applied for over 30 years to elucidate macromolecular dynamics. This review mainly focuses on the methodological aspects of PCA and related methods and their applications for investigating protein dynamics.
Collapse
|
6
|
The effects of cardiolipin on the structural dynamics of the mitochondrial ADP/ATP carrier in its cytosol-open state. J Lipid Res 2022; 63:100227. [PMID: 35569528 PMCID: PMC9189224 DOI: 10.1016/j.jlr.2022.100227] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Revised: 05/09/2022] [Accepted: 05/10/2022] [Indexed: 11/21/2022] Open
|
7
|
Pasquadibisceglie A, Polticelli F. Structural determinants of ligands recognition by the human mitochondrial basic amino acids transporter SLC25A29. Insights from molecular dynamics simulations of the c-state. Comput Struct Biotechnol J 2021; 19:5600-5612. [PMID: 34849194 PMCID: PMC8598871 DOI: 10.1016/j.csbj.2021.10.007] [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: 05/28/2021] [Revised: 09/16/2021] [Accepted: 10/04/2021] [Indexed: 11/26/2022] Open
Abstract
In mitochondria, metabolic processes require the trafficking of solutes and organic molecules, such as amino acids. This task is accomplished by the Mitochondrial Carrier Family members (also known as SLC25), among which the SLC25A29 is responsible for the translocation of basic amino acids. In this regard, nitric oxide levels originated by the arginine mitochondrial catabolism have been shown to strongly affect cancer cells' metabolic status. Furthermore, the metabolic disease saccharopinuria has been linked to a mitochondrial dysregulation caused by a toxic intermediate of the lysine catabolism. In both cases, a reduction of the activity of SLC25A29 has been shown to ameliorate these pathological conditions. However, no detailed structural data are available on SLC25A29. In the present work, molecular modelling, docking and dynamics simulations have been employed to analyse the structural determinants of ligands recognition by SLC25A29 in the c-state. Results confirm and reinforce earlier predictions that Asn73, Arg160 and Glu161, and Arg257 represent the ligand contact points I, II, and III, respectively, and that Arg160, Trp204 and Arg257 form a stable interaction, likely critical for ligand binding and translocation. These results are discussed in view of the experimental data available for SLC25A29 and other homologous carriers of the same family.
Collapse
Affiliation(s)
| | - Fabio Polticelli
- Department of Sciences, Roma Tre University, 00146 Rome, Italy.,National Institute of Nuclear Physics, Roma Tre Section, 00146 Rome, Italy
| |
Collapse
|
8
|
Montalvo-Acosta JJ, Kunji ERS, Ruprecht JJ, Dehez F, Chipot C. Structure, substrate binding, and symmetry of the mitochondrial ADP/ATP carrier in its matrix-open state. Biophys J 2021; 120:5187-5195. [PMID: 34748764 DOI: 10.1016/j.bpj.2021.11.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 10/14/2021] [Accepted: 11/02/2021] [Indexed: 01/28/2023] Open
Abstract
The mitochondrial ADP/ATP carrier (AAC) performs the first and last step in oxidative phosphorylation by exchanging ADP and ATP across the mitochondrial inner membrane. Its optimal function has been shown to be dependent on cardiolipins (CLs), unique phospholipids located almost exclusively in the mitochondrial membrane. In addition, AAC exhibits an enthralling threefold pseudosymmetry, a unique feature of members of the SLC25 family. Recently, its conformation poised for binding of ATP was solved by x-ray crystallography referred to as the matrix state. Binding of the substrate leads to conformational changes that export of ATP to the mitochondrial intermembrane space. In this contribution, we investigate the influence of CLs on the structure, substrate-binding properties, and structural symmetry of the matrix state, employing microsecond-scale molecular dynamics simulations. Our findings demonstrate that CLs play a minor stabilizing role on the AAC structure. The interdomain salt bridges and hydrogen bonds forming the cytoplasmic network and tyrosine braces, which ensure the integrity of the global AAC scaffold, highly benefit from the presence of CLs. Under these conditions, the carrier is found to be organized in a more compact structure in its interior, as revealed by analyses of the electrostatic potential, measure of the AAC cavity aperture, and the substrate-binding assays. Introducing a convenient structure-based symmetry metric, we quantified the structural threefold pseudosymmetry of AAC, not only for the crystallographic structure, but also for conformational states of the carrier explored in the molecular dynamics simulations. Our results suggest that CLs moderately contribute to preserve the pseudosymmetric structure of AAC.
Collapse
Affiliation(s)
- Joel José Montalvo-Acosta
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, Vandœuvre-lès-Nancy cedex, France
| | - Edmund R S Kunji
- Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Cambridge, United Kingdom
| | - Jonathan J Ruprecht
- Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Cambridge, United Kingdom
| | - François Dehez
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, Vandœuvre-lès-Nancy cedex, France.
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, Vandœuvre-lès-Nancy cedex, France; Theoretical and Computational Biophysics Group, Beckman Institute for Advanced Science and Technology, Urbana, Illinois; Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois.
| |
Collapse
|
9
|
Pasquadibisceglie A, Polticelli F. Computational studies of the mitochondrial carrier family SLC25. Present status and future perspectives. BIO-ALGORITHMS AND MED-SYSTEMS 2021. [DOI: 10.1515/bams-2021-0018] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Abstract
The members of the mitochondrial carrier family, also known as solute carrier family 25 (SLC25), are transmembrane proteins involved in the translocation of a plethora of small molecules between the mitochondrial intermembrane space and the matrix. These transporters are characterized by three homologous domains structure and a transport mechanism that involves the transition between different conformations. Mutations in regions critical for these transporters’ function often cause several diseases, given the crucial role of these proteins in the mitochondrial homeostasis. Experimental studies can be problematic in the case of membrane proteins, in particular concerning the characterization of the structure–function relationships. For this reason, computational methods are often applied in order to develop new hypotheses or to support/explain experimental evidence. Here the computational analyses carried out on the SLC25 members are reviewed, describing the main techniques used and the outcome in terms of improved knowledge of the transport mechanism. Potential future applications on this protein family of more recent and advanced in silico methods are also suggested.
Collapse
Affiliation(s)
| | - Fabio Polticelli
- Department of Sciences , Roma Tre University , Rome , Italy
- National Institute of Nuclear Physics, Roma Tre Section , Rome , Italy
| |
Collapse
|
10
|
Tanimoto S, Tamura K, Hayashi S, Yoshida N, Nakano H. A computational method to simulate global conformational changes of proteins induced by cosolvent. J Comput Chem 2021; 42:552-563. [PMID: 33433010 DOI: 10.1002/jcc.26481] [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: 09/27/2020] [Revised: 12/09/2020] [Accepted: 12/29/2020] [Indexed: 12/14/2022]
Abstract
A computational method to investigate the global conformational change of a protein is proposed by combining the linear response path following (LRPF) method and three-dimensional reference interaction site model (3D-RISM) theory, which is referred to as the LRPF/3D-RISM method. The proposed method makes it possible to efficiently simulate protein conformational changes caused by either solutions of varying concentrations or the presence of cosolvent species by taking advantage of the LRPF and 3D-RISM. The proposed method is applied to the urea-induced denaturation of ubiquitin. The LRPF/3D-RISM trajectories successfully simulate the early stage of the denaturation process within the simulation time of 300 ns, whereas no significant structural change is observed even in the 1 μs standard MD simulation. The obtained LRPF/3D-RISM trajectories reproduce the mechanism of the urea denaturation of ubiquitin reported in previous studies, and demonstrate the high efficiency of the method.
Collapse
Affiliation(s)
- Shoichi Tanimoto
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Koichi Tamura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shigehiko Hayashi
- Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| | - Haruyuki Nakano
- Department of Chemistry, Graduate School of Science, Kyushu University, Fukuoka, Japan
| |
Collapse
|
11
|
Mao X, Yao S, Yi Q, Xu ZM, Cang X. Function-related asymmetry of the specific cardiolipin binding sites on the mitochondrial ADP/ATP carrier. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2021; 1863:183466. [PMID: 32871114 DOI: 10.1016/j.bbamem.2020.183466] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2020] [Revised: 08/16/2020] [Accepted: 08/24/2020] [Indexed: 01/04/2023]
Abstract
The ADP/ATP carrier (AAC) transports matrix ATP and cytosolic ADP across the inner mitochondrial membrane (IMM). It is well known that cardiolipin (CL) plays an important role in regulating the function of AAC, yet the underlying mechanism still remains elusive. AAC is composed of three homologous domains, and three specific CL binding sites are located at the domain-domain interfaces near the matrix side. Here we report an in-depth investigation on the dynamic properties of the bound CL within the three specific sites through all-atom molecular dynamics simulations of up to 13 μs in total. Our results highlight the importance of the basic and polar residues in CL binding. The basic residues from the linker helix and/or the [Y/W/F][K/R]G motif enable the bound CL to form an intra-domain binding mode, and the canonical inter-domain binding mode only forms when these basic residues are occupied by an additional phospholipid. Of special significance, differences in the basic and polar residues lead to remarkable asymmetry among the three specific CL binding sites. We found that the bound CL at the interface of domains 2 and 3 predominantly adopts inter-domain binding mode, while CLs at the other two sites have much more intra-domain populations. This is consistent with the asymmetric crystal structure of the matrix state (m-state) AAC which implies an asymmetric transport mechanism. The dynamic equilibrium between the inter-domain and intra-domain binding modes observed in our simulations could be highly important for the bound CLs to adapt to the movements during state transitions.
Collapse
Affiliation(s)
- Xiaoting Mao
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310052, China; Institute of Genetics, Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Shihao Yao
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310052, China; Institute of Genetics, Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Qiuzi Yi
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310052, China; Institute of Genetics, Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Zhe-Ming Xu
- Department of Urology, The Children's Hospital, Zhejiang University School of Medicine, National Clinical Research Center for Child Health and National Children's Regional Medical Center, Hangzhou, Zhejiang 310052, China
| | - Xiaohui Cang
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310052, China; Institute of Genetics, Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China.
| |
Collapse
|
12
|
Škulj S, Brkljača Z, Vazdar M. Molecular Dynamics Simulations of the Elusive Matrix‐Open State of Mitochondrial ADP/ATP Carrier. Isr J Chem 2020. [DOI: 10.1002/ijch.202000011] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sanja Škulj
- Division of Organic Chemistry and BiochemistryRuđer Bošković Institute Bijenička 54 HR-10000 Zagreb Croatia
| | - Zlatko Brkljača
- Division of Organic Chemistry and BiochemistryRuđer Bošković Institute Bijenička 54 HR-10000 Zagreb Croatia
| | - Mario Vazdar
- Division of Organic Chemistry and BiochemistryRuđer Bošković Institute Bijenička 54 HR-10000 Zagreb Croatia
| |
Collapse
|
13
|
Tamura K, Sugita Y. Free Energy Analysis of a Conformational Change of Heme ABC Transporter BhuUV-T. J Phys Chem Lett 2020; 11:2824-2829. [PMID: 32202796 PMCID: PMC10961826 DOI: 10.1021/acs.jpclett.0c00547] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The heme ATP-binding cassette (ABC) transporter BhuUV-T of bacterial pathogen Burkholderia cenocepacia is required to transport heme across the inner cell membrane. The current hypothesis is that the binding of two ATPs to the nucleotide-binding domains of the transporter drives the initial steps of the transport cycle in which the empty transport sites are reoriented from the cytosol to the periplasm. Molecular details are missing because the structure of a key occluded intermediate remains hypothetical. Here we perform molecular simulations to analyze the free energy surface (FES) of the first step of the reorientation, namely the transition from an open inward-facing (IF) transport site to an occluded (Occ) conformation. We have modeled the latter structure in silico in a previous study. A simple annealing procedure removes residual bias originating from non-equilibrium targeted molecular dynamics. The calculated FES reveals the role of the ATPs in inducing the IF → Occ conformational change and validates the modeled Occ conformation.
Collapse
Affiliation(s)
- Koichi Tamura
- Computational Biophysics
Research
Team, RIKEN Center for Computational Science, 6-7-1 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Computational Biophysics
Research
Team, RIKEN Center for Computational Science, 6-7-1 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Laboratory for Biomolecular Function
Simulation, RIKEN Center for Biosystems
Dynamics Research, 6-7-1
minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
14
|
Tamura K, Sugimoto H, Shiro Y, Sugita Y. Chemo-Mechanical Coupling in the Transport Cycle of a Heme ABC Transporter. J Phys Chem B 2019; 123:7270-7281. [PMID: 31362510 DOI: 10.1021/acs.jpcb.9b04356] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The heme importer from pathogenic bacteria is a member of the ATP-binding cassette (ABC) transporter family, which uses the energy of ATP-binding and hydrolysis for extensive conformational changes. Previous studies have indicated that conformational changes after heme translocation are triggered by ATP-binding to nucleotide binding domains (NBDs) and then, in turn, induce conformational transitions of the transmembrane domains (TMDs). In this study, we applied a template-based iterative all-atom molecular dynamics (MD) simulation to predict the ATP-bound outward-facing conformation of the Burkholderia cenocepacia heme importer BhuUV-T. The resulting model showed a stable conformation of the TMD with the cytoplasmic gate in the closed state and the periplasmic gate in the open state. Furthermore, targeted MD simulation predicted the intermediate structure of an occluded form (Occ) with bound ATP, in which both ends of the heme translocation channel are closed. The MD simulation of the predicted Occ revealed that Ser147 on the ABC signature motifs (LSGG[Q/E]) of NBDs occasionally flips and loses the active conformation required for ATP-hydrolysis. The flipping motion was found to be coupled to the inter-NBD distance. Our results highlight the functional significance of the signature motif of ABC transporters in regulation of ATPase and chemo-mechanical coupling mechanism.
Collapse
Affiliation(s)
- Koichi Tamura
- Computational Biophysics Research Team , RIKEN Center for Computational Science , 6-7-1 minatojima-Minamimachi, Chuo-ku , Kobe , Hyogo 650-0047 , Japan
| | - Hiroshi Sugimoto
- Graduate School of Life Science , University of Hyogo , 3-2-1 Kouto, Kamigori , Ako , Hyogo 678-1297 , Japan.,Synchrotron Radiation Life Science Instrumentation Team , RIKEN SPring-8 Center , 1-1-1 Kouto , Sayo , Hyogo 679-5148 , Japan
| | - Yoshitsugu Shiro
- Graduate School of Life Science , University of Hyogo , 3-2-1 Kouto, Kamigori , Ako , Hyogo 678-1297 , Japan
| | - Yuji Sugita
- Computational Biophysics Research Team , RIKEN Center for Computational Science , 6-7-1 minatojima-Minamimachi, Chuo-ku , Kobe , Hyogo 650-0047 , Japan.,Theoretical Molecular Science Laboratory , RIKEN Cluster for Pioneering Research , 2-1 Hirosawa , Wako , Saitama 351-0198 , Japan.,Laboratory for Biomolecular Function Simulation , RIKEN Center for Biosystems Dynamics Research , 6-7-1 minatojima-Minamimachi, Chuo-ku , Kobe , Hyogo 650-0047 , Japan
| |
Collapse
|
15
|
Yi Q, Li Q, Yao S, Chen Y, Guan MX, Cang X. Molecular dynamics simulations on apo ADP/ATP carrier shed new lights on the featured motif of the mitochondrial carriers. Mitochondrion 2019; 47:94-102. [PMID: 31129042 DOI: 10.1016/j.mito.2019.05.006] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2018] [Revised: 04/30/2019] [Accepted: 05/20/2019] [Indexed: 02/04/2023]
Abstract
The ADP/ATP carrier (AAC) is a transporter responsible for the equal molar exchange of cytosolic ADP and ATP synthesized within mitochondrial matrix across the mitochondrial membrane. Its primary structure consists of three homologous repeats, and each repeat contains a conserved motif that is shared by all members of the mitochondrial carrier family (MCF). Although these MCF motif residues cluster together in the crystal structure of AAC, detailed analyses on the interactions among the motif residues are still limited. In the present study, all-atom molecular dynamics (MD) simulations of up to 10 μs have been carried out on AAC, and interactions and structural dynamics of the MCF motif residues have been specifically investigated. Our simulations have revealed: i) a very asymmetrical electrostatic network at the bottom of the pocket of apo AAC, ii) the asymmetrical interactions between the Pro kink region and the [YWF][KR] G motif in three repeats, iii) the role of the conserved Arg residues in stabilizing the C-ends of the odd-numbered helices, iv) the structural change of the [YWF][KR] G motif and its potential involvement in substrate translocation process. Our results highlight the asymmetry of the MCF residues in the three repeats, which might contribute to the ability of the carriers to transport the asymmetrical substrates. Our observations provide microscopic basis for further research on the translocation mechanism of mitochondrial carriers.
Collapse
Affiliation(s)
- Qiuzi Yi
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, China; Institute of Genetics, Zhejiang University and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Qiang Li
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, China; Institute of Genetics, Zhejiang University and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Shihao Yao
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, China; Institute of Genetics, Zhejiang University and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Ye Chen
- School of Information and Electric Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221008, China
| | - Min-Xin Guan
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, China; Institute of Genetics, Zhejiang University and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China
| | - Xiaohui Cang
- Division of Medical Genetics and Genomics, The Children's Hospital, Zhejiang University School of Medicine, China; Institute of Genetics, Zhejiang University and Department of Genetics, Zhejiang University School of Medicine, Hangzhou, Zhejiang 310058, China.
| |
Collapse
|
16
|
Chipot C, Dehez F, Schnell JR, Zitzmann N, Pebay-Peyroula E, Catoire LJ, Miroux B, Kunji ERS, Veglia G, Cross TA, Schanda P. Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies. Chem Rev 2018; 118:3559-3607. [PMID: 29488756 PMCID: PMC5896743 DOI: 10.1021/acs.chemrev.7b00570] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Indexed: 12/25/2022]
Abstract
Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.
Collapse
Affiliation(s)
- Christophe Chipot
- SRSMC, UMR 7019 Université de Lorraine CNRS, Vandoeuvre-les-Nancy F-54500, France
- Laboratoire
International Associé CNRS and University of Illinois at Urbana−Champaign, Vandoeuvre-les-Nancy F-54506, France
- Department
of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - François Dehez
- SRSMC, UMR 7019 Université de Lorraine CNRS, Vandoeuvre-les-Nancy F-54500, France
- Laboratoire
International Associé CNRS and University of Illinois at Urbana−Champaign, Vandoeuvre-les-Nancy F-54506, France
| | - Jason R. Schnell
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Nicole Zitzmann
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | | | - Laurent J. Catoire
- Laboratory
of Biology and Physico-Chemistry of Membrane Proteins, Institut de Biologie Physico-Chimique (IBPC), UMR
7099 CNRS, Paris 75005, France
- University
Paris Diderot, Paris 75005, France
- PSL
Research University, Paris 75005, France
| | - Bruno Miroux
- Laboratory
of Biology and Physico-Chemistry of Membrane Proteins, Institut de Biologie Physico-Chimique (IBPC), UMR
7099 CNRS, Paris 75005, France
- University
Paris Diderot, Paris 75005, France
- PSL
Research University, Paris 75005, France
| | - Edmund R. S. Kunji
- Medical
Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge CB2 0XY, United Kingdom
| | - Gianluigi Veglia
- Department
of Biochemistry, Molecular Biology, and Biophysics, and Department
of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Timothy A. Cross
- National
High Magnetic Field Laboratory, Florida
State University, Tallahassee, Florida 32310, United States
| | - Paul Schanda
- Université
Grenoble Alpes, CEA, CNRS, IBS, Grenoble F-38000, France
| |
Collapse
|