scholarly journals Assessment of Software Methods for Estimating Protein-Protein Relative Binding Affinities

2020 ◽  
Author(s):  
Tawny R. Gonzalez ◽  
Kyle P. Martin ◽  
Jonathan E. Barnes ◽  
Jagdish Suresh Patel ◽  
F. Marty Ytreberg
Keyword(s):  
Molecular Dynamics ◽  
Computational Cost ◽  
Structural Features ◽  
New Drugs ◽  
Binding Affinities ◽  
Evolutionary Mechanisms ◽  
Relative Binding ◽  
Amino Acid Mutations ◽  
Dynamics Simulations ◽  
The Impact

AbstractA growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < −0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

scholarly journals Assessment of software methods for estimating protein-protein relative binding affinities

PLoS ONE ◽  
2020 ◽  
Vol 15 (12) ◽  
pp. e0240573
Author(s):  
Tawny R. Gonzalez ◽  
Kyle P. Martin ◽  
Jonathan E. Barnes ◽  
Jagdish Suresh Patel ◽  
F. Marty Ytreberg
Keyword(s):  
Molecular Dynamics ◽  
Computational Cost ◽  
Structural Features ◽  
New Drugs ◽  
Binding Affinities ◽  
Evolutionary Mechanisms ◽  
Relative Binding ◽  
Amino Acid Mutations ◽  
Dynamics Simulations ◽  
The Impact

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83–98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

scholarly journals Clustering Molecular Dynamics Trajectories for Optimizing Docking Experiments

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Renata De Paris ◽  
Christian V. Quevedo ◽  
Duncan D. Ruiz ◽  
Osmar Norberto de Souza ◽  
Rodrigo C. Barros
Keyword(s):  
Molecular Dynamics ◽  
Substrate Binding ◽  
Clustering Algorithms ◽  
Computational Cost ◽  
Structural Features ◽  
Dynamics Simulation ◽  
Rational Drug ◽  
Protein Receptors ◽  
Binding Cavity ◽  
Dynamics Simulations

Molecular dynamics simulations of protein receptors have become an attractive tool for rational drug discovery. However, the high computational cost of employing molecular dynamics trajectories in virtual screening of large repositories threats the feasibility of this task. Computational intelligence techniques have been applied in this context, with the ultimate goal of reducing the overall computational cost so the task can become feasible. Particularly, clustering algorithms have been widely used as a means to reduce the dimensionality of molecular dynamics trajectories. In this paper, we develop a novel methodology for clustering entire trajectories using structural features from the substrate-binding cavity of the receptor in order to optimize docking experiments on a cloud-based environment. The resulting partition was selected based on three clustering validity criteria, and it was further validated by analyzing the interactions between 20 ligands and a fully flexible receptor (FFR) model containing a 20 ns molecular dynamics simulation trajectory. Our proposed methodology shows that taking into account features of the substrate-binding cavity as input for thek-means algorithm is a promising technique for accurately selecting ensembles of representative structures tailored to a specific ligand.

Parameterization of Divalent Cations for Coarse-Grained Simulations

2020 ◽  
Author(s):  
Florencia Klein ◽  
Daniela Cáceres-Rojas ◽  
Monica Carrasco ◽  
Juan Carlos Tapia ◽  
Julio Caballero ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Metal Ions ◽  
Molecular Dynamics Simulations ◽  
Divalent Cations ◽  
Computational Cost ◽  
Data Bank ◽  
Coarse Grained ◽  
Interaction Parameters ◽  
Dynamics Simulations ◽  
Dynamical Description

<p>Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges to achieve an accurate structural and dynamical description of many biological assemblies. This is particularly the case for coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or simply not available for the vast majority of CG- force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations. Using this simple approach, we provide a set of interaction parameters for Calcium, Magnesium, and Zinc ions, which cover more than 80% of the metal-bound structures reported on the PDB. Simulations performed using the SIRAH force field on several proteins and DNA systems show that using the present approach it is possible to obtain non-bonded interaction parameters that obviate the use of topological constraints. </p>

Active Learning Accelerates Ab Initio Molecular Dynamics on Pericyclic Reactive Energy Surfaces

2020 ◽  
Author(s):  
Shi Jun Ang ◽  
Wujie Wang ◽  
Daniel Schwalbe-Koda ◽  
Simon Axelrod ◽  
Rafael Gomez-Bombarelli
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Transition State ◽  
Potential Energy Surface ◽  
Ab Initio ◽  
Energy Surface ◽  
Computational Cost ◽  
Reactive Trajectories ◽  
Dynamics Simulations ◽  
Energy Surfaces

<div>Modeling dynamical effects in chemical reactions, such as post-transition state bifurcation, requires <i>ab initio</i> molecular dynamics simulations due to the breakdown of simpler static models like transition state theory. However, these simulations tend to be restricted to lower-accuracy electronic structure methods and scarce sampling because of their high computational cost. Here, we report the use of statistical learning to accelerate reactive molecular dynamics simulations by combining high-throughput ab initio calculations, graph-convolution interatomic potentials and active learning. This pipeline was demonstrated on an ambimodal trispericyclic reaction involving 8,8-dicyanoheptafulvene and 6,6-dimethylfulvene. With a dataset size of approximately</div><div>31,000 M062X/def2-SVP quantum mechanical calculations, the computational cost of exploring the reactive potential energy surface was reduced by an order of magnitude. Thousands of virtually costless picosecond-long reactive trajectories suggest that post-transition state bifurcation plays a minor role for the reaction in vacuum. Furthermore, a transfer-learning strategy effectively upgraded the potential energy surface to higher</div><div>levels of theory ((SMD-)M06-2X/def2-TZVPD in vacuum and three other solvents, as well as the more accurate DLPNO-DSD-PBEP86 D3BJ/def2-TZVPD) using about 10% additional calculations for each surface. Since the larger basis set and the dynamic correlation capture intramolecular non-covalent interactions more accurately, they uncover longer lifetimes for the charge-separated intermediate on the more accurate potential energy surfaces. The character of the intermediate switches from entropic to thermodynamic upon including implicit solvation effects, with lifetimes increasing with solvent polarity. Analysis of 2,000 reactive trajectories on the chloroform PES shows a qualitative agreement with the experimentally-reported periselectivity for this reaction. This overall approach is broadly applicable and opens a door to the study of dynamical effects in larger, previously-intractable reactive systems.</div>

Quantitative Ranking of Ligand Binding Kinetics with a Multiscale Milestoning Simulation Approach

2018 ◽  
Author(s):  
Benjamin R. Jagger ◽  
Christoper T. Lee ◽  
Rommie Amaro
Keyword(s):  
Molecular Dynamics ◽  
Drug Discovery ◽  
Small Molecule ◽  
Brownian Dynamics ◽  
Model Simulation ◽  
Computational Cost ◽  
Lead Selection ◽  
Delta G ◽  
Dynamics Simulations

<p>The ranking of small molecule binders by their kinetic (kon and koff) and thermodynamic (delta G) properties can be a valuable metric for lead selection and optimization in a drug discovery campaign, as these quantities are often indicators of in vivo efficacy. Efficient and accurate predictions of these quantities can aid the in drug discovery effort, acting as a screening step. We have previously described a hybrid molecular dynamics, Brownian dynamics, and milestoning model, Simulation Enabled Estimation of Kinetic Rates (SEEKR), that can predict kon’s, koff’s, and G’s. Here we demonstrate the effectiveness of this approach for ranking a series of seven small molecule compounds for the model system, -cyclodextrin, based on predicted kon’s and koff’s. We compare our results using SEEKR to experimentally determined rates as well as rates calculated using long-timescale molecular dynamics simulations and show that SEEKR can effectively rank the compounds by koff and G with reduced computational cost. We also provide a discussion of convergence properties and sensitivities of calculations with SEEKR to establish “best practices” for its future use.</p>

Design of Inhibitors for Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH) Enzyme of Leishmania mexicana

2020 ◽  
Vol 16 (6) ◽  
pp. 784-795
Author(s):  
Krisnna M.A. Alves ◽  
Fábio José Bonfim Cardoso ◽  
Kathia M. Honorio ◽  
Fábio A. de Molfetta
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Point Of View ◽  
New Drugs ◽  
Leishmania Mexicana ◽  
Receptor Complexes ◽  
Starting Point ◽  
Dynamics Simulations ◽  
Selection Of

Background:: Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease. Objective:: The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana). Methods: A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations. Results:: Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity. Conclusion:: he use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.

Molecular Dynamics Study of Ion Impact Phenomena and Compressive Stress in Thin Films

1993 ◽  
Vol 317 ◽  
Author(s):  
N.A. Marks ◽  
P. Guan ◽  
D.R. Mckenzie ◽  
B.A. PailThorpe
Keyword(s):  
Molecular Dynamics ◽  
Three Dimensional ◽  
Carbon Films ◽  
Loss Mechanism ◽  
Ion Energy ◽  
Lennard Jones ◽  
Impact Phenomena ◽  
Ion Impact ◽  
Dynamics Simulations ◽  
The Impact

ABSTRACTMolecular dynamics simulations of nickel and carbon have been used to study the phenomena due to ion impact. The nickel and carbon interactions were described using the Lennard-Jones and Stillinger-Weber potentials respectively. The phenomena occurring after the impact of 100 e V to 1 keV ions were studied in the nickel simulations, which were both two and three-dimensional. Supersonic focussed collision sequences (or focusons) were observed, and associated with these focusons were unexpected sonic bow waves, which were a major energy loss mechanism for the focuson. A number of 2D carbon films were grown and the stress in the films as a function of incident ion energy was Measured. With increasing energy the stress changed from tensile to compressive and reached a maximum around 50 eV, in agreement with experiment.

Cluster, Surface and Bulk Properties of ZnCd Binary Alloys: Molecular-Dynamics Simulations

2005 ◽  
Vol 502 ◽  
pp. 51-56 ◽  
Author(s):  
Sakir Erkoc
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Density Functional ◽  
Structural Features ◽  
Bulk Materials ◽  
Bulk Properties ◽  
Atom Clusters ◽  
Dynamics Simulations ◽  
Mixed Clusters ◽  
Stable Structures

The structural and electronic properties of isolated neutral ZnmCdn clusters for m+n £ 3 have been investigated by performing density functional theory calculations at B3LYP level. The optimum geometries, vibrational frequencies, electronic structures, and the possible dissosiation channels of the clusters considered have been obtained. An empirical many-body potential energy function (PEF), which comprices two- and three-body atomic interactions, has been developed to investigate the structural features and energetics of ZnmCdn (m+n=3,4) microclusters. The most stable structures were found to be triangular for the three-atom clusters and tetrahedral for the four-atom clusters. On the other hand, the structural features and energetics of Znn-mCdm (n=7,8) microclusters, and Zn50, Cd50, Zn25Cd25, Zn12Cd38, and Zn38Cd12 nanoparticles have been investigated by performing molecular-dynamics computer simulations using the developed PEF. The most stable structures were found to be compact and three-dimensional for all elemental and mixed clusters. An interesting structural feature of the mixed clusters is that Zn and Cd atoms do not mix in mixed clusters, they come together almost without mixing. Surface and bulk properties of Zn, Cd, and ZnCd systems have been investigated too by performing molecular-dynamics simulations using the developed PEF. Surface reconstruction and multilayer relaxation on clean surfaces, adatom on surface, substitutional atom on surface and bulk materials, and vacancy on surface and bulk materials have been studied extensively.

Sign in / Sign up

Export Citation Format

Share Document