1
|
Kamsri B, Kamsri P, Punkvang A, Chimprasit A, Saparpakorn P, Hannongbua S, Spencer J, Oliveira ASF, Mulholland AJ, Pungpo P. Signal Propagation in the ATPase Domain of Mycobacterium tuberculosis DNA Gyrase from Dynamical-Nonequilibrium Molecular Dynamics Simulations. Biochemistry 2024; 63:1493-1504. [PMID: 38742407 PMCID: PMC11154950 DOI: 10.1021/acs.biochem.4c00161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Revised: 04/26/2024] [Accepted: 05/01/2024] [Indexed: 05/16/2024]
Abstract
DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A2B2) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.
Collapse
Affiliation(s)
- Bundit Kamsri
- Department
of Chemistry and Center of Excellence for Innovation in Chemistry,
Faculty of Science, Ubon Ratchathani University, Ubon Ratchathani 34190, Thailand
| | - Pharit Kamsri
- Division
of Chemistry, Faculty of Science, Nakhon
Phanom University, Nakhon
Phanom 48000, Thailand
| | - Auradee Punkvang
- Division
of Chemistry, Faculty of Science, Nakhon
Phanom University, Nakhon
Phanom 48000, Thailand
| | - Aunlika Chimprasit
- Department
of Chemistry, Faculty of Science, Kasetsart
University, Bangkok 10900, Thailand
| | | | - Supa Hannongbua
- Department
of Chemistry, Faculty of Science, Kasetsart
University, Bangkok 10900, Thailand
| | - James Spencer
- School
of Cellular and Molecular Medicine, Biomedical Sciences Building, University of Bristol, Bristol BS8 1TD, U.K.
| | - A. Sofia F. Oliveira
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Adrian J. Mulholland
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Pornpan Pungpo
- Department
of Chemistry and Center of Excellence for Innovation in Chemistry,
Faculty of Science, Ubon Ratchathani University, Ubon Ratchathani 34190, Thailand
| |
Collapse
|
2
|
Castelli M, Marchetti F, Osuna S, F. Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations. J Am Chem Soc 2024; 146:901-919. [PMID: 38116743 PMCID: PMC10785808 DOI: 10.1021/jacs.3c11396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/06/2023] [Accepted: 12/07/2023] [Indexed: 12/21/2023]
Abstract
Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.
Collapse
Affiliation(s)
- Matteo Castelli
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Filippo Marchetti
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
- INSTM, via G. Giusti 9, 50121 Florence, Italy
- E4
Computer Engineering, via Martiri delle libertà 66, 42019 Scandiano (RE), Italy
| | - Sílvia Osuna
- Institut
de Química Computacional i Catàlisi (IQCC) and Departament
de Química, Universitat de Girona, Girona, Catalonia E-17071, Spain
- ICREA, Barcelona, Catalonia E-08010, Spain
| | - A. Sofia F. Oliveira
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Adrian J. Mulholland
- Centre for
Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K.
| | - Stefano A. Serapian
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| | - Giorgio Colombo
- Department
of Chemistry, University of Pavia, viale T. Taramelli 12, 27100 Pavia, Italy
| |
Collapse
|
3
|
Soya N, Xu H, Roldan A, Yang Z, Ye H, Jiang F, Premchandar A, Veit G, Cole SPC, Kappes J, Hegedüs T, Lukacs GL. Folding correctors can restore CFTR posttranslational folding landscape by allosteric domain-domain coupling. Nat Commun 2023; 14:6868. [PMID: 37891162 PMCID: PMC10611759 DOI: 10.1038/s41467-023-42586-8] [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: 11/16/2022] [Accepted: 10/16/2023] [Indexed: 10/29/2023] Open
Abstract
The folding/misfolding and pharmacological rescue of multidomain ATP-binding cassette (ABC) C-subfamily transporters, essential for organismal health, remain incompletely understood. The ABCC transporters core consists of two nucleotide binding domains (NBD1,2) and transmembrane domains (TMD1,2). Using molecular dynamic simulations, biochemical and hydrogen deuterium exchange approaches, we show that the mutational uncoupling or stabilization of NBD1-TMD1/2 interfaces can compromise or facilitate the CFTR(ABCC7)-, MRP1(ABCC1)-, and ABCC6-transporters posttranslational coupled domain-folding in the endoplasmic reticulum. Allosteric or orthosteric binding of VX-809 and/or VX-445 folding correctors to TMD1/2 can rescue kinetically trapped CFTR posttranslational folding intermediates of cystic fibrosis (CF) mutants of NBD1 or TMD1 by global rewiring inter-domain allosteric-networks. We propose that dynamic allosteric domain-domain communications not only regulate ABCC-transporters function but are indispensable to tune the folding landscape of their posttranslational intermediates. These allosteric networks can be compromised by CF-mutations, and reinstated by correctors, offering a framework for mechanistic understanding of ABCC-transporters (mis)folding.
Collapse
Affiliation(s)
- Naoto Soya
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Haijin Xu
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Ariel Roldan
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Zhengrong Yang
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Haoxin Ye
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Fan Jiang
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Aiswarya Premchandar
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Guido Veit
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada
| | - Susan P C Cole
- Division of Cancer Biology and Genetics, Department of Pathology and Molecular Medicine, Queen's University Cancer Research Institute, Kingston, ON, Canada
| | - John Kappes
- Heersink School of Medicine, University of Alabama School of Medicine, Birmingham, AL, USA
| | - Tamás Hegedüs
- Department of Biophysics and Radiation Biology, Semmelweis University, 1085, Budapest, Hungary
- ELKH-SE Biophysical Virology Research Group, Eötvös Loránd Research Network, Budapest, Hungary
| | - Gergely L Lukacs
- Department of Physiology and Biochemistry, McGill University, Montréal, QC, Canada.
| |
Collapse
|
4
|
Soya N, Xu H, Roldan A, Yang Z, Ye H, Jiang F, Premchandar A, Veit G, Cole SPC, Kappes J, Hegedus T, Lukacs GL. Folding correctors can restore CFTR posttranslational folding landscape by allosteric domain-domain coupling. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.19.563107. [PMID: 37905074 PMCID: PMC10614980 DOI: 10.1101/2023.10.19.563107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The folding/misfolding and pharmacological rescue of multidomain ATP-binding cassette (ABC) C-subfamily transporters, essential for organismal health, remain incompletely understood. The ABCC transporters core consists of two nucleotide binding domains (NBD1,2) and transmembrane domains (TMD1,2). Using molecular dynamic simulations, biochemical and hydrogen deuterium exchange approaches, we show that the mutational uncoupling or stabilization of NBD1-TMD1/2 interfaces can compromise or facilitate the CFTR(ABCC7)-, MRP1(ABCC1)-, and ABCC6-transporters posttranslational coupled domain-folding in the endoplasmic reticulum. Allosteric or orthosteric binding of VX-809 and/or VX-445 folding correctors to TMD1/2 can rescue kinetically trapped CFTR post-translational folding intermediates of cystic fibrosis (CF) mutants of NBD1 or TMD1 by global rewiring inter-domain allosteric-networks. We propose that dynamic allosteric domain-domain communications not only regulate ABCC-transporters function but are indispensable to tune the folding landscape of their post-translational intermediates. These allosteric networks can be compromised by CF-mutations, and reinstated by correctors, offering a framework for mechanistic understanding of ABCC-transporters (mis)folding. One-Sentence Summary Allosteric interdomain communication and its modulation are critical determinants of ABCC-transporters post-translational conformational biogenesis, misfolding, and pharmacological rescue.
Collapse
|
5
|
Oliveira ASF, Shoemark DK, Davidson AD, Berger I, Schaffitzel C, Mulholland AJ. SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid. J Mol Cell Biol 2023; 15:mjad021. [PMID: 36990513 PMCID: PMC10563148 DOI: 10.1093/jmcb/mjad021] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 11/07/2022] [Accepted: 03/28/2023] [Indexed: 03/31/2023] Open
Abstract
The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.
Collapse
Affiliation(s)
- A Sofia F Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | | | - Andrew D Davidson
- School of Cellular and Molecular Medicine, University of Bristol, University Walk, Bristol BS8 1TD, UK
| | - Imre Berger
- School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
- School of Chemistry, Max Planck Bristol Centre for Minimal Biology, Bristol BS8 1TS, UK
| | | | - Adrian J Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
6
|
Chan HT, Oliveira ASF, Schofield CJ, Mulholland AJ, Duarte F. Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease. JACS AU 2023; 3:1767-1774. [PMID: 37384148 PMCID: PMC10262681 DOI: 10.1021/jacsau.3c00185] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 05/19/2023] [Accepted: 05/19/2023] [Indexed: 06/30/2023]
Abstract
The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.
Collapse
Affiliation(s)
- H. T.
Henry Chan
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - A. Sofia F. Oliveira
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
- School
of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Christopher J. Schofield
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - Adrian J. Mulholland
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
| | - Fernanda Duarte
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| |
Collapse
|
7
|
Agajanian S, Alshahrani M, Bai F, Tao P, Verkhivker GM. Exploring and Learning the Universe of Protein Allostery Using Artificial Intelligence Augmented Biophysical and Computational Approaches. J Chem Inf Model 2023; 63:1413-1428. [PMID: 36827465 PMCID: PMC11162550 DOI: 10.1021/acs.jcim.2c01634] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/26/2023]
Abstract
Allosteric mechanisms are commonly employed regulatory tools used by proteins to orchestrate complex biochemical processes and control communications in cells. The quantitative understanding and characterization of allosteric molecular events are among major challenges in modern biology and require integration of innovative computational experimental approaches to obtain atomistic-level knowledge of the allosteric states, interactions, and dynamic conformational landscapes. The growing body of computational and experimental studies empowered by emerging artificial intelligence (AI) technologies has opened up new paradigms for exploring and learning the universe of protein allostery from first principles. In this review we analyze recent developments in high-throughput deep mutational scanning of allosteric protein functions; applications and latest adaptations of Alpha-fold structural prediction methods for studies of protein dynamics and allostery; new frontiers in integrating machine learning and enhanced sampling techniques for characterization of allostery; and recent advances in structural biology approaches for studies of allosteric systems. We also highlight recent computational and experimental studies of the SARS-CoV-2 spike (S) proteins revealing an important and often hidden role of allosteric regulation driving functional conformational changes, binding interactions with the host receptor, and mutational escape mechanisms of S proteins which are critical for viral infection. We conclude with a summary and outlook of future directions suggesting that AI-augmented biophysical and computer simulation approaches are beginning to transform studies of protein allostery toward systematic characterization of allosteric landscapes, hidden allosteric states, and mechanisms which may bring about a new revolution in molecular biology and drug discovery.
Collapse
Affiliation(s)
- Steve Agajanian
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Mohammed Alshahrani
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
| | - Fang Bai
- Shanghai Institute for Advanced Immunochemical Studies, School of Life Science and Technology and Information Science and Technology, Shanghai Tech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Peng Tao
- Department of Chemistry, Center for Research Computing, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, Texas 75205, United States
| | - Gennady M Verkhivker
- Keck Center for Science and Engineering, Graduate Program in Computational and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| |
Collapse
|
8
|
McDonald EF, Woods H, Smith ST, Kim M, Schoeder CT, Plate L, Meiler J. Structural Comparative Modeling of Multi-Domain F508del CFTR. Biomolecules 2022; 12:biom12030471. [PMID: 35327663 PMCID: PMC8946492 DOI: 10.3390/biom12030471] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 03/15/2022] [Accepted: 03/17/2022] [Indexed: 12/07/2022] Open
Abstract
Cystic fibrosis (CF) is a rare genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), an epithelial anion channel expressed in several vital organs. Absence of functional CFTR results in imbalanced osmotic equilibrium and subsequent mucus build up in the lungs-which increases the risk of infection and eventually causes death. CFTR is an ATP-binding cassette (ABC) transporter family protein composed of two transmembrane domains (TMDs), two nucleotide binding domains (NBDs), and an unstructured regulatory domain. The most prevalent patient mutation is the deletion of F508 (F508del), making F508del CFTR the primary target for current FDA approved CF therapies. However, no experimental multi-domain F508del CFTR structure has been determined and few studies have modeled F508del using multi-domain WT CFTR structures. Here, we used cryo-EM density data and Rosetta comparative modeling (RosettaCM) to compare a F508del model with published experimental data on CFTR NBD1 thermodynamics. We then apply this modeling method to generate multi-domain WT and F508del CFTR structural models. These models demonstrate the destabilizing effects of F508del on NBD1 and the NBD1/TMD interface in both the inactive and active conformation of CFTR. Furthermore, we modeled F508del/R1070W and F508del bound to the CFTR corrector VX-809. Our models reveal the stabilizing effects of VX-809 on multi-domain models of F508del CFTR and pave the way for rational design of additional drugs that target F508del CFTR for treatment of CF.
Collapse
Affiliation(s)
- Eli Fritz McDonald
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA; (E.F.M.); (C.T.S.); (L.P.)
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA; (H.W.); (S.T.S.)
| | - Hope Woods
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA; (H.W.); (S.T.S.)
- Program in Chemical and Physical Biology, Vanderbilt University, Nashville, TN 37235, USA;
| | - Shannon T. Smith
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA; (H.W.); (S.T.S.)
- Program in Chemical and Physical Biology, Vanderbilt University, Nashville, TN 37235, USA;
| | - Minsoo Kim
- Program in Chemical and Physical Biology, Vanderbilt University, Nashville, TN 37235, USA;
| | - Clara T. Schoeder
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA; (E.F.M.); (C.T.S.); (L.P.)
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA; (H.W.); (S.T.S.)
- Leipzig Medical School, Leipzig University, 04109 Leipzig, Germany
| | - Lars Plate
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA; (E.F.M.); (C.T.S.); (L.P.)
- Department of Biological Sciences, Vanderbilt University, Nashville, TN 37235, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA; (E.F.M.); (C.T.S.); (L.P.)
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37235, USA; (H.W.); (S.T.S.)
- Leipzig Medical School, Leipzig University, 04109 Leipzig, Germany
- Department of Pharmacology, Vanderbilt University, Nashville, TN 37235, USA
- Institute for Drug Discovery, Leipzig University, 04109 Leipzig, Germany
- Correspondence: ; Tel.: +1-(615)-936-2211
| |
Collapse
|
9
|
Prins S, Corradi V, Sheppard DN, Tieleman DP, Vergani P. Can two wrongs make a right? F508del-CFTR ion channel rescue by second-site mutations in its transmembrane domains. J Biol Chem 2022; 298:101615. [PMID: 35065958 PMCID: PMC8861112 DOI: 10.1016/j.jbc.2022.101615] [Citation(s) in RCA: 4] [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: 09/01/2021] [Revised: 01/13/2022] [Accepted: 01/17/2022] [Indexed: 11/20/2022] Open
Abstract
Deletion of phenylalanine 508 (F508del) in the cystic fibrosis transmembrane conductance regulator (CFTR) anion channel is the most common cause of cystic fibrosis. The F508 residue is located on nucleotide-binding domain 1 (NBD1) in contact with the cytosolic extensions of the transmembrane helices, in particular intracellular loop 4 (ICL4). To investigate how absence of F508 at this interface impacts the CFTR protein, we carried out a mutagenesis scan of ICL4 by introducing second-site mutations at 11 positions in cis with F508del. Using an image-based fluorescence assay, we measured how each mutation affected membrane proximity and ion-channel function. The scan strongly validated the effectiveness of R1070W at rescuing F508del defects. Molecular dynamics simulations highlighted two features characterizing the ICL4/NBD1 interface of F508del/R1070W-CFTR: flexibility, with frequent transient formation of interdomain hydrogen bonds, and loosely stacked aromatic sidechains (F1068, R1070W, and F1074, mimicking F1068, F508, and F1074 in WT CFTR). F508del-CFTR displayed a distorted aromatic stack, with F1068 displaced toward the space vacated by F508, while in F508del/R1070F-CFTR, which largely retained F508del defects, R1070F could not form hydrogen bonds and the interface was less flexible. Other ICL4 second-site mutations which partially rescued F508del-CFTR included F1068M and F1074M. Methionine side chains allow hydrophobic interactions without the steric rigidity of aromatic rings, possibly conferring flexibility to accommodate the absence of F508 and retain a dynamic interface. These studies highlight how both hydrophobic interactions and conformational flexibility might be important at the ICL4/NBD1 interface, suggesting possible structural underpinnings of F508del-induced dysfunction.
Collapse
Affiliation(s)
- Stella Prins
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, UK
| | - Valentina Corradi
- Department of Biological Sciences, Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - David N Sheppard
- School of Physiology, Pharmacology and Neuroscience, University of Bristol, Bristol, UK
| | - D Peter Tieleman
- Department of Biological Sciences, Centre for Molecular Simulation, University of Calgary, Calgary, Alberta, Canada
| | - Paola Vergani
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, UK.
| |
Collapse
|
10
|
Sofia F. Oliveira A, Shoemark DK, Avila Ibarra A, Davidson AD, Berger I, Schaffitzel C, Mulholland AJ. The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour. Comput Struct Biotechnol J 2021; 20:139-147. [PMID: 34934478 PMCID: PMC8670790 DOI: 10.1016/j.csbj.2021.12.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 12/08/2021] [Accepted: 12/08/2021] [Indexed: 12/12/2022] Open
Abstract
The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.
Collapse
Key Words
- ACE2, angiotensin-converting 2 enzyme
- CD, connector domain
- CH, central helix
- FA, fatty acid
- FP, fusion peptide
- FPPR, fusion-peptide proximal region
- HR1, heptad repeat 1
- LA, Linoleic acid
- MD, Molecular dynamics
- MERS, middle east respiratory syndrome
- NTD, N-terminal domain
- RBD, receptor binding domain
- RBM, receptor-binding motif
- RMB, receptor binding motif
- SARS, severe acute respiratory syndrome
- SARS-CoV-2, severe acute respiratory syndrome 2
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
| | - Deborah K. Shoemark
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Amaurys Avila Ibarra
- Research Software Engineering, Advanced Computing Research Centre, University of Bristol, Bristol BS1 5QD, UK
| | - Andrew D. Davidson
- School of Cellular and Molecular Medicine, University of Bristol, University Walk, Bristol BS8 1TD, UK
| | - Imre Berger
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
- Max Planck Bristol Centre for Minimal Biology, Cantock's Close, Bristol BS8 1TS, UK
| | - Christiane Schaffitzel
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol BS8 1TQ, UK
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Adrian J. Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, UK
| |
Collapse
|
11
|
Abreu B, Cruz C, Oliveira ASF, Soares CM. ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK 2E importer. Sci Rep 2021; 11:10591. [PMID: 34012037 PMCID: PMC8134467 DOI: 10.1038/s41598-021-89556-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2020] [Accepted: 04/27/2021] [Indexed: 02/03/2023] Open
Abstract
ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK2E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK2E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.
Collapse
Affiliation(s)
- Bárbara Abreu
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - Carlos Cruz
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| | - A. Sofia F. Oliveira
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal ,grid.5337.20000 0004 1936 7603School of Biochemistry and Centre for Computational Chemistry, University of Bristol, Bristol, UK
| | - Cláudio M. Soares
- grid.10772.330000000121511713ITQB NOVA, Instituto de Tecnologia Química E Biológica António Xavier, Universidade Nova de Lisboa, Oeiras, Portugal
| |
Collapse
|
12
|
Galdadas I, Qu S, Oliveira ASF, Olehnovics E, Mack AR, Mojica MF, Agarwal PK, Tooke CL, Gervasio FL, Spencer J, Bonomo RA, Mulholland AJ, Haider S. Allosteric communication in class A β-lactamases occurs via cooperative coupling of loop dynamics. eLife 2021; 10:e66567. [PMID: 33755013 PMCID: PMC8060031 DOI: 10.7554/elife.66567] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 03/19/2021] [Indexed: 12/16/2022] Open
Abstract
Understanding allostery in enzymes and tools to identify it offer promising alternative strategies to inhibitor development. Through a combination of equilibrium and nonequilibrium molecular dynamics simulations, we identify allosteric effects and communication pathways in two prototypical class A β-lactamases, TEM-1 and KPC-2, which are important determinants of antibiotic resistance. The nonequilibrium simulations reveal pathways of communication operating over distances of 30 Å or more. Propagation of the signal occurs through cooperative coupling of loop dynamics. Notably, 50% or more of clinically relevant amino acid substitutions map onto the identified signal transduction pathways. This suggests that clinically important variation may affect, or be driven by, differences in allosteric behavior, providing a mechanism by which amino acid substitutions may affect the relationship between spectrum of activity, catalytic turnover, and potential allosteric behavior in this clinically important enzyme family. Simulations of the type presented here will help in identifying and analyzing such differences.
Collapse
Affiliation(s)
- Ioannis Galdadas
- University College London, Department of ChemistryLondonUnited Kingdom
| | - Shen Qu
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| | - Ana Sofia F Oliveira
- University of Bristol, Centre for Computational Chemistry, School of ChemistryBristolUnited Kingdom
| | - Edgar Olehnovics
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| | - Andrew R Mack
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Molecular Biology and MicrobiologyClevelandUnited States
| | - Maria F Mojica
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Infectious Diseases, School of MedicineClevelandUnited States
| | - Pratul K Agarwal
- Department of Physiological Sciences and High-Performance Computing Center, Oklahoma State UniversityStillwaterUnited States
| | - Catherine L Tooke
- University of Bristol, School of Cellular and Molecular MedicineBristolUnited Kingdom
| | - Francesco Luigi Gervasio
- University College London, Department of ChemistryLondonUnited Kingdom
- University College London, Institute of Structural and Molecular BiologyLondonUnited Kingdom
- University of Geneva, Pharmaceutical SciencesGenevaSwitzerland
| | - James Spencer
- University of Bristol, School of Cellular and Molecular MedicineBristolUnited Kingdom
| | - Robert A Bonomo
- Veterans Affairs Northeast Ohio Healthcare System, Research ServiceClevelandUnited States
- Case Western Reserve University, Department of Molecular Biology and MicrobiologyClevelandUnited States
- Case Western Reserve University, Department of Infectious Diseases, School of MedicineClevelandUnited States
- Case Western Reserve University, Department of BiochemistryClevelandUnited States
- Case Western Reserve University, Department of PharmacologyClevelandUnited States
- Case Western Reserve University, Department of Proteomics and BioinformaticsClevelandUnited States
- CWRU-Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES)ClevelandUnited States
| | - Adrian J Mulholland
- University of Bristol, Centre for Computational Chemistry, School of ChemistryBristolUnited Kingdom
| | - Shozeb Haider
- University College London School of Pharmacy, Pharmaceutical and Biological ChemistryLondonUnited Kingdom
| |
Collapse
|
13
|
Oliveira ASF, Ciccotti G, Haider S, Mulholland AJ. Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:144. [PMID: 34720710 PMCID: PMC8549953 DOI: 10.1140/epjb/s10051-021-00157-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/05/2021] [Indexed: 05/03/2023]
Abstract
ABSTRACT A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation. Besides providing a detailed picture of the system's dynamic structural response to an external perturbation, this approach also has the advantage that the statistical significance of the response can be assessed. The D-NEMD approach has been used recently to identify a general mechanism of inter-domain signal propagation in nicotinic acetylcholine receptors, and allosteric effects in β -lactamase enzymes, for example. It complements equilibrium MD and is a very promising approach to identifying and analysing allosteric effects. Here, we review the D-NEMD approach and its application to biomolecular systems, including transporters, receptors, and enzymes.
Collapse
Affiliation(s)
- A. Sofia F. Oliveira
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
- BrisSynBio, Life Sciences Building, Tyndall Avenue, Bristol, BS8 1TQ UK
| | - Giovanni Ciccotti
- Institute for Applied Computing “Mauro Picone” (IAC), CNR, Via dei Taurini 19, 00185 Rome, Italy
- School of Physics, University College of Dublin, UCD-Belfield, Dublin 4, Ireland
- Università di Roma La Sapienza, Ple. A. Moro 5, 00185 Rome, Italy
| | - Shozeb Haider
- School of Pharmacy, University College London, London, WC1N 1AX UK
| | - Adrian J. Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, BS8 1TS UK
| |
Collapse
|