1
|
Badaoui M, Buigues PJ, Berta D, Mandana GM, Gu H, Földes T, Dickson CJ, Hornak V, Kato M, Molteni C, Parsons S, Rosta E. Combined Free-Energy Calculation and Machine Learning Methods for Understanding Ligand Unbinding Kinetics. J Chem Theory Comput 2022; 18:2543-2555. [PMID: 35195418 PMCID: PMC9097281 DOI: 10.1021/acs.jctc.1c00924] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
![]()
The
determination of drug residence times, which define the time
an inhibitor is in complex with its target, is a fundamental part
of the drug discovery process. Synthesis and experimental measurements
of kinetic rate constants are, however, expensive and time consuming.
In this work, we aimed to obtain drug residence times computationally.
Furthermore, we propose a novel algorithm to identify molecular design
objectives based on ligand unbinding kinetics. We designed an enhanced
sampling technique to accurately predict the free-energy profiles
of the ligand unbinding process, focusing on the free-energy barrier
for unbinding. Our method first identifies unbinding paths determining
a corresponding set of internal coordinates (ICs) that form contacts
between the protein and the ligand; it then iteratively updates these
interactions during a series of biased molecular dynamics (MD) simulations
to reveal the ICs that are important for the whole of the unbinding
process. Subsequently, we performed finite-temperature string simulations
to obtain the free-energy barrier for unbinding using the set of ICs
as a complex reaction coordinate. Importantly, we also aimed to enable
the further design of drugs focusing on improved residence times.
To this end, we developed a supervised machine learning (ML) approach
with inputs from unbiased “downhill” trajectories initiated
near the transition state (TS) ensemble of the string unbinding path.
We demonstrate that our ML method can identify key ligand–protein
interactions driving the system through the TS. Some of the most important
drugs for cancer treatment are kinase inhibitors. One of these kinase
targets is cyclin-dependent kinase 2 (CDK2), an appealing target for
anticancer drug development. Here, we tested our method using two
different CDK2 inhibitors for the potential further development of
these compounds. We compared the free-energy barriers obtained from
our calculations with those observed in available experimental data.
We highlighted important interactions at the distal ends of the ligands
that can be targeted for improved residence times. Our method provides
a new tool to determine unbinding rates and to identify key structural
features of the inhibitors that can be used as starting points for
novel design strategies in drug discovery.
Collapse
Affiliation(s)
- Magd Badaoui
- Department of Chemistry, King's College London, London SE1 1DB, United Kingdom.,Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Pedro J Buigues
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Dénes Berta
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Gaurav M Mandana
- Department of Chemistry, King's College London, London SE1 1DB, United Kingdom
| | - Hankang Gu
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Tamás Földes
- Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Callum J Dickson
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Viktor Hornak
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Mitsunori Kato
- Computer-Aided Drug Discovery, Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Carla Molteni
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Simon Parsons
- School of Computer Science, University of Lincoln, Lincoln LN6 7TS, United Kingdom
| | - Edina Rosta
- Department of Chemistry, King's College London, London SE1 1DB, United Kingdom.,Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
2
|
Jones PE, Pérez-Segura C, Bryer AJ, Perilla JR, Hadden-Perilla JA. Molecular dynamics of the viral life cycle: progress and prospects. Curr Opin Virol 2021; 50:128-138. [PMID: 34464843 PMCID: PMC8651149 DOI: 10.1016/j.coviro.2021.08.003] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Revised: 08/09/2021] [Accepted: 08/09/2021] [Indexed: 01/29/2023]
Abstract
Molecular dynamics (MD) simulations across spatiotemporal resolutions are widely applied to study viruses and represent the central technique uniting the field of computational virology. We discuss the progress of MD in elucidating the dynamics of the viral life cycle, including the status of modeling intact extracellular virions and leveraging advanced simulations to mimic active life cycle processes. We further remark on the prospects of MD for continued contributions to the basic science characterization of viruses, especially given the increasing availability of high-quality experimental data and supercomputing power. Overall, integrative computational methods that are closely guided by experiments are unmatched in the level of detail they provide, enabling-now and in the future-new discoveries relevant to thwarting viral infection.
Collapse
Affiliation(s)
- Peter Eugene Jones
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Carolina Pérez-Segura
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Alexander J Bryer
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Jodi A Hadden-Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States.
| |
Collapse
|
3
|
Abstract
The protein HIV Reverse Transcriptase (HIV RT) synthesizes a DNA strand according to a template. During the synthesis, the polymerase slides on the double stranded DNA to allow the entry of a new nucleotide to the active site. We use Molecular Dynamics simulations to estimate the free energy profile and the time scale of the DNA-protein's relative displacement in the complex's closed state. We illustrate that the presence of the catalytic magnesium slows down the process. Upon removing the catalytic magnesium ion, the process is rapid and significantly faster than reopening the active site in preparation for the new substrate. We speculate that magnesium regulates DNA translocation. The magnesium locks the DNA into a specific orientation during the chemical addition of the nucleotide. The release of Mg2+ eases DNA sliding and the acceptance of a new substrate.
Collapse
Affiliation(s)
- Hao Wang
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Ron Elber
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States.,Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
4
|
Alberini G, Benfenati F, Maragliano L. Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:1354-1367. [PMID: 33570938 PMCID: PMC8023575 DOI: 10.1021/acs.jcim.0c01407] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Activation of voltage-gated
ion channels is regulated by conformational
changes of the voltage sensor domains (VSDs), four water- and ion-impermeable
modules peripheral to the central, permeable pore domain. Anomalous
currents, defined as ω-currents, have been recorded in response
to mutations of residues on the VSD S4 helix and associated with ion
fluxes through the VSDs. In humans, gene defects in the potassium
channel Kv7.2 result in a broad range of epileptic disorders, from
benign neonatal seizures to severe epileptic encephalopathies. Experimental
evidence suggests that the R207Q mutation in S4, associated with peripheral
nerve hyperexcitability, induces ω-currents at depolarized potentials,
but the fine structural details are still elusive. In this work, we
use atom-detailed molecular dynamics simulations and a refined model
structure of the Kv7.2 VSD in the active conformation in a membrane/water
environment to study the effect of R207Q and four additional mutations
of proven clinical importance. Our results demonstrate that the R207Q
mutant shows the most pronounced increase of hydration in the internal
VSD cavity, a feature favoring the occurrence of ω-currents.
Free energy and kinetics calculations of sodium permeation through
the native and mutated VSD indicate as more favorable the formation
of a cationic current in the latter. Overall, our simulations establish
a mechanistic linkage between genetic variations and their physiological
outcome, by providing a computational description that includes both
thermodynamic and kinetic features of ion permeation associated with
ω-currents.
Collapse
Affiliation(s)
- Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Experimental Medicine, Università degli Studi di Genova, Viale Benedetto XV, 3, 16132 Genova, Italy
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi, 10, 16132 Genova, Italy
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| |
Collapse
|
5
|
Elber R, Fathizadeh A, Ma P, Wang H. Modeling molecular kinetics with Milestoning. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1512] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ron Elber
- Department of Chemistry, The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Arman Fathizadeh
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| | - Piao Ma
- Department of Chemistry University of Texas at Austin Austin Texas USA
| | - Hao Wang
- The Oden Institute for Computational Engineering and Sciences University of Texas at Austin Austin Texas USA
| |
Collapse
|
6
|
Dangerfield TL, Johnson KA. Optimized incorporation of an unnatural fluorescent amino acid affords measurement of conformational dynamics governing high-fidelity DNA replication. J Biol Chem 2020; 295:17265-17280. [PMID: 33020184 PMCID: PMC7863912 DOI: 10.1074/jbc.ra120.015557] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Revised: 09/29/2020] [Indexed: 12/20/2022] Open
Abstract
DNA polymerase from bacteriophage T7 undergoes large, substrate-induced conformational changes that are thought to account for high replication fidelity, but prior studies were adversely affected by mutations required to construct a Cys-lite variant needed for site-specific fluorescence labeling. Here we have optimized the direct incorporation of a fluorescent un-natural amino acid, (7-hydroxy-4-coumarin-yl)-ethylglycine, using orthogonal amber suppression machinery in Escherichia coli MS methods verify that the unnatural amino acid is only incorporated at one position with minimal background. We show that the single fluorophore provides a signal to detect nucleotide-induced conformational changes through equilibrium and stopped-flow kinetic measurements of correct nucleotide binding and incorporation. Pre-steady-state chemical quench methods show that the kinetics and fidelity of DNA replication catalyzed by the labeled enzyme are largely unaffected by the unnatural amino acid. These advances enable rigorous analysis to establish the kinetic and mechanistic basis for high-fidelity DNA replication.
Collapse
Affiliation(s)
- Tyler L Dangerfield
- Institute for Cellular and Molecular Biology, Department of Molecular Biosciences, University of Texas, Austin, Texas, USA
| | - Kenneth A Johnson
- Institute for Cellular and Molecular Biology, Department of Molecular Biosciences, University of Texas, Austin, Texas, USA.
| |
Collapse
|