1
|
Zhao M, Yu W, MacKerell AD. Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms. J Phys Chem B 2024; 128:7362-7375. [PMID: 39031121 PMCID: PMC11294009 DOI: 10.1021/acs.jpcb.4c03045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2024]
Abstract
In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.
Collapse
Affiliation(s)
- Mingtian Zhao
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| |
Collapse
|
2
|
Goel H, Hazel A, Yu W, Jo S, MacKerell AD. Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design. NEW J CHEM 2022; 46:919-932. [PMID: 35210743 PMCID: PMC8863107 DOI: 10.1039/d1nj04028f] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Site Identification by Ligand Competitive Saturation (SILCS) is a molecular simulation approach that uses diverse small solutes in aqueous solution to obtain functional group affinity patterns of a protein or other macromolecule. This involves employing a combined Grand Canonical Monte Carlo (GCMC)-molecular dynamics (MD) method to sample the full 3D space of the protein, including deep binding pockets and interior cavities from which functional group free energy maps (FragMaps) are obtained. The information content in the maps, which include contributions from protein flexibilty and both protein and functional group desolvation contributions, can be used in many aspects of the drug discovery process. These include identification of novel ligand binding pockets, including allosteric sites, pharmacophore modeling, prediction of relative protein-ligand binding affinities for database screening and lead optimization efforts, evaluation of protein-protein interactions as well as in the formulation of biologics-based drugs including monoclonal antibodies. The present article summarizes the various tools developed in the context of the SILCS methodology and their utility in computer-aided drug design (CADD) applications, showing how the SILCS toolset can improve the drug-development process on a number of fronts with respect to both accuracy and throughput representing a new avenue of CADD applications.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Anthony Hazel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States
| | - Sunhwan Jo
- SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20, Penn St. Baltimore, Maryland 21201, United States., SilcsBio LLC, 1100 Wicomico St. Suite 323, Baltimore, MD, 21230, United States.,, Tel: 410-706-7442, Fax: 410-706-5017
| |
Collapse
|
3
|
Yanagisawa K, Moriwaki Y, Terada T, Shimizu K. EXPRORER: Rational Cosolvent Set Construction Method for Cosolvent Molecular Dynamics Using Large-Scale Computation. J Chem Inf Model 2021; 61:2744-2753. [PMID: 34061535 DOI: 10.1021/acs.jcim.1c00134] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Cosolvent molecular dynamics (CMD) simulations involve an MD simulation of a protein in the presence of explicit water molecules mixed with cosolvent molecules to perform hotspot detection, binding site identification, and binding energy estimation, while other existing methods (e.g., MixMD, SILCS, and MDmix) utilize small molecules that represent functional groups of compounds. However, the cosolvent selections employed in these methods differ and there are only a few cosolvents that are commonly used in these methods. In this study, we proposed a systematic method for constructing a set of cosolvents for drug discovery, termed the EXtended PRObes set construction by REpresentative Retrieval (EXPRORER). First, we extracted typical substructures from FDA-approved drugs, generated 138 cosolvent structures, and for each cosolvent molecule, we conducted CMD simulations to generate a spatial probability distribution map of cosolvent atoms (PMAP). Analyses of PMAP similarity revealed that a cosolvent pair with a PMAP similarity greater than 0.70-0.75 shared similar structural features. We present a method for the construction of a cosolvent subset that satisfies a similarity threshold for all cosolvents, and we tested the constructed sets for four proteins. To our knowledge, this is the first study to include a systematic proposal for cosolvent set construction, and thus, the EXPRORER cosolvents will provide deeper insights into ligand binding sites of various proteins.
Collapse
Affiliation(s)
- Keisuke Yanagisawa
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Tokyo 152-8550, Japan.,Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Yoshitaka Moriwaki
- Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Tohru Terada
- Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan.,Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo 113-8657, Japan.,Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| | - Kentaro Shimizu
- Department of Biotechnology, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan.,Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo 113-8657, Japan.,Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo 113-8657, Japan
| |
Collapse
|
4
|
Nunes RS, Vila-Viçosa D, Costa PJ. Halogen Bonding: An Underestimated Player in Membrane–Ligand Interactions. J Am Chem Soc 2021; 143:4253-4267. [DOI: 10.1021/jacs.0c12470] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Rafael Santana Nunes
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
- Centro de Química Estrutural, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
| | - Diogo Vila-Viçosa
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| | - Paulo J. Costa
- BioISI - Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, Campo Grande, C8 bdg, 1749-016 Lisboa, Portugal
| |
Collapse
|
5
|
Tan YS, Verma CS. Straightforward Incorporation of Multiple Ligand Types into Molecular Dynamics Simulations for Efficient Binding Site Detection and Characterization. J Chem Theory Comput 2020; 16:6633-6644. [PMID: 32810406 DOI: 10.1021/acs.jctc.0c00405] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Binding site identification and characterization is an important initial step in structure-based drug design. To account for the effects of protein flexibility and solvation, several cosolvent molecular dynamics (MD) simulation methods that incorporate small organic molecules into the protein's solvent box to probe for binding sites have been developed. However, most of these methods are highly inefficient, as they allow for the use of only one probe type at a time, which means that multiple sets of simulations have to be performed to map different types of binding sites. The high probe concentrations used in some of these methods also necessitate the use of artificial repulsive forces to prevent the probes from aggregating. Here, we present multiple-ligand-mapping MD (mLMMD), a method that incorporates multiple types of probes for simultaneous and efficient mapping of different types of binding sites without the need for introduction of artificial forces that may cause unintended mapping artifacts. We validate the method on a diverse set of 10 proteins and show that the mLMMD probes are able to reliably identify hydrophobic, hydrogen-bonding, charged, and cryptic binding sites in all of the test cases. Our results also highlight the potential utility of mLMMD for virtual screening and rational drug design.
Collapse
Affiliation(s)
- Yaw Sing Tan
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Chandra S Verma
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671.,Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543.,School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551
| |
Collapse
|
6
|
Dong TG, Peng H, He XF, Wang X, Gao J. Hybrid Molecular Dynamics for Elucidating Cooperativity Between Halogen Bond and Water Molecules During the Interaction of p53-Y220C and the PhiKan5196 Complex. Front Chem 2020; 8:344. [PMID: 32457871 PMCID: PMC7221198 DOI: 10.3389/fchem.2020.00344] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 04/02/2020] [Indexed: 12/15/2022] Open
Abstract
The cooperativity between hydrogen and halogen bonds plays an important role in rational drug design. However, mimicking the dynamic cooperation between these bonds is a challenging issue, which has impeded the development of the halogen bond force field. In this study, the Y220C–PhiKan5196 complex of p53 protein was adopted as a model, and the functions of three water molecules that formed hydrogen bonds with halogen atoms were analyzed by the simulation method governed by the hybrid quantum mechanical/molecular mechanical molecular dynamics. A comparison with the water-free model revealed that the strength of the halogen bond in the complex was consistently stronger. This confirmed that the water molecules formed weak hydrogen bonds with the halogen atom and cooperated with the halogen atom to enhance the halogen bond. Further, it was discovered that the roles of the three water molecules were not the same. Therefore, the results obtained herein can facilitate a rational drug design. Further, this work emphasizes on the fact that, in addition to protein pockets and ligands, the role of voids should also be considered with regard to the water molecules surrounding them.
Collapse
Affiliation(s)
- Tian-Ge Dong
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Hui Peng
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xue-Feng He
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xiaocong Wang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| |
Collapse
|
7
|
Ustach VD, Lakkaraju SK, Jo S, Yu W, Jiang W, MacKerell AD. Optimization and Evaluation of Site-Identification by Ligand Competitive Saturation (SILCS) as a Tool for Target-Based Ligand Optimization. J Chem Inf Model 2019; 59:3018-3035. [PMID: 31034213 PMCID: PMC6597307 DOI: 10.1021/acs.jcim.9b00210] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Chemical fragment cosolvent sampling techniques have become a versatile tool in ligand-protein binding prediction. Site-identification by ligand competitive saturation (SILCS) is one such method that maps the distribution of chemical fragments on a protein as free energy fields called FragMaps. Ligands are then simulated via Monte Carlo techniques in the field of the FragMaps (SILCS-MC) to predict their binding conformations and relative affinities for the target protein. Application of SILCS-MC using a number of different scoring schemes and MC sampling protocols against multiple protein targets was undertaken to evaluate and optimize the predictive capability of the method. Seven protein targets and 551 ligands with broad chemical variability were used to evaluate and optimize the model to maximize Pearson's correlation coefficient, Pearlman's predictive index, correct relative binding affinity, and root-mean-square error versus the absolute experimental binding affinities. Across the protein-ligand sets, the relative affinities of the ligands were predicted correctly an average of 69% of the time for the highest overall SILCS protocol. Training the FragMap weighting factors using a Bayesian machine learning (ML) algorithm led to an increase to an average 75% relative correct affinity predictions. Furthermore, once the optimal protocol is identified for a specific protein-ligand system average predictabilities of 76% are achieved. The ML algorithm is successful with small training sets of data (30 or more compounds) due to the use of physically correct FragMap weights as priors. Notably, the 76% correct relative prediction rate is similar to or better than free energy perturbation methods that are significantly computationally more expensive than SILCS. The results further support the utility of SILCS as a powerful and computationally accessible tool to support lead optimization and development in drug discovery.
Collapse
Affiliation(s)
- Vincent D. Ustach
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| | | | - Sunhwan Jo
- SilcsBio, LLC, 8 Market Place, Suite 300, Baltimore, MD 21202
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| | - Wenjuan Jiang
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
- SilcsBio, LLC, 8 Market Place, Suite 300, Baltimore, MD 21202
| |
Collapse
|