1
|
Pant S, Dehghani-Ghahnaviyeh S, Trebesch N, Rasouli A, Chen T, Kapoor K, Wen PC, Tajkhorshid E. Dissecting Large-Scale Structural Transitions in Membrane Transporters Using Advanced Simulation Technologies. J Phys Chem B 2025; 129:3703-3719. [PMID: 40100959 DOI: 10.1021/acs.jpcb.5c00104] [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/20/2025]
Abstract
Membrane transporters are integral membrane proteins that act as gatekeepers of the cell, controlling fundamental processes such as recruitment of nutrients and expulsion of waste material. At a basic level, transporters operate using the "alternating access model," in which transported substances are accessible from only one side of the membrane at a time. This model usually involves large-scale structural changes in the transporter, which often cannot be captured using unbiased, conventional molecular simulation techniques. In this article, we provide an overview of some of the major simulation techniques that have been applied to characterize the structural dynamics and energetics involved in the transition of membrane transporters between their functional states. After briefly introducing each technique, we discuss some of their advantages and limitations and provide some recent examples of their application to membrane transporters.
Collapse
Affiliation(s)
- Shashank Pant
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Sepehr Dehghani-Ghahnaviyeh
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Noah Trebesch
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Ali Rasouli
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Tianle Chen
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Karan Kapoor
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Po-Chao Wen
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, Department of Biochemistry, and Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois 61801-3028, United States
| |
Collapse
|
2
|
Chan MC, Alfawaz Y, Paul A, Shukla D. Molecular insights into the elevator-type mechanism of the cyanobacterial bicarbonate transporter BicA. Biophys J 2025; 124:379-392. [PMID: 39674889 PMCID: PMC11788499 DOI: 10.1016/j.bpj.2024.12.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2022] [Revised: 03/17/2024] [Accepted: 12/10/2024] [Indexed: 12/17/2024] Open
Abstract
Cyanobacteria are responsible for up to 80% of aquatic carbon dioxide fixation and have evolved a specialized carbon concentrating mechanism to increase photosynthetic yield. As such, cyanobacteria are attractive targets for synthetic biology and engineering approaches to address the demands of global energy security, food production, and climate change for an increasing world's population. The bicarbonate transporter BicA is a sodium-dependent, low-affinity, high-flux bicarbonate symporter expressed in the plasma membrane of cyanobacteria. Despite extensive biochemical characterization of BicA, including the resolution of the BicA crystal structure, the dynamic understanding of the bicarbonate transport mechanism remains elusive. To this end, we have collected over 1 ms of all-atom molecular dynamics simulation data of the BicA dimer to elucidate the structural rearrangements involved in the substrate transport process. We further characterized the energetics of the transition of BicA protomers and investigated potential mutations that are shown to decrease the free energy barrier of conformational transitions. In all, our study illuminates a detailed mechanistic understanding of the conformational dynamics of bicarbonate transporters and provides atomistic insights to engineering these transporters for enhanced photosynthetic production.
Collapse
Affiliation(s)
- Matthew C Chan
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois
| | - Yazeed Alfawaz
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois
| | - Arnav Paul
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, Illinois
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, Illinois; Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, Illinois; Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, Illinois; Department of Plant Biology, University of Illinois Urbana-Champaign, Urbana, Illinois; Department of Bioengineering, University of Illinois Urbana-Champaign, Urbana, Illinois.
| |
Collapse
|
3
|
Liu B, Boysen JG, Unarta IC, Du X, Li Y, Huang X. Exploring transition states of protein conformational changes via out-of-distribution detection in the hyperspherical latent space. Nat Commun 2025; 16:349. [PMID: 39753544 PMCID: PMC11699157 DOI: 10.1038/s41467-024-55228-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2024] [Accepted: 12/05/2024] [Indexed: 01/06/2025] Open
Abstract
Identifying transitional states is crucial for understanding protein conformational changes that underlie numerous biological processes. Markov state models (MSMs), built from Molecular Dynamics (MD) simulations, capture these dynamics through transitions among metastable conformational states, and have demonstrated success in studying protein conformational changes. However, MSMs face challenges in identifying transition states, as they partition MD conformations into discrete metastable states (or free energy minima), lacking description of transition states located at the free energy barriers. Here, we introduce Transition State identification via Dispersion and vAriational principle Regularized neural networks (TS-DAR), a deep learning framework inspired by out-of-distribution (OOD) detection in trustworthy artificial intelligence (AI). TS-DAR offers an end-to-end pipeline that can simultaneously detect all transition states between multiple free minima from MD simulations using the regularized hyperspherical embeddings in latent space. The key insight of TS-DAR lies in treating transition state structures as OOD data, recognizing that they are sparsely populated and exhibit a distributional shift from metastable states. We demonstrate the power of TS-DAR by applying it to a 2D potential, alanine dipeptide, and the translocation of a DNA motor protein on DNA, where it outperforms previous methods in identifying transition states.
Collapse
Affiliation(s)
- Bojun Liu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Data Science Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Jordan G Boysen
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Ilona Christy Unarta
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Data Science Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuefeng Du
- Department of Computer Sciences, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Yixuan Li
- Data Science Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Department of Computer Sciences, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA.
- Data Science Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA.
| |
Collapse
|
4
|
Dmitrieva N, Gholami S, Alleva C, Carloni P, Alfonso-Prieto M, Fahlke C. Transport mechanism of DgoT, a bacterial homolog of SLC17 organic anion transporters. EMBO J 2024; 43:6740-6765. [PMID: 39455803 PMCID: PMC11649914 DOI: 10.1038/s44318-024-00279-y] [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: 03/19/2024] [Revised: 10/01/2024] [Accepted: 10/07/2024] [Indexed: 10/28/2024] Open
Abstract
The solute carrier 17 (SLC17) family contains anion transporters that accumulate neurotransmitters in secretory vesicles, remove carboxylated monosaccharides from lysosomes, or extrude organic anions from the kidneys and liver. We combined classical molecular dynamics simulations, Markov state modeling and hybrid first principles quantum mechanical/classical mechanical (QM/MM) simulations with experimental approaches to describe the transport mechanisms of a model bacterial protein, the D-galactonate transporter DgoT, at atomic resolution. We found that protonation of D46 and E133 precedes galactonate binding and that substrate binding induces closure of the extracellular gate, with the conserved R47 coupling substrate binding to transmembrane helix movement. After isomerization to an inward-facing conformation, deprotonation of E133 and subsequent proton transfer from D46 to E133 opens the intracellular gate and permits galactonate dissociation either in its unprotonated form or after proton transfer from E133. After release of the second proton, apo DgoT returns to the outward-facing conformation. Our results provide a framework to understand how various SLC17 transport functions with distinct transport stoichiometries can be attained through subtle variations in proton and substrate binding/unbinding.
Collapse
Affiliation(s)
- Natalia Dmitrieva
- Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, Forschungszentrum Jülich, 52425, Jülich, Germany
| | - Samira Gholami
- Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, Forschungszentrum Jülich, 52425, Jülich, Germany
- Institute for Advanced Simulation (IAS-5) and Institute of Neuroscience and Medicine (INM-9), Computational Biomedicine, Forschungszentrum Jülich, 52425, Jülich, Germany
| | - Claudia Alleva
- Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, Forschungszentrum Jülich, 52425, Jülich, Germany
- Department of Biochemistry and Biophysics and Science for Life Laboratory, Stockholm University, Stockholm, Sweden
| | - Paolo Carloni
- Institute for Advanced Simulation (IAS-5) and Institute of Neuroscience and Medicine (INM-9), Computational Biomedicine, Forschungszentrum Jülich, 52425, Jülich, Germany
- JARA-HPC, Forschungszentrum Jülich, 54245, Jülich, Germany
- Department of Physics, RWTH Aachen University, 52056, Aachen, Germany
- JARA Institute Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich, 52425, Jülich, Germany
| | - Mercedes Alfonso-Prieto
- Institute for Advanced Simulation (IAS-5) and Institute of Neuroscience and Medicine (INM-9), Computational Biomedicine, Forschungszentrum Jülich, 52425, Jülich, Germany
- Cécile and Oskar Vogt Institute for Brain Research, University Hospital Düsseldorf, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Christoph Fahlke
- Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, Forschungszentrum Jülich, 52425, Jülich, Germany.
| |
Collapse
|
5
|
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
|
6
|
Weigle AT, Shukla D. The Arabidopsis AtSWEET13 transporter discriminates sugars by selective facial and positional substrate recognition. Commun Biol 2024; 7:764. [PMID: 38914639 PMCID: PMC11196581 DOI: 10.1038/s42003-024-06291-6] [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/26/2022] [Accepted: 05/03/2024] [Indexed: 06/26/2024] Open
Abstract
Transporters are targeted by endogenous metabolites and exogenous molecules to reach cellular destinations, but it is generally not understood how different substrate classes exploit the same transporter's mechanism. Any disclosure of plasticity in transporter mechanism when treated with different substrates becomes critical for developing general selectivity principles in membrane transport catalysis. Using extensive molecular dynamics simulations with an enhanced sampling approach, we select the Arabidopsis sugar transporter AtSWEET13 as a model system to identify the basis for glucose versus sucrose molecular recognition and transport. Here we find that AtSWEET13 chemical selectivity originates from a conserved substrate facial selectivity demonstrated when committing alternate access, despite mono-/di-saccharides experiencing differing degrees of conformational and positional freedom throughout other stages of transport. However, substrate interactions with structural hallmarks associated with known functional annotations can help reinforce selective preferences in molecular transport.
Collapse
Affiliation(s)
- Austin T Weigle
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - Diwakar Shukla
- Department of Chemical & Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
| |
Collapse
|
7
|
Jaiswal N, Kumar A. Modulators of Candida albicans Membrane Drug Transporters: A Lucrative Portfolio for the Development of Effective Antifungals. Mol Biotechnol 2024; 66:960-974. [PMID: 38206530 DOI: 10.1007/s12033-023-01017-1] [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/28/2023] [Accepted: 12/01/2023] [Indexed: 01/12/2024]
Abstract
The escalating prevalence of membrane drug transporters and drug efflux pumps in pathogenic yeast like Candida albicans necessitates a comprehensive understanding of their roles in MDR. The overexpression of drug transporter families, ABC and MFS, implicated in MDR through drug efflux and poses a significant challenge in the diagnosis and treatment of fungal infection. Various mechanisms have been proposed for MDR; however, the upregulation of ABC and MFS superfamily transporters is most noticeable in MDR. The direct inhibition of these transporters seems an efficient strategy to overcome this problem. The goal of the article is to present an overview of the prospect of utilizing these modulators of C. albicans drug transports as effective antifungal molecules against MDR addressing a critical gap in the field. The review tries to address to prevent drug extrusion by modulating the expression of drug transporters of C. albicans. The review discussed the progress in identifying potent, selective, and non-toxic modulators of these transporters to develop some effective antifungals and overcome MDR. We reviewed major studies in this area and found that recent work has shifted toward the exploration of natural compounds as potential modulators to restore drug sensitivity in MDR fungal cells. The focus of this review is to survey and interpret current research information on modulators of C. albicans drug transporters from natural sources emphasizing those compounds that are potent antifungal agents.
Collapse
Affiliation(s)
- Neha Jaiswal
- Department of Biotechnology, National Institute of Technology, Raipur, CG, 492010, India
| | - Awanish Kumar
- Department of Biotechnology, National Institute of Technology, Raipur, CG, 492010, India.
| |
Collapse
|
8
|
Wu Y, Cao S, Qiu Y, Huang X. Tutorial on how to build non-Markovian dynamic models from molecular dynamics simulations for studying protein conformational changes. J Chem Phys 2024; 160:121501. [PMID: 38516972 PMCID: PMC10964226 DOI: 10.1063/5.0189429] [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: 11/28/2023] [Accepted: 02/20/2024] [Indexed: 03/23/2024] Open
Abstract
Protein conformational changes play crucial roles in their biological functions. In recent years, the Markov State Model (MSM) constructed from extensive Molecular Dynamics (MD) simulations has emerged as a powerful tool for modeling complex protein conformational changes. In MSMs, dynamics are modeled as a sequence of Markovian transitions among metastable conformational states at discrete time intervals (called lag time). A major challenge for MSMs is that the lag time must be long enough to allow transitions among states to become memoryless (or Markovian). However, this lag time is constrained by the length of individual MD simulations available to track these transitions. To address this challenge, we have recently developed Generalized Master Equation (GME)-based approaches, encoding non-Markovian dynamics using a time-dependent memory kernel. In this Tutorial, we introduce the theory behind two recently developed GME-based non-Markovian dynamic models: the quasi-Markov State Model (qMSM) and the Integrative Generalized Master Equation (IGME). We subsequently outline the procedures for constructing these models and provide a step-by-step tutorial on applying qMSM and IGME to study two peptide systems: alanine dipeptide and villin headpiece. This Tutorial is available at https://github.com/xuhuihuang/GME_tutorials. The protocols detailed in this Tutorial aim to be accessible for non-experts interested in studying the biomolecular dynamics using these non-Markovian dynamic models.
Collapse
Affiliation(s)
- Yue Wu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Author to whom correspondence should be addressed:
| |
Collapse
|
9
|
Ohnuki J, Jaunet-Lahary T, Yamashita A, Okazaki KI. Accelerated Molecular Dynamics and AlphaFold Uncover a Missing Conformational State of Transporter Protein OxlT. J Phys Chem Lett 2024; 15:725-732. [PMID: 38215403 DOI: 10.1021/acs.jpclett.3c03052] [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: 01/14/2024]
Abstract
Transporter proteins change their conformations to carry their substrate across the cell membrane. The conformational dynamics is vital to understanding the transport function. We have studied the oxalate transporter (OxlT), an oxalate:formate antiporter from Oxalobacter formigenes, significant in avoiding kidney stone formation. The atomic structure of OxlT has been recently solved in the outward-open and occluded states. However, the inward-open conformation is still missing, hindering a complete understanding of the transporter. Here, we performed a Gaussian accelerated molecular dynamics simulation to sample the extensive conformational space of OxlT and successfully predicted the inward-open conformation where cytoplasmic substrate formate binding was preferred over oxalate binding. We also identified critical interactions for the inward-open conformation. The results were complemented by an AlphaFold2 structure prediction. Although AlphaFold2 solely predicted OxlT in the outward-open conformation, mutation of the identified critical residues made it partly predict the inward-open conformation, identifying possible state-shifting mutations.
Collapse
Affiliation(s)
- Jun Ohnuki
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8585, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi 444-8585, Japan
| | - Titouan Jaunet-Lahary
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8585, Japan
| | - Atsuko Yamashita
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama 700-8530, Japan
| | - Kei-Ichi Okazaki
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8585, Japan
- Graduate Institute for Advanced Studies, SOKENDAI, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
10
|
Kleiman DE, Nadeem H, Shukla D. Adaptive Sampling Methods for Molecular Dynamics in the Era of Machine Learning. J Phys Chem B 2023; 127:10669-10681. [PMID: 38081185 DOI: 10.1021/acs.jpcb.3c04843] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Molecular dynamics (MD) simulations are fundamental computational tools for the study of proteins and their free energy landscapes. However, sampling protein conformational changes through MD simulations is challenging due to the relatively long time scales of these processes. Many enhanced sampling approaches have emerged to tackle this problem, including biased sampling and path-sampling methods. In this Perspective, we focus on adaptive sampling algorithms. These techniques differ from other approaches because the thermodynamic ensemble is preserved and the sampling is enhanced solely by restarting MD trajectories at particularly chosen seeds rather than introducing biasing forces. We begin our treatment with an overview of theoretically transparent methods, where we discuss principles and guidelines for adaptive sampling. Then, we present a brief summary of select methods that have been applied to realistic systems in the past. Finally, we discuss recent advances in adaptive sampling methodology powered by deep learning techniques, as well as their shortcomings.
Collapse
Affiliation(s)
- Diego E Kleiman
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Hassan Nadeem
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Diwakar Shukla
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
11
|
Chon NL, Schultz NJ, Zheng H, Lin H. Anion Pathways in the NarK Nitrate/Nitrite Exchanger. J Chem Inf Model 2023; 63:5142-5152. [PMID: 37585651 PMCID: PMC10482320 DOI: 10.1021/acs.jcim.3c00295] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Indexed: 08/18/2023]
Abstract
NarK nitrate/nitrite antiporter imports nitrate (a mineral form of the essential element nitrogen) into the cell and exports nitrite (a metabolite that can be toxic in high concentrations) out of the cell. However, many details about its operational mechanism remain poorly understood. In this work, we performed steered molecular dynamics simulations of anion translocations and quantum-chemistry model calculations of the binding sites to study the wild-type NarK protein and its R89K mutant. Our results shed light on the importance of the two strictly conserved binding-site arginine residues (R89 and R305) and two glycine-rich signature motifs (G164-M176 and G408-F419) in anion movement through the pore. We also observe conformational changes of the protein during anion migration. For the R89K mutant, our quantum calculations reveal a competition for a proton between the anion (especially nitrite) and lysine, which can potentially slow down or even trap the anion in the pore. Our findings provide a possible explanation for the striking experimental finding that the arginine-to-lysine mutation, despite preserving the charge, impedes or abolishes anion transport in such mutants of NarK and other similar nitrate/nitrite exchangers.
Collapse
Affiliation(s)
- Nara Lee Chon
- Department
of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Natalie Jean Schultz
- Department
of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| | - Hongjin Zheng
- Department
of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, Colorado 80045, United States
| | - Hai Lin
- Department
of Chemistry, University of Colorado Denver, Denver, Colorado 80217, United States
| |
Collapse
|
12
|
Qiu Y, O’Connor MS, Xue M, Liu B, Huang X. An Efficient Path Classification Algorithm Based on Variational Autoencoder to Identify Metastable Path Channels for Complex Conformational Changes. J Chem Theory Comput 2023; 19:4728-4742. [PMID: 37382437 PMCID: PMC11042546 DOI: 10.1021/acs.jctc.3c00318] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Conformational changes (i.e., dynamic transitions between pairs of conformational states) play important roles in many chemical and biological processes. Constructing the Markov state model (MSM) from extensive molecular dynamics (MD) simulations is an effective approach to dissect the mechanism of conformational changes. When combined with transition path theory (TPT), MSM can be applied to elucidate the ensemble of kinetic pathways connecting pairs of conformational states. However, the application of TPT to analyze complex conformational changes often results in a vast number of kinetic pathways with comparable fluxes. This obstacle is particularly pronounced in heterogeneous self-assembly and aggregation processes. The large number of kinetic pathways makes it challenging to comprehend the molecular mechanisms underlying conformational changes of interest. To address this challenge, we have developed a path classification algorithm named latent-space path clustering (LPC) that efficiently lumps parallel kinetic pathways into distinct metastable path channels, making them easier to comprehend. In our algorithm, MD conformations are first projected onto a low-dimensional space containing a small set of collective variables (CVs) by time-structure-based independent component analysis (tICA) with kinetic mapping. Then, MSM and TPT are constructed to obtain the ensemble of pathways, and a deep learning architecture named the variational autoencoder (VAE) is used to learn the spatial distributions of kinetic pathways in the continuous CV space. Based on the trained VAE model, the TPT-generated ensemble of kinetic pathways can be embedded into a latent space, where the classification becomes clear. We show that LPC can efficiently and accurately identify the metastable path channels in three systems: a 2D potential, the aggregation of two hydrophobic particles in water, and the folding of the Fip35 WW domain. Using the 2D potential, we further demonstrate that our LPC algorithm outperforms the previous path-lumping algorithms by making substantially fewer incorrect assignments of individual pathways to four path channels. We expect that LPC can be widely applied to identify the dominant kinetic pathways underlying complex conformational changes.
Collapse
Affiliation(s)
- Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Michael S. O’Connor
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Mingyi Xue
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Bojun Liu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| |
Collapse
|
13
|
Dominic AJ, Cao S, Montoya-Castillo A, Huang X. Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently. J Am Chem Soc 2023; 145:9916-9927. [PMID: 37104720 DOI: 10.1021/jacs.3c01095] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.
Collapse
Affiliation(s)
- Anthony J Dominic
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | | | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
14
|
Sauve S, Williamson J, Polasa A, Moradi M. Ins and Outs of Rocker Switch Mechanism in Major Facilitator Superfamily of Transporters. MEMBRANES 2023; 13:membranes13050462. [PMID: 37233523 DOI: 10.3390/membranes13050462] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Revised: 04/20/2023] [Accepted: 04/23/2023] [Indexed: 05/27/2023]
Abstract
The major facilitator superfamily (MFS) of transporters consists of three classes of membrane transporters: symporters, uniporters, and antiporters. Despite such diverse functions, MFS transporters are believed to undergo similar conformational changes within their distinct transport cycles, known as the rocker-switch mechanism. While the similarities between conformational changes are noteworthy, the differences are also important since they could potentially explain the distinct functions of symporters, uniporters, and antiporters of the MFS superfamily. We reviewed a variety of experimental and computational structural data on a select number of antiporters, symporters, and uniporters from the MFS family to compare the similarities and differences of the conformational dynamics of three different classes of transporters.
Collapse
Affiliation(s)
- Stephanie Sauve
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| | - Joseph Williamson
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| | - Adithya Polasa
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| | - Mahmoud Moradi
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, AR 72701, USA
| |
Collapse
|
15
|
Jaunet-Lahary T, Shimamura T, Hayashi M, Nomura N, Hirasawa K, Shimizu T, Yamashita M, Tsutsumi N, Suehiro Y, Kojima K, Sudo Y, Tamura T, Iwanari H, Hamakubo T, Iwata S, Okazaki KI, Hirai T, Yamashita A. Structure and mechanism of oxalate transporter OxlT in an oxalate-degrading bacterium in the gut microbiota. Nat Commun 2023; 14:1730. [PMID: 37012268 PMCID: PMC10070484 DOI: 10.1038/s41467-023-36883-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2021] [Accepted: 02/20/2023] [Indexed: 04/05/2023] Open
Abstract
An oxalate-degrading bacterium in the gut microbiota absorbs food-derived oxalate to use this as a carbon and energy source, thereby reducing the risk of kidney stone formation in host animals. The bacterial oxalate transporter OxlT selectively uptakes oxalate from the gut to bacterial cells with a strict discrimination from other nutrient carboxylates. Here, we present crystal structures of oxalate-bound and ligand-free OxlT in two distinct conformations, occluded and outward-facing states. The ligand-binding pocket contains basic residues that form salt bridges with oxalate while preventing the conformational switch to the occluded state without an acidic substrate. The occluded pocket can accommodate oxalate but not larger dicarboxylates, such as metabolic intermediates. The permeation pathways from the pocket are completely blocked by extensive interdomain interactions, which can be opened solely by a flip of a single side chain neighbouring the substrate. This study shows the structural basis underlying metabolic interactions enabling favourable symbiosis.
Collapse
Affiliation(s)
- Titouan Jaunet-Lahary
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, 444-8585, Japan
| | - Tatsuro Shimamura
- Graduate School of Medicine, Kyoto University, Kyoto, 606-8501, Japan.
| | - Masahiro Hayashi
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Norimichi Nomura
- Graduate School of Medicine, Kyoto University, Kyoto, 606-8501, Japan
| | - Kouta Hirasawa
- Graduate School of Medicine, Kyoto University, Kyoto, 606-8501, Japan
| | | | | | - Naotaka Tsutsumi
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
- School of Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Yuta Suehiro
- School of Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Keiichi Kojima
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Yuki Sudo
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Takashi Tamura
- Graduate School of Environmental and Life Sciences, Okayama University, Okayama, 700-8530, Japan
| | - Hiroko Iwanari
- Research Center for Advanced Science and Technology, The University of Tokyo, Tokyo, 153-8904, Japan
| | - Takao Hamakubo
- Research Center for Advanced Science and Technology, The University of Tokyo, Tokyo, 153-8904, Japan
| | - So Iwata
- Graduate School of Medicine, Kyoto University, Kyoto, 606-8501, Japan
| | - Kei-Ichi Okazaki
- Research Center for Computational Science, Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, 444-8585, Japan.
| | | | - Atsuko Yamashita
- Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan.
- RIKEN SPring-8 Center, Sayo, 679-5148, Japan.
- School of Pharmaceutical Sciences, Okayama University, Okayama, 700-8530, Japan.
| |
Collapse
|
16
|
Unarta IC, Goonetilleke EC, Wang D, Huang X. Nucleotide addition and cleavage by RNA polymerase II: Coordination of two catalytic reactions using a single active site. J Biol Chem 2022; 299:102844. [PMID: 36581202 PMCID: PMC9860460 DOI: 10.1016/j.jbc.2022.102844] [Citation(s) in RCA: 6] [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: 07/03/2022] [Revised: 12/19/2022] [Accepted: 12/22/2022] [Indexed: 12/28/2022] Open
Abstract
RNA polymerase II (Pol II) incorporates complementary ribonucleotides into the growing RNA chain one at a time via the nucleotide addition cycle. The nucleotide addition cycle, however, is prone to misincorporation of noncomplementary nucleotides. Thus, to ensure transcriptional fidelity, Pol II backtracks and then cleaves the misincorporated nucleotides. These two reverse reactions, nucleotide addition and cleavage, are catalyzed in the same active site of Pol II, which is different from DNA polymerases or other endonucleases. Recently, substantial progress has been made to understand how Pol II effectively performs its dual role in the same active site. Our review highlights these recent studies and provides an overall model of the catalytic mechanisms of Pol II. In particular, RNA extension follows the two-metal-ion mechanism, and several Pol II residues play important roles to facilitate the catalysis. In sharp contrast, the cleavage reaction is independent of any Pol II residues. Interestingly, Pol II relies on its residues to recognize the misincorporated nucleotides during the backtracking process, prior to cleavage. In this way, Pol II efficiently compartmentalizes its two distinct catalytic functions using the same active site. Lastly, we also discuss a new perspective on the potential third Mg2+ in the nucleotide addition and intrinsic cleavage reactions.
Collapse
Affiliation(s)
- Ilona Christy Unarta
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Eshani C Goonetilleke
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Dong Wang
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, California, USA; Department of Cellular and Molecular Medicine, School of Medicine, University of California, San Diego, La Jolla, California, USA; Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California, USA.
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin, USA.
| |
Collapse
|
17
|
Kleiman DE, Shukla D. Multiagent Reinforcement Learning-Based Adaptive Sampling for Conformational Dynamics of Proteins. J Chem Theory Comput 2022; 18:5422-5434. [PMID: 36044642 DOI: 10.1021/acs.jctc.2c00683] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Machine learning is increasingly applied to improve the efficiency and accuracy of molecular dynamics (MD) simulations. Although the growth of distributed computer clusters has allowed researchers to obtain higher amounts of data, unbiased MD simulations have difficulty sampling rare states, even under massively parallel adaptive sampling schemes. To address this issue, several algorithms inspired by reinforcement learning (RL) have arisen to promote exploration of the slow collective variables (CVs) of complex systems. Nonetheless, most of these algorithms are not well-suited to leverage the information gained by simultaneously sampling a system from different initial states (e.g., a protein in different conformations associated with distinct functional states). To fill this gap, we propose two algorithms inspired by multiagent RL that extend the functionality of closely related techniques (REAP and TSLC) to situations where the sampling can be accelerated by learning from different regions of the energy landscape through coordinated agents. Essentially, the algorithms work by remembering which agent discovered each conformation and sharing this information with others at the action-space discretization step. A stakes function is introduced to modulate how different agents sense rewards from discovered states of the system. The consequences are three-fold: (i) agents learn to prioritize CVs using only relevant data, (ii) redundant exploration is reduced, and (iii) agents that obtain higher stakes are assigned more actions. We compare our algorithm with other adaptive sampling techniques (least counts, REAP, TSLC, and AdaptiveBandit) to show and rationalize the gain in performance.
Collapse
Affiliation(s)
- Diego E Kleiman
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Diwakar Shukla
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.,Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
18
|
Integration of machine learning with computational structural biology of plants. Biochem J 2022; 479:921-928. [PMID: 35484946 DOI: 10.1042/bcj20200942] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 04/01/2022] [Accepted: 04/06/2022] [Indexed: 11/17/2022]
Abstract
Computational structural biology of proteins has developed rapidly in recent decades with the development of new computational tools and the advancement of computing hardware. However, while these techniques have widely been used to make advancements in human medicine, these methods have seen less utilization in the plant sciences. In the last several years, machine learning methods have gained popularity in computational structural biology. These methods have enabled the development of new tools which are able to address the major challenges that have hampered the wide adoption of the computational structural biology of plants. This perspective examines the remaining challenges in computational structural biology and how the development of machine learning techniques enables more in-depth computational structural biology of plants.
Collapse
|
19
|
Chan MC, Procko E, Shukla D. Structural Rearrangement of the Serotonin Transporter Intracellular Gate Induced by Thr276 Phosphorylation. ACS Chem Neurosci 2022; 13:933-945. [PMID: 35258286 DOI: 10.1021/acschemneuro.1c00714] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
The reuptake of the neurotransmitter serotonin from the synaptic cleft by the serotonin transporter, SERT, is essential for proper neurological signaling. Biochemical studies have shown that Thr276 of transmembrane helix 5 is a site of PKG-mediated SERT phosphorylation, which has been proposed to shift the SERT conformational equilibria to promote inward-facing states, thus enhancing 5-HT transport. Recent structural and simulation studies have provided insights into the conformation transitions during substrate transport but have not shed light on SERT regulation via post-translational modifications. Using molecular dynamics simulations and Markov state models, we investigate how Thr276 phosphorylation impacts the SERT mechanism and its role in enhancing transporter stability and function. Our simulations show that Thr276 phosphorylation alters the hydrogen-bonding network involving residues on transmembrane helix 5. This in turn decreases the free energy barriers for SERT to transition to the inward-facing state, thus facilitating 5-HT import. The results provide atomistic insights into in vivo SERT regulation and can be extended to other pharmacologically important transporters in the solute carrier family.
Collapse
Affiliation(s)
- Matthew C. Chan
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Erik Procko
- Department of Biochemistry, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Neuroscience Program, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Cancer Center at Illinois, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Cancer Center at Illinois, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- National Center for Supercomputing Applications, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
- Department of Bioengineering, University of Illinois Urbana−Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
20
|
Konovalov K, Unarta IC, Cao S, Goonetilleke EC, Huang X. Markov State Models to Study the Functional Dynamics of Proteins in the Wake of Machine Learning. JACS AU 2021; 1:1330-1341. [PMID: 34604842 PMCID: PMC8479766 DOI: 10.1021/jacsau.1c00254] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Indexed: 05/19/2023]
Abstract
Markov state models (MSMs) based on molecular dynamics (MD) simulations are routinely employed to study protein folding, however, their application to functional conformational changes of biomolecules is still limited. In the past few years, the field of computational chemistry has experienced a surge of advancements stemming from machine learning algorithms, and MSMs have not been left out. Unlike global processes, such as protein folding, the application of MSMs to functional conformational changes is challenging because they mostly consist of localized structural transitions. Therefore, it is critical to properly select a subset of structural features that can describe the slowest dynamics of these functional conformational changes. To address this challenge, we recommend several automatic feature selection methods such as Spectral-OASIS. To identify states in MSMs, the chosen features can be subject to dimensionality reduction methods such as TICA or deep learning based VAMPNets to project MD conformations onto a few collective variables for subsequent clustering. Another challenge for the application of MSMs to the study of functional conformational changes is the ability to comprehend their biophysical mechanisms, as MSMs built for these processes often require a large number of states. We recommend the recently developed quasi-MSMs (qMSMs) to address this issue. Compared to MSMs, qMSMs encode the non-Markovian dynamics via the generalized master equation and can significantly reduce the number of states. As a result, qMSMs can be built with a handful of states to facilitate the interpretation of functional conformational changes. In the wake of machine learning, we believe that the rapid advancement in the MSM methodology will lead to their wider application in studying functional conformational changes of biomolecules.
Collapse
Affiliation(s)
- Kirill
A. Konovalov
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Ilona Christy Unarta
- Department
of Chemical and Biological Engineering, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Siqin Cao
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Eshani C. Goonetilleke
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Xuhui Huang
- Department
of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Department
of Chemical and Biological Engineering, The Hong Kong University of Science and Technology, Kowloon, Hong Kong
- Hong
Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| |
Collapse
|