1
|
Strahan J, Lorpaiboon C, Weare J, Dinner AR. BAD-NEUS: Rapidly converging trajectory stratification. J Chem Phys 2024; 161:084109. [PMID: 39185846 PMCID: PMC11349377 DOI: 10.1063/5.0215975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2024] [Accepted: 07/25/2024] [Indexed: 08/27/2024] Open
Abstract
An issue for molecular dynamics simulations is that events of interest often involve timescales that are much longer than the simulation time step, which is set by the fastest timescales of the model. Because of this timescale separation, direct simulation of many events is prohibitively computationally costly. This issue can be overcome by aggregating information from many relatively short simulations that sample segments of trajectories involving events of interest. This is the strategy of Markov state models (MSMs) and related approaches, but such methods suffer from approximation error because the variables defining the states generally do not capture the dynamics fully. By contrast, once converged, the weighted ensemble (WE) method aggregates information from trajectory segments so as to yield unbiased estimates of both thermodynamic and kinetic statistics. Unfortunately, errors decay no faster than unbiased simulation in WE as originally formulated and commonly deployed. Here, we introduce a theoretical framework for describing WE that shows that the introduction of an approximate stationary distribution on top of the stratification, as in nonequilibrium umbrella sampling (NEUS), accelerates convergence. Building on ideas from MSMs and related methods, we generalize the NEUS approach in such a way that the approximation error can be reduced systematically. We show that the improved algorithm can decrease the simulation time required to achieve the desired precision by orders of magnitude.
Collapse
Affiliation(s)
- John Strahan
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Chatipat Lorpaiboon
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - Jonathan Weare
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
| | - Aaron R. Dinner
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
2
|
Wang J, Miao Y. Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. J Chem Theory Comput 2024; 20:5829-5841. [PMID: 39002136 DOI: 10.1021/acs.jctc.4c00502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
3
|
Wang J, Miao Y. Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.06.592668. [PMID: 38766067 PMCID: PMC11100592 DOI: 10.1101/2024.05.06.592668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| |
Collapse
|
4
|
Xu C, Zhang X, Zhao L, Verkhivker GM, Bai F. Accurate Characterization of Binding Kinetics and Allosteric Mechanisms for the HSP90 Chaperone Inhibitors Using AI-Augmented Integrative Biophysical Studies. JACS AU 2024; 4:1632-1645. [PMID: 38665669 PMCID: PMC11040708 DOI: 10.1021/jacsau.4c00123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 04/28/2024]
Abstract
The binding kinetics of drugs to their targets are gradually being recognized as a crucial indicator of the efficacy of drugs in vivo, leading to the development of various computational methods for predicting the binding kinetics in recent years. However, compared with the prediction of binding affinity, the underlying structure and dynamic determinants of binding kinetics are more complicated. Efficient and accurate methods for predicting binding kinetics are still lacking. In this study, quantitative structure-kinetics relationship (QSKR) models were developed using 132 inhibitors targeting the ATP binding domain of heat shock protein 90α (HSP90α) to predict the dissociation rate constant (koff), enabling a direct assessment of the drug-target residence time. These models demonstrated good predictive performance, where hydrophobic and hydrogen bond interactions significantly influence the koff prediction. In subsequent applications, our models were used to assist in the discovery of new inhibitors for the N-terminal domain of HSP90α (N-HSP90α), demonstrating predictive capabilities on an experimental validation set with a new scaffold. In X-ray crystallography experiments, the loop-middle conformation of apo N-HSP90α was observed for the first time (previously, the loop-middle conformation had only been observed in holo-N-HSP90α structures). Interestingly, we observed different conformations of apo N-HSP90α simultaneously in an asymmetric unit, which was also observed in a holo-N-HSP90α structure, suggesting an equilibrium of conformations between different states in solution, which could be one of the determinants affecting the binding kinetics of the ligand. Different ligands can undergo conformational selection or alter the equilibrium of conformations, inducing conformational rearrangements and resulting in different effects on binding kinetics. We then used molecular dynamics simulations to describe conformational changes of apo N-HSP90α in different conformational states. In summary, the study of the binding kinetics and molecular mechanisms of N-HSP90α provides valuable information for the development of more targeted therapeutic approaches.
Collapse
Affiliation(s)
- Chao Xu
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Xianglei Zhang
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Lianghao Zhao
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Graduate Program in Computational
and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Fang Bai
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
- School
of Information Science and Technology, ShanghaiTech
University, 393 Middle Huaxia Road, Shanghai 201210, China
- Shanghai
Clinical Research and Trial Center, Shanghai 201210, China
| |
Collapse
|
5
|
Bose S, Lotz SD, Deb I, Shuck M, Lee KSS, Dickson A. How Robust Is the Ligand Binding Transition State? J Am Chem Soc 2023; 145:25318-25331. [PMID: 37943667 PMCID: PMC11059145 DOI: 10.1021/jacs.3c08940] [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: 11/12/2023]
Abstract
For many drug targets, it has been shown that the kinetics of drug binding (e.g., on rate and off rate) is more predictive of drug efficacy than thermodynamic quantities alone. This motivates the development of predictive computational models that can be used to optimize compounds on the basis of their kinetics. The structural details underpinning these computational models are found not only in the bound state but also in the short-lived ligand binding transition states. Although transition states cannot be directly observed experimentally due to their extremely short lifetimes, recent successes have demonstrated that modeling the ligand binding transition state is possible with the help of enhanced sampling molecular dynamics methods. Previously, we generated unbinding paths for an inhibitor of soluble epoxide hydrolase (sEH) with a residence time of 11 min. Here, we computationally modeled unbinding events with the weighted ensemble method REVO (resampling of ensembles by variation optimization) for five additional inhibitors of sEH with residence times ranging from 14.25 to 31.75 min, with average prediction accuracy within an order of magnitude. The unbinding ensembles are analyzed in detail, focusing on features of the ligand binding transition state ensembles (TSEs). We find that ligands with similar bound poses can show significant differences in their ligand binding TSEs, in terms of their spatial distribution and protein-ligand interactions. However, we also find similarities across the TSEs when examining more general features such as ligand degrees of freedom. Together these findings show significant challenges for rational, kinetics-based drug design.
Collapse
Affiliation(s)
- Samik Bose
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Samuel D Lotz
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Indrajit Deb
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Megan Shuck
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Kin Sing Stephen Lee
- Department of Pharmacology and Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Institute of Integrative Toxicology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Alex Dickson
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
6
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python Package for the Analysis of Simulations of Intrinsically Disordered Proteins. J Chem Theory Comput 2023; 19:5609-5620. [PMID: 37463458 PMCID: PMC11188088 DOI: 10.1021/acs.jctc.3c00190] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite simulation analysis of unfolded regions of proteins (SOURSOP), an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
Affiliation(s)
- Jared M. Lalmansingh
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex T. Keeley
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana-Champaign, IL 61801, USA
| | - Kiersten M. Ruff
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Rohit V. Pappu
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO 63130, USA
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
| | - Alex S. Holehouse
- Center for Biomolecular Condensates, Washington University in St. Louis, St. Louis, MO 63130, USA
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO, 63110, USA
| |
Collapse
|
7
|
Babin KM, Karim JA, Gordon PH, Lennon J, Dickson A, Pioszak AA. Adrenomedullin 2/intermedin is a slow off-rate, long-acting endogenous agonist of the adrenomedullin 2 G protein-coupled receptor. J Biol Chem 2023:104785. [PMID: 37146967 DOI: 10.1016/j.jbc.2023.104785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 04/20/2023] [Accepted: 05/01/2023] [Indexed: 05/07/2023] Open
Abstract
Adrenomedullin 2/intermedin (AM2/IMD), adrenomedullin (AM), and calcitonin gene-related peptide (CGRP) have signaling functions in the cardiovascular, lymphatic, and nervous systems by activating three heterodimeric receptors comprised of the class B GPCR CLR and a RAMP1, -2, or -3 modulatory subunit. CGRP and AM prefer the RAMP1 and RAMP2/3 complexes, respectively, whereas AM2/IMD is thought to be relatively non-selective. Accordingly, AM2/IMD exhibits overlapping actions with CGRP and AM, so the rationale for this third agonist for the CLR-RAMP complexes is unclear. Here, we report that AM2/IMD is kinetically selective for CLR-RAMP3, known as the AM2R, and we define the structural basis for its distinct kinetics. In live cell biosensor assays, AM2/IMD-AM2R elicited substantially longer duration cAMP signaling than the eight other peptide-receptor combinations. AM2/IMD and AM bound the AM2R with similar equilibrium affinities, but AM2/IMD had a much slower off-rate and longer receptor residence time, thus explaining its prolonged signaling capacity. Peptide and receptor chimeras and mutagenesis were used to map the regions responsible for the distinct binding and signaling kinetics to the AM2/IMD mid-region and the RAMP3 extracellular domain (ECD). Molecular dynamics simulations revealed how the former forms stable interactions at the CLR ECD-transmembrane domain interface and how the latter augments the CLR ECD binding pocket to anchor the AM2/IMD C-terminus. These two strong binding components only combine in the AM2R. Our findings uncover AM2/IMD-AM2R as a cognate pair with unique temporal features, reveal how AM2/IMD and RAMP3 collaborate to shape CLR signaling, and have significant implications for AM2/IMD biology.
Collapse
Affiliation(s)
- Katie M Babin
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - Jordan A Karim
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - Peyton H Gordon
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - James Lennon
- Departments of Biochemistry and Molecular Biology and Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824
| | - Alex Dickson
- Departments of Biochemistry and Molecular Biology and Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824.
| | - Augen A Pioszak
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104.
| |
Collapse
|
8
|
Roussey NM, Dickson A. Quality over quantity: Sampling high probability rare events with the weighted ensemble algorithm. J Comput Chem 2023; 44:935-947. [PMID: 36510846 PMCID: PMC10164457 DOI: 10.1002/jcc.27054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 10/27/2022] [Accepted: 11/27/2022] [Indexed: 12/15/2022]
Abstract
The prediction of (un)binding rates and free energies is of great significance to the drug design process. Although many enhanced sampling algorithms and approaches have been developed, there is not yet a reliable workflow to predict these quantities. Previously we have shown that free energies and transition rates can be calculated by directly simulating the binding and unbinding processes with our variant of the WE algorithm "Resampling of Ensembles by Variation Optimization", or "REVO". Here, we calculate binding free energies retrospectively for three SAMPL6 host-guest systems and prospectively for a SAMPL9 system to test a modification of REVO that restricts its cloning behavior in quasi-unbound states. Specifically, trajectories cannot clone if they meet a physical requirement that represents a high likelihood of unbinding, which in the case of this work is a center-of-mass to center-of-mass distance. The overall effect of this change was difficult to predict, as it results in fewer unbinding events each of which with a much higher statistical weight. For all four systems tested, this new strategy produced either more accurate unbinding free energies or more consistent results between simulations than the standard REVO algorithm. This approach is highly flexible, and any feature of interest for a system can be used to determine cloning eligibility. These findings thus constitute an important improvement in the calculation of transition rates and binding free energies with the weighted ensemble method.
Collapse
Affiliation(s)
- Nicole M Roussey
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA
| | - Alex Dickson
- Department of Biochemistry and Molecular Biology, Department of Computational Mathematics, Science, and Engineering, Michigan State University, East Lansing, Michigan, USA
| |
Collapse
|
9
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
10
|
Lalmansingh JM, Keeley AT, Ruff KM, Pappu RV, Holehouse AS. SOURSOP: A Python package for the analysis of simulations of intrinsically disordered proteins. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.16.528879. [PMID: 36824878 PMCID: PMC9949127 DOI: 10.1101/2023.02.16.528879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
Abstract
Conformational heterogeneity is a defining hallmark of intrinsically disordered proteins and protein regions (IDRs). The functions of IDRs and the emergent cellular phenotypes they control are associated with sequence-specific conformational ensembles. Simulations of conformational ensembles that are based on atomistic and coarse-grained models are routinely used to uncover the sequence-specific interactions that may contribute to IDR functions. These simulations are performed either independently or in conjunction with data from experiments. Functionally relevant features of IDRs can span a range of length scales. Extracting these features requires analysis routines that quantify a range of properties. Here, we describe a new analysis suite SOURSOP, an object-oriented and open-source toolkit designed for the analysis of simulated conformational ensembles of IDRs. SOURSOP implements several analysis routines motivated by principles in polymer physics, offering a unique collection of simple-to-use functions to characterize IDR ensembles. As an extendable framework, SOURSOP supports the development and implementation of new analysis routines that can be easily packaged and shared.
Collapse
|
11
|
Ligand Gaussian Accelerated Molecular Dynamics 2 (LiGaMD2): Improved Calculations of Ligand Binding Thermodynamics and Kinetics with Closed Protein Pocket. J Chem Theory Comput 2023; 19:733-745. [PMID: 36706316 DOI: 10.1021/acs.jctc.2c01194] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Ligand binding thermodynamics and kinetics are critical parameters for drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics from molecular simulations due to limited simulation timescales. Protein dynamics, especially in the ligand binding pocket, often plays an important role in ligand binding. Based on our previously developed Ligand Gaussian accelerated molecular dynamics (LiGaMD), here we present LiGaMD2 in which a selective boost potential was applied to both the ligand and protein residues in the binding pocket to improve sampling of ligand binding and dissociation. To validate the performance of LiGaMD2, the T4 lysozyme (T4L) mutants with open and closed pockets bound by different ligands were chosen as model systems. LiGaMD2 could efficiently capture repetitive ligand dissociation and binding within microsecond simulations of all T4L systems. The obtained ligand binding kinetic rates and free energies agreed well with available experimental values and previous modeling results. Therefore, LiGaMD2 provides an improved approach to sample opening of closed protein pockets for ligand dissociation and binding, thereby allowing for efficient calculations of ligand binding thermodynamics and kinetics.
Collapse
|
12
|
Babin KM, Karim JA, Gordon PH, Lennon J, Dickson A, Pioszak AA. Adrenomedullin 2/intermedin is a slow off-rate, long-acting endogenous agonist of the adrenomedullin 2 G protein-coupled receptor. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.13.523955. [PMID: 36711519 PMCID: PMC9882245 DOI: 10.1101/2023.01.13.523955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
The signaling peptides adrenomedullin 2/intermedin (AM2/IMD), adrenomedullin (AM), and CGRP have overlapping and distinct functions in the cardiovascular, lymphatic, and nervous systems by activating three shared receptors comprised of the class B GPCR CLR in complex with a RAMP1, -2, or -3 modulatory subunit. Here, we report that AM2/IMD, which is thought to be a non-selective agonist, is kinetically selective for CLR-RAMP3, known as the AM 2 R. AM2/IMD-AM 2 R elicited substantially longer duration cAMP signaling than the eight other peptide-receptor combinations due to AM2/IMD slow off-rate binding kinetics. The regions responsible for the slow off-rate were mapped to the AM2/IMD mid-region and the RAMP3 extracellular domain. MD simulations revealed how these bestow enhanced stability to the complex. Our results uncover AM2/IMD-AM 2 R as a cognate pair with unique temporal features, define the mechanism of kinetic selectivity, and explain how AM2/IMD and RAMP3 collaborate to shape the signaling output of a clinically important GPCR.
Collapse
Affiliation(s)
- Katie M. Babin
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - Jordan A. Karim
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - Peyton H. Gordon
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
| | - James Lennon
- Departments of Biochemistry and Molecular Biology and Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824
| | - Alex Dickson
- Departments of Biochemistry and Molecular Biology and Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824
| | - Augen A. Pioszak
- Department of Biochemistry and Molecular Biology, University of Oklahoma Health Sciences Center, Oklahoma City, OK 73104
- Lead contact
| |
Collapse
|
13
|
Aristoff D, Copperman J, Simpson G, Webber RJ, Zuckerman DM. Weighted ensemble: Recent mathematical developments. J Chem Phys 2023; 158:014108. [PMID: 36610976 PMCID: PMC9822651 DOI: 10.1063/5.0110873] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2022] [Accepted: 12/04/2022] [Indexed: 12/12/2022] Open
Abstract
Weighted ensemble (WE) is an enhanced sampling method based on periodically replicating and pruning trajectories generated in parallel. WE has grown increasingly popular for computational biochemistry problems due, in part, to improved hardware and accessible software implementations. Algorithmic and analytical improvements have played an important role, and progress has accelerated in recent years. Here, we discuss and elaborate on the WE method from a mathematical perspective, highlighting recent results that enhance the computational efficiency. The mathematical theory reveals a new strategy for optimizing trajectory management that approaches the best possible variance while generalizing to systems of arbitrary dimension.
Collapse
Affiliation(s)
- D. Aristoff
- Mathematics, Colorado State University, Fort Collins, CO 80521 USA
| | - J. Copperman
- Biomedical Engineering, Oregon Health and Science University, Portland, OR 97239 USA
| | - G. Simpson
- Mathematics, Drexel University, Philadelphia, Pennsylvania 19104 USA
| | - R. J. Webber
- Computing and Mathematical Sciences, California Institute of Technology, Pasadena, California 91125 USA
| | - D. M. Zuckerman
- Biomedical Engineering, Oregon Health and Science University, Portland, OR 97239 USA
| |
Collapse
|
14
|
Wang KW, Lee J, Zhang H, Suh D, Im W. CHARMM-GUI Implicit Solvent Modeler for Various Generalized Born Models in Different Simulation Programs. J Phys Chem B 2022; 126:7354-7364. [PMID: 36117287 DOI: 10.1021/acs.jpcb.2c05294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Implicit solvent models are widely used because they are advantageous to speed up simulations by drastically decreasing the number of solvent degrees of freedom, which allows one to achieve long simulation time scales for large system sizes. CHARMM-GUI, a web-based platform, has been developed to support the setup of complex multicomponent molecular systems and prepare input files. This study describes an Implicit Solvent Modeler (ISM) in CHARMM-GUI for various generalized Born (GB) implicit solvent simulations in different molecular dynamics programs such as AMBER, CHARMM, GENESIS, NAMD, OpenMM, and Tinker. The GB models available in ISM include GB-HCT, GB-OBC, GB-neck, GBMV, and GBSW with the CHARMM and Amber force fields for protein, DNA, RNA, glycan, and ligand systems. Using the system and input files generated by ISM, implicit solvent simulations of protein, DNA, and RNA systems produce similar results for different simulation packages with the same input information. Protein-ligand systems are also considered to further validate the systems and input files generated by ISM. Simple ligand root-mean-square deviation (RMSD) and molecular mechanics generalized Born surface area (MM/GBSA) calculations show that the performance of implicit simulations is better than docking and can be used for early stage ligand screening. These reasonable results indicate that ISM is a useful and reliable tool to provide various implicit solvent simulation applications.
Collapse
Affiliation(s)
- Kye Won Wang
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jumin Lee
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Han Zhang
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Donghyuk Suh
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
15
|
Vymětal J, Vondrášek J. Iterative Landmark-Based Umbrella Sampling (ILBUS) Protocol for Sampling of Conformational Space of Biomolecules. J Chem Inf Model 2022; 62:4783-4798. [PMID: 36122323 DOI: 10.1021/acs.jcim.2c00370] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Computer simulations of biomolecules such as molecular dynamics often suffer from insufficient sampling. Due to limited computational resources, insufficient sampling prevents obtaining proper equilibrium distributions of observed properties. To deal with this problem, we proposed a simulation protocol for efficient resampling of collected off-equilibrium trajectories. These trajectories are utilized for the initial mapping of the conformational space, which is later properly resampled by the introduced Iterative Landmark-Based Umbrella Sampling (ILBUS) method. Reconstruction of static equilibrium properties is achieved by the multistate Bennett acceptance ratio (MBAR) method, which enables efficient use of simulated data. The ILBUS protocol is geometry-based and does not demand any additional collective variable or a dimensional-reduction technique. The only requirement is a set of suitably spaced reference conformations, which serve as landmarks in the mapped conformational space. Additionally, the ILBUS protocol encompasses an iterative process that optimizes the force constant used in the umbrella sampling simulation. Such tuning is an inherent feature of the protocol and does not need to be performed by the user in advance. Furthermore, even the simulations with suboptimal force constants can be used in estimates by MBAR. We demonstrate the feasibility and the performance of this approach in the study of the conformational landscape of the alanine dipeptide, met-enkephalin, and adenylate kinase.
Collapse
Affiliation(s)
- Jiří Vymětal
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 542/2, 160 00 Praha 6, Czech Republic
| | - Jiří Vondrášek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 542/2, 160 00 Praha 6, Czech Republic
| |
Collapse
|
16
|
Ahmad K, Rizzi A, Capelli R, Mandelli D, Lyu W, Carloni P. Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective. Front Mol Biosci 2022; 9:899805. [PMID: 35755817 PMCID: PMC9216551 DOI: 10.3389/fmolb.2022.899805] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Accepted: 05/09/2022] [Indexed: 12/12/2022] Open
Abstract
The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.
Collapse
Affiliation(s)
- Katya Ahmad
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Andrea Rizzi
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Atomistic Simulations, Istituto Italiano di Tecnologia, Genova, Italy
| | - Riccardo Capelli
- Department of Applied Science and Technology (DISAT), Politecnico di Torino, Torino, Italy
| | - Davide Mandelli
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Wenping Lyu
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen, China
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, China
| | - Paolo Carloni
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich, Jülich, Germany
| |
Collapse
|
17
|
Gilboa T, Ogata AF, Reilly C, Walt DR. Single-molecule studies reveal method for tuning the heterogeneous activity of Alkaline Phosphatase. Biophys J 2022; 121:2027-2034. [DOI: 10.1016/j.bpj.2022.05.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2021] [Revised: 03/09/2022] [Accepted: 05/04/2022] [Indexed: 11/28/2022] Open
|
18
|
Roussey NM, Dickson A. Local Ion Densities can Influence Transition Paths of Molecular Binding. Front Mol Biosci 2022; 9:858316. [PMID: 35558558 PMCID: PMC9086317 DOI: 10.3389/fmolb.2022.858316] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Accepted: 04/01/2022] [Indexed: 11/22/2022] Open
Abstract
Improper reaction coordinates can pose significant problems for path-based binding free energy calculations. Particularly, omission of long timescale motions can lead to over-estimation of the energetic barriers between the bound and unbound states. Many methods exist to construct the optimal reaction coordinate using a pre-defined basis set of features. Although simulations are typically conducted in explicit solvent, the solvent atoms are often excluded by these feature sets—resulting in little being known about their role in reaction coordinates, and ultimately, their role in determining (un)binding rates and free energies. In this work, analysis is done on an extensive set of host-guest unbinding trajectories, working to characterize differences between high and low probability unbinding trajectories with a focus on solvent-based features, including host-ion interactions, guest-ion interactions and location-dependent ion densities. We find that differences in ion densities as well as guest-ion interactions strongly correlate with differences in the probabilities of reactive paths that are used to determine free energies of (un)binding and play a significant role in the unbinding process.
Collapse
Affiliation(s)
- Nicole M. Roussey
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, United States
| | - Alex Dickson
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, United States
- Department of Computational Mathematics, Science, and Engineering, Michigan State University, East Lansing, MI, United States
- *Correspondence: Alex Dickson,
| |
Collapse
|
19
|
Russo JD, Zhang S, Leung JMG, Bogetti AT, Thompson JP, DeGrave AJ, Torrillo PA, Pratt AJ, Wong KF, Xia J, Copperman J, Adelman JL, Zwier MC, LeBard DN, Zuckerman DM, Chong LT. WESTPA 2.0: High-Performance Upgrades for Weighted Ensemble Simulations and Analysis of Longer-Timescale Applications. J Chem Theory Comput 2022; 18:638-649. [PMID: 35043623 PMCID: PMC8825686 DOI: 10.1021/acs.jctc.1c01154] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The weighted ensemble (WE) family of methods is one of several statistical mechanics-based path sampling strategies that can provide estimates of key observables (rate constants and pathways) using a fraction of the time required by direct simulation methods such as molecular dynamics or discrete-state stochastic algorithms. WE methods oversee numerous parallel trajectories using intermittent overhead operations at fixed time intervals, enabling facile interoperability with any dynamics engine. Here, we report on the major upgrades to the WESTPA software package, an open-source, high-performance framework that implements both basic and recently developed WE methods. These upgrades offer substantial improvements over traditional WE methods. The key features of the new WESTPA 2.0 software enhance the efficiency and ease of use: an adaptive binning scheme for more efficient surmounting of large free energy barriers, streamlined handling of large simulation data sets, exponentially improved analysis of kinetics, and developer-friendly tools for creating new WE methods, including a Python API and resampler module for implementing both binned and "binless" WE strategies.
Collapse
Affiliation(s)
- John D Russo
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - She Zhang
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Jeremy M G Leung
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Anthony T Bogetti
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Jeff P Thompson
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Alex J DeGrave
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Paul A Torrillo
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - A J Pratt
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kim F Wong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Junchao Xia
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Jeremy Copperman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - Joshua L Adelman
- Department of Biological Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Matthew C Zwier
- Department of Chemistry, Drake University, Des Moines, Iowa 50311-4505, United States
| | - David N LeBard
- OpenEye Scientific, Santa Fe, New Mexico 87508, United States
| | - Daniel M Zuckerman
- Department of Biomedical Engineering, Oregon Health and Science University, Portland, Oregon 97239-3098, United States
| | - Lillian T Chong
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
20
|
Dixon T, Lotz SD, Dickson A. Creating Maps of the Ligand Binding Landscape for Kinetics-Based Drug Discovery. Methods Mol Biol 2022; 2385:325-334. [PMID: 34888727 DOI: 10.1007/978-1-0716-1767-0_15] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Simulations of ligand-protein interactions can be very useful for drug design and to gain biological insight. Full pathways of ligand-protein binding can be used to get information about ligand binding transition states, which form the rate-limiting step of the binding and release processes. However, these simulations are typically limited by the presence of large energy barriers that separate stable poses of interest. Here we describe a simulation protocol for exploring and analyzing landscapes of ligand-protein interactions that makes use of molecular docking, enhanced molecular simulation with the weighted ensemble algorithm, and network analysis. It can be accomplished using a modest cluster of graphics processing units and freely accessible software. This protocol focuses on the construction and analysis of a network model of ligand binding poses and provides links to resources that describe the other steps in more detail. The end result of this protocol is a map of the ligand-protein binding landscape that identifies transition states of the ligand binding pathway, as well as alternative bound poses that could be stabilized with modifications to the ligand.
Collapse
Affiliation(s)
- Tom Dixon
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI, USA
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI, USA
| | - Samuel D Lotz
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI, USA
- Roivant Sciences, New York, NY, USA
| | - Alex Dickson
- Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, MI, USA.
- Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI, USA.
| |
Collapse
|
21
|
Gomez D, Peña Ccoa WJ, Singh Y, Rojas E, Hocky GM. Molecular Paradigms for Biological Mechanosensing. J Phys Chem B 2021; 125:12115-12124. [PMID: 34709040 DOI: 10.1021/acs.jpcb.1c06330] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many proteins in living cells are subject to mechanical forces, which can be generated internally by molecular machines, or externally, e.g., by pressure gradients. In general, these forces fall in the piconewton range, which is similar in magnitude to forces experienced by a molecule due to thermal fluctuations. While we would naively expect such moderate forces to produce only minimal changes, a wide variety of "mechanosensing" proteins have evolved with functions that are responsive to forces in this regime. The goal of this article is to provide a physical chemistry perspective on protein-based molecular mechanosensing paradigms used in living systems, and how these paradigms can be explored using novel computational methods.
Collapse
Affiliation(s)
- David Gomez
- Department of Biology, New York University, New York, New York 10003, United States.,Department of Chemistry, New York University, New York, New York 10003, United States
| | - Willmor J Peña Ccoa
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Yuvraj Singh
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Enrique Rojas
- Department of Biology, New York University, New York, New York 10003, United States
| | - Glen M Hocky
- Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
22
|
Nunes-Alves A, Kokh DB, Wade RC. Ligand unbinding mechanisms and kinetics for T4 lysozyme mutants from τRAMD simulations. Curr Res Struct Biol 2021; 3:106-111. [PMID: 34235490 PMCID: PMC8244441 DOI: 10.1016/j.crstbi.2021.04.001] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2021] [Revised: 04/20/2021] [Accepted: 04/25/2021] [Indexed: 01/14/2023] Open
Abstract
The protein-ligand residence time, τ, influences molecular function in biological networks and has been recognized as an important determinant of drug efficacy. To predict τ, computational methods must overcome the problem that τ often exceeds the timescales accessible to conventional molecular dynamics (MD) simulation. Here, we apply the τ-Random Acceleration Molecular Dynamics (τRAMD) method to a set of kinetically characterized complexes of T4 lysozyme mutants with small, engineered binding cavities. τRAMD yields relative ligand dissociation rates in good accordance with experiments across this diverse set of complexes that differ with regard to measurement temperature, ligand identity, protein mutation and binding cavity. τRAMD thereby allows a comprehensive characterization of the ligand egress routes and determinants of τ. Although ligand dissociation by multiple egress routes is observed, we find that egress via the predominant route determines the value of τ. We also find that the presence of a greater number of metastable states along egress pathways leads to slower protein-ligand dissociation. These physical insights could be exploited in the rational optimization of the kinetic properties of drug candidates. Relative residence times are computed for T4 lysozyme mutant-ligand complexes. τ-Random Acceleration Molecular Dynamics provide efficient sampling of unbinding. Computed dissociation rates show good agreement with all available measured values. Ligand egress via the predominant route determines the value of the residence time. The presence of metastable states along egress pathways slows down dissociation.
Collapse
Affiliation(s)
- Ariane Nunes-Alves
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120, Heidelberg, Germany
| | - Daria B Kokh
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany.,Center for Molecular Biology (ZMBH), DKFZ-ZMBH Alliance, Heidelberg University, Im Neuenheimer Feld 282, 69120, Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Im Neuenheimer Feld 205, Heidelberg, Germany
| |
Collapse
|