51
|
Yang JF, Wang F, Chen YZ, Hao GF, Yang GF. LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor. Brief Bioinform 2019; 21:2206-2218. [DOI: 10.1093/bib/bbz141] [Citation(s) in RCA: 69] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2019] [Revised: 10/12/2019] [Accepted: 10/14/2019] [Indexed: 12/14/2022] Open
Abstract
Abstract
Protein dynamics is central to all biological processes, including signal transduction, cellular regulation and biological catalysis. Among them, in-depth exploration of ligand-driven protein dynamics contributes to an optimal understanding of protein function, which is particularly relevant to drug discovery. Hence, a wide range of computational tools have been designed to investigate the important dynamic information in proteins. However, performing and analyzing protein dynamics is still challenging due to the complicated operation steps, giving rise to great difficulty, especially for nonexperts. Moreover, there is a lack of web protocol to provide online facility to investigate and visualize ligand-driven protein dynamics. To this end, in this study, we integrated several bioinformatic tools to develop a protocol, named Ligand and Receptor Molecular Dynamics (LARMD, http://chemyang.ccnu.edu.cn/ccb/server/LARMD/ and http://agroda.gzu.edu.cn:9999/ccb/server/LARMD/), for profiling ligand-driven protein dynamics. To be specific, estrogen receptor (ER) was used as a case to reveal ERβ-selective mechanism, which plays a vital role in the treatment of inflammatory diseases and many types of cancers in clinical practice. Two different residues (Ile373/Met421 and Met336/Leu384) in the pocket of ERβ/ERα were the significant determinants for selectivity, especially Met336 of ERβ. The helix H8, helix H11 and H7-H8 loop influenced the migration of selective agonist (WAY-244). These computational results were consistent with the experimental results. Therefore, LARMD provides a user-friendly online protocol to study the dynamic property of protein and to design new ligand or site-directed mutagenesis.
Collapse
Affiliation(s)
- Jing-Fang Yang
- Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P.R.China
- International Joint Research Center for Intelligent Biosensor Technology and Health, Central China Normal University,Wuhan, 430079, China
| | - Fan Wang
- Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P.R.China
- International Joint Research Center for Intelligent Biosensor Technology and Health, Central China Normal University,Wuhan, 430079, China
| | - Yu-Zong Chen
- Department of Pharmacy, National University of Singapore, 18 Science Drive 4, Singapore 117543
| | - Ge-Fei Hao
- Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P.R.China
- International Joint Research Center for Intelligent Biosensor Technology and Health, Central China Normal University,Wuhan, 430079, China
- State Key Laboratory Breeding Base of Green Pesticide and Agricultural Bioengineering, Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Research and Development Center for Fine Chemicals, Guizhou University, Guiyang 550025, P. R. China
| | - Guang-Fu Yang
- Key Laboratory of Pesticide & Chemical Biology, Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, P.R.China
- International Joint Research Center for Intelligent Biosensor Technology and Health, Central China Normal University,Wuhan, 430079, China
- Collaborative Innovation Center of Chemical Science and Engineering, Tianjing 300072, P.R.China
| |
Collapse
|
52
|
Tian J, Liu F, Xu Z, Shi J, Liang T, Zhang Y, Da LT. Regulatory Role of One Critical Catalytic Loop of Polypeptide N-Acetyl-Galactosaminyltransferase-2 in Substrate Binding and Catalysis during Mucin-Type O-Glycosylation. ACS Catal 2019. [DOI: 10.1021/acscatal.9b03782] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Jiaqi Tian
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Feng Liu
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Zhijue Xu
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Jingjing Shi
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Tao Liang
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Yan Zhang
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| | - Lin-Tai Da
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
| |
Collapse
|
53
|
Li D, Ji B. Protein conformational transitions coupling with ligand interactions: Simulations from molecules to medicine. MEDICINE IN NOVEL TECHNOLOGY AND DEVICES 2019. [DOI: 10.1016/j.medntd.2019.100026] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|
54
|
Encounter complexes and hidden poses of kinase-inhibitor binding on the free-energy landscape. Proc Natl Acad Sci U S A 2019; 116:18404-18409. [PMID: 31451651 PMCID: PMC6744929 DOI: 10.1073/pnas.1904707116] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Modern drug discovery increasingly focuses on the drug-target binding kinetics which depend on drug (un)binding pathways. The conventional molecular dynamics simulation can observe only a few binding events even using the fastest supercomputer. Here, we develop 2D gREST/REUS simulation with enhanced flexibility of the ligand and the protein binding site. Simulation (43 μs in total) applied to an inhibitor binding to c-Src kinase covers 100 binding and unbinding events. On the statistically converged free-energy landscapes, we succeed in predicting the X-ray binding structure, including water positions. Furthermore, we characterize hidden semibound poses and transient encounter complexes on the free-energy landscapes. Regulatory residues distant from the catalytic core are responsible for the initial inhibitor uptake and regulation of subsequent bindings, which was unresolved by experiments. Stabilizing/blocking of either the semibound poses or the encounter complexes can be an effective strategy to optimize drug-target residence time.
Collapse
|
55
|
Vergara R, Romero‐Romero S, Velázquez‐López I, Espinoza‐Pérez G, Rodríguez‐Hernández A, Pulido NO, Sosa‐Peinado A, Rodríguez‐Romero A, Fernández‐Velasco DA. The interplay of protein–ligand and water‐mediated interactions shape affinity and selectivity in the LAO binding protein. FEBS J 2019; 287:763-782. [DOI: 10.1111/febs.15019] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Revised: 06/25/2019] [Accepted: 07/24/2019] [Indexed: 12/16/2022]
Affiliation(s)
- Renan Vergara
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| | - Sergio Romero‐Romero
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| | - Isabel Velázquez‐López
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| | - Georgina Espinoza‐Pérez
- Laboratorio de Química de Biomacromoléculas 3, Departamento de Química de Biomacromoléculas, Instituto de Química Universidad Nacional Autónoma de México Ciudad de México México
| | - Annia Rodríguez‐Hernández
- Laboratorio de Química de Biomacromoléculas 3, Departamento de Química de Biomacromoléculas, Instituto de Química Universidad Nacional Autónoma de México Ciudad de México México
| | - Nancy O. Pulido
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| | - Alejandro Sosa‐Peinado
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| | - Adela Rodríguez‐Romero
- Laboratorio de Química de Biomacromoléculas 3, Departamento de Química de Biomacromoléculas, Instituto de Química Universidad Nacional Autónoma de México Ciudad de México México
| | - Daniel Alejandro Fernández‐Velasco
- Laboratorio de Fisicoquímica e Ingeniería de Proteínas, Departamento de Bioquímica, Facultad de Medicina Universidad Nacional Autónoma de México Ciudad de México México
| |
Collapse
|
56
|
Sumino A, Sumikama T, Uchihashi T, Oiki S. High-speed AFM reveals accelerated binding of agitoxin-2 to a K + channel by induced fit. SCIENCE ADVANCES 2019; 5:eaax0495. [PMID: 31281899 PMCID: PMC6609221 DOI: 10.1126/sciadv.aax0495] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2019] [Accepted: 05/23/2019] [Indexed: 06/09/2023]
Abstract
Agitoxin-2 (AgTx2) from scorpion venom is a potent blocker of K+ channels. The docking model has been elucidated, but it remains unclear whether binding dynamics are described by a two-state model (AgTx2-bound and AgTx2-unbound) or a more complicated mechanism, such as induced fit or conformational selection. Here, we observed the binding dynamics of AgTx2 to the KcsA channel using high-speed atomic force microscopy. From images of repeated binding and dissociation of AgTx2 to the channel, single-molecule kinetic analyses revealed that the affinity of the channel for AgTx2 increased during persistent binding and decreased during persistent dissociation. We propose a four-state model, including high- and low-affinity states of the channel, with relevant rate constants. An induced-fit pathway was dominant and accelerated binding by 400 times. This is the first analytical imaging of scorpion toxin binding in real time, which is applicable to various biological dynamics including channel ligands, DNA-modifier proteins, and antigen-antibody complexes.
Collapse
Affiliation(s)
- A. Sumino
- Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kanazawa 920-1192, Japan
- Institute for Frontier Science Initiative, Kanazawa University, Kanazawa 920-1192, Japan
- PRESTO, Japan Science and Technology Agency (JST), Saitama 332-0012, Japan
- Department of Molecular Physiology and Biophysics, Faculty of Medical Sciences, University of Fukui, Fukui 910-1193, Japan
| | - T. Sumikama
- Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kanazawa 920-1192, Japan
- Department of Molecular Physiology and Biophysics, Faculty of Medical Sciences, University of Fukui, Fukui 910-1193, Japan
| | - T. Uchihashi
- Department of Physics and Structural Biology Research Center, Nagoya University, Chikusa-ku, Nagoya 464-8602, Japan
- Exploratory Research Center on Life and Living Systems, National Institute of Natural Sciences, 5-1 Higashiyama, Myodaiji, Okazaki, Aichi 444-8787, Japan
| | - S. Oiki
- Department of Molecular Physiology and Biophysics, Faculty of Medical Sciences, University of Fukui, Fukui 910-1193, Japan
| |
Collapse
|
57
|
Egawa T, Callender R. General mathematical formula for near equilibrium relaxation kinetics of basic enzyme reactions and its applications to find conformational selection steps. Math Biosci 2019; 313:61-70. [PMID: 30935841 DOI: 10.1016/j.mbs.2019.03.007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 03/27/2019] [Accepted: 03/27/2019] [Indexed: 10/27/2022]
Abstract
A general mathematical formula of basic enzyme reactions was derived with nearly no dependence on conditions nor assumptions on relaxation kinetic processes near equilibrium in a simple single-substrate-single-product enzyme reaction. The new formula gives precise relationships between the rate constants of the elementary reaction steps and the apparent relaxation rate constant, rather than the initial velocity that is generally used to determine enzymatic parameters according to the Michaelis-Menten theory. The present formula is shown to be complementary to the Michaelis-Menten formulae in a sense that the initial velocity and the relaxation rate constant data together could determine the enzyme-substrate dissociation constant KES, which has been usually conditionally approximated by the Michaelis constant KM within the framework of the Michaelis-Menten formulae. We also describe relaxation kinetics of enzyme reactions that include the conformational selection processes, in which only one enzymatic conformer among a conformational ensemble can bind with either the substrate or product. The present mathematical approaches, together with numerical computation analyses, suggested that the presence of conformational selection steps in enzymatic reactions can be experimentally detected simply by enzymatic assays with catalytic amounts of enzyme.
Collapse
Affiliation(s)
- Tsuyoshi Egawa
- Department of Biochemistry, Albert Einstein College of Medicine, United States.
| | - Robert Callender
- Department of Biochemistry, Albert Einstein College of Medicine, United States.
| |
Collapse
|
58
|
Zhu L, Sheong FK, Cao S, Liu S, Unarta IC, Huang X. TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromolecules. J Chem Phys 2019; 150:124105. [DOI: 10.1063/1.5082633] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen, Guangdong 518172, China
| | - Fu Kit Sheong
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Siqin Cao
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Song Liu
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Ilona C. Unarta
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Bioengineering Program, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- HKUST-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China
| |
Collapse
|
59
|
Harada R, Shigeta Y. How low-resolution structural data predict the conformational changes of a protein: a study on data-driven molecular dynamics simulations. Phys Chem Chem Phys 2019; 20:17790-17798. [PMID: 29922770 DOI: 10.1039/c8cp02246a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Parallel cascade selection molecular dynamics (PaCS-MD) is a conformational sampling method for generating transition pathways between a given reactant and a product. PaCS-MD repeats the following two steps: (1) selections of initial structures relevant to transitions and (2) their conformational resampling. When selecting the initial structures, several measures are utilized to identify their potential to undergo transitions. In the present study, low-resolution structural data obtained from small angle scattering (SAXS) and cryo-electron microscopy (EM) are adopted as the measures in PaCS-MD to promote the conformational transitions of proteins, which is defined as SAXS-/EM-driven targeted PaCS-MD. By selecting the essential structures that have high correlations with the low-resolution structural data, the SAXS-/EM-driven targeted PaCS-MD identifies a set of transition pathways between the reactant and the product. As a demonstration, the present method successfully predicted the open-closed transition pathway of the lysine-, arginine-, ornithine-binding protein with a ns-order simulation time, indicating that the data-driven PaCS-MD simulation might work to promote the conformational transitions of proteins efficiently.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Ibaraki 305-8577, Japan.
| | | |
Collapse
|
60
|
Swenson DWH, Prinz JH, Noe F, Chodera JD, Bolhuis PG. OpenPathSampling: A Python Framework for Path Sampling Simulations. 2. Building and Customizing Path Ensembles and Sample Schemes. J Chem Theory Comput 2019; 15:837-856. [PMID: 30359525 PMCID: PMC6374748 DOI: 10.1021/acs.jctc.8b00627] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2018] [Indexed: 11/28/2022]
Abstract
The OpenPathSampling (OPS) package provides an easy-to-use framework to apply transition path sampling methodologies to complex molecular systems with a minimum of effort. Yet, the extensibility of OPS allows for the exploration of new path sampling algorithms by building on a variety of basic operations. In a companion paper [ Swenson et al. J. Chem. Theory Comput. 2018 , 10.1021/acs.jctc.8b00626 ] we introduced the basic concepts and the structure of the OPS package, and how it can be employed to perform standard transition path sampling and (replica exchange) transition interface sampling. In this paper, we elaborate on two theoretical developments that went into the design of OPS. The first development relates to the construction of path ensembles, the what is being sampled. We introduce a novel set-based notation for the path ensemble, which provides an alternative paradigm for constructing path ensembles and allows building arbitrarily complex path ensembles from fundamental ones. The second fundamental development is the structure for the customization of Monte Carlo procedures; how path ensembles are being sampled. We describe in detail the OPS objects that implement this approach to customization, the MoveScheme and the PathMover, and provide tools to create and manipulate these objects. We illustrate both the path ensemble building and sampling scheme customization with several examples. OPS thus facilitates both standard path sampling application in complex systems as well as the development of new path sampling methodology, beyond the default.
Collapse
Affiliation(s)
- David W. H. Swenson
- van’t
Hoff Institute for Molecular Sciences, University
of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
- Computational
and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Jan-Hendrik Prinz
- Computational
and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
- Department
of Mathematics and Computer Science, Arnimallee 6, Freie Universität Berlin, 14195 Berlin, Germany
| | - Frank Noe
- Department
of Mathematics and Computer Science, Arnimallee 6, Freie Universität Berlin, 14195 Berlin, Germany
| | - John D. Chodera
- Computational
and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Peter G. Bolhuis
- van’t
Hoff Institute for Molecular Sciences, University
of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
61
|
Wang D, Weng J, Wang W. An unconventional ligand‐binding mechanism of substrate‐binding proteins: MD simulation and Markov state model analysis of BtuF. J Comput Chem 2019; 40:1440-1448. [DOI: 10.1002/jcc.25798] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2018] [Revised: 12/22/2018] [Accepted: 01/28/2019] [Indexed: 01/01/2023]
Affiliation(s)
- Dongdong Wang
- Department of Chemistry, Institutes of Biomedical Sciences and Multiscale Research Institute of Complex System Fudan University Shanghai 200438 People's Republic of China
| | - Jingwei Weng
- Department of Chemistry, Institutes of Biomedical Sciences and Multiscale Research Institute of Complex System Fudan University Shanghai 200438 People's Republic of China
| | - Wenning Wang
- Department of Chemistry, Institutes of Biomedical Sciences and Multiscale Research Institute of Complex System Fudan University Shanghai 200438 People's Republic of China
| |
Collapse
|
62
|
On identifying collective displacements in apo-proteins that reveal eventual binding pathways. PLoS Comput Biol 2019; 15:e1006665. [PMID: 30645590 PMCID: PMC6333327 DOI: 10.1371/journal.pcbi.1006665] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2018] [Accepted: 11/23/2018] [Indexed: 01/19/2023] Open
Abstract
Binding of small molecules to proteins often involves large conformational changes in the latter, which open up pathways to the binding site. Observing and pinpointing these rare events in large scale, all-atom, computations of specific protein-ligand complexes, is expensive and to a great extent serendipitous. Further, relevant collective variables which characterise specific binding or un-binding scenarios are still difficult to identify despite the large body of work on the subject. Here, we show that possible primary and secondary binding pathways can be discovered from short simulations of the apo-protein without waiting for an actual binding event to occur. We use a projection formalism, introduced earlier to study deformation in solids, to analyse local atomic displacements into two mutually orthogonal subspaces—those which are “affine” i.e. expressible as a homogeneous deformation of the native structure, and those which are not. The susceptibility to non-affine displacements among the various residues in the apo- protein is then shown to correlate with typical binding pathways and sites crucial for allosteric modifications. We validate our observation with all-atom computations of three proteins, T4-Lysozyme, Src kinase and Cytochrome P450. Designing drugs which target specific proteins involved in diseases consumes a lot of time and effort in the pharmaceutical industry. In recent times, in silico design of drugs using all-atom molecular modelling has started to provide crucial inputs. Even so, discovery of binding pathways of small molecules both at the primary binding site, as well as sites for allosteric control, is time consuming and often fortuitous. We provide here a framework within which critical conformational changes likely to occur during binding are quantified from statistical analysis of configurations of proteins in their apo, or inactive form, greatly simplifying identification of target residues. We illustrate this idea by analysing ligand binding pathways for three proteins T4- Lysozyme, P450 and Src kinase, which are active respectively in the immune system, metabolism and cancer.
Collapse
|
63
|
Swenson DWH, Prinz JH, Noe F, Chodera JD, Bolhuis PG. OpenPathSampling: A Python Framework for Path Sampling Simulations. 1. Basics. J Chem Theory Comput 2018; 15:813-836. [PMID: 30336030 PMCID: PMC6374749 DOI: 10.1021/acs.jctc.8b00626] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
![]()
Transition
path sampling techniques allow molecular dynamics simulations of complex
systems to focus on rare dynamical events, providing
insight into mechanisms and the ability to calculate rates inaccessible
by ordinary dynamics simulations. While path sampling algorithms are
conceptually as simple as importance sampling Monte Carlo, the technical
complexity of their implementation has kept these techniques out of
reach of the broad community. Here, we introduce an easy-to-use Python
framework called OpenPathSampling (OPS) that facilitates path sampling
for (bio)molecular systems with minimal effort and yet is still extensible.
Interfaces to OpenMM and an internal dynamics engine for simple models
are provided in the initial release, but new molecular simulation
packages can easily be added. Multiple ready-to-use transition path
sampling methodologies are implemented, including standard transition
path sampling (TPS) between reactant and product states and transition
interface sampling (TIS) and its replica exchange variant (RETIS),
as well as recent multistate and multiset extensions of transition
interface sampling (MSTIS, MISTIS). In addition, tools are provided
to facilitate the implementation of new path sampling schemes built
on basic path sampling components. In this paper, we give an overview
of the design of this framework and illustrate the simplicity of applying
the available path sampling algorithms to a variety of benchmark problems.
Collapse
Affiliation(s)
- David W H Swenson
- van 't Hoff Institute for Molecular Sciences , University of Amsterdam , P.O. Box 94157, 1090 GD Amsterdam , The Netherlands.,Computational and Systems Biology Program , Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center , New York , New York 10065 , United States
| | - Jan-Hendrik Prinz
- Computational and Systems Biology Program , Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center , New York , New York 10065 , United States.,Department of Mathematics and Computer Science, Arnimallee 6 , Freie Universität Berlin , 14195 Berlin , Germany
| | - Frank Noe
- Department of Mathematics and Computer Science, Arnimallee 6 , Freie Universität Berlin , 14195 Berlin , Germany
| | - John D Chodera
- Computational and Systems Biology Program , Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center , New York , New York 10065 , United States
| | - Peter G Bolhuis
- van 't Hoff Institute for Molecular Sciences , University of Amsterdam , P.O. Box 94157, 1090 GD Amsterdam , The Netherlands
| |
Collapse
|
64
|
Konovalov KA, Wang W, Huang X. Conformational selection turns on phenylalanine hydroxylase. J Biol Chem 2018; 293:19544-19545. [PMID: 30578407 DOI: 10.1074/jbc.h118.006676] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Phenylalanine hydroxylase catalyzes a critical step in the phenylalanine catabolic pathway, and impairment of the human enzyme is linked to phenylketonuria. Phenylalanine is also a positive allosteric regulator of the enzyme, and the allosteric binding site has been determined by crystallography. However, the allosteric activation mechanism remains unclear. Using large-scale simulations to explore how phenylalanine binds to the regulatory site, Ge et al. discovered gating motions of the protein that suggest a conformational selection mechanism.
Collapse
Affiliation(s)
- Kirill A Konovalov
- From the Department of Chemistry, Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Wei Wang
- From the Department of Chemistry, Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- From the Department of Chemistry, Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| |
Collapse
|
65
|
Yang W, Riley BT, Lei X, Porebski BT, Kass I, Buckle AM, McGowan S. Mapping the Pathway and Dynamics of Bestatin Inhibition of the
Plasmodium falciparum
M1 Aminopeptidase
Pf
A‐M1. ChemMedChem 2018; 13:2504-2513. [DOI: 10.1002/cmdc.201800563] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2018] [Revised: 10/04/2018] [Indexed: 12/26/2022]
Affiliation(s)
- Wei Yang
- Infection and Immunity Program, Department of Microbiology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
| | - Blake T. Riley
- Infection and Immunity Program, Department of Biochemistry and Molecular Biology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
| | - Xiangyun Lei
- School of Chemical and Biomolecular EngineeringGeorgia Institute of Technology Atlanta GA USA
| | - Benjamin T. Porebski
- Infection and Immunity Program, Department of Biochemistry and Molecular Biology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
- Current address: Medical Research Council Laboratory of Molecular Biology Francis Crick Avenue Cambridge CB2 0QH UK
| | - Itamar Kass
- Infection and Immunity Program, Department of Biochemistry and Molecular Biology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
- Victorian Life Sciences Computation CentreMonash University Clayton VIC 3800 Australia
| | - Ashley M. Buckle
- Infection and Immunity Program, Department of Biochemistry and Molecular Biology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
| | - Sheena McGowan
- Infection and Immunity Program, Department of Microbiology, Biomedicine Discovery InstituteMonash University Clayton VIC 3800 Australia
| |
Collapse
|
66
|
Menzer WM, Li C, Sun W, Xie B, Minh DDL. Simple Entropy Terms for End-Point Binding Free Energy Calculations. J Chem Theory Comput 2018; 14:6035-6049. [PMID: 30296084 DOI: 10.1021/acs.jctc.8b00418] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We introduce a number of computationally inexpensive modifications to the MM/PBSA and MM/GBSA estimators for binding free energies, which are based on average receptor-ligand interaction energies in simulations of a noncovalent complex, to improve the treatment of entropy: second- and higher-order terms in a cumulant expansion and a confining potential on ligand external degrees of freedom. We also consider a filter for snapshots where ligands have drifted from the initial binding pose. The variations were tested on six sets of systems for which binding modes and free energies have previously been experimentally determined. For some data sets, none of the tested estimators led to results significantly correlated with measured free energies. In data sets with nontrivial correlation, a ligand RMSD cutoff of 3 Å and a second-order truncation of the cumulant expansion was found to be comparable or better than the average interaction energy by several statistical metrics.
Collapse
|
67
|
Sun X, Singh S, Blumer KJ, Bowman GR. Simulation of spontaneous G protein activation reveals a new intermediate driving GDP unbinding. eLife 2018; 7:e38465. [PMID: 30289386 PMCID: PMC6224197 DOI: 10.7554/elife.38465] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Accepted: 10/04/2018] [Indexed: 12/12/2022] Open
Abstract
Activation of heterotrimeric G proteins is a key step in many signaling cascades. However, a complete mechanism for this process, which requires allosteric communication between binding sites that are ~30 Å apart, remains elusive. We construct an atomically detailed model of G protein activation by combining three powerful computational methods: metadynamics, Markov state models (MSMs), and CARDS analysis of correlated motions. We uncover a mechanism that is consistent with a wide variety of structural and biochemical data. Surprisingly, the rate-limiting step for GDP release correlates with tilting rather than translation of the GPCR-binding helix 5. β-Strands 1 - 3 and helix 1 emerge as hubs in the allosteric network that links conformational changes in the GPCR-binding site to disordering of the distal nucleotide-binding site and consequent GDP release. Our approach and insights provide foundations for understanding disease-implicated G protein mutants, illuminating slow events in allosteric networks, and examining unbinding processes with slow off-rates.
Collapse
Affiliation(s)
- Xianqiang Sun
- Department of Biochemistry and Molecular BiophysicsWashington University School of MedicineMissouriUnited States
| | - Sukrit Singh
- Department of Biochemistry and Molecular BiophysicsWashington University School of MedicineMissouriUnited States
| | - Kendall J Blumer
- Department of Cell Biology and PhysiologyWashington University School of MedicineMissouriUnited States
| | - Gregory R Bowman
- Department of Biochemistry and Molecular BiophysicsWashington University School of MedicineMissouriUnited States
- Center for Biological Systems EngineeringWashington University School of MedicineMissouriUnited States
| |
Collapse
|
68
|
Huggins DJ, Biggin PC, Dämgen MA, Essex JW, Harris SA, Henchman RH, Khalid S, Kuzmanic A, Laughton CA, Michel J, Mulholland AJ, Rosta E, Sansom MSP, van der Kamp MW. Biomolecular simulations: From dynamics and mechanisms to computational assays of biological activity. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1393] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- David J. Huggins
- TCM Group, Cavendish Laboratory University of Cambridge Cambridge UK
- Unilever Centre, Department of Chemistry University of Cambridge Cambridge UK
- Department of Physiology and Biophysics Weill Cornell Medical College New York NY
| | | | - Marc A. Dämgen
- Department of Biochemistry University of Oxford Oxford UK
| | - Jonathan W. Essex
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | - Sarah A. Harris
- School of Physics and Astronomy University of Leeds Leeds UK
- Astbury Centre for Structural and Molecular Biology University of Leeds Leeds UK
| | - Richard H. Henchman
- Manchester Institute of Biotechnology The University of Manchester Manchester UK
- School of Chemistry The University of Manchester Oxford UK
| | - Syma Khalid
- School of Chemistry University of Southampton Southampton UK
- Institute for Life Sciences University of Southampton Southampton UK
| | | | - Charles A. Laughton
- School of Pharmacy University of Nottingham Nottingham UK
- Centre for Biomolecular Sciences University of Nottingham Nottingham UK
| | - Julien Michel
- EaStCHEM school of Chemistry University of Edinburgh Edinburgh UK
| | - Adrian J. Mulholland
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
| | - Edina Rosta
- Department of Chemistry King's College London London UK
| | | | - Marc W. van der Kamp
- Centre of Computational Chemistry, School of Chemistry University of Bristol Bristol UK
- School of Biochemistry, Biomedical Sciences Building University of Bristol Bristol UK
| |
Collapse
|
69
|
Wehmeyer C, Noé F. Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics. J Chem Phys 2018; 148:241703. [PMID: 29960344 DOI: 10.1063/1.5011399] [Citation(s) in RCA: 170] [Impact Index Per Article: 28.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Inspired by the success of deep learning techniques in the physical and chemical sciences, we apply a modification of an autoencoder type deep neural network to the task of dimension reduction of molecular dynamics data. We can show that our time-lagged autoencoder reliably finds low-dimensional embeddings for high-dimensional feature spaces which capture the slow dynamics of the underlying stochastic processes-beyond the capabilities of linear dimension reduction techniques.
Collapse
Affiliation(s)
- Christoph Wehmeyer
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
70
|
Dibak M, Del Razo MJ, De Sancho D, Schütte C, Noé F. MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations. J Chem Phys 2018; 148:214107. [PMID: 29884049 DOI: 10.1063/1.5020294] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long time-scale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large length scales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time and length scales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step toward MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B ⇌ C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported.
Collapse
Affiliation(s)
- Manuel Dibak
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Mauricio J Del Razo
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - David De Sancho
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), and Donostia International Physics Center (DIPC), P.K. 1072, 20080 Donostia, Euskadi, Spain
| | - Christof Schütte
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
71
|
Zeng X, Li ZW, Zheng X, Zhu L, Sun ZY, Lu ZY, Huang X. Improving the productivity of monodisperse polyhedral cages by the rational design of kinetic self-assembly pathways. Phys Chem Chem Phys 2018; 20:10030-10037. [PMID: 29620122 DOI: 10.1039/c8cp00522b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Hollow polyhedral cages hold great potential for application in nanotechnological and biomedical fields. Understanding the formation mechanism of these self-assembled structures could provide guidance for the rational design of the desired polyhedral cages. Here, by constructing kinetic network models from extensive coarse-grained molecular dynamics simulations, we elucidated the formation mechanism of the dodecahedral cage, which is formed by the self-assembly of patchy particles. We found that the dodecahedral cage is formed through increasing the aggregate size followed by structure rearrangement. Based on this mechanistic understanding, we improved the productivity of the dodecahedral cage through the rational design of the patch arrangement of patchy particles, which promotes the structural rearrangement process. Our results demonstrate that it should be a feasible strategy to achieve the rational design of the desired nanostructures via the kinetic analysis. We anticipate that this methodology could be extended to other self-assembly systems for the fabrication of functional nanomaterials.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | | | | | | | | | | | | |
Collapse
|
72
|
Daskalakis V. Protein-protein interactions within photosystem II under photoprotection: the synergy between CP29 minor antenna, subunit S (PsbS) and zeaxanthin at all-atom resolution. Phys Chem Chem Phys 2018; 20:11843-11855. [PMID: 29658553 DOI: 10.1039/c8cp01226a] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The assembly and disassembly of protein complexes within cells are crucial life-sustaining processes. In photosystem II (PSII) of higher plants, there is a delicate yet obscure balance between light harvesting and photo-protection under fluctuating light conditions, that involves protein-protein complexes. Recent breakthroughs in molecular dynamics (MD) simulations are combined with new approaches herein to provide structural and energetic insight into such a complex between the CP29 minor antenna and the PSII subunit S (PsbS). The microscopic model involves extensive sampling of bound and dissociated states at atomic resolution in the presence of photo-protective zeaxanthin (Zea), and reveals well defined protein-protein cross-sections. The complex is placed within PSII, and macroscopic connections are emerging (PsbS-CP29-CP24-CP47) along the energy transfer pathways from the antenna to the PSII core. These connections explain macroscopic observations in the literature, while the previously obscured atomic scale details are now revealed. The implications of these findings are discussed in the context of the Non-Photochemical Quenching (NPQ) of chlorophyll fluorescence, the down-regulatory mechanism of photosynthesis, that enables the protection of PSII against excess excitation load. Zea is found at the PsbS-CP29 cross-section and a pH-dependent equilibrium between PsbS dimer/monomers and the PsbS-CP29 dissociation/association is identified as the target for engineering tolerant plants with increased crop and biomass yields. Finally, the new MD based approaches can be used to probe protein-protein interactions in general, and the PSII structure provided can initiate large scale molecular simulations of the photosynthetic apparatus, under NPQ conditions.
Collapse
Affiliation(s)
- Vangelis Daskalakis
- Department of Environmental Science and Technology, Cyprus University of Technology (CUT), 30 Archbishop Kyprianou Str., 3603, Limassol, Cyprus.
| |
Collapse
|
73
|
Banda-Vázquez J, Shanmugaratnam S, Rodríguez-Sotres R, Torres-Larios A, Höcker B, Sosa-Peinado A. Redesign of LAOBP to bind novel l-amino acid ligands. Protein Sci 2018. [PMID: 29524280 DOI: 10.1002/pro.3403] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Computational protein design is still a challenge for advancing structure-function relationships. While recent advances in this field are promising, more information for genuine predictions is needed. Here, we discuss different approaches applied to install novel glutamine (Gln) binding into the Lysine/Arginine/Ornithine binding protein (LAOBP) from Salmonella typhimurium. We studied the ligand binding behavior of two mutants: a binding pocket grafting design based on a structural superposition of LAOBP to the Gln binding protein QBP from Escherichia coli and a design based on statistical coupled positions. The latter showed the ability to bind Gln even though the protein was not very stable. Comparison of both approaches highlighted a nonconservative shared point mutation between LAOBP_graft and LAOBP_sca. This context dependent L117K mutation in LAOBP turned out to be sufficient for introducing Gln binding, as confirmed by different experimental techniques. Moreover, the crystal structure of LAOBP_L117K in complex with its ligand is reported.
Collapse
Affiliation(s)
| | - Sooruban Shanmugaratnam
- Max Planck Institute for Developmental Biology, Tübingen, Germany.,Universität Bayreuth, Bayreuth, Germany
| | | | | | - Birte Höcker
- Max Planck Institute for Developmental Biology, Tübingen, Germany.,Universität Bayreuth, Bayreuth, Germany
| | | |
Collapse
|
74
|
Zeng X, Zhu L, Zheng X, Cecchini M, Huang X. Harnessing complexity in molecular self-assembly using computer simulations. Phys Chem Chem Phys 2018; 20:6767-6776. [PMID: 29479585 DOI: 10.1039/c7cp06181a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In molecular self-assembly, hundreds of thousands of freely-diffusing molecules associate to form ordered and functional architectures in the absence of an actuator. This intriguing phenomenon plays a critical role in biology and has become a powerful tool for the fabrication of advanced nanomaterials. Due to the limited spatial and temporal resolutions of current experimental techniques, computer simulations offer a complementary strategy to explore self-assembly with atomic resolution. Here, we review recent computational studies focusing on both thermodynamic and kinetic aspects. As we shall see, thermodynamic approaches based on modeling and statistical mechanics offer initial guidelines to design nanostructures with modest computational effort. Computationally more intensive analyses based on molecular dynamics simulations and kinetic network models (KNMs) reach beyond it, opening the door to the rational design of self-assembly pathways. Current limitations of these methodologies are discussed. We anticipate that the synergistic use of thermodynamic and kinetic analyses based on computer simulations will provide an important contribution to the de novo design of self-assembly.
Collapse
Affiliation(s)
- Xiangze Zeng
- Department of Chemistry, Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong.
| | | | | | | | | |
Collapse
|
75
|
McCluskey K, Carlos Penedo J. An integrated perspective on RNA aptamer ligand-recognition models: clearing muddy waters. Phys Chem Chem Phys 2018; 19:6921-6932. [PMID: 28225108 DOI: 10.1039/c6cp08798a] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Riboswitches are short RNA motifs that sensitively and selectively bind cognate ligands to modulate gene expression. Like protein receptor-ligand pairs, their binding dynamics are traditionally categorized as following one of two paradigmatic mechanisms: conformational selection and induced fit. In conformational selection, ligand binding stabilizes a particular state already present in the receptor's dynamic ensemble. In induced fit, ligand-receptor interactions enable the system to overcome the energetic barrier into a previously inaccessible state. In this article, we question whether a polarized division of RNA binding mechanisms truly meets the conceptual needs of the field. We will review the history behind this classification of RNA-ligand interactions, and the way induced fit in particular has been rehabilitated by single-molecule studies of RNA aptamers. We will highlight several recent results from single-molecule experimental studies of riboswitches that reveal gaps or even contradictions between common definitions of the two terms, and we will conclude by proposing a more robust framework that considers the range of RNA behaviors unveiled in recent years as a reality to be described, rather than an increasingly unwieldy set of exceptions to the traditional models.
Collapse
Affiliation(s)
- K McCluskey
- Department of Physics and Astronomy, University of St. Andrews, St. Andrews, KY16 9SS, UK.
| | - J Carlos Penedo
- Department of Physics and Astronomy, University of St. Andrews, St. Andrews, KY16 9SS, UK. and Biomolecular Sciences Research Complex, University of St. Andrews, St. Andrews, KY16 9SS, UK.
| |
Collapse
|
76
|
Zhu L, Sheong FK, Zeng X, Huang X. Elucidation of the conformational dynamics of multi-body systems by construction of Markov state models. Phys Chem Chem Phys 2018; 18:30228-30235. [PMID: 27314275 DOI: 10.1039/c6cp02545e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Constructing Markov State Models (MSMs) based on short molecular dynamics simulations is a powerful computational technique to complement experiments in predicting long-time kinetics of biomolecular processes at atomic resolution. Even though the MSM approach has been widely applied to study one-body processes such as protein folding and enzyme conformational changes, the majority of biological processes, e.g. protein-ligand recognition, signal transduction, and protein aggregation, essentially involve multiple entities. Here we review the attempts at constructing MSMs for multi-body systems, point out the challenges therein and discuss recent algorithmic progresses that alleviate these challenges. In particular, we describe an automatic kinetics based partitioning method that achieves optimal definition of the conformational states in a multi-body system, and discuss a novel maximum-likelihood approach that efficiently estimates the slow uphill kinetics utilizing pre-computed equilibrium populations of all states. We expect that these new algorithms and their combinations may boost investigations of important multi-body biological processes via the efficient construction of MSMs.
Collapse
Affiliation(s)
- Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China.
| | - Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China. and Centre of Systems Biology and Human Health, School of Science and Institute for Advance Study, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| |
Collapse
|
77
|
Unarta IC, Zhu L, Tse CKM, Cheung PPH, Yu J, Huang X. Molecular mechanisms of RNA polymerase II transcription elongation elucidated by kinetic network models. Curr Opin Struct Biol 2018; 49:54-62. [PMID: 29414512 DOI: 10.1016/j.sbi.2018.01.002] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2017] [Revised: 12/22/2017] [Accepted: 01/02/2018] [Indexed: 12/30/2022]
Abstract
Transcription elongation cycle (TEC) of RNA polymerase II (Pol II) is a process of adding a nucleoside triphosphate to the growing messenger RNA chain. Due to the long timescale events in Pol II TEC, an advanced computational technique, such as Markov State Model (MSM), is needed to provide atomistic mechanism and reaction rates. The combination of MSM and experimental results can be used to build a kinetic network model (KNM) of the whole TEC. This review provides a brief protocol to build MSM and KNM of the whole TEC, along with the latest findings of MSM and other computational studies of Pol II TEC. Lastly, we offer a perspective on potentially using a sequence dependent KNM to predict genome-wide transcription error.
Collapse
Affiliation(s)
- Ilona Christy Unarta
- Bioengineering Graduate Program, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, Hong Kong
| | - Lizhe Zhu
- Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, Hong Kong; Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Carmen Ka Man Tse
- Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, Hong Kong; Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Peter Pak-Hang Cheung
- Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, Hong Kong
| | - Jin Yu
- Beijing Computational Science Research Center, Beijing 100084, China
| | - Xuhui Huang
- Center of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & Reconstruction, Hong Kong; Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong; HKUST-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China.
| |
Collapse
|
78
|
Zhou Y, Hussain M, Kuang G, Zhang J, Tu Y. Mechanistic insights into peptide and ligand binding of the ATAD2-bromodomain via atomistic simulations disclosing a role of induced fit and conformational selection. Phys Chem Chem Phys 2018; 20:23222-23232. [DOI: 10.1039/c8cp03860k] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Atomistic simulations of the ATAD2-bromodomain disclose a role of induced fit and conformational selection upon ligand and peptide binding.
Collapse
Affiliation(s)
- Yang Zhou
- Department of Theoretical Chemistry and Biology
- KTH Royal Institute of Technology
- AlbaNova University Center
- Stockholm
- Sweden
| | - Muzammal Hussain
- Guangdong Provincial Key Laboratory of Biocomputing
- Institute of Chemical Biology
- Guangzhou Institutes of Biomedicine and Health
- Chinese Academy of Sciences
- Guangzhou 510530
| | - Guanglin Kuang
- Department of Theoretical Chemistry and Biology
- KTH Royal Institute of Technology
- AlbaNova University Center
- Stockholm
- Sweden
| | - Jiancun Zhang
- Guangdong Provincial Key Laboratory of Biocomputing
- Institute of Chemical Biology
- Guangzhou Institutes of Biomedicine and Health
- Chinese Academy of Sciences
- Guangzhou 510530
| | - Yaoquan Tu
- Department of Theoretical Chemistry and Biology
- KTH Royal Institute of Technology
- AlbaNova University Center
- Stockholm
- Sweden
| |
Collapse
|
79
|
Cortes-Hernandez P, Domínguez-Ramírez L. Role of cis-trans proline isomerization in the function of pathogenic enterobacterial Periplasmic Binding Proteins. PLoS One 2017; 12:e0188935. [PMID: 29190818 PMCID: PMC5708682 DOI: 10.1371/journal.pone.0188935] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2017] [Accepted: 11/15/2017] [Indexed: 12/31/2022] Open
Abstract
Periplasmic Binding Proteins (PBPs) trap nutrients for their internalization into bacteria by ABC transporters. Ligand binding triggers PBP closure by bringing its two domains together like a Venus flytrap. The atomic determinants that control PBP opening and closure for nutrient capture and release are not known, although it is proposed that opening and ligand release occur while in contact with the ABC transporter for concurrent substrate translocation. In this paper we evaluated the effect of the isomerization of a conserved proline, located near the binding site, on the propensity of PBPs to open and close. ArgT/LAO from Salmonella typhimurium and HisJ from Escherichia coli were studied through molecular mechanics at two different temperatures: 300 and 323 K. Eight microseconds were simulated per protein to analyze protein opening and closure in the absence of the ABC transporter. We show that when the studied proline is in trans, closed empty LAO and HisJ can open. In contrast, with the proline in cis, opening transitions were much less frequent and characterized by smaller changes. The proline in trans also renders the open trap prone to close over a ligand. Our data suggest that the isomerization of this conserved proline modulates the PBP mechanism: the proline in trans allows the exploration of conformational space to produce trap opening and closure, while in cis it restricts PBP movement and could limit ligand release until in productive contact with the ABC transporter. This is the first time that a proline isomerization has been related to the control of a large conformational change like the PBP flytrap mechanism.
Collapse
Affiliation(s)
- Paulina Cortes-Hernandez
- Centro de Investigacion Biomedica de Oriente (CIBIOR), Instituto Mexicano del Seguro Social (IMSS), Metepec, Puebla, Mexico
| | - Lenin Domínguez-Ramírez
- Chemical and Biological Sciences, Universidad de las Américas Puebla (UDLAP), Cholula, Puebla, Mexico
- * E-mail:
| |
Collapse
|
80
|
Ma H, Li A, Gao K. Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme. ACS OMEGA 2017; 2:8414-8420. [PMID: 30023582 PMCID: PMC6045336 DOI: 10.1021/acsomega.7b01414] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/22/2017] [Accepted: 11/15/2017] [Indexed: 05/30/2023]
Abstract
Human carbonic anhydrase II (HCA II) is an enzyme that catalyzes the reversible hydration of CO2 into bicarbonate (HCO3-) and a proton (H+) as well as other reactions at an extremely high rate. This enzyme plays fundamental roles in human physiology/pathology, such as controlling the pH level in cells and so on. However, the binding mechanism between apo-HCA II and CO2 or other ligands as well as related conformational changes remains poorly understood, and atomic investigation into it could promote our understanding of related internal physiological/pathological mechanisms. In this study, long-time atomic molecular dynamics simulations as well as the clustering and free-energy analysis were performed to reveal the dynamics of apo-HCA II as well as the mechanism upon ligand binding. Our simulations indicate that the crystallographic B-factors considerably underestimate the loop dynamics: multiple conformations can be adopted by loops 1 and 2, especially for loop 1 because loop 1 is one side of the binding pocket, and its left-to-right movement can compress or extend the binding pocket, leading to one inactive (closed) state, three intermediate (semiopen) states, and one active (open) state; CO2 cannot get into the binding pocket of the inactive state but can get into those of intermediate and active states. The coexistence of multiple conformational states proposes a possible conformational selection model for the binding mechanism between apo-HCA II and CO2 or other ligands, revising our previous view of its functional mechanism of conformational change upon ligand binding and offering valuable structural insights into the workings of HCA II.
Collapse
Affiliation(s)
- Huishu Ma
- Institute of Biophysics and Department
of Physics, Central China Normal University, Wuhan 430079, P. R. China
| | - Anbang Li
- Institute of Biophysics and Department
of Physics, Central China Normal University, Wuhan 430079, P. R. China
| | - Kaifu Gao
- Institute of Biophysics and Department
of Physics, Central China Normal University, Wuhan 430079, P. R. China
| |
Collapse
|
81
|
Paul F, Wehmeyer C, Abualrous ET, Wu H, Crabtree MD, Schöneberg J, Clarke J, Freund C, Weikl TR, Noé F. Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations. Nat Commun 2017; 8:1095. [PMID: 29062047 PMCID: PMC5653669 DOI: 10.1038/s41467-017-01163-6] [Citation(s) in RCA: 102] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2017] [Accepted: 08/22/2017] [Indexed: 11/10/2022] Open
Abstract
Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25-109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.
Collapse
Affiliation(s)
- Fabian Paul
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA
- Max Planck Institute of Colloids and Interfaces, Department of Theory and Bio-Systems, 14476, Potsdam, Germany
| | - Christoph Wehmeyer
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA
| | - Esam T Abualrous
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA
| | - Hao Wu
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA
| | - Michael D Crabtree
- Department of Chemistry, University of Cambridge, Cambridge, CB2 1EW, UK
| | - Johannes Schöneberg
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA
| | - Jane Clarke
- Department of Chemistry, University of Cambridge, Cambridge, CB2 1EW, UK
| | - Christian Freund
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Thielallee 63, 14195, Berlin, Germany
| | - Thomas R Weikl
- Max Planck Institute of Colloids and Interfaces, Department of Theory and Bio-Systems, 14476, Potsdam, Germany
| | - Frank Noé
- Department of Molecular and Cell Biology and California Institute for Quantitative Biosciences, University of California, Berkeley, CA, 94720, USA.
| |
Collapse
|
82
|
Wang W, Cao S, Zhu L, Huang X. Constructing Markov State Models to elucidate the functional conformational changes of complex biomolecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1343] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Wei Wang
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Siqin Cao
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Lizhe Zhu
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
| | - Xuhui Huang
- Department of ChemistryThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Center of Systems Biology and Human HealthThe Hong Kong University of Science and Technology Kowloon Hong Kong
- Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration & ReconstructionThe Hong Kong University of Science and Technology Kowloon Hong Kong
- HKUST‐Shenzhen Research Institute Shenzhen China
| |
Collapse
|
83
|
Fonin AV, Golikova AD, Zvereva IA, D'Auria S, Staiano M, Uversky VN, Kuznetsova IM, Turoverov KK. Osmolyte-Like Stabilizing Effects of Low GdnHCl Concentrations on d-Glucose/d-Galactose-Binding Protein. Int J Mol Sci 2017; 18:E2008. [PMID: 28925982 PMCID: PMC5618657 DOI: 10.3390/ijms18092008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 09/14/2017] [Accepted: 09/16/2017] [Indexed: 11/16/2022] Open
Abstract
The ability of d-glucose/d-galactose-binding protein (GGBP) to reversibly interact with its ligands, glucose and galactose, makes this protein an attractive candidate for sensing elements of glucose biosensors. This potential is largely responsible for attracting researchers to study the conformational properties of this protein. Previously, we showed that an increase in the fluorescence intensity of the fluorescent dye 6-bromoacetyl-2-dimetylaminonaphtalene (BADAN) is linked to the holo-form of the GGBP/H152C mutant in solutions containing sub-denaturing concentrations of guanidine hydrochloride (GdnHCl). It was hypothesized that low GdnHCl concentrations might lead to compaction of the protein, thereby facilitating ligand binding. In this work, we utilize BADAN fluorescence spectroscopy, intrinsic protein UV fluorescence spectroscopy, and isothermal titration calorimetry (ITC) to show that the sub-denaturing GdnHCl concentrations possess osmolyte-like stabilizing effects on the structural dynamics, conformational stability, and functional activity of GGBP/H152C and the wild type of this protein (wtGGBP). Our data are consistent with the model where low GdnHCl concentrations promote a shift in the dynamic distribution of the protein molecules toward a conformational ensemble enriched in molecules with a tighter structure and a more closed conformation. This promotes the increase in the configurational complementarity between the protein and glucose molecules that leads to the increase in glucose affinity in both GGBP/H152C and wtGGBP.
Collapse
Affiliation(s)
- Alexander V Fonin
- Institute of Cytology of the Russian Academy of Sciences, Laboratory of Structural Dynamics, Stability and Folding of Proteins, Tikhoretsky av. 4, 194064 St. Petersburg, Russia.
| | - Alexandra D Golikova
- Saint Petersburg State University, Universitetskaya nab. 7/9, 199034 St. Petersburg, Russia.
| | - Irina A Zvereva
- Saint Petersburg State University, Universitetskaya nab. 7/9, 199034 St. Petersburg, Russia.
| | - Sabato D'Auria
- CNR, Institute of Food Science, via Roma 64, 83100 Avellino, Italy.
| | - Maria Staiano
- CNR, Institute of Food Science, via Roma 64, 83100 Avellino, Italy.
| | - Vladimir N Uversky
- Department of Molecular Medicine and Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, 12901 Bruce B. Downs Blvd., Tampa, FL 33612, USA.
| | - Irina M Kuznetsova
- Institute of Cytology of the Russian Academy of Sciences, Laboratory of Structural Dynamics, Stability and Folding of Proteins, Tikhoretsky av. 4, 194064 St. Petersburg, Russia.
| | - Konstantin K Turoverov
- Institute of Cytology of the Russian Academy of Sciences, Laboratory of Structural Dynamics, Stability and Folding of Proteins, Tikhoretsky av. 4, 194064 St. Petersburg, Russia.
- Department of Biophysics, Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya av. 29, 195251 St. Petersburg, Russia.
| |
Collapse
|
84
|
Abramyan AM, Stolzenberg S, Li Z, Loland CJ, Noé F, Shi L. The Isomeric Preference of an Atypical Dopamine Transporter Inhibitor Contributes to Its Selection of the Transporter Conformation. ACS Chem Neurosci 2017; 8:1735-1746. [PMID: 28441487 DOI: 10.1021/acschemneuro.7b00094] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Cocaine, a widely abused psychostimulant, inhibits the dopamine transporter (DAT) by trapping the protein in an outward-open conformation, whereas atypical DAT inhibitors such as benztropine have low abuse liability and prefer less outward-open conformations. Here, we use a spectrum of computational modeling and simulation approaches to obtain the underlying molecular mechanism in atomistic detail. Interestingly, our quantum mechanical calculations and molecular dynamics (MD) simulations suggest that a benztropine derivative JHW007 prefers a different stereoisomeric conformation of tropane in binding to DAT compared to that of a cocaine derivative, CFT. To further investigate the different inhibition mechanisms of DAT, we carried out MD simulations in combination with Markov state modeling analysis of wild-type and Y156F DAT in the absence of any ligand or the presence of CFT or JHW007. Our results indicate that the Y156F mutation and CFT shift the conformational equilibrium toward an outward-open conformation, whereas JHW007 prefers an inward-occluded conformation. Our findings reveal the mechanistic details of DAT inhibition by JHW007 at the atomistic level, which provide clues for rational design of atypical inhibitors.
Collapse
Affiliation(s)
- Ara M. Abramyan
- Computational
Chemistry and Molecular Biophysics Unit, Molecular Targets and Medications
Discovery Branch, NIH/NIDA/IRP, Baltimore, Maryland 21224, United States
| | - Sebastian Stolzenberg
- Computational
Molecular Biology group, Institute for Mathematics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Zheng Li
- Department
of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| | - Claus J. Loland
- Molecular
Neuropharmacology Group, Department of Neuroscience and Pharmacology,
The Faculty of Health Sciences, The Panum Institute, University of Copenhagen, DK-2200 Copenhagen, Denmark
| | - Frank Noé
- Computational
Molecular Biology group, Institute for Mathematics, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| | - Lei Shi
- Computational
Chemistry and Molecular Biophysics Unit, Molecular Targets and Medications
Discovery Branch, NIH/NIDA/IRP, Baltimore, Maryland 21224, United States
- Department
of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| |
Collapse
|
85
|
Harada R, Shigeta Y. An assessment of optimal time scale of conformational resampling for parallel cascade selection molecular dynamics. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1362696] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, Tsukuba, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, Tsukuba, Japan
| |
Collapse
|
86
|
Meng L, Sheong FK, Zeng X, Zhu L, Huang X. Path lumping: An efficient algorithm to identify metastable path channels for conformational dynamics of multi-body systems. J Chem Phys 2017; 147:044112. [DOI: 10.1063/1.4995558] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Luming Meng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Fu Kit Sheong
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xiangze Zeng
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Lizhe Zhu
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Center of Systems Biology and Human Health, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
- HKUST-Shenzhen Research Institute, Hi-Tech Park, Nanshan, Shenzhen 518057, China
| |
Collapse
|
87
|
Verma R, Mitchell-Koch K. In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function. Catalysts 2017; 7:212. [PMID: 30464857 PMCID: PMC6241538 DOI: 10.3390/catal7070212] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.
Collapse
Affiliation(s)
- Rajni Verma
- Department of Chemistry, McKinley Hall, Wichita State University, 1845 Fairmount, Wichita, KS 67260-0051, USA
| | - Katie Mitchell-Koch
- Department of Chemistry, McKinley Hall, Wichita State University, 1845 Fairmount, Wichita, KS 67260-0051, USA
| |
Collapse
|
88
|
Ge Y, Kier BL, Andersen NH, Voelz VA. Computational and Experimental Evaluation of Designed β-Cap Hairpins Using Molecular Simulations and Kinetic Network Models. J Chem Inf Model 2017; 57:1609-1620. [PMID: 28614661 DOI: 10.1021/acs.jcim.7b00132] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Molecular simulation has been used to model the detailed folding properties of peptides, yet prospective computational peptide design by such approaches remains challenging and nontrivial. To test the accuracy of simulation-based hairpin design, we characterized the folding properties of a series of so-called β-cap hairpin peptides designed to mimic a conserved hairpin of LapD, a bacterial intracellular signaling protein, both experimentally by NMR spectroscopy and computationally by implicit-solvent replica-exchange molecular dynamics using three different AMBER force fields (ff96, ff99sb-ildn, and ff99sb-ildn-NMR). A unique challenge presented by these designs is the presence of both a terminal Trp-Trp capping motif and a conserved GWxQ motif in the hairpin turn required for binding to LapG. Consistent with previous studies, we found AMBER ff96 to be the most accurate when used with the OBC GBSA implicit solvent model, despite its known bias toward β-sheet conformations when used in explicit-solvent simulations. To gain microscopic insight into the folding landscape of the hairpin designs, we additionally performed parallel simulations on the Folding@home distributed computing platform using AMBER ff99sb-ildn-NMR with TIP3P explicit solvent. Markov state models (MSMs) built from trajectory data reveal a number of non-native interactions between Trp and other amino acid side chains, creating potential problems in achieving well-folded hairpin structures in solution.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Chemistry, Temple University , Philadelphia, Pennsylvania 19122, United States
| | - Brandon L Kier
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Niels H Andersen
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Vincent A Voelz
- Department of Chemistry, Temple University , Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
89
|
Complete protein-protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modelling. Nat Chem 2017; 9:1005-1011. [PMID: 28937668 DOI: 10.1038/nchem.2785] [Citation(s) in RCA: 237] [Impact Index Per Article: 33.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2016] [Accepted: 04/21/2017] [Indexed: 12/23/2022]
Abstract
Protein-protein association is fundamental to many life processes. However, a microscopic model describing the structures and kinetics during association and dissociation is lacking on account of the long lifetimes of associated states, which have prevented efficient sampling by direct molecular dynamics (MD) simulations. Here we demonstrate protein-protein association and dissociation in atomistic resolution for the ribonuclease barnase and its inhibitor barstar by combining adaptive high-throughput MD simulations and hidden Markov modelling. The model reveals experimentally consistent intermediate structures, energetics and kinetics on timescales from microseconds to hours. A variety of flexibly attached intermediates and misbound states funnel down to a transition state and a native basin consisting of the loosely bound near-native state and the tightly bound crystallographic state. These results offer a deeper level of insight into macromolecular recognition and our approach opens the door for understanding and manipulating a wide range of macromolecular association processes.
Collapse
|
90
|
Thayer KM, Lakhani B, Beveridge DL. Molecular Dynamics-Markov State Model of Protein Ligand Binding and Allostery in CRIB-PDZ: Conformational Selection and Induced Fit. J Phys Chem B 2017; 121:5509-5514. [PMID: 28489401 DOI: 10.1021/acs.jpcb.7b02083] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Conformational selection and induced fit are well-known contributors to ligand binding and allosteric effects in proteins. Molecular dynamics (MD) simulations now enable the theoretical study of protein-ligand binding in terms of ensembles of interconverting microstates and the population shifts characteristic of "dynamical allostery." Here we investigate protein-ligand binding and allostery based on a Markov state model (MSM) with states and rates obtained from all-atom MD simulations. As an exemplary case, we consider the single domain protein par-6 PDZ with and without ligand and allosteric effector. This is one of the smallest proteins in which allostery has been experimentally observed. In spite of the increased complexity intrinsic to a statistical ensemble perspective, we find that conformational selection and induced fit mechanisms can be readily identified in the analysis. In the nonallosteric pathway, MD-MSM shows that PDZ binds ligand via conformational selection. However, the allosteric pathway requires an activation step that involves a conformational change induced by the allosteric effector Cdc42. Once in the allosterically activated state, we find that ligand binding can proceed by conformational selection. Our MD-MSM model predicts that allostery in this and possibly other systems involves both induced fit and conformational selection, not just one or the other.
Collapse
Affiliation(s)
- Kelly M Thayer
- Departments of Chemistry, ‡Molecular Biology & Biochemistry, §Molecular Biophysics Program, and ∥Department of Computer Science, Wesleyan University , Middletown, Connecticut 06459, United States of America
| | - Bharat Lakhani
- Departments of Chemistry, ‡Molecular Biology & Biochemistry, §Molecular Biophysics Program, and ∥Department of Computer Science, Wesleyan University , Middletown, Connecticut 06459, United States of America
| | - David L Beveridge
- Departments of Chemistry, ‡Molecular Biology & Biochemistry, §Molecular Biophysics Program, and ∥Department of Computer Science, Wesleyan University , Middletown, Connecticut 06459, United States of America
| |
Collapse
|
91
|
Singh S, Bowman GR. Quantifying Allosteric Communication via Both Concerted Structural Changes and Conformational Disorder with CARDS. J Chem Theory Comput 2017; 13:1509-1517. [PMID: 28282132 PMCID: PMC5934993 DOI: 10.1021/acs.jctc.6b01181] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Allosteric (i.e., long-range) communication within proteins is crucial for many biological processes, such as the activation of signaling cascades in response to specific stimuli. However, the physical basis for this communication remains unclear. Existing computational methods for identifying allostery focus on the role of concerted structural changes, but recent experimental work demonstrates that disorder is also an important factor. Here, we introduce the Correlation of All Rotameric and Dynamical States (CARDS) framework for quantifying correlations between both the structure and disorder of different regions of a protein. To quantify disorder, we draw inspiration from methods for quantifying "dynamic heterogeneity" from chemical physics to classify segments of a dihedral's time evolution as being in either ordered or disordered regimes. To demonstrate the utility of this approach, we apply CARDS to the Catabolite Activator Protein (CAP), a transcriptional activator that is regulated by Cyclic Adenosine MonoPhosphate (cAMP) binding. We find that CARDS captures allosteric communication between the two cAMP-Binding Domains (CBDs). Importantly, CARDS reveals that this coupling is dominated by disorder-mediated correlations, consistent with NMR experiments that establish allosteric coupling between the CBDs occurs without a concerted structural change. CARDS also recapitulates an enhanced role for disorder in the communication between the DNA-Binding Domains (DBDs) and CBDs in the S62F variant of CAP. Finally, we demonstrate that using CARDS to find communication hotspots identifies regions of CAP that are in allosteric communication without foreknowledge of their identities. Therefore, we expect CARDS to be of great utility for both understanding and predicting allostery.
Collapse
Affiliation(s)
- Sukrit Singh
- Department of Biochemistry and Molecular Biophysics, Washington University in St. Louis, St. Louis, MO
| | - Gregory R. Bowman
- Department of Biochemistry and Molecular Biophysics, Washington University in St. Louis, St. Louis, MO
- Center for Biological Systems Engineering, Washington University in St. Louis, St. Louis, MO
| |
Collapse
|
92
|
Diamantis P, Unke OT, Meuwly M. Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N. PLoS Comput Biol 2017; 13:e1005450. [PMID: 28358830 PMCID: PMC5391117 DOI: 10.1371/journal.pcbi.1005450] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2016] [Revised: 04/13/2017] [Accepted: 03/13/2017] [Indexed: 11/18/2022] Open
Abstract
In heme proteins, the efficient transport of ligands such as NO or O2 to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it. Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.
Collapse
Affiliation(s)
| | - Oliver T. Unke
- Department of Chemistry, University of Basel, Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Basel, Switzerland
- * E-mail:
| |
Collapse
|
93
|
Gao K, Zhao Y. A Network of Conformational Transitions in the Apo Form of NDM-1 Enzyme Revealed by MD Simulation and a Markov State Model. J Phys Chem B 2017; 121:2952-2960. [PMID: 28319394 DOI: 10.1021/acs.jpcb.7b00062] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
New Delhi metallo-β-lactamase-1 (NDM-1) is a novel β-lactamase enzyme that confers enteric bacteria with nearly complete resistance to all β-lactam antibiotics, so it raises a formidable and global threat to human health. However, the binding mechanism between apo-NDM-1 and antibiotics as well as related conformational changes remains poorly understood, which largely hinders the overcoming of its antibiotic resistance. In our study, long-time conventional molecular dynamics simulation and Markov state models were applied to reveal both the dynamical and conformational landscape of apo-NDM-1: the MD simulation demonstrates that loop L3, which is responsible for antibiotic binding, is the most flexible and undergoes dramatic conformational changes; moreover, the Markov state model built from the simulation maps four metastable states including open, semiopen, and closed conformations of loop L3 as well as frequent transitions between the states. Our findings propose a possible conformational selection model for the binding mechanism between apo-NDM-1 and antibiotics, which facilitates the design of novel inhibitors and antibiotics.
Collapse
Affiliation(s)
- Kaifu Gao
- Institute of Biophysics and Department of Physics, Central China Normal University , Wuhan 430079, P. R. China
| | - Yunjie Zhao
- Institute of Biophysics and Department of Physics, Central China Normal University , Wuhan 430079, P. R. China
| |
Collapse
|
94
|
Bernetti M, Cavalli A, Mollica L. Protein-ligand (un)binding kinetics as a new paradigm for drug discovery at the crossroad between experiments and modelling. MEDCHEMCOMM 2017; 8:534-550. [PMID: 30108770 PMCID: PMC6072069 DOI: 10.1039/c6md00581k] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Accepted: 01/25/2017] [Indexed: 12/14/2022]
Abstract
In the last three decades, protein and nucleic acid structure determination and comprehension of the mechanisms, leading to their physiological and pathological functions, have become a cornerstone of biomedical sciences. A deep understanding of the principles governing the fates of cells and tissue at the molecular level has been gained over the years, offering a solid basis for the rational design of drugs aimed at the pharmacological treatment of numerous diseases. Historically, affinity indicators (i.e. Kd and IC50/EC50) have been assumed to be valid indicators of the in vivo efficacy of a drug. However, recent studies pointed out that the kinetics of the drug-receptor binding process could be as important or even more important than affinity in determining the drug efficacy. This eventually led to a growing interest in the characterisation and prediction of the rate constants of protein-ligand association and dissociation. For instance, a drug with a longer residence time can kinetically select a given receptor over another, even if the affinity for both receptors is comparable, thus increasing its therapeutic index. Therefore, understanding the molecular features underlying binding and unbinding processes is of central interest towards the rational control of drug binding kinetics. In this review, we report the theoretical framework behind protein-ligand association and highlight the latest advances in the experimental and computational approaches exploited to investigate the binding kinetics.
Collapse
Affiliation(s)
- M Bernetti
- Department of Pharmacy and Biotechnology , University of Bologna , via Belmeloro 6 , 40126 Bologna , Italy
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| | - A Cavalli
- Department of Pharmacy and Biotechnology , University of Bologna , via Belmeloro 6 , 40126 Bologna , Italy
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| | - L Mollica
- CompuNet , Istituto Italiano di Tecnologia , via Morego 30 , 16163 Genova , Italy .
| |
Collapse
|
95
|
Moritsugu K, Terada T, Kidera A. Free-Energy Landscape of Protein–Ligand Interactions Coupled with Protein Structural Changes. J Phys Chem B 2017; 121:731-740. [DOI: 10.1021/acs.jpcb.6b11696] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kei Moritsugu
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehirocho, Tsurumi, Yokohama 230-0045, Japan
| | - Tohru Terada
- Graduate
School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo, Tokyo 113-8657, Japan
| | - Akinori Kidera
- Graduate
School of Medical Life Science, Yokohama City University, 1-7-29
Suehirocho, Tsurumi, Yokohama 230-0045, Japan
| |
Collapse
|
96
|
Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches. Methods Mol Biol 2016. [PMID: 27924488 DOI: 10.1007/978-1-4939-6563-2_18] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.
Collapse
|
97
|
Kar RK, Kharerin H, Padinhateeri R, Bhat PJ. Multiple Conformations of Gal3 Protein Drive the Galactose-Induced Allosteric Activation of the GAL Genetic Switch of Saccharomyces cerevisiae. J Mol Biol 2016; 429:158-176. [PMID: 27913116 DOI: 10.1016/j.jmb.2016.11.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2016] [Revised: 11/04/2016] [Accepted: 11/04/2016] [Indexed: 01/01/2023]
Abstract
Gal3p is an allosteric monomeric protein that activates the GAL genetic switch of Saccharomyces cerevisiae in response to galactose. Expression of constitutive mutant of Gal3p or overexpression of wild-type Gal3p activates the GAL switch in the absence of galactose. These data suggest that Gal3p exists as an ensemble of active and inactive conformations. Structural data have indicated that Gal3p exists in open (inactive) and closed (active) conformations. However, a mutant of Gal3p that predominantly exists in inactive conformation and is yet capable of responding to galactose has not been isolated. To understand the mechanism of allosteric transition, we have isolated a triple mutant of Gal3p with V273I, T404A, and N450D substitutions, which, upon overexpression, fails to activate the GAL switch on its own but activates the switch in response to galactose. Overexpression of Gal3p mutants with single or double mutations in any of the three combinations failed to exhibit the behavior of the triple mutant. Molecular dynamics analysis of the wild-type and the triple mutant along with two previously reported constitutive mutants suggests that the wild-type Gal3p may also exist in super-open conformation. Furthermore, our results suggest that the dynamics of residue F237 situated in the hydrophobic pocket located in the hinge region drives the transition between different conformations. Based on this study, we suggest that conformational selection mechanism is the driving force in the allosteric transition of Gal3p, which may have implications in other signaling pathways involving monomeric proteins.
Collapse
Affiliation(s)
- Rajesh Kumar Kar
- Department of Biosciences and Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India
| | - Hungyo Kharerin
- Department of Biosciences and Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India
| | - Ranjith Padinhateeri
- Department of Biosciences and Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India
| | - Paike Jayadeva Bhat
- Department of Biosciences and Bioengineering, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India.
| |
Collapse
|
98
|
Feng Y, Zhang L, Wu S, Liu Z, Gao X, Zhang X, Liu M, Liu J, Huang X, Wang W. Conformational Dynamics of apo-GlnBP Revealed by Experimental and Computational Analysis. Angew Chem Int Ed Engl 2016. [DOI: 10.1002/ange.201606613] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Yitao Feng
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Lu Zhang
- Department of Chemistry; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
| | - Shaowen Wu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Zhijun Liu
- National Center for Protein Science; Shanghai Institute of Biochemistry and Cell Biology; Chinese Academy of Sciences; Shanghai China
| | - Xin Gao
- King Abdullah University of Science and Technology (KAUST); Computational Bioscience Research Center (CBRC); Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division; Thuwal 23955 Saudi Arabia
| | - Xu Zhang
- Key Laboratory of Magnetic and Resonance in Biological Systems; State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics; Centre for Magnetic Resonance; Wuhan Institute of Physics and Mathematics; Chinese Academy of Sciences; Wuhan China
| | - Maili Liu
- Key Laboratory of Magnetic and Resonance in Biological Systems; State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics; Centre for Magnetic Resonance; Wuhan Institute of Physics and Mathematics; Chinese Academy of Sciences; Wuhan China
| | - Jianwei Liu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Xuhui Huang
- Department of Chemistry; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
- Division of Biomedical Engineering; Center of Systems Biology and Human Health; Institute for Advance Study and School of Science; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
| | - Wenning Wang
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| |
Collapse
|
99
|
Feng Y, Zhang L, Wu S, Liu Z, Gao X, Zhang X, Liu M, Liu J, Huang X, Wang W. Conformational Dynamics of apo-GlnBP Revealed by Experimental and Computational Analysis. Angew Chem Int Ed Engl 2016; 55:13990-13994. [PMID: 27730716 DOI: 10.1002/anie.201606613] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Revised: 08/28/2016] [Indexed: 01/01/2023]
Affiliation(s)
- Yitao Feng
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Lu Zhang
- Department of Chemistry; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
| | - Shaowen Wu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Zhijun Liu
- National Center for Protein Science; Shanghai Institute of Biochemistry and Cell Biology; Chinese Academy of Sciences; Shanghai China
| | - Xin Gao
- King Abdullah University of Science and Technology (KAUST); Computational Bioscience Research Center (CBRC); Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division; Thuwal 23955 Saudi Arabia
| | - Xu Zhang
- Key Laboratory of Magnetic and Resonance in Biological Systems; State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics; Centre for Magnetic Resonance; Wuhan Institute of Physics and Mathematics; Chinese Academy of Sciences; Wuhan China
| | - Maili Liu
- Key Laboratory of Magnetic and Resonance in Biological Systems; State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics; Centre for Magnetic Resonance; Wuhan Institute of Physics and Mathematics; Chinese Academy of Sciences; Wuhan China
| | - Jianwei Liu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| | - Xuhui Huang
- Department of Chemistry; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
- Division of Biomedical Engineering; Center of Systems Biology and Human Health; Institute for Advance Study and School of Science; The Hong Kong University of Science and Technology; Clear Water Bay Kowloon Hong Kong
| | - Wenning Wang
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials; Department of Chemistry, and Institutes of Biomedical Sciences; Fudan University; Shanghai P.R. China
| |
Collapse
|
100
|
Hart KM, Ho CMW, Dutta S, Gross ML, Bowman GR. Modelling proteins' hidden conformations to predict antibiotic resistance. Nat Commun 2016; 7:12965. [PMID: 27708258 PMCID: PMC5477488 DOI: 10.1038/ncomms12965] [Citation(s) in RCA: 94] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2016] [Accepted: 08/19/2016] [Indexed: 11/30/2022] Open
Abstract
TEM β-lactamase confers bacteria with resistance to many antibiotics and rapidly evolves activity against new drugs. However, functional changes are not easily explained by differences in crystal structures. We employ Markov state models to identify hidden conformations and explore their role in determining TEM's specificity. We integrate these models with existing drug-design tools to create a new technique, called Boltzmann docking, which better predicts TEM specificity by accounting for conformational heterogeneity. Using our MSMs, we identify hidden states whose populations correlate with activity against cefotaxime. To experimentally detect our predicted hidden states, we use rapid mass spectrometric footprinting and confirm our models' prediction that increased cefotaxime activity correlates with reduced Ω-loop flexibility. Finally, we design novel variants to stabilize the hidden cefotaximase states, and find their populations predict activity against cefotaxime in vitro and in vivo. Therefore, we expect this framework to have numerous applications in drug and protein design.
Collapse
Affiliation(s)
- Kathryn M. Hart
- Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, 660 South Euclid Avenue, St Louis, Missouri 63110, USA
| | - Chris M. W. Ho
- Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, 660 South Euclid Avenue, St Louis, Missouri 63110, USA
| | - Supratik Dutta
- Department of Chemistry, Washington University in St Louis, One Brookings Drive, St Louis, Missouri 63130, USA
| | - Michael L. Gross
- Department of Chemistry, Washington University in St Louis, One Brookings Drive, St Louis, Missouri 63130, USA
| | - Gregory R. Bowman
- Department of Biochemistry & Molecular Biophysics, Washington University School of Medicine, 660 South Euclid Avenue, St Louis, Missouri 63110, USA
- Department of Biomedical Engineering, and Center for Biological Systems Engineering, Washington University in St Louis, One Brookings Drive, St Louis, Missouri 63130, USA
| |
Collapse
|