1
|
Skoptsova AA, Geronikaki A, Novichikhina NP, Sulimov AV, Ilin IS, Sulimov VB, Bykov GA, Podoplelova NA, Pyankov OV, Shikhaliev KS. Design, Synthesis, and Evaluation of New Hybrid Derivatives of 5,6-Dihydro-4 H-pyrrolo[3,2,1- ij]quinolin-2(1 H)-one as Potential Dual Inhibitors of Blood Coagulation Factors Xa and XIa. Molecules 2024; 29:373. [PMID: 38257286 PMCID: PMC10818416 DOI: 10.3390/molecules29020373] [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: 10/30/2023] [Revised: 12/29/2023] [Accepted: 01/08/2024] [Indexed: 01/24/2024] Open
Abstract
Cardiovascular diseases caused by blood coagulation system disorders are one of the leading causes of morbidity and mortality in the world. Research shows that blood clotting factors are involved in these thrombotic processes. Among them, factor Xa occupies a key position in the blood coagulation cascade. Another coagulation factor, XIa, is also a promising target because its inhibition can suppress thrombosis with a limited contribution to normal hemostasis. In this regard, the development of dual inhibitors as new generation anticoagulants is an urgent problem. Here we report the synthesis and evaluation of novel potential dual inhibitors of coagulation factors Xa and XIa. Based on the principles of molecular design, we selected a series of compounds that combine in their structure fragments of pyrrolo[3,2,1-ij]quinolin-2-one and thiazole, connected through a hydrazine linker. The production of new hybrid molecules was carried out using a two-stage method. The reaction of 5,6-dihydropyrrolo[3,2,1-ij]quinoline-1,2-diones with thiosemicarbazide gave the corresponding hydrazinocarbothioamides. The reaction of the latter with DMAD led to the target methyl 2-(4-oxo-2-(2-(2-oxo-5,6-dihydro-4H-pyrrolo[3,2,1-ij]quinolin-1(2H)-ylidene)hydrazineyl)thiazol-5(4H)-ylidene)acetates in high yields. In vitro testing of the synthesized molecules revealed that ten of them showed high inhibition values for both the coagulation factors Xa and XIa, and the IC50 value for some compounds was also assessed. The resulting structures were also tested for their ability to inhibit thrombin.
Collapse
Affiliation(s)
- Anna A. Skoptsova
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia; (A.A.S.); (N.P.N.)
| | - Athina Geronikaki
- School of Pharmacy, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
| | - Nadezhda P. Novichikhina
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia; (A.A.S.); (N.P.N.)
| | - Alexey V. Sulimov
- Research Computing Center, Lomonosov Moscow State University, 119992 Moscow, Russia; (A.V.S.); (I.S.I.); (V.B.S.)
| | - Ivan S. Ilin
- Research Computing Center, Lomonosov Moscow State University, 119992 Moscow, Russia; (A.V.S.); (I.S.I.); (V.B.S.)
| | - Vladimir B. Sulimov
- Research Computing Center, Lomonosov Moscow State University, 119992 Moscow, Russia; (A.V.S.); (I.S.I.); (V.B.S.)
| | - Georgii A. Bykov
- Department of Biophysics at the Faculty of Physics, Lomonosov Moscow State University, 119992 Moscow, Russia;
| | | | - Oleg V. Pyankov
- State Research Center of Virology and Biotechnology “Vector”, 630559 Koltsovo, Russia;
| | - Khidmet S. Shikhaliev
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia; (A.A.S.); (N.P.N.)
| |
Collapse
|
2
|
Sulimov A, Ilin I, Kutov D, Shikhaliev K, Shcherbakov D, Pyankov O, Stolpovskaya N, Medvedeva S, Sulimov V. New Chemicals Suppressing SARS-CoV-2 Replication in Cell Culture. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27175732. [PMID: 36080498 PMCID: PMC9457583 DOI: 10.3390/molecules27175732] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Revised: 08/30/2022] [Accepted: 08/30/2022] [Indexed: 01/18/2023]
Abstract
Candidates to being inhibitors of the main protease (Mpro) of SARS-CoV-2 were selected from the database of Voronezh State University using molecular modeling. The database contained approximately 19,000 compounds represented by more than 41,000 ligand conformers. These ligands were docked into Mpro using the SOL docking program. For one thousand ligands with best values of the SOL score, the protein–ligand binding enthalpy was calculated by the PM7 quantum-chemical method with the COSMO solvent model. Using the SOL score and the calculated protein–ligand binding enthalpies, eighteen compounds were selected for the experiments. Several of these inhibitors suppressed the replication of the coronavirus in cell culture, and we used the best three among them in the search for chemical analogs. Selection among analogs using the same procedure followed by experiments led to identification of seven inhibitors of the SARS-CoV-2 replication in cell culture with EC50 values at the micromolar level. The identified inhibitors belong to three chemical classes. The three inhibitors, 4,4-dimethyldithioquinoline derivatives, inhibit SARS-CoV-2 replication in Vero E6 cell culture just as effectively as the best published non-covalent inhibitors, and show low cytotoxicity. These results open up a possibility to develop antiviral drugs against the SARS-CoV-2 coronavirus.
Collapse
Affiliation(s)
- Alexey Sulimov
- Dimonta Ltd., 15 Nagornaya Str., Bldg 8, 117186 Moscow, Russia
- Research Computing Center, Lomonosov Moscow State University, Leninskie Gory, 1, Building 4, 119234 Moscow, Russia
| | - Ivan Ilin
- Dimonta Ltd., 15 Nagornaya Str., Bldg 8, 117186 Moscow, Russia
- Research Computing Center, Lomonosov Moscow State University, Leninskie Gory, 1, Building 4, 119234 Moscow, Russia
| | - Danil Kutov
- Dimonta Ltd., 15 Nagornaya Str., Bldg 8, 117186 Moscow, Russia
- Research Computing Center, Lomonosov Moscow State University, Leninskie Gory, 1, Building 4, 119234 Moscow, Russia
- Correspondence: (D.K.); (V.S.)
| | - Khidmet Shikhaliev
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia
| | - Dmitriy Shcherbakov
- State Research Centre of Virology and Biotechnology “Vector”, 630559 Koltsovo, Russia
| | - Oleg Pyankov
- State Research Centre of Virology and Biotechnology “Vector”, 630559 Koltsovo, Russia
| | - Nadezhda Stolpovskaya
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia
| | - Svetlana Medvedeva
- Department of Organic Chemistry, Faculty of Chemistry, Voronezh State University, 1 Universitetskaya Sq., 394018 Voronezh, Russia
| | - Vladimir Sulimov
- Dimonta Ltd., 15 Nagornaya Str., Bldg 8, 117186 Moscow, Russia
- Research Computing Center, Lomonosov Moscow State University, Leninskie Gory, 1, Building 4, 119234 Moscow, Russia
- Correspondence: (D.K.); (V.S.)
| |
Collapse
|
3
|
Sulimov A, Kutov D, Ilin I, Sulimov V. Quantum-Chemical Quasi-Docking for Molecular Dynamics Calculations. NANOMATERIALS (BASEL, SWITZERLAND) 2022; 12:274. [PMID: 35055291 PMCID: PMC8781293 DOI: 10.3390/nano12020274] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 01/11/2022] [Accepted: 01/13/2022] [Indexed: 01/14/2023]
Abstract
The quantum quasi-docking procedure is used to compare the docking accuracies of two quantum-chemical semiempirical methods, namely, PM6-D3H4X and PM7. Quantum quasi-docking is an approximation to quantum docking. In quantum docking, it is necessary to search directly for the global minimum of the energy of the protein-ligand complex calculated by the quantum-chemical method. In quantum quasi-docking, firstly, we look for a wide spectrum of low-energy minima, calculated using the MMFF94 force field, and secondly, we recalculate the energies of all these minima using the quantum-chemical method, and among these recalculated energies we determine the lowest energy and the corresponding ligand position. Both PM6-D3H4X and PM7 are novel methods that describe well-dispersion interactions, hydrogen and halogen bonds. The PM6-D3H4X and PM7 methods are used with the COSMO implicit solvent model as it is implemented in the MOPAC program. The comparison is made for 25 high quality protein-ligand complexes. Firstly, the docking positioning accuracies have been compared, and we demonstrated that PM7+COSMO provides better positioning accuracy than PM6-D3H4X. Secondly, we found that PM7+COSMO demonstrates a much higher correlation between the calculated and measured protein-ligand binding enthalpies than PM6-D3H4X. For future quantum docking PM7+COSMO is preferable, but the COSMO model must be improved.
Collapse
Affiliation(s)
- Alexey Sulimov
- Dimonta, Ltd., 117186 Moscow, Russia; (A.S.); (I.I.)
- Research Computer Center, Lomonosov Moscow State University, 119992 Moscow, Russia
| | - Danil Kutov
- Dimonta, Ltd., 117186 Moscow, Russia; (A.S.); (I.I.)
- Research Computer Center, Lomonosov Moscow State University, 119992 Moscow, Russia
| | - Ivan Ilin
- Dimonta, Ltd., 117186 Moscow, Russia; (A.S.); (I.I.)
- Research Computer Center, Lomonosov Moscow State University, 119992 Moscow, Russia
| | - Vladimir Sulimov
- Dimonta, Ltd., 117186 Moscow, Russia; (A.S.); (I.I.)
- Research Computer Center, Lomonosov Moscow State University, 119992 Moscow, Russia
| |
Collapse
|
4
|
Abstract
Molecular Docking is used to positioning the computer-generated 3D structure of small
ligands into a receptor structure in a variety of orientations, conformations and positions. This
method is useful in drug discovery and medicinal chemistry providing insights into molecular
recognition. Docking has become an integral part of Computer-Aided Drug Design and Discovery
(CADDD). Traditional docking methods suffer from limitations of semi-flexible or static treatment
of targets and ligand. Over the last decade, advances in the field of computational, proteomics and
genomics have also led to the development of different docking methods which incorporate
protein-ligand flexibility and their different binding conformations. Receptor flexibility accounts
for more accurate binding pose predictions and a more rational depiction of protein binding
interactions with the ligand. Protein flexibility has been included by generating protein ensembles
or by dynamic docking methods. Dynamic docking considers solvation, entropic effects and also
fully explores the drug-receptor binding and recognition from both energetic and mechanistic point
of view. Though in the fast-paced drug discovery program, dynamic docking is computationally
expensive but is being progressively used for screening of large compound libraries to identify the
potential drugs. In this review, a quick introduction is presented to the available docking methods
and their application and limitations in drug discovery.
Collapse
Affiliation(s)
- Ritu Jakhar
- Center for Bioinformatics, Maharshi Dayanand University, Rohtak, India
| | - Mehak Dangi
- Center for Bioinformatics, Maharshi Dayanand University, Rohtak, India
| | - Alka Khichi
- Center for Bioinformatics, Maharshi Dayanand University, Rohtak, India
| | | |
Collapse
|
5
|
Sulimov VB, Kutov DC, Sulimov AV. Advances in Docking. Curr Med Chem 2020; 26:7555-7580. [PMID: 30182836 DOI: 10.2174/0929867325666180904115000] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Revised: 07/04/2018] [Accepted: 07/06/2018] [Indexed: 11/22/2022]
Abstract
BACKGROUND Design of small molecules which are able to bind to the protein responsible for a disease is the key step of the entire process of the new medicine discovery. Atomistic computer modeling can significantly improve effectiveness of such design. The accurate calculation of the free energy of binding a small molecule (a ligand) to the target protein is the most important problem of such modeling. Docking is one of the most popular molecular modeling methods for finding ligand binding poses in the target protein and calculating the protein-ligand binding energy. This energy is used for finding the most active compounds for the given target protein. This short review aims to give a concise description of distinctive features of docking programs focusing on computation methods and approximations influencing their accuracy. METHODS This review is based on the peer-reviewed research literature including author's own publications. The main features of several representative docking programs are briefly described focusing on their characteristics influencing docking accuracy: force fields, energy calculations, solvent models, algorithms of the best ligand pose search, global and local optimizations, ligand and target protein flexibility, and the simplifications made for the docking accelerating. Apart from other recent reviews focused mainly on the performance of different docking programs, in this work, an attempt is made to extract the most important functional characteristics defining the docking accuracy. Also a roadmap for increasing the docking accuracy is proposed. This is based on the new generation of docking programs which have been realized recently. These programs and respective new global optimization algorithms are described shortly. RESULTS Several popular conventional docking programs are considered. Their search of the best ligand pose is based explicitly or implicitly on the global optimization problem. Several algorithms are used to solve this problem, and among them, the heuristic genetic algorithm is distinguished by its popularity and an elaborate design. All conventional docking programs for their acceleration use the preliminary calculated grids of protein-ligand interaction potentials or preferable points of protein and ligand conjugation. These approaches and commonly used fitting parameters restrict strongly the docking accuracy. Solvent is considered in exceedingly simplified approaches in the course of the global optimization and the search for the best ligand poses. More accurate approaches on the base of implicit solvent models are used frequently for more careful binding energy calculations after docking. The new generation of docking programs are developed recently. They find the spectrum of low energy minima of a protein-ligand complex including the global minimum. These programs should be more accurate because they do not use a preliminary calculated grid of protein-ligand interaction potentials and other simplifications, the energy of any conformation of the molecular system is calculated in the frame of a given force field and there are no fitting parameters. A new docking algorithm is developed and fulfilled specially for the new docking programs. This algorithm allows docking a flexible ligand into a flexible protein with several dozen mobile atoms on the base of the global energy minimum search. Such docking results in improving the accuracy of ligand positioning in the docking process. The adequate choice of the method of molecular energy calculations also results in the better docking positioning accuracy. An advancement in the application of quantum chemistry methods to docking and scoring is revealed. CONCLUSION The findings of this review confirm the great demand in docking programs for discovery of new medicine substances with the help of molecular modeling. New trends in docking programs design are revealed. These trends are focused on the increase of the docking accuracy at the expense of more accurate molecular energy calculations without any fitting parameters, including quantum-chemical methods and implicit solvent models, and by using new global optimization algorithms which make it possible to treat flexibility of ligands and mobility of protein atoms simultaneously. Finally, it is shown that all the necessary prerequisites for increasing the docking accuracy can be accomplished in practice.
Collapse
Affiliation(s)
- Vladimir B Sulimov
- Dimonta, Ltd., Nagornaya Street 15, Building 8, 117186 Moscow, Russian Federation.,Research Computer Center, Moscow State University, Leninskie Gory 1, Building 4, 119991 Moscow, Russian Federation
| | - Danil C Kutov
- Dimonta, Ltd., Nagornaya Street 15, Building 8, 117186 Moscow, Russian Federation.,Research Computer Center, Moscow State University, Leninskie Gory 1, Building 4, 119991 Moscow, Russian Federation
| | - Alexey V Sulimov
- Dimonta, Ltd., Nagornaya Street 15, Building 8, 117186 Moscow, Russian Federation.,Research Computer Center, Moscow State University, Leninskie Gory 1, Building 4, 119991 Moscow, Russian Federation
| |
Collapse
|
6
|
Izadi S, Harris RC, Fenley MO, Onufriev AV. Accuracy Comparison of Generalized Born Models in the Calculation of Electrostatic Binding Free Energies. J Chem Theory Comput 2018; 14:1656-1670. [PMID: 29378399 DOI: 10.1021/acs.jctc.7b00886] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The need for accurate yet efficient representation of the aqueous environment in biomolecular modeling has led to the development of a variety of generalized Born (GB) implicit solvent models. While many studies have focused on the accuracy of available GB models in predicting solvation free energies, a systematic assessment of the quality of these models in binding free energy calculations, crucial for rational drug design, has not been undertaken. Here, we evaluate the accuracies of eight common GB flavors (GB-HCT, GB-OBC, GB-neck2, GBNSR6, GBSW, GBMV1, GBMV2, and GBMV3), available in major molecular dynamics packages, in predicting the electrostatic binding free energies ( ΔΔ Gel) for a diverse set of 60 biomolecular complexes belonging to four main classes: protein-protein, protein-drug, RNA-peptide, and small complexes. The GB flavors are examined in terms of their ability to reproduce the results from the Poisson-Boltzmann (PB) model, commonly used as accuracy reference in this context. We show that the agreement with the PB of ΔΔ Gel estimates varies widely between different GB models and also across different types of biomolecular complexes, with R2 correlations ranging from 0.3772 to 0.9986. A surface-based "R6" GB model recently implemented in AMBER shows the closest overall agreement with reference PB ( R2 = 0.9949, RMSD = 8.75 kcal/mol). The RNA-peptide and protein-drug complex sets appear to be most challenging for all but one model, as indicated by the large deviations from the PB in ΔΔ Gel. Small neutral complexes present the least challenge for most of the GB models tested. The quantitative demonstration of the strengths and weaknesses of the GB models across the diverse complex types provided here can be used as a guide for practical computations and future development efforts.
Collapse
Affiliation(s)
- Saeed Izadi
- Early Stage Pharmaceutical Development , Genentech Inc. , 1 DNA Way , South San Francisco , California 94080 , United States
| | - Robert C Harris
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Marcia O Fenley
- Institute of Molecular Biophysics , Florida State University , Tallahassee , Florida 32306-3408 , United States
| | | |
Collapse
|
7
|
Brieg M, Setzler J, Albert S, Wenzel W. Generalized Born implicit solvent models for small molecule hydration free energies. Phys Chem Chem Phys 2018; 19:1677-1685. [PMID: 27995260 DOI: 10.1039/c6cp07347f] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Hydration free energy estimation of small molecules from all-atom simulations was widely investigated in recent years, as it provides an essential test of molecular force fields and our understanding of solvation effects. While explicit solvent representations result in highly accurate models, they also require extensive sampling due to the high number of solvent degrees of freedom. Implicit solvent models, such as those based on the generalized Born model for electrostatic solvation effects and a solvent accessible surface area term for nonpolar contributions (GBSA), significantly reduce the number of degrees of freedom and the computational cost to estimate hydration free energies. However, a recent survey revealed a gap in the accuracy between explicit TIP3P solvent estimates and those computed with many common GBSA models. Here we address this shortcoming by providing a thorough comparison of the performance of three implicit solvent models with different nonpolar contributions and a generalized Born term to estimate experimental hydration free energies. Starting with a minimal set of only ten atom types, we demonstrate that a nonpolar term with atom type dependent surface tension coefficients in combination with an accurate generalized Born term and fully optimized parameters performs best in estimating hydration free energies, even yielding comparable results to the explicit TIP3P water model. Analysis of our results provides evidence that the asymmetric behavior of water around oppositely charged atoms is one of the main sources of error for two of the three implicit solvent models. Explicitly accounting for this effect in the parameterization reduces the corresponding errors, suggesting this as a general strategy for improving implicit solvent models. The findings presented here will help to improve the existing generalized Born based implicit solvent models implemented in state-of-the-art molecular simulation packages.
Collapse
Affiliation(s)
- Martin Brieg
- Steinbuch Centre for Computing (SCC), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany and Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany.
| | - Julia Setzler
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany.
| | - Steffen Albert
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany.
| | - Wolfgang Wenzel
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany.
| |
Collapse
|
8
|
Search for approaches to improving the calculation accuracy of the protein—ligand binding energy by docking. Russ Chem Bull 2018. [DOI: 10.1007/s11172-017-1966-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
9
|
Onufriev AV, Izadi S. Water models for biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1347] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics; Virginia Tech; Blacksburg VA USA
- Department of Computer Science; Virginia Tech; Blacksburg VA USA
- Center for Soft Matter and Biological Physics; Virginia Tech; Blacksburg VA USA
| | - Saeed Izadi
- Early Stage Pharmaceutical Development; Genentech Inc.; South San Francisco, CA USA
| |
Collapse
|
10
|
Forouzesh N, Izadi S, Onufriev AV. Grid-Based Surface Generalized Born Model for Calculation of Electrostatic Binding Free Energies. J Chem Inf Model 2017; 57:2505-2513. [DOI: 10.1021/acs.jcim.7b00192] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Saeed Izadi
- Early Stage Pharmaceutical
Development, Genentech Inc., 1 DNA
Way, South San Francisco, California 94080, United States
| | - Alexey V. Onufriev
- Center
for Soft Matter and Biological Physics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|
11
|
Sulimov AV, Zheltkov DA, Oferkin IV, Kutov DC, Katkova EV, Tyrtyshnikov EE, Sulimov VB. Evaluation of the novel algorithm of flexible ligand docking with moveable target-protein atoms. Comput Struct Biotechnol J 2017; 15:275-285. [PMID: 28377797 PMCID: PMC5367798 DOI: 10.1016/j.csbj.2017.02.004] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2016] [Accepted: 02/28/2017] [Indexed: 11/28/2022] Open
Abstract
We present the novel docking algorithm based on the Tensor Train decomposition and the TT-Cross global optimization. The algorithm is applied to the docking problem with flexible ligand and moveable protein atoms. The energy of the protein-ligand complex is calculated in the frame of the MMFF94 force field in vacuum. The grid of precalculated energy potentials of probe ligand atoms in the field of the target protein atoms is not used. The energy of the protein-ligand complex for any given configuration is computed directly with the MMFF94 force field without any fitting parameters. The conformation space of the system coordinates is formed by translations and rotations of the ligand as a whole, by the ligand torsions and also by Cartesian coordinates of the selected target protein atoms. Mobility of protein and ligand atoms is taken into account in the docking process simultaneously and equally. The algorithm is realized in the novel parallel docking SOL-P program and results of its performance for a set of 30 protein-ligand complexes are presented. Dependence of the docking positioning accuracy is investigated as a function of parameters of the docking algorithm and the number of protein moveable atoms. It is shown that mobility of the protein atoms improves docking positioning accuracy. The SOL-P program is able to perform docking of a flexible ligand into the active site of the target protein with several dozens of protein moveable atoms: the native crystallized ligand pose is correctly found as the global energy minimum in the search space with 157 dimensions using 4700 CPU ∗ h at the Lomonosov supercomputer.
Collapse
Affiliation(s)
- Alexey V Sulimov
- Dimonta, Ltd, Nagornaya Street 15, Bldg. 8, Moscow 117186, Russia; Research Computer Center, Moscow State University, Leninskie Gory 1, Bldg. 4, Moscow 119992, Russia
| | - Dmitry A Zheltkov
- Faculty of Computational Mathematics and Cybernetics of Lomonosov Moscow State University, Leninskie Gory 1, Bldg. 52, Moscow 119992, Russia
| | - Igor V Oferkin
- Dimonta, Ltd, Nagornaya Street 15, Bldg. 8, Moscow 117186, Russia
| | - Danil C Kutov
- Dimonta, Ltd, Nagornaya Street 15, Bldg. 8, Moscow 117186, Russia; Research Computer Center, Moscow State University, Leninskie Gory 1, Bldg. 4, Moscow 119992, Russia
| | - Ekaterina V Katkova
- Dimonta, Ltd, Nagornaya Street 15, Bldg. 8, Moscow 117186, Russia; Research Computer Center, Moscow State University, Leninskie Gory 1, Bldg. 4, Moscow 119992, Russia
| | - Eugene E Tyrtyshnikov
- Faculty of Computational Mathematics and Cybernetics of Lomonosov Moscow State University, Leninskie Gory 1, Bldg. 52, Moscow 119992, Russia; Institute of Numerical Mathematics of Russian Academy of Sciences, Gubkin Street 8, Moscow, 119333, Russia
| | - Vladimir B Sulimov
- Dimonta, Ltd, Nagornaya Street 15, Bldg. 8, Moscow 117186, Russia; Research Computer Center, Moscow State University, Leninskie Gory 1, Bldg. 4, Moscow 119992, Russia
| |
Collapse
|
12
|
Katkova EV, Onufriev AV, Aguilar B, Sulimov VB. Accuracy comparison of several common implicit solvent models and their implementations in the context of protein-ligand binding. J Mol Graph Model 2016; 72:70-80. [PMID: 28064081 DOI: 10.1016/j.jmgm.2016.12.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Revised: 12/07/2016] [Accepted: 12/15/2016] [Indexed: 11/18/2022]
Abstract
In this study several commonly used implicit solvent models are compared with respect to their accuracy of estimating solvation energies of small molecules and proteins, as well as desolvation penalty in protein-ligand binding. The test set consists of 19 small proteins, 104 small molecules, and 15 protein-ligand complexes. We compared predicted hydration energies of small molecules with their experimental values; the results of the solvation and desolvation energy calculations for small molecules, proteins and protein-ligand complexes in water were also compared with Thermodynamic Integration calculations based on TIP3P water model and Amber12 force field. The following implicit solvent (water) models considered here are: PCM (Polarized Continuum Model implemented in DISOLV and MCBHSOLV programs), GB (Generalized Born method implemented in DISOLV program, S-GB, and GBNSR6 stand-alone version), COSMO (COnductor-like Screening Model implemented in the DISOLV program and the MOPAC package) and the Poisson-Boltzmann model (implemented in the APBS program). Different parameterizations of the molecules were examined: we compared MMFF94 force field, Amber12 force field and the quantum-chemical semi-empirical PM7 method implemented in the MOPAC package. For small molecules, all of the implicit solvent models tested here yield high correlation coefficients (0.87-0.93) between the calculated solvation energies and the experimental values of hydration energies. For small molecules high correlation (0.82-0.97) with the explicit solvent energies is seen as well. On the other hand, estimated protein solvation energies and protein-ligand binding desolvation energies show substantial discrepancy (up to 10kcal/mol) with the explicit solvent reference. The correlation of polar protein solvation energies and protein-ligand desolvation energies with the corresponding explicit solvent results is 0.65-0.99 and 0.76-0.96 respectively, though this difference in correlations is caused more by different parameterization and less by methods and indicates the need for further improvement of implicit solvent models parameterization. Within the same parameterization, various implicit methods give practically the same correlation with results obtained in explicit solvent model for ligands and proteins: e.g. correlation values of polar ligand solvation energies and the corresponding energies in the frame of explicit solvent were 0.953-0.966 for the APBS program, the GBNSR6 program and all models used in the DISOLV program. The DISOLV program proved to be on a par with the other used programs in the case of proteins and ligands solvation energy calculation. However, the solution of the Poisson-Boltzmann equation (APBS program) and Generalized Born method (implemented in the GBNSR6 program) proved to be the most accurate in calculating the desolvation energies of complexes.
Collapse
Affiliation(s)
- E V Katkova
- Dimonta, Ltd., Nagornaya Street 15, Bldg 8, Moscow, 117186, Russia; Research Computer Center, Lomonosov Moscow State University, Leninskie Gory 1,Bldg 4, Moscow, 119992, Russia.
| | - A V Onufriev
- Departments of Computer Science and Physics, Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, VA, USA
| | - B Aguilar
- Institute for Systems Biology, Seattle, WA, USA
| | - V B Sulimov
- Dimonta, Ltd., Nagornaya Street 15, Bldg 8, Moscow, 117186, Russia; Research Computer Center, Lomonosov Moscow State University, Leninskie Gory 1,Bldg 4, Moscow, 119992, Russia
| |
Collapse
|
13
|
Grigoriev FV, Sulimov VB. Simple combined explicit/implicit water model. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2016.1218488] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
14
|
Evaluation of Docking Target Functions by the Comprehensive Investigation of Protein-Ligand Energy Minima. Adv Bioinformatics 2015; 2015:126858. [PMID: 26693223 PMCID: PMC4674582 DOI: 10.1155/2015/126858] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2015] [Revised: 10/28/2015] [Accepted: 11/04/2015] [Indexed: 12/19/2022] Open
Abstract
The adequate choice of the docking target function impacts the accuracy of the ligand positioning as well as the accuracy of the protein-ligand binding energy calculation. To evaluate a docking target function we compared positions of its minima with the experimentally known pose of the ligand in the protein active site. We evaluated five docking target functions based on either the MMFF94 force field or the PM7 quantum-chemical method with or without implicit solvent models: PCM, COSMO, and SGB. Each function was tested on the same set of 16 protein-ligand complexes. For exhaustive low-energy minima search the novel MPI parallelized docking program FLM and large supercomputer resources were used. Protein-ligand binding energies calculated using low-energy minima were compared with experimental values. It was demonstrated that the docking target function on the base of the MMFF94 force field in vacuo can be used for discovery of native or near native ligand positions by finding the low-energy local minima spectrum of the target function. The importance of solute-solvent interaction for the correct ligand positioning is demonstrated. It is shown that docking accuracy can be improved by replacement of the MMFF94 force field by the new semiempirical quantum-chemical PM7 method.
Collapse
|
15
|
Application of Molecular Modeling to Development of New Factor Xa Inhibitors. BIOMED RESEARCH INTERNATIONAL 2015; 2015:120802. [PMID: 26484350 PMCID: PMC4592935 DOI: 10.1155/2015/120802] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2015] [Revised: 08/12/2015] [Accepted: 08/20/2015] [Indexed: 12/30/2022]
Abstract
In consequence of the key role of factor Xa in the clotting cascade and absence of its activity in the processes that do not affect coagulation, this protein is an attractive target for development of new blood coagulation inhibitors. Factor Xa is more effective and convenient target for creation of anticoagulants than thrombin, inhibition of which may cause some side effects. This study is aimed at finding new inhibitors of factor Xa by molecular computer modeling including docking SOL and postdocking optimization DISCORE programs. After validation of molecular modeling methods on well-known factor Xa inhibitors the virtual screening of NCI Diversity and Voronezh State University databases of ready-made low molecular weight species has been carried out. Seventeen compounds selected on the basis of modeling results have been tested experimentally in vitro. It has been found that 12 of them showed activity against factor Xa (IC50 = 1.8-40 μM). Based on analysis of the results, the new original compound was synthesized and experimentally verified. It shows activity against factor Xa with IC50 value of 0.7 μM.
Collapse
|
16
|
Izadi S, Aguilar B, Onufriev AV. Protein-Ligand Electrostatic Binding Free Energies from Explicit and Implicit Solvation. J Chem Theory Comput 2015; 11:4450-9. [PMID: 26575935 PMCID: PMC5217485 DOI: 10.1021/acs.jctc.5b00483] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Accurate yet efficient computational models of solvent environment are central for most calculations that rely on atomistic modeling, such as prediction of protein-ligand binding affinities. In this study, we evaluate the accuracy of a recently developed generalized Born implicit solvent model, GBNSR6 (Aguilar et al. J. Chem. Theory Comput. 2010, 6, 3613-3639), in estimating the electrostatic solvation free energies (ΔG(pol)) and binding free energies (ΔΔG(pol)) for small protein-ligand complexes. We also compare estimates based on three different explicit solvent models (TIP3P, TIP4PEw, and OPC). The two main findings are as follows. First, the deviation (RMSD = 7.04 kcal/mol) of GBNSR6 binding affinities from commonly used TIP3P reference values is comparable to the deviations between explicit models themselves, e.g. TIP4PEw vs TIP3P (RMSD = 5.30 kcal/mol). A simple uniform adjustment of the atomic radii by a single scaling factor reduces the RMS deviation of GBNSR6 from TIP3P to within the above "error margin" - differences between ΔΔG(pol) estimated by different common explicit solvent models. The simple radii scaling virtually eliminates the systematic deviation (ΔΔG(pol)) between GBNSR6 and two out of the three explicit water models and significantly reduces the deviation from the third explicit model. Second, the differences between electrostatic binding energy estimates from different explicit models is disturbingly large; for example, the deviation between TIP4PEw and TIP3P estimates of ΔΔG(pol) values can be up to ∼50% or ∼9 kcal/mol, which is significantly larger than the "chemical accuracy" goal of ∼1 kcal/mol. The absolute ΔG(pol) calculated with different explicit models could differ by tens of kcal/mol. These discrepancies point to unacceptably high sensitivity of binding affinity estimates to the choice of common explicit water models. The absence of a clear "gold standard" among these models strengthens the case for the use of accurate implicit solvation models for binding energetics, which may be orders of magnitude faster.
Collapse
Affiliation(s)
- Saeed Izadi
- Department of Biomedical Engineering and Mechanics, Department of Computer Science, and Departments of Computer Science and Physics, Virginia Tech , Blacksburg, Virginia 24060, United States
| | - Boris Aguilar
- Department of Biomedical Engineering and Mechanics, Department of Computer Science, and Departments of Computer Science and Physics, Virginia Tech , Blacksburg, Virginia 24060, United States
| | - Alexey V Onufriev
- Department of Biomedical Engineering and Mechanics, Department of Computer Science, and Departments of Computer Science and Physics, Virginia Tech , Blacksburg, Virginia 24060, United States
| |
Collapse
|
17
|
Setzler J, Seith C, Brieg M, Wenzel W. SLIM: an improved generalized Born implicit membrane model. J Comput Chem 2015; 35:2027-39. [PMID: 25243932 DOI: 10.1002/jcc.23717] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2013] [Revised: 07/02/2014] [Accepted: 07/28/2014] [Indexed: 12/24/2022]
Abstract
In most implicit continuum models, membranes are represented as heterogeneous dielectric environments, but their treatment within computationally efficient generalized Born (GB) models is challenging. Despite several previous attempts, an adequate description of multiple dielectric regions in implicit GB-based membrane models that reproduce the qualitative and quantitative features of Poisson-Boltzmann (PB) electrostatics remains an unmet prerequisite of qualitatively correct implicit membrane models. A novel scheme (SLIM) to decompose one environment consisting of multiple dielectric regions into a sum of multiple environments consisting only of two dielectric regions each is proposed to solve this issue. These simpler environments can be treated with established GB methods. This approach captures qualitative features of PB electrostatic that are not present in previous models. Simulations of three membrane proteins demonstrate that this model correctly reproduces known properties of these proteins in agreement with experimental or other computational studies.
Collapse
Affiliation(s)
- Julia Setzler
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021, Karlsruhe, Germany
| | | | | | | |
Collapse
|
18
|
Decherchi S, Masetti M, Vyalov I, Rocchia W. Implicit solvent methods for free energy estimation. Eur J Med Chem 2014; 91:27-42. [PMID: 25193298 DOI: 10.1016/j.ejmech.2014.08.064] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2014] [Revised: 08/21/2014] [Accepted: 08/23/2014] [Indexed: 12/12/2022]
Abstract
Solvation is a fundamental contribution in many biological processes and especially in molecular binding. Its estimation can be performed by means of several computational approaches. The aim of this review is to give an overview of existing theories and methods to estimate solvent effects giving a specific focus on the category of implicit solvent models and their use in Molecular Dynamics. In many of these models, the solvent is considered as a continuum homogenous medium, while the solute can be represented at the atomic detail and at different levels of theory. Despite their degree of approximation, implicit methods are still widely employed due to their trade-off between accuracy and efficiency. Their derivation is rooted in the statistical mechanics and integral equations disciplines, some of the related details being provided here. Finally, methods that combine implicit solvent models and molecular dynamics simulation, are briefly described.
Collapse
Affiliation(s)
- Sergio Decherchi
- CONCEPT Lab, D3 Computation, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Via Belmeloro 6, 40126 Bologna, Italy
| | - Ivan Vyalov
- CONCEPT Lab, D3 Computation, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Walter Rocchia
- CONCEPT Lab, D3 Computation, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy.
| |
Collapse
|
19
|
Onufriev AV, Aguilar B. Accuracy of continuum electrostatic calculations based on three common dielectric boundary definitions. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2014; 13. [PMID: 26236064 DOI: 10.1142/s0219633614400069] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We investigate the influence of three common definitions of the solute/solvent dielectric boundary (DB) on the accuracy of the electrostatic solvation energy ΔGel computed within the Poisson Boltzmann and the generalized Born models of implicit solvation. The test structures include small molecules, peptides and small proteins; explicit solvent ΔGel are used as accuracy reference. For common atomic radii sets BONDI, PARSE (and ZAP9 for small molecules) the use of van der Waals (vdW) DB results, on average, in considerably larger errors in ΔGel than the molecular surface (MS) DB. The optimal probe radius ρw for which the MS DB yields the most accurate ΔGel varies considerably between structure types. The solvent accessible surface (SAS) DB becomes optimal at ρw ~ 0.2 Å (exact value is sensitive to the structure and atomic radii), at which point the average accuracy of ΔGel is comparable to that of the MS-based boundary. The geometric equivalence of SAS to vdW surface based on the same atomic radii uniformly increased by ρw gives the corresponding optimal vdW DB. For small molecules, the optimal vdW DB based on BONDI + 0.2 Å radii can yield ΔGel estimates at least as accurate as those based on the optimal MS DB. Also, in small molecules, pairwise charge-charge interactions computed with the optimal vdW DB are virtually equal to those computed with the MS DB, suggesting that in this case the two boundaries are practically equivalent by the electrostatic energy criteria. In structures other than small molecules, the optimal vdW and MS dielectric boundaries are not equivalent: the respective pairwise electrostatic interactions in the presence of solvent can differ by up to 5 kcal/mol for individual atomic pairs in small proteins, even when the total ΔGel are equal. For small proteins, the average decrease in pairwise electrostatic interactions resulting from the switch from optimal MS to optimal vdW DB definition can be mimicked within the MS DB definition by doubling of the solute dielectric constant. However, the use of the higher interior dielectric does not eliminate the large individual deviations between pairwise interactions computed within the two DB definitions. It is argued that while the MS based definition of the dielectric boundary is more physically correct in some types of practical calculations, the choice is not so clear in some other common scenarios.
Collapse
Affiliation(s)
- Alexey V Onufriev
- Department of Computer Science and Department of Physics, Virginia Tech, Blacksburg, VA 24060, and Department of Computer Science, Virginia Tech, Blacksburg, VA 24060
| | - Boris Aguilar
- Department of Computer Science and Department of Physics, Virginia Tech, Blacksburg, VA 24060, and Department of Computer Science, Virginia Tech, Blacksburg, VA 24060
| |
Collapse
|
20
|
Application of molecular modeling to urokinase inhibitors development. BIOMED RESEARCH INTERNATIONAL 2014; 2014:625176. [PMID: 24967388 PMCID: PMC4055159 DOI: 10.1155/2014/625176] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/12/2013] [Accepted: 04/22/2014] [Indexed: 01/01/2023]
Abstract
Urokinase-type plasminogen activator (uPA) plays an important role in the regulation of diverse physiologic and pathologic processes. Experimental research has shown that elevated uPA expression is associated with cancer progression, metastasis, and shortened survival in patients, whereas suppression of proteolytic activity of uPA leads to evident decrease of metastasis. Therefore, uPA has been considered as a promising molecular target for development of anticancer drugs. The present study sets out to develop the new selective uPA inhibitors using computer-aided structural based drug design methods. Investigation involves the following stages: computer modeling of the protein active site, development and validation of computer molecular modeling methods: docking (SOL program), postprocessing (DISCORE program), direct generalized docking (FLM program), and the application of the quantum chemical calculations (MOPAC package), search of uPA inhibitors among molecules from databases of ready-made compounds to find new uPA inhibitors, and design of new chemical structures and their optimization and experimental examination. On the basis of known uPA inhibitors and modeling results, 18 new compounds have been designed, calculated using programs mentioned above, synthesized, and tested in vitro. Eight of them display inhibitory activity and two of them display activity about 10 μM.
Collapse
|
21
|
TANG XIAOCHUAN, DUAN YONG. VERIFICATION OF THE GENERALIZED BORN MODEL AT SHORT DISTANCES. J MECH MED BIOL 2014. [DOI: 10.1142/s0219519413400204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The generalized Born (GB) model, one of the implicit solvent models, is widely applied in molecular dynamics (MD) simulations as a simple description of the solvation effect. In the GB model, an empirical function called the Still's formula, with the algorithmic simplicity, is utilized to calculate the solvation energy due to the polarization, termed as ΔG pol . Applications of the GB model have exhibited reasonable accuracy and high computational efficiency. However, there is still room for improvements. Most of the attempts to improve the GB model focus on optimizing effective Born radii. Contrarily, limited researches have been performed to improve the feasibility of the Still's formula. In this paper, analytical methods was applied to investigate the validity of the Still's formula at short distance. Taking advantage of the toroidal coordinates and Mehler–Fock transform, the analytical solutions of the GB model at short distances was derived explicitly for the first time. Additionally, the solvation energy was numerically computed using proper algorithms based on the analytical solutions and compared with ΔG pol calculated in the GB model. With the analysis on the deficiencies of the Still's formula at short distances, potential methods to improve the validity of the GB model were discussed.
Collapse
Affiliation(s)
- XIAOCHUAN TANG
- 4335 Genome and Biomedical Sciences Facility, 451 Health Sciences Drive, Davis, California 95616, USA
- Applied Science Graduate Program and Genome Center, University of California Davis, One Shields Avenue, Davis, California 95616, USA
| | - YONG DUAN
- Applied Science Graduate Program and Genome Center, University of California Davis, One Shields Avenue, Davis, California 95616, USA
| |
Collapse
|
22
|
Sulimov AV, Kutov DC, Oferkin IV, Katkova EV, Sulimov VB. Application of the Docking Program SOL for CSAR Benchmark. J Chem Inf Model 2013; 53:1946-56. [DOI: 10.1021/ci400094h] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alexey V. Sulimov
- Research Computer Center, Moscow State University, Leninskie Gory 1, bldg 4,
Moscow 119992, Russia
- Dimonta, Ltd., Nagornaya Street 15, bldg 8, Moscow
117186, Russia
| | - Danil C. Kutov
- Research Computer Center, Moscow State University, Leninskie Gory 1, bldg 4,
Moscow 119992, Russia
| | - Igor V. Oferkin
- Research Computer Center, Moscow State University, Leninskie Gory 1, bldg 4,
Moscow 119992, Russia
- Dimonta, Ltd., Nagornaya Street 15, bldg 8, Moscow
117186, Russia
| | - Ekaterina V. Katkova
- Research Computer Center, Moscow State University, Leninskie Gory 1, bldg 4,
Moscow 119992, Russia
- Dimonta, Ltd., Nagornaya Street 15, bldg 8, Moscow
117186, Russia
| | - Vladimir B. Sulimov
- Research Computer Center, Moscow State University, Leninskie Gory 1, bldg 4,
Moscow 119992, Russia
- Dimonta, Ltd., Nagornaya Street 15, bldg 8, Moscow
117186, Russia
| |
Collapse
|
23
|
Parameterization of the hydration free energy computations for organic solutes in the framework of the implicit solvent model with the nonuniform dielectric function. COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2012.12.027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
24
|
Brieg M, Wenzel W. PowerBorn: A Barnes-Hut Tree Implementation for Accurate and Efficient Born Radii Computation. J Chem Theory Comput 2013; 9:1489-98. [PMID: 26587611 DOI: 10.1021/ct300870s] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Implicit solvent models are one of the standard tools in computational biophysics. While Poisson-Boltzmann methods offer highly accurate results within this framework, generalized Born models have been used due to their higher computational efficiency in many (bio)molecular simulations, where computational power is a limiting factor. In recent years, there have been remarkable advances to reduce some deficiencies in the generalized Born models. On the other hand, these advances come at an increased computational cost that contrasts the reasons for choosing generalized Born models over Poisson-Boltzmann methods. To address this performance issue, we present a new algorithm for Born radii computation, one performance critical part in the evaluation of generalized Born models, which is based on a Barnes-Hut tree code scheme. We show that an implementation of this algorithm provides accurate Born radii and polar solvation free energies in comparison to Poisson-Boltzmann computations, while delivering up to an order of magnitude better performance over existing, similarly accurate methods. The C++ implementation of this algorithm will be available at http://www.int.kit.edu/nanosim/ .
Collapse
Affiliation(s)
- Martin Brieg
- Steinbuch Centre for Computing (SCC), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), P.O. Box 3640, 76021 Karlsruhe, Germany
| |
Collapse
|
25
|
|
26
|
Lange AW, Herbert JM. Improving Generalized Born Models by Exploiting Connections to Polarizable Continuum Models. I. An Improved Effective Coulomb Operator. J Chem Theory Comput 2012; 8:1999-2011. [DOI: 10.1021/ct300111m] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Adrian W. Lange
- Department
of Chemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department
of Chemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
27
|
Fogolari F, Corazza A, Yarra V, Jalaru A, Viglino P, Esposito G. Bluues: a program for the analysis of the electrostatic properties of proteins based on generalized Born radii. BMC Bioinformatics 2012; 13 Suppl 4:S18. [PMID: 22536964 PMCID: PMC3434445 DOI: 10.1186/1471-2105-13-s4-s18] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
BACKGROUND The Poisson-Boltzmann (PB) equation and its linear approximation have been widely used to describe biomolecular electrostatics. Generalized Born (GB) models offer a convenient computational approximation for the more fundamental approach based on the Poisson-Boltzmann equation, and allows estimation of pairwise contributions to electrostatic effects in the molecular context. RESULTS We have implemented in a single program most common analyses of the electrostatic properties of proteins. The program first computes generalized Born radii, via a surface integral and then it uses generalized Born radii (using a finite radius test particle) to perform electrostaic analyses. In particular the ouput of the program entails, depending on user's requirement: 1) the generalized Born radius of each atom; 2) the electrostatic solvation free energy; 3) the electrostatic forces on each atom (currently in a developmental stage); 4) the pH-dependent properties (total charge and pH-dependent free energy of folding in the pH range -2 to 18; 5) the pKa of all ionizable groups; 6) the electrostatic potential at the surface of the molecule; 7) the electrostatic potential in a volume surrounding the molecule; CONCLUSIONS Although at the expense of limited flexibility the program provides most common analyses with requirement of a single input file in PQR format. The results obtained are comparable to those obtained using state-of-the-art Poisson-Boltzmann solvers. A Linux executable with example input and output files is provided as supplementary material.
Collapse
Affiliation(s)
- Federico Fogolari
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
- Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, Roma 00136, Italy
| | - Alessandra Corazza
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
- Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, Roma 00136, Italy
| | - Vijaylakshmi Yarra
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
| | - Anusha Jalaru
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
| | - Paolo Viglino
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
- Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, Roma 00136, Italy
| | - Gennaro Esposito
- Dipartimento di Scienze Mediche e Biologiche. Università di Udine, Piazzale Kolbe, 4, Udine 33100, Italy
- Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, Roma 00136, Italy
| |
Collapse
|
28
|
Xu Z, Cai W. Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations. SIAM REVIEW. SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS 2011; 53:683-720. [PMID: 23745011 PMCID: PMC3671632 DOI: 10.1137/090774288] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
We review recent developments of fast analytical methods for macroscopic electrostatic calculations in biological applications, including the Poisson-Boltzmann (PB) and the generalized Born models for electrostatic solvation energy. The focus is on analytical approaches for hybrid solvation models, especially the image charge method for a spherical cavity, and also the generalized Born theory as an approximation to the PB model. This review places much emphasis on the mathematical details behind these methods.
Collapse
Affiliation(s)
- Zhenli Xu
- Department of Mathematics and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, People's Republic of China, and Department of Mathematics and Statistics, University of North Carolina, Charlotte, NC 28223 ()
| | - Wei Cai
- Department of Mathematics and Statistics, University of North Carolina, Charlotte, NC 28223 (), and Beijing International Center for Mathematical Research, Beijing, People's Republic of China, 100871
| |
Collapse
|
29
|
Kupervasser OY, Kikot’ IP. Enlarged surface meshes and normalization conditions for columns and rows of matrices in the COSMO method. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2011. [DOI: 10.1134/s199079311109020x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
30
|
Bardhan JP, Knepley MG. Mathematical analysis of the boundary-integral based electrostatics estimation approximation for molecular solvation: exact results for spherical inclusions. J Chem Phys 2011; 135:124107. [PMID: 21974512 DOI: 10.1063/1.3641485] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We analyze the mathematically rigorous BIBEE (boundary-integral based electrostatics estimation) approximation of the mixed-dielectric continuum model of molecular electrostatics, using the analytically solvable case of a spherical solute containing an arbitrary charge distribution. Our analysis, which builds on Kirkwood's solution using spherical harmonics, clarifies important aspects of the approximation and its relationship to generalized Born models. First, our results suggest a new perspective for analyzing fast electrostatic models: the separation of variables between material properties (the dielectric constants) and geometry (the solute dielectric boundary and charge distribution). Second, we find that the eigenfunctions of the reaction-potential operator are exactly preserved in the BIBEE model for the sphere, which supports the use of this approximation for analyzing charge-charge interactions in molecular binding. Third, a comparison of BIBEE to the recent GBε theory suggests a modified BIBEE model capable of predicting electrostatic solvation free energies to within 4% of a full numerical Poisson calculation. This modified model leads to a projection-framework understanding of BIBEE and suggests opportunities for future improvements.
Collapse
Affiliation(s)
- Jaydeep P Bardhan
- Department of Molecular Biophysics and Physiology, Rush University Medical Center, Chicago, Illinois 60612, USA.
| | | |
Collapse
|
31
|
Abstract
We have developed a treecode-based O(N log N) algorithm for the generalized Born (GB) implicit solvation model. Our treecode-based GB (tGB) is based on the GBr6 [J. Phys. Chem. B 111, 3055 (2007)], an analytical GB method with a pairwise descreening approximation for the R6 volume integral expression. The algorithm is composed of a cutoff scheme for the effective Born radii calculation, and a treecode implementation of the GB charge-charge pair interactions. Test results demonstrate that the tGB algorithm can reproduce the vdW surface based Poisson solvation energy with an average relative error less than 0.6% while providing an almost linear-scaling calculation for a representative set of 25 proteins with different sizes (from 2815 atoms to 65456 atoms). For a typical system of 10k atoms, the tGB calculation is three times faster than the direct summation as implemented in the original GBr6 model. Thus, our tGB method provides an efficient way for performing implicit solvent GB simulations of larger biomolecular systems at longer time scales.
Collapse
Affiliation(s)
- Zhenli Xu
- Department of Mathematics, and Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China.
| | | | | |
Collapse
|
32
|
Sinauridze EI, Romanov AN, Gribkova IV, Kondakova OA, Surov SS, Gorbatenko AS, Butylin AA, Monakov MY, Bogolyubov AA, Kuznetsov YV, Sulimov VB, Ataullakhanov FI. New synthetic thrombin inhibitors: molecular design and experimental verification. PLoS One 2011; 6:e19969. [PMID: 21603576 PMCID: PMC3095642 DOI: 10.1371/journal.pone.0019969] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2010] [Accepted: 04/19/2011] [Indexed: 12/04/2022] Open
Abstract
Background The development of new anticoagulants is an important goal for the improvement of thromboses treatments. Objectives The design, synthesis and experimental testing of new safe and effective small molecule direct thrombin inhibitors for intravenous administration. Methods Computer-aided molecular design of new thrombin inhibitors was performed using our original docking program SOL, which is based on the genetic algorithm of global energy minimization in the framework of a Merck Molecular Force Field. This program takes into account the effects of solvent. The designed molecules with the best scoring functions (calculated binding energies) were synthesized and their thrombin inhibitory activity evaluated experimentally in vitro using a chromogenic substrate in a buffer system and using a thrombin generation test in isolated plasma and in vivo using the newly developed model of hemodilution-induced hypercoagulation in rats. The acute toxicities of the most promising new thrombin inhibitors were evaluated in mice, and their stabilities in aqueous solutions were measured. Results New compounds that are both effective direct thrombin inhibitors (the best KI was <1 nM) and strong anticoagulants in plasma (an IC50 in the thrombin generation assay of approximately 100 nM) were discovered. These compounds contain one of the following new residues as the basic fragment: isothiuronium, 4-aminopyridinium, or 2-aminothiazolinium. LD50 values for the best new inhibitors ranged from 166.7 to >1111.1 mg/kg. A plasma-substituting solution supplemented with one of the new inhibitors prevented hypercoagulation in the rat model of hemodilution-induced hypercoagulation. Activities of the best new inhibitors in physiological saline (1 µM solutions) were stable after sterilization by autoclaving, and the inhibitors remained stable at long-term storage over more than 1.5 years at room temperature and at 4°C. Conclusions The high efficacy, stability and low acute toxicity reveal that the inhibitors that were developed may be promising for potential medical applications.
Collapse
Affiliation(s)
- Elena I Sinauridze
- Laboratory of Physical Biochemistry, National Research Center for Hematology, Russian Academy of Medical Sciences, Moscow, Russia.
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
33
|
Fedichev PO, Getmantsev EG, Menshikov LI. O(NlogN) Continuous electrostatics method for fast calculation of solvation energies of biomolecules. J Comput Chem 2010; 32:1368-76. [DOI: 10.1002/jcc.21719] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2010] [Revised: 10/20/2010] [Accepted: 10/21/2010] [Indexed: 11/07/2022]
|
34
|
Aguilar B, Shadrach R, Onufriev AV. Reducing the Secondary Structure Bias in the Generalized Born Model via R6 Effective Radii. J Chem Theory Comput 2010. [DOI: 10.1021/ct100392h] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Boris Aguilar
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24060, United States, Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States, and Departments of Computer Science and Physics, Virginia Tech, Blacksburg, Virginia 24060, United States
| | - Richard Shadrach
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24060, United States, Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States, and Departments of Computer Science and Physics, Virginia Tech, Blacksburg, Virginia 24060, United States
| | - Alexey V. Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg, Virginia 24060, United States, Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States, and Departments of Computer Science and Physics, Virginia Tech, Blacksburg, Virginia 24060, United States
| |
Collapse
|
35
|
Bardhan JP. Interpreting the Coulomb-field approximation for generalized-Born electrostatics using boundary-integral equation theory. J Chem Phys 2009; 129:144105. [PMID: 19045132 DOI: 10.1063/1.2987409] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The importance of molecular electrostatic interactions in aqueous solution has motivated extensive research into physical models and numerical methods for their estimation. The computational costs associated with simulations that include many explicit water molecules have driven the development of implicit-solvent models, with generalized-Born (GB) models among the most popular of these. In this paper, we analyze a boundary-integral equation interpretation for the Coulomb-field approximation (CFA), which plays a central role in most GB models. This interpretation offers new insights into the nature of the CFA, which traditionally has been assessed using only a single point charge in the solute. The boundary-integral interpretation of the CFA allows the use of multiple point charges, or even continuous charge distributions, leading naturally to methods that eliminate the interpolation inaccuracies associated with the Still equation. This approach, which we call boundary-integral-based electrostatic estimation by the CFA (BIBEE/CFA), is most accurate when the molecular charge distribution generates a smooth normal displacement field at the solute-solvent boundary, and CFA-based GB methods perform similarly. Conversely, both methods are least accurate for charge distributions that give rise to rapidly varying or highly localized normal displacement fields. Supporting this analysis are comparisons of the reaction-potential matrices calculated using GB methods and boundary-element-method (BEM) simulations. An approximation similar to BIBEE/CFA exhibits complementary behavior, with superior accuracy for charge distributions that generate rapidly varying normal fields and poorer accuracy for distributions that produce smooth fields. This approximation, BIBEE by preconditioning (BIBEE/P), essentially generates initial guesses for preconditioned Krylov-subspace iterative BEMs. Thus, iterative refinement of the BIBEE/P results recovers the BEM solution; excellent agreement is obtained in only a few iterations. The boundary-integral-equation framework may also provide a means to derive rigorous results explaining how the empirical correction terms in many modern GB models significantly improve accuracy despite their simple analytical forms.
Collapse
Affiliation(s)
- Jaydeep P Bardhan
- Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA.
| |
Collapse
|
36
|
Grigoriev FV, Gabin SN, Romanov AN, Sulimov VB. Computation of the contribution from the cavity effect to protein-ligand binding free energy. J Phys Chem B 2008; 112:15355-60. [PMID: 18991438 DOI: 10.1021/jp8041439] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present results of the investigation of the cavity creation/annihilation effect in view of formation of the protein-ligand (PL) complexes. The protein and ligand were considered as rigid structures. The change of the cavity creation/annihilation free energy DeltaG(cav) was calculated for three PL complexes using the thermodynamic integration procedure with the original algorithm for growing the interaction potential between the cavity and the water molecules. The thermodynamic cycle consists of two stages, annihilation of the cavity of the ligand for the unbound state and its creation at the active site of the protein (bound state). It was revealed that for all complexes under investigation, the values of DeltaG(cav) are negative and favorable for binding. The main contribution to DeltaG(cav) appears due to the annihilation of the cavity of the ligand. All computations were made using the parallel version of CAVE code, elaborated in our preceding work.
Collapse
Affiliation(s)
- F V Grigoriev
- Research Computing Center, M.V. Lomonosov Moscow State University, Russia, 119992 Moscow, Vorobjovi Gory, MGU NIVC.
| | | | | | | |
Collapse
|
37
|
Mongan J, Svrcek-Seiler WA, Onufriev A. Analysis of integral expressions for effective Born radii. J Chem Phys 2008; 127:185101. [PMID: 18020664 DOI: 10.1063/1.2783847] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Generalized Born (GB) models provide a computationally efficient means of representing the electrostatic effects of solvent and are widely used, especially in molecular dynamics (MD). Accurate and facile computation of the effective Born radii is a key for the performance of GB models. Here, we examine a simple integral prescription, R6, based on the exact solution of the Poisson-Boltzmann (PB) equation for a perfect sphere. Numerical tests on 22 molecules representing a variety of structural classes show that R6 may be more accurate than the more complex integral-based approaches such as GBMV2. At the same time, R6 is computationally less demanding. Fundamental limitations of current integration-based methods for calculating effective radii, including R6, are explored and the deviations from the numerical PB results are correlated with specific topological and geometrical features of the molecular surface. A small systematic bias observed in the R6-based radii can be removed with a single, transferable constant offset; when the resulting effective radii are used in the "classical" (Still et al.'s) GB formula to compute the electrostatic solvation free energy, the average deviation from the PB reference is no greater than when the "perfect" (PB-based) effective radii are used. This deviation is also appreciably smaller than the uncertainty of the PB reference itself, as estimated by comparison to explicit solvent.
Collapse
Affiliation(s)
- John Mongan
- Bioinformatics Program, Medical Scientist Training Program, Center for Theoretical Biological Physics, and Department of Chemistry and Biochemistry, UC San Diego, La Jolla, California 92093-0365, USA
| | | | | |
Collapse
|
38
|
Mongan J, Simmerling C, McCammon JA, Case DA, Onufriev A. Generalized Born model with a simple, robust molecular volume correction. J Chem Theory Comput 2006; 3:156-169. [PMID: 21072141 DOI: 10.1021/ct600085e] [Citation(s) in RCA: 295] [Impact Index Per Article: 16.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Generalized Born (GB) models provide a computationally efficient means of representing the electrostatic effects of solvent and are widely used, especially in molecular dynamics (MD). A class of particularly fast GB models is based on integration over an interior volume approximated as a pairwise union of atom spheres-effectively, the interior is defined by a van der Waals rather than Lee-Richards molecular surface. The approximation is computationally efficient, but if uncorrected, allows for high dielectric (water) regions smaller than a water molecule between atoms, leading to decreased accuracy. Here, an earlier pairwise GB model is extended by a simple analytic correction term that largely alleviates the problem by correctly describing the solvent-excluded volume of each pair of atoms. The correction term introduces a free energy barrier to the separation of non-bonded atoms. This free energy barrier is seen in explicit solvent and Lee-Richards molecular surface implicit solvent calculations, but has been absent from earlier pairwise GB models. When used in MD, the correction term yields protein hydrogen bond length distributions and polypeptide conformational ensembles that are in better agreement with explicit solvent results than earlier pairwise models. The robustness and simplicity of the correction preserves the efficiency of the pairwise GB models while making them a better approximation to reality.
Collapse
Affiliation(s)
- John Mongan
- Bioinformatics Program, Medical Scientist Training Program, Center for Theoretical Biological Physics, UC San Diego, La Jolla, CA 92093-0365
| | | | | | | | | |
Collapse
|
39
|
Zhu K, Pincus DL, Zhao S, Friesner RA. Long loop prediction using the protein local optimization program. Proteins 2006; 65:438-52. [PMID: 16927380 DOI: 10.1002/prot.21040] [Citation(s) in RCA: 104] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We have developed an improved sampling algorithm and energy model for protein loop prediction, the combination of which has yielded the first methodology capable of achieving good results for the prediction of loop backbone conformations of 11 residue length or greater. Applied to our newly constructed test suite of 104 loops ranging from 11 to 13 residues, our method obtains average/median global backbone root-mean-square deviations (RMSDs) to the native structure (superimposing the body of the protein, not the loop itself) of 1.00/0.62 A for 11 residue loops, 1.15/0.60 A for 12 residue loops, and 1.25/0.76 A for 13 residue loops. Sampling errors are virtually eliminated, while energy errors leading to large backbone RMSDs are very infrequent compared to any previously reported efforts, including our own previous study. We attribute this success to both an improved sampling algorithm and, more critically, the inclusion of a hydrophobic term, which appears to approximately fix a major flaw in SGB solvation model that we have been employing. A discussion of these results in the context of the general question of the accuracy of continuum solvation models is presented.
Collapse
Affiliation(s)
- Kai Zhu
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | | | | |
Collapse
|
40
|
Ruvinsky AM, Kozintsev AV. New and fast statistical-thermodynamic method for computation of protein-ligand binding entropy substantially improves docking accuracy. J Comput Chem 2005; 26:1089-95. [PMID: 15929088 DOI: 10.1002/jcc.20246] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We present a novel method to estimate the contributions of translational and rotational entropy to protein-ligand binding affinity. The method is based on estimates of the configurational integral through the sizes of clusters obtained from multiple docking positions. Cluster sizes are defined as the intervals of variation of center of ligand mass and Euler angles in the cluster. Then we suggest a method to consider the entropy of torsional motions. We validate the suggested methods on a set of 135 PDB protein-ligand complexes by comparing the averaged root-mean square deviations (RMSD) of the top-scored ligand docked positions, accounting and not accounting for entropy contributions, relative to the experimentally determined positions. We demonstrate that the method increases docking accuracy by 10-21% when used in conjunction with the AutoDock docking program, thus reducing the percent of incorrectly docked ligands by 1.4-fold to four-fold, so that in some cases the percent of ligands correctly docked to within an RMSD of 2 A is above 90%. We show that the suggested method to account for entropy of relative motions is identical to the method based on the Monte Carlo integration over intervals of variation of center of ligand mass and Euler angles in the cluster.
Collapse
Affiliation(s)
- A M Ruvinsky
- Force Field Laboratory, Algodign, LLC, B. Sadovaya, 8, 103379, Moscow, Russia.
| | | |
Collapse
|