1
|
Sisk TR, Robustelli P. Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model. Proc Natl Acad Sci U S A 2024; 121:e2313360121. [PMID: 38294935 PMCID: PMC10861926 DOI: 10.1073/pnas.2313360121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 11/22/2023] [Indexed: 02/02/2024] Open
Abstract
A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.
Collapse
Affiliation(s)
- Thomas R. Sisk
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| | - Paul Robustelli
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| |
Collapse
|
2
|
Raddi RM, Voelz VA. Markov State Model of Solvent Features Reveals Water Dynamics in Protein-Peptide Binding. J Phys Chem B 2023; 127:10682-10690. [PMID: 38078851 DOI: 10.1021/acs.jpcb.3c04775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
In this work, we investigate the role of solvent in the binding reaction of the p53 transactivation domain (TAD) peptide to its receptor MDM2. Previously, our group generated 831 μs of explicit-solvent aggregate molecular simulation trajectory data for the MDM2-p53 peptide binding reaction using large-scale distributed computing and subsequently built a Markov State Model (MSM) of the binding reaction (Zhou et al. 2017). Here, we perform a tICA analysis and construct an MSM with similar hyperparameters while using only solvent-based structural features. We find a remarkably similar landscape but accelerated implied timescales for the slowest motions. The solvent shells contributing most to the first tICA eigenvector are those centered on Lys24 and Thr18 of the p53 TAD peptide in the range of 3-6 Å. Important solvent shells were visualized to reveal solvation and desolvation transitions along the peptide-protein binding trajectories. Our results provide a solvent-centric view of the hydrophobic effect in action for a realistic peptide-protein binding scenario.
Collapse
Affiliation(s)
- Robert M Raddi
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Vincent A Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
3
|
Hofmann H. All over or overall - Do we understand allostery? Curr Opin Struct Biol 2023; 83:102724. [PMID: 37898005 DOI: 10.1016/j.sbi.2023.102724] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 09/27/2023] [Accepted: 09/28/2023] [Indexed: 10/30/2023]
Abstract
Allostery is probably the most important concept in the regulation of cellular processes. Models to explain allostery are plenty. Each sheds light on different aspects but their entirety conveys an ambiguous feeling of comprehension and disappointment. Here, I discuss the most popular allostery models, their roots, similarities, and limitations. All of them are thermodynamic models. Naturally this bears a certain degree of redundancy, which forms the center of this review. After sixty years, many questions remain unanswered, mainly because our human longing for causality as base for understanding is not satisfied by thermodynamics alone. A description of allostery in terms of pathways, i.e., as a temporal chain of events, has been-, and still is-, a missing piece of the puzzle.
Collapse
Affiliation(s)
- Hagen Hofmann
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Herzl St. 234, 76100 Rehovot, Israel.
| |
Collapse
|
4
|
Sisk T, Robustelli P. Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.21.550103. [PMID: 37546728 PMCID: PMC10401938 DOI: 10.1101/2023.07.21.550103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.
Collapse
Affiliation(s)
- Thomas Sisk
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755
| | - Paul Robustelli
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755
| |
Collapse
|
5
|
Qiu Y, O’Connor MS, Xue M, Liu B, Huang X. An Efficient Path Classification Algorithm Based on Variational Autoencoder to Identify Metastable Path Channels for Complex Conformational Changes. J Chem Theory Comput 2023; 19:4728-4742. [PMID: 37382437 PMCID: PMC11042546 DOI: 10.1021/acs.jctc.3c00318] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Conformational changes (i.e., dynamic transitions between pairs of conformational states) play important roles in many chemical and biological processes. Constructing the Markov state model (MSM) from extensive molecular dynamics (MD) simulations is an effective approach to dissect the mechanism of conformational changes. When combined with transition path theory (TPT), MSM can be applied to elucidate the ensemble of kinetic pathways connecting pairs of conformational states. However, the application of TPT to analyze complex conformational changes often results in a vast number of kinetic pathways with comparable fluxes. This obstacle is particularly pronounced in heterogeneous self-assembly and aggregation processes. The large number of kinetic pathways makes it challenging to comprehend the molecular mechanisms underlying conformational changes of interest. To address this challenge, we have developed a path classification algorithm named latent-space path clustering (LPC) that efficiently lumps parallel kinetic pathways into distinct metastable path channels, making them easier to comprehend. In our algorithm, MD conformations are first projected onto a low-dimensional space containing a small set of collective variables (CVs) by time-structure-based independent component analysis (tICA) with kinetic mapping. Then, MSM and TPT are constructed to obtain the ensemble of pathways, and a deep learning architecture named the variational autoencoder (VAE) is used to learn the spatial distributions of kinetic pathways in the continuous CV space. Based on the trained VAE model, the TPT-generated ensemble of kinetic pathways can be embedded into a latent space, where the classification becomes clear. We show that LPC can efficiently and accurately identify the metastable path channels in three systems: a 2D potential, the aggregation of two hydrophobic particles in water, and the folding of the Fip35 WW domain. Using the 2D potential, we further demonstrate that our LPC algorithm outperforms the previous path-lumping algorithms by making substantially fewer incorrect assignments of individual pathways to four path channels. We expect that LPC can be widely applied to identify the dominant kinetic pathways underlying complex conformational changes.
Collapse
Affiliation(s)
- Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Michael S. O’Connor
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Mingyi Xue
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Bojun Liu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| |
Collapse
|
6
|
Asadollahi K, Scott DJ, Gooley PR. NMR applications to GPCR recognition by peptide ligands. Curr Opin Pharmacol 2023; 70:102366. [PMID: 37003111 DOI: 10.1016/j.coph.2023.102366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 01/30/2023] [Accepted: 02/11/2023] [Indexed: 04/03/2023]
Abstract
Peptides form the largest group of ligands that modulate the activity of more than 120 different GPCRs. Among which linear disordered peptide ligands usually undergo significant conformational changes upon binding that is essential for receptor recognition and activation. Conformational selection and induced fit are the extreme mechanisms of coupled folding and binding that can be distinguished by analysis of binding pathways by methods that include NMR. However, the large size of GPCRs in membrane-mimetic environments limits NMR applications. In this review, we highlight advances in the field that can be adopted to address coupled folding and binding of peptide ligands to their cognate receptors.
Collapse
Affiliation(s)
- Kazem Asadollahi
- Department of Biochemistry and Pharmacology, University of Melbourne, Parkville, Victoria, 3010, Australia; Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, 3010, Australia; The Florey Institute of Neuroscience and Mental Health, Parkville, Victoria, 3010, Australia
| | - Daniel J Scott
- Department of Biochemistry and Pharmacology, University of Melbourne, Parkville, Victoria, 3010, Australia; The Florey Institute of Neuroscience and Mental Health, Parkville, Victoria, 3010, Australia
| | - Paul R Gooley
- Department of Biochemistry and Pharmacology, University of Melbourne, Parkville, Victoria, 3010, Australia; Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, 3010, Australia.
| |
Collapse
|
7
|
A litmus test for classifying recognition mechanisms of transiently binding proteins. Nat Commun 2022; 13:3792. [PMID: 35778416 PMCID: PMC9249894 DOI: 10.1038/s41467-022-31374-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 06/15/2022] [Indexed: 11/17/2022] Open
Abstract
Partner recognition in protein binding is critical for all biological functions, and yet, delineating its mechanism is challenging, especially when recognition happens within microseconds. We present a theoretical and experimental framework based on straight-forward nuclear magnetic resonance relaxation dispersion measurements to investigate protein binding mechanisms on sub-millisecond timescales, which are beyond the reach of standard rapid-mixing experiments. This framework predicts that conformational selection prevails on ubiquitin’s paradigmatic interaction with an SH3 (Src-homology 3) domain. By contrast, the SH3 domain recognizes ubiquitin in a two-state binding process. Subsequent molecular dynamics simulations and Markov state modeling reveal that the ubiquitin conformation selected for binding exhibits a characteristically extended C-terminus. Our framework is robust and expandable for implementation in other binding scenarios with the potential to show that conformational selection might be the design principle of the hubs in protein interaction networks. The authors provide a litmus test for the recognition mechanism of transiently binding proteins based on nuclear magnetic resonance and find a conformational selection binding mechanism through concentration-dependent kinetics of ubiquitin and SH3.
Collapse
|
8
|
Antoszewski A, Lorpaiboon C, Strahan J, Dinner AR. Kinetics of Phenol Escape from the Insulin R 6 Hexamer. J Phys Chem B 2021; 125:11637-11649. [PMID: 34648712 DOI: 10.1021/acs.jpcb.1c06544] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Therapeutic preparations of insulin often contain phenolic molecules, which can impact both pharmacokinetics and shelf life. Thus, understanding the interactions of insulin and phenolic molecules can aid in designing improved therapeutics. In this study, we use molecular dynamics to investigate phenol release from the insulin hexamer. Leveraging recent advances in methods for analyzing molecular dynamics data, we expand on existing simulation studies to identify and quantitatively characterize six phenol binding/unbinding pathways for wild-type and A10 Ile → Val and B13 Glu → Gln mutant insulins. A number of these pathways involve large-scale opening of the primary escape channel, suggesting that the hexamer is much more dynamic than previously appreciated. We show that phenol unbinding is a multipathway process, with no single pathway representing more than 50% of the reactive current and all pathways representing at least 10%. We use the mutant simulations to show how the contributions of specific pathways can be rationally manipulated. Predicting the net effects of mutations is more challenging because the kinetics depend on all of the pathways, demanding quantitatively accurate simulations and experiments.
Collapse
Affiliation(s)
- Adam Antoszewski
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chatipat Lorpaiboon
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - John Strahan
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States.,James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States.,Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
9
|
Glielmo A, Husic BE, Rodriguez A, Clementi C, Noé F, Laio A. Unsupervised Learning Methods for Molecular Simulation Data. Chem Rev 2021; 121:9722-9758. [PMID: 33945269 PMCID: PMC8391792 DOI: 10.1021/acs.chemrev.0c01195] [Citation(s) in RCA: 116] [Impact Index Per Article: 38.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Indexed: 12/21/2022]
Abstract
Unsupervised learning is becoming an essential tool to analyze the increasingly large amounts of data produced by atomistic and molecular simulations, in material science, solid state physics, biophysics, and biochemistry. In this Review, we provide a comprehensive overview of the methods of unsupervised learning that have been most commonly used to investigate simulation data and indicate likely directions for further developments in the field. In particular, we discuss feature representation of molecular systems and present state-of-the-art algorithms of dimensionality reduction, density estimation, and clustering, and kinetic models. We divide our discussion into self-contained sections, each discussing a specific method. In each section, we briefly touch upon the mathematical and algorithmic foundations of the method, highlight its strengths and limitations, and describe the specific ways in which it has been used-or can be used-to analyze molecular simulation data.
Collapse
Affiliation(s)
- Aldo Glielmo
- International
School for Advanced Studies (SISSA) 34014 Trieste, Italy
| | - Brooke E. Husic
- Freie
Universität Berlin, Department of Mathematics
and Computer Science, 14195 Berlin, Germany
| | - Alex Rodriguez
- International Centre for Theoretical
Physics (ICTP), Condensed Matter and Statistical
Physics Section, 34100 Trieste, Italy
| | - Cecilia Clementi
- Freie
Universität Berlin, Department for
Physics, 14195 Berlin, Germany
- Rice
University Houston, Department of Chemistry, Houston, Texas 77005, United States
| | - Frank Noé
- Freie
Universität Berlin, Department of Mathematics
and Computer Science, 14195 Berlin, Germany
- Freie
Universität Berlin, Department for
Physics, 14195 Berlin, Germany
- Rice
University Houston, Department of Chemistry, Houston, Texas 77005, United States
| | - Alessandro Laio
- International
School for Advanced Studies (SISSA) 34014 Trieste, Italy
- International Centre for Theoretical
Physics (ICTP), Condensed Matter and Statistical
Physics Section, 34100 Trieste, Italy
| |
Collapse
|
10
|
Clerc I, Sagar A, Barducci A, Sibille N, Bernadó P, Cortés J. The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective. Comput Struct Biotechnol J 2021; 19:3817-3828. [PMID: 34285781 PMCID: PMC8273358 DOI: 10.1016/j.csbj.2021.06.031] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Revised: 06/17/2021] [Accepted: 06/21/2021] [Indexed: 01/15/2023] Open
Abstract
Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.
Collapse
Affiliation(s)
- Ilinka Clerc
- LAAS-CNRS, Université de Toulouse, CNRS, Toulouse, France
| | - Amin Sagar
- Centre de Biochimie Structurale, INSERM, CNRS, Université de Montpellier, France
| | - Alessandro Barducci
- Centre de Biochimie Structurale, INSERM, CNRS, Université de Montpellier, France
| | - Nathalie Sibille
- Centre de Biochimie Structurale, INSERM, CNRS, Université de Montpellier, France
| | - Pau Bernadó
- Centre de Biochimie Structurale, INSERM, CNRS, Université de Montpellier, France
| | - Juan Cortés
- LAAS-CNRS, Université de Toulouse, CNRS, Toulouse, France
| |
Collapse
|
11
|
Ge Y, Zhang S, Erdelyi M, Voelz VA. Solution-State Preorganization of Cyclic β-Hairpin Ligands Determines Binding Mechanism and Affinities for MDM2. J Chem Inf Model 2021; 61:2353-2367. [PMID: 33905247 PMCID: PMC9960209 DOI: 10.1021/acs.jcim.1c00029] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Understanding mechanisms of protein folding and binding is crucial to designing their molecular function. Molecular dynamics (MD) simulations and Markov state model (MSM) approaches provide a powerful way to understand complex conformational change that occurs over long time scales. Such dynamics are important for the design of therapeutic peptidomimetic ligands, whose affinity and binding mechanism are dictated by a combination of folding and binding. To examine the role of preorganization in peptide binding to protein targets, we performed massively parallel explicit-solvent MD simulations of cyclic β-hairpin ligands designed to mimic the p53 transactivation domain and competitively bind mouse double minute 2 homologue (MDM2). Disrupting the MDM2-p53 interaction is a therapeutic strategy to prevent degradation of the p53 tumor suppressor in cancer cells. MSM analysis of over 3 ms of aggregate trajectory data enabled us to build a detailed mechanistic model of coupled folding and binding of four cyclic peptides which we compare to experimental binding affinities and rates. The results show a striking relationship between the relative preorganization of each ligand in solution and its affinity for MDM2. Specifically, changes in peptide conformational populations predicted by the MSMs suggest that entropy loss upon binding is the main factor influencing affinity. The MSMs also enable detailed examination of non-native interactions which lead to misfolded states and comparison of structural ensembles with experimental NMR measurements. In contrast to an MSM study of p53 transactivation domain (TAD) binding to MDM2, MSMs of cyclic β-hairpin binding show a conformational selection mechanism. Finally, we make progress toward predicting accurate off rates of cyclic peptides using multiensemble Markov models (MEMMs) constructed from unbiased and biased simulated trajectories.
Collapse
Affiliation(s)
- Yunui Ge
- Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| | - Si Zhang
- Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| | - Mate Erdelyi
- Department of Chemistry - BMC, Uppsala University, SE-75123 Uppsala, Sweden
| | - Vincent A. Voelz
- Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| |
Collapse
|
12
|
Freitas FC, Ferreira PHB, Favaro DC, Oliveira RJD. Shedding Light on the Inhibitory Mechanisms of SARS-CoV-1/CoV-2 Spike Proteins by ACE2-Designed Peptides. J Chem Inf Model 2021; 61:1226-1243. [PMID: 33619962 PMCID: PMC7931628 DOI: 10.1021/acs.jcim.0c01320] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Indexed: 01/07/2023]
Abstract
Angiotensin-converting enzyme 2 (ACE2) is the host cellular receptor that locks onto the surface spike protein of the 2002 SARS coronavirus (SARS-CoV-1) and of the novel, highly transmissible and deadly 2019 SARS-CoV-2, responsible for the COVID-19 pandemic. One strategy to avoid the virus infection is to design peptides by extracting the human ACE2 peptidase domain α1-helix, which would bind to the coronavirus surface protein, preventing the virus entry into the host cells. The natural α1-helix peptide has a stronger affinity to SARS-CoV-2 than to SARS-CoV-1. Another peptide was designed by joining α1 with the second portion of ACE2 that is far in the peptidase sequence yet grafted in the spike protein interface with ACE2. Previous studies have shown that, among several α1-based peptides, the hybrid peptidic scaffold is the one with the highest/strongest affinity for SARS-CoV-1, which is comparable to the full-length ACE2 affinity. In this work, binding and folding dynamics of the natural and designed ACE2-based peptides were simulated by the well-known coarse-grained structure-based model, with the computed thermodynamic quantities correlating with the experimental binding affinity data. Furthermore, theoretical kinetic analysis of native contact formation revealed the distinction between these processes in the presence of the different binding partners SARS-CoV-1 and SARS-CoV-2 spike domains. Additionally, our results indicate the existence of a two-state folding mechanism for the designed peptide en route to bind to the spike proteins, in contrast to a downhill mechanism for the natural α1-helix peptides. The presented low-cost simulation protocol demonstrated its efficiency in evaluating binding affinities and identifying the mechanisms involved in the neutralization of spike-ACE2 interaction by designed peptides. Finally, the protocol can be used as a computer-based screening of more potent designed peptides by experimentalists searching for new therapeutics against COVID-19.
Collapse
Affiliation(s)
- Frederico Campos Freitas
- Laboratório de Biofísica Teórica,
Departamento de Física, Instituto de Ciências Exatas, Naturais
e Educação, Universidade Federal do Triângulo
Mineiro, Uberaba, MG 38064-200, Brazil
| | - Paulo Henrique Borges Ferreira
- Laboratório de Biofísica Teórica,
Departamento de Física, Instituto de Ciências Exatas, Naturais
e Educação, Universidade Federal do Triângulo
Mineiro, Uberaba, MG 38064-200, Brazil
| | - Denize Cristina Favaro
- Departamento de Química Orgânica,
Instituto de Química, Universidade Estadual de
Campinas, São Paulo, SP 13083-970, Brazil
| | - Ronaldo Junio de Oliveira
- Laboratório de Biofísica Teórica,
Departamento de Física, Instituto de Ciências Exatas, Naturais
e Educação, Universidade Federal do Triângulo
Mineiro, Uberaba, MG 38064-200, Brazil
| |
Collapse
|
13
|
Wang W. Recent advances in atomic molecular dynamics simulation of intrinsically disordered proteins. Phys Chem Chem Phys 2021; 23:777-784. [PMID: 33355572 DOI: 10.1039/d0cp05818a] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Intrinsically disordered proteins (IDPs) play important roles in cellular functions. The inherent structural heterogeneity of IDPs makes the high-resolution experimental characterization of IDPs extremely difficult. Molecular dynamics (MD) simulation could provide the atomic-level description of the structural and dynamic properties of IDPs. This perspective reviews the recent progress in atomic MD simulation studies of IDPs, including the development of force fields and sampling methods, as well as applications in IDP-involved protein-protein interactions. The employment of large-scale simulations and advanced sampling techniques allows more accurate estimation of the thermodynamics and kinetics of IDP-mediated protein interactions, and the holistic landscape of the binding process of IDPs is emerging.
Collapse
Affiliation(s)
- Wenning Wang
- Department of Chemistry, Multiscale Research Institute of Complex Systems and Institute of Biomedical Sciences, Fudan University, Shanghai 200438, China.
| |
Collapse
|
14
|
Abstract
Molecular dynamics simulations can now routinely access the microsecond timescale, making feasible direct sampling of ligand association events. While Markov State Model (MSM) approaches offer a useful framework for analyzing such trajectory data to gain insight into binding mechanisms, accurate modeling of ligand association pathways and kinetics must be done carefully. We describe methods and good practices for constructing MSMs of ligand binding from unbiased trajectory data and discuss how to use time-lagged independent component analysis (tICA) to build informative models, using as an example recent simulation work to model the binding of phenylalanine to the regulatory ACT domain dimer of phenylalanine hydroxylase. We describe a variety of methods for estimating association rates from MSMs and discuss how to distinguish between conformational selection and induced-fit mechanisms using MSMs. In addition, we review some examples of MSMs constructed to elucidate the mechanisms by which p53 transactivation domain (TAD) and related peptides bind the oncoprotein MDM2.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Vincent A Voelz
- Department of Chemistry, Temple University, Philadelphia, PA, USA.
| |
Collapse
|
15
|
Zou Y, Yang S, Sanders JN, Li W, Yu P, Wang H, Tang Z, Liu W, Houk KN. Computational Investigation of the Mechanism of Diels-Alderase PyrI4. J Am Chem Soc 2020; 142:20232-20239. [PMID: 33190496 DOI: 10.1021/jacs.0c10813] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We studied the mechanisms of activation and stereoselectivity of a monofunctional Diels-Alderase (PyrI4)-catalyzed intramolecular Diels-Alder reaction that leads to formation of the key spiro-tetramate moiety in the biosynthesis of the pyrroindomycin family of natural products. Key activation effects of PyrI4 include acid catalysis and an induced-fit mechanism that cooperate with the unique "lid" feature of PyrI4 to stabilize the Diels-Alder transition state. PyrI4 enhances the intrinsic Diels-Alder stereoselectivity of the substrate and leads to stereospecific formation of the product.
Collapse
Affiliation(s)
- Yike Zou
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Song Yang
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Jacob N Sanders
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Wei Li
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Peiyuan Yu
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Hongbo Wang
- State Key Laboratory of Bioorganic and Natural Products Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, 345 Lingling Road, Shanghai 200032, China
| | - Zhijun Tang
- State Key Laboratory of Bioorganic and Natural Products Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, 345 Lingling Road, Shanghai 200032, China
| | - Wen Liu
- State Key Laboratory of Bioorganic and Natural Products Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, 345 Lingling Road, Shanghai 200032, China
| | - K N Houk
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| |
Collapse
|
16
|
Antoszewski A, Feng CJ, Vani BP, Thiede EH, Hong L, Weare J, Tokmakoff A, Dinner AR. Insulin Dissociates by Diverse Mechanisms of Coupled Unfolding and Unbinding. J Phys Chem B 2020; 124:5571-5587. [PMID: 32515958 PMCID: PMC7774804 DOI: 10.1021/acs.jpcb.0c03521] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The protein hormone insulin exists in various oligomeric forms, and a key step in binding its cellular receptor is dissociation of the dimer. This dissociation process and its corresponding association process have come to serve as paradigms of coupled (un)folding and (un)binding more generally. Despite its fundamental and practical importance, the mechanism of insulin dimer dissociation remains poorly understood. Here, we use molecular dynamics simulations, leveraging recent developments in umbrella sampling, to characterize the energetic and structural features of dissociation in unprecedented detail. We find that the dissociation is inherently multipathway with limiting behaviors corresponding to conformational selection and induced fit, the two prototypical mechanisms of coupled folding and binding. Along one limiting path, the dissociation leads to detachment of the C-terminal segment of the insulin B chain from the protein core, a feature believed to be essential for receptor binding. We simulate IR spectroscopy experiments to aid in interpreting current experiments and identify sites where isotopic labeling can be most effective for distinguishing the contributions of the limiting mechanisms.
Collapse
Affiliation(s)
- Adam Antoszewski
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chi-Jui Feng
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Bodhi P Vani
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Erik H Thiede
- Department of Computer Science, The University of Chicago, Chicago, Illinois 60637, United States
- Department of Statistics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Lu Hong
- Graduate Program in Biophysical Sciences, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jonathan Weare
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
| | - Andrei Tokmakoff
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
- Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
17
|
Robustelli P, Piana S, Shaw DE. Mechanism of Coupled Folding-upon-Binding of an Intrinsically Disordered Protein. J Am Chem Soc 2020; 142:11092-11101. [DOI: 10.1021/jacs.0c03217] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Paul Robustelli
- D. E. Shaw Research, New York, New York 10036, United States
| | - Stefano Piana
- D. E. Shaw Research, New York, New York 10036, United States
| | - David E. Shaw
- D. E. Shaw Research, New York, New York 10036, United States
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, New York 10032, United States
| |
Collapse
|
18
|
Yang J, Zeng Y, Liu Y, Gao M, Liu S, Su Z, Huang Y. Electrostatic interactions in molecular recognition of intrinsically disordered proteins. J Biomol Struct Dyn 2019; 38:4883-4894. [DOI: 10.1080/07391102.2019.1692073] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Jing Yang
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Yifan Zeng
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Yunfei Liu
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Meng Gao
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Sen Liu
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Zhengding Su
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| | - Yongqi Huang
- Key Laboratory of Industrial Fermentation (Ministry of Education), Hubei University of Technology, Wuhan, China
- Hubei Key Laboratory of Industrial Microbiology, Hubei University of Technology, Wuhan, China
| |
Collapse
|
19
|
Sun B, Vaughan D, Tikunova S, Creamer TP, Davis JP, Kekenes-Huskey PM. Calmodulin-Calcineurin Interaction beyond the Calmodulin-Binding Region Contributes to Calcineurin Activation. Biochemistry 2019; 58:4070-4085. [PMID: 31483613 DOI: 10.1021/acs.biochem.9b00626] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Calcineurin (CaN) is a calcium-dependent phosphatase involved in numerous signaling pathways. Its activation is in part driven by the binding of calmodulin (CaM) to a CaM recognition region (CaMBR) within CaN's regulatory domain (RD). However, secondary interactions between CaM and the CaN RD may be necessary to fully activate CaN. Specifically, it is established that the CaN RD folds upon CaM binding and a region C-terminal to CaMBR, the "distal helix", assumes an α-helix fold and contributes to activation [Dunlap, T. B., et al. (2013) Biochemistry 52, 8643-8651]. We hypothesized in that previous study that this distal helix can bind CaM in a region distinct from the canonical CaMBR. To test this hypothesis, we utilized molecular simulations, including replica-exchange molecular dynamics, protein-protein docking, and computational mutagenesis, to determine potential distal helix-binding sites on CaM's surface. We isolated a potential binding site on CaM (site D) that facilitates moderate-affinity interprotein interactions and predicted that mutation of site D residues K30 and G40 on CaM would weaken CaN distal helix binding. We experimentally confirmed that two variants (K30E and G40D) indicate weaker binding of a phosphate substrate p-nitrophenyl phosphate to the CaN catalytic site by a phosphatase assay. This weakened substrate affinity is consistent with competitive binding of the CaN autoinhibition domain to the catalytic site, which we suggest is due to the weakened distal helix-CaM interactions. This study therefore suggests a novel mechanism for CaM regulation of CaN that may extend to other CaM targets.
Collapse
Affiliation(s)
- Bin Sun
- Department of Chemistry , University of Kentucky , Lexington , Kentucky 40506 , United States
| | - Darin Vaughan
- Department of Chemistry , University of Kentucky , Lexington , Kentucky 40506 , United States
| | - Svetlana Tikunova
- Department of Physiology and Cell Biology , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Trevor P Creamer
- Center for Structural Biology and Department of Molecular & Cellular Biochemistry , University of Kentucky , Lexington , Kentucky 40536 , United States
| | - Jonathan P Davis
- Department of Physiology and Cell Biology , The Ohio State University , Columbus , Ohio 43210 , United States
| | - P M Kekenes-Huskey
- Department of Chemical and Materials Engineering , University of Kentucky , Lexington , Kentucky 40506 , United States.,Department of Cell and Molecular Physiology , Loyola University Chicago , Maywood , Illinois 60153 , United States
| |
Collapse
|
20
|
Ge Y, Borne E, Stewart S, Hansen MR, Arturo EC, Jaffe EK, Voelz VA. Simulations of the regulatory ACT domain of human phenylalanine hydroxylase (PAH) unveil its mechanism of phenylalanine binding. J Biol Chem 2018; 293:19532-19543. [PMID: 30287685 DOI: 10.1074/jbc.ra118.004909] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Revised: 09/17/2018] [Indexed: 12/20/2022] Open
Abstract
Phenylalanine hydroxylase (PAH) regulates phenylalanine (Phe) levels in mammals to prevent neurotoxicity resulting from high Phe concentrations as observed in genetic disorders leading to hyperphenylalaninemia and phenylketonuria. PAH senses elevated Phe concentrations by transient allosteric Phe binding to a protein-protein interface between ACT domains of different subunits in a PAH tetramer. This interface is present in an activated PAH (A-PAH) tetramer and absent in a resting-state PAH (RS-PAH) tetramer. To investigate this allosteric sensing mechanism, here we used the GROMACS molecular dynamics simulation suite on the Folding@home computing platform to perform extensive molecular simulations and Markov state model (MSM) analysis of Phe binding to ACT domain dimers. These simulations strongly implicated a conformational selection mechanism for Phe association with ACT domain dimers and revealed protein motions that act as a gating mechanism for Phe binding. The MSMs also illuminate a highly mobile hairpin loop, consistent with experimental findings also presented here that the PAH variant L72W does not shift the PAH structural equilibrium toward the activated state. Finally, simulations of ACT domain monomers are presented, in which spontaneous transitions between resting-state and activated conformations are observed, also consistent with a mechanism of conformational selection. These mechanistic details provide detailed insight into the regulation of PAH activation and provide testable hypotheses for the development of new allosteric effectors to correct structural and functional defects in PAH.
Collapse
Affiliation(s)
- Yunhui Ge
- From the Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122
| | - Elias Borne
- Fox Chase Cancer Center, Temple University Health System, Philadelphia, Pennsylvania 19111, and
| | - Shannon Stewart
- Fox Chase Cancer Center, Temple University Health System, Philadelphia, Pennsylvania 19111, and
| | - Michael R Hansen
- Fox Chase Cancer Center, Temple University Health System, Philadelphia, Pennsylvania 19111, and
| | - Emilia C Arturo
- Fox Chase Cancer Center, Temple University Health System, Philadelphia, Pennsylvania 19111, and.,Drexel University College of Medicine, Philadelphia, Pennsylvania 19129
| | - Eileen K Jaffe
- Fox Chase Cancer Center, Temple University Health System, Philadelphia, Pennsylvania 19111, and
| | - Vincent A Voelz
- From the Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122,
| |
Collapse
|
21
|
Salmaso V, Moro S. Bridging Molecular Docking to Molecular Dynamics in Exploring Ligand-Protein Recognition Process: An Overview. Front Pharmacol 2018; 9:923. [PMID: 30186166 PMCID: PMC6113859 DOI: 10.3389/fphar.2018.00923] [Citation(s) in RCA: 311] [Impact Index Per Article: 51.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Accepted: 07/26/2018] [Indexed: 12/22/2022] Open
Abstract
Computational techniques have been applied in the drug discovery pipeline since the 1980s. Given the low computational resources of the time, the first molecular modeling strategies relied on a rigid view of the ligand-target binding process. During the years, the evolution of hardware technologies has gradually allowed simulating the dynamic nature of the binding event. In this work, we present an overview of the evolution of structure-based drug discovery techniques in the study of ligand-target recognition phenomenon, going from the static molecular docking toward enhanced molecular dynamics strategies.
Collapse
Affiliation(s)
- Veronica Salmaso
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy
| | - Stefano Moro
- Molecular Modeling Section, Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Padova, Italy
| |
Collapse
|