1
|
Kalathingal M, Rhee YM. Thermodynamic consequences of stapling side-chains on a peptide ligand using a lactam-bridge: A theoretical study on anti-angiogenic peptides targeting VEGF. Proteins 2024; 92:959-974. [PMID: 38602129 DOI: 10.1002/prot.26692] [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: 01/19/2024] [Revised: 03/28/2024] [Accepted: 03/29/2024] [Indexed: 04/12/2024]
Abstract
Peptides are promising therapeutic agents for various biological targets due to their high efficacy and low toxicity, and the design of peptide ligands with high binding affinity to the target of interest is of utmost importance in peptide-based drug design. Introducing a conformational constraint to a flexible peptide ligand using a side-chain lactam-bridge is a convenient and efficient method to improve its binding affinity to the target. However, in general, such a small structural modification to a flexible ligand made with the intent of lowering the configurational entropic penalty for binding may have unintended consequences in different components of the binding enthalpy and entropy, including the configurational entropy component, which are still not clearly understood. Toward probing this, we examine different components of the binding enthalpy and entropy as well as the underlying structure and dynamics, for a side-chain lactam-bridged peptide inhibitor and its flexible analog forming complexes with vascular endothelial growth factor (VEGF), using all-atom molecular dynamics simulations. It is found that introducing a side-chain lactam-bridge constraint into the flexible peptide analog led to a gain in configurational entropy change but losses in solvation entropy, solute internal energy, and solvation energy changes upon binding, pinpointing the opportunities and challenges in drug design. The present study features an interplay between configurational and solvation entropy changes, as well as the one between binding enthalpy and entropy, in ligand-target binding upon imposing a conformational constraint into a flexible ligand.
Collapse
Affiliation(s)
- Mahroof Kalathingal
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| |
Collapse
|
2
|
Ahad Hossain M, Sultana S, Alanazi MM, Hadni H, Bhat AR, Hasan I, Kawsar SM. In vitro antimicrobial, anticancer evaluation, and in silico studies of mannopyranoside analogs against bacterial and fungal proteins: Acylation leads to improved antimicrobial activity. Saudi Pharm J 2024; 32:102093. [PMID: 38737807 PMCID: PMC11087236 DOI: 10.1016/j.jsps.2024.102093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2024] [Accepted: 04/28/2024] [Indexed: 05/14/2024] Open
Abstract
Carbohydrate analogs are an important, well-established class of clinically useful medicinal agents that exhibit potent antimicrobial activity. Thus, we explored the various therapeutic potential of methyl α-D-mannopyranoside (MαDM) analogs, including their ability to synthesize and assess their antibacterial, antifungal, and anticancer properties; additionally, molecular docking, molecular dynamics simulation, and ADMET analysis were performed. The structure of the synthesized MαDM analogs was ascertained by spectroscopic techniques and physicochemical and elemental analysis. In vitro antimicrobial activity was assessed and revealed significant inhibitory effects, particularly against gram-negative bacteria along with the prediction of activity spectra for substances (PASS). Concurrently, MαDM analogs showed good results against antifungal pathogens and exhibited promising anticancer effects in vitro, demonstrating dose-dependent cytotoxicity against Ehrlich ascites carcinoma (EAC) cancer cells while sparing normal cells from compound 5, with an IC50 of 4511.65 µg/mL according to the MTT colorimetric assay. A structure-activity relationship (SAR) study revealed that hexose combined with the acyl chains of decanoyl (C-10) and benzenesulfonyl (C6H5SO2-) had synergistic effects on the bacteria and fungi that were examined. Molecular docking was performed against the Escherichia coli (6KZV) and Candida albicans (1EAG) proteins to acquire insights into the molecular interactions underlying the observed biological activities. The docking results were further supported by 100 ns molecular dynamics simulations, which provided a dynamic view of the stability and flexibility of complexes involving MαDM and its targets. In addition, ADMET analysis was used to evaluate the toxicological and pharmacokinetic profiles. Owing to their promising drug-like properties, these MαDM analogs exhibit potential as prospective therapeutic candidates for future development.
Collapse
Affiliation(s)
- Md. Ahad Hossain
- Laboratory of Carbohydrate and Nucleoside Chemistry (LCNC), Department of Chemistry, Faculty of Science, University of Chittagong, Chittagong 4331, Bangladesh
| | - Shahin Sultana
- Laboratory of Carbohydrate and Nucleoside Chemistry (LCNC), Department of Chemistry, Faculty of Science, University of Chittagong, Chittagong 4331, Bangladesh
| | - Mohammed M. Alanazi
- Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, P.O. Box 2457, Riyadh 11451, Saudi Arabia
| | - Hanine Hadni
- LIMAS, Faculty of Sciences Dhar El Mahraz, Sidi Mohamed Ben Abdellah University, Fez, Morocco
| | - Ajmal R. Bhat
- Department of Chemistry, RTM Nagpur University, Nagpur 440033, India
| | - Imtiaj Hasan
- Department of Biochemistry and Molecular Biology, Faculty of Science, University of Rajshahi, Rajshahi 6205, Bangladesh
| | - Sarkar M.A. Kawsar
- Laboratory of Carbohydrate and Nucleoside Chemistry (LCNC), Department of Chemistry, Faculty of Science, University of Chittagong, Chittagong 4331, Bangladesh
| |
Collapse
|
3
|
Hoogstraten CA, Koenderink JB, van Straaten CE, Scheer-Weijers T, Smeitink JAM, Schirris TJJ, Russel FGM. Pyruvate dehydrogenase is a potential mitochondrial off-target for gentamicin based on in silico predictions and in vitro inhibition studies. Toxicol In Vitro 2024; 95:105740. [PMID: 38036072 DOI: 10.1016/j.tiv.2023.105740] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Revised: 11/08/2023] [Accepted: 11/22/2023] [Indexed: 12/02/2023]
Abstract
During the drug development process, organ toxicity leads to an estimated failure of one-third of novel chemical entities. Drug-induced toxicity is increasingly associated with mitochondrial dysfunction, but identifying the underlying molecular mechanisms remains a challenge. Computational modeling techniques have proven to be a good tool in searching for drug off-targets. Here, we aimed to identify mitochondrial off-targets of the nephrotoxic drugs tenofovir and gentamicin using different in silico approaches (KRIPO, ProBis and PDID). Dihydroorotate dehydrogenase (DHODH) and pyruvate dehydrogenase (PDH) were predicted as potential novel off-target sites for tenofovir and gentamicin, respectively. The predicted targets were evaluated in vitro, using (colorimetric) enzymatic activity measurements. Tenofovir did not inhibit DHODH activity, while gentamicin potently reduced PDH activity. In conclusion, the use of in silico methods appeared a valuable approach in predicting PDH as a mitochondrial off-target of gentamicin. Further research is required to investigate the contribution of PDH inhibition to overall renal toxicity of gentamicin.
Collapse
Affiliation(s)
- Charlotte A Hoogstraten
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands; Radboud Center for Mitochondrial Medicine, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands
| | - Jan B Koenderink
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands
| | - Carolijn E van Straaten
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands
| | - Tom Scheer-Weijers
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands
| | - Jan A M Smeitink
- Radboud Center for Mitochondrial Medicine, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands; Department of Pediatrics, Amalia Children's Hospital, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands; Khondrion BV, Nijmegen 6525 EX, the Netherlands
| | - Tom J J Schirris
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands; Radboud Center for Mitochondrial Medicine, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands
| | - Frans G M Russel
- Division of Pharmacology and Toxicology, Department of Pharmacy, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands; Radboud Center for Mitochondrial Medicine, Radboud University Medical Center, Nijmegen 6500 HB, the Netherlands.
| |
Collapse
|
4
|
Das S, Merz KM. Molecular Gas-Phase Conformational Ensembles. J Chem Inf Model 2024; 64:749-760. [PMID: 38206321 DOI: 10.1021/acs.jcim.3c01309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.
Collapse
Affiliation(s)
- Susanta Das
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| | - Kenneth M Merz
- Department of Chemistry, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824, United States
| |
Collapse
|
5
|
Shakour N, Hoseinpoor S, Rajabian F, Azimi SG, Iranshahi M, Sadeghi-Aliabadi H, Hadizadeh F. Discovery of non-peptide GLP-1r natural agonists for enhancing coronary safety in type 2 diabetes patients. J Biomol Struct Dyn 2024:1-18. [PMID: 38165453 DOI: 10.1080/07391102.2023.2298734] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 12/17/2023] [Indexed: 01/03/2024]
Abstract
This study explores the computational discovery of non-peptide agonists targeting the Glucagon-Like Peptide-1 Receptor (GLP-1R) to enhance the safety of major coronary outcomes in individuals affected by Type 2 Diabetes. The objective is to identify novel compounds that can activate the GLP-1R pathway without the limitations associated with peptide agonists. Type 2 diabetes mellitus (T2DM) is associated with an increased risk of cardiovascular disease (CVD) and mortality, which is attributed to the accumulation of fat in organs, including the heart. Glucagon-like peptide-1 receptor agonists (GLP-1RAs) are frequently used to manage T2DM and could potentially offer cardiovascular benefits. Therefore, this study examines non-peptide agonists of GLP-1R to improve coronary safety in type 2 diabetes patients. After rigorous assessments, two standout candidates were identified, with natural compound 12 emerging as the most promising. This study represents a notable advancement in enhancing the management of coronary outcomes among individuals with type 2 diabetes. The computational methodology employed successfully pinpointed potential GLP-1R natural agonists, providing optimism for the development of safer and more effective therapeutic interventions. Although computational methodologies have provided crucial insights, realizing the full potential of these compounds requires extensive experimental investigations, crucial in advancing therapeutic strategies for this critical patient population.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Neda Shakour
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
- Student Research Committee, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Saeideh Hoseinpoor
- Department of Biochemistry and Biophysics, Faculty of Sciences, Mashhad Branch, Islamic Azad University, Mashhad, Iran
- International UNESCO Center for Health-Related Basic Sciences and Human Nutrition, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Fatemeh Rajabian
- Department of Pharmacodynamics and Toxicology, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Sabikeh G Azimi
- Department of Chemistry, Faculty of Sciences, University of Birjand, Birjand, Iran
| | - Mehrdad Iranshahi
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Hojjat Sadeghi-Aliabadi
- Department of Pharmaceutical Chemistry, Isfahan University of Medical Sciences, Isfahan, Iran
| | - Farzin Hadizadeh
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| |
Collapse
|
6
|
Baillif B, Cole J, Giangreco I, McCabe P, Bender A. Applying atomistic neural networks to bias conformer ensembles towards bioactive-like conformations. J Cheminform 2023; 15:124. [PMID: 38129933 PMCID: PMC10740246 DOI: 10.1186/s13321-023-00794-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 12/10/2023] [Indexed: 12/23/2023] Open
Abstract
Identifying bioactive conformations of small molecules is an essential process for virtual screening applications relying on three-dimensional structure such as molecular docking. For most small molecules, conformer generators retrieve at least one bioactive-like conformation, with an atomic root-mean-square deviation (ARMSD) lower than 1 Å, among the set of low-energy conformers generated. However, there is currently no general method to prioritise these likely target-bound conformations in the ensemble. In this work, we trained atomistic neural networks (AtNNs) on 3D information of generated conformers of a curated subset of PDBbind ligands to predict the ARMSD to their closest bioactive conformation, and evaluated the early enrichment of bioactive-like conformations when ranking conformers by AtNN prediction. AtNN ranking was compared with bioactivity-unaware baselines such as ascending Sage force field energy ranking, and a slower bioactivity-based baseline ranking by ascending Torsion Fingerprint Deviation to the Maximum Common Substructure to the most similar molecule in the training set (TFD2SimRefMCS). On test sets from random ligand splits of PDBbind, ranking conformers using ComENet, the AtNN encoding the most 3D information, leads to early enrichment of bioactive-like conformations with a median BEDROC of 0.29 ± 0.02, outperforming the best bioactivity-unaware Sage energy ranking baseline (median BEDROC of 0.18 ± 0.02), and performing on a par with the bioactivity-based TFD2SimRefMCS baseline (median BEDROC of 0.31 ± 0.02). The improved performance of the AtNN and TFD2SimRefMCS baseline is mostly observed on test set ligands that bind proteins similar to proteins observed in the training set. On a more challenging subset of flexible molecules, the bioactivity-unaware baselines showed median BEDROCs up to 0.02, while AtNNs and TFD2SimRefMCS showed median BEDROCs between 0.09 and 0.13. When performing rigid ligand re-docking of PDBbind ligands with GOLD using the 1% top-ranked conformers, ComENet ranked conformers showed a higher successful docking rate than bioactivity-unaware baselines, with a rate of 0.48 ± 0.02 compared to CSD probability baseline with a rate of 0.39 ± 0.02. Similarly, on a pharmacophore searching experiment, selecting the 20% top-ranked conformers ranked by ComENet showed higher hit rate compared to baselines. Hence, the approach presented here uses AtNNs successfully to focus conformer ensembles towards bioactive-like conformations, representing an opportunity to reduce computational expense in virtual screening applications on known targets that require input conformations.
Collapse
Affiliation(s)
- Benoit Baillif
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1EW, UK
| | - Jason Cole
- Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge, CB2 1EZ, UK
| | - Ilenia Giangreco
- Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge, CB2 1EZ, UK
- Exscientia plc, The Schrödinger Building, Oxford Science Park, Oxford, OX4 4GE, UK
| | - Patrick McCabe
- Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge, CB2 1EZ, UK
| | - Andreas Bender
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge, CB2 1EW, UK.
| |
Collapse
|
7
|
Han F, Hao D, He X, Wang L, Niu T, Wang J. Distribution of Bound Conformations in Conformational Ensembles for X-ray Ligands Predicted by the ANI-2X Machine Learning Potential. J Chem Inf Model 2023; 63:6608-6618. [PMID: 37899502 PMCID: PMC10647024 DOI: 10.1021/acs.jcim.3c01350] [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] [Received: 08/22/2023] [Revised: 10/11/2023] [Accepted: 10/13/2023] [Indexed: 10/31/2023]
Abstract
In this study, we systematically studied the energy distribution of bioactive conformations of small molecular ligands in their conformational ensembles using ANI-2X, a machine learning potential, in conjunction with one of our recently developed geometry optimization algorithms, known as a conjugate gradient with backtracking line search (CG-BS). We first evaluated the combination of these methods (ANI-2X/CG-BS) using two molecule sets. For the 231-molecule set, ab initio calculations were performed at both the ωB97X/6-31G(d) and B3LYP-D3BJ/DZVP levels for accuracy comparison, while for the 8,992-molecule set, ab initio calculations were carried out at the B3LYP-D3BJ/DZVP level. For each molecule in the two molecular sets, up to 10 conformations were generated, which diminish the influence of individual outliers on the performance evaluation. Encouraged by the performance of ANI-2x/CG-BS in these evaluations, we calculated the energy distributions using ANI-2x/CG-BS for more than 27,000 ligands in the protein data bank (PDB). Each ligand has at least one conformation bound to a biological molecule, and this ligand conformation is labeled as a bound conformation. Besides the bound conformations, up to 200 conformations were generated using OpenEye's Omega2 software (https://docs.eyesopen.com/applications/ omega/) for each conformation. We performed a statistical analysis of how the bound conformation energies are distributed in the ensembles for 17,197 PDB ligands that have their bound conformation energies within the energy ranges of the Omega2-generated conformation ensembles. We found that half of the ligands have their relative conformation energy lower than 2.91 kcal/mol for the bound conformations in comparison with the global conformations, and about 90% of the bound conformations are within 10 kcal/mol above the global conformation energies. This information is useful to guide the construction of libraries for shape-based virtual screening and to improve the docking algorithm to efficiently sample bound conformations.
Collapse
Affiliation(s)
- Fengyang Han
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Dongxiao Hao
- School
of Electronics and Information Engineering, Ankang University, Ankang 725000, China
| | - Xibing He
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Luxuan Wang
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Taoyu Niu
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Junmei Wang
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
8
|
Jain AN, Brueckner AC, Jorge C, Cleves AE, Khandelwal P, Cortes JC, Mueller L. Complex peptide macrocycle optimization: combining NMR restraints with conformational analysis to guide structure-based and ligand-based design. J Comput Aided Mol Des 2023; 37:519-535. [PMID: 37535171 PMCID: PMC10505130 DOI: 10.1007/s10822-023-00524-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/20/2023] [Indexed: 08/04/2023]
Abstract
Systematic optimization of large macrocyclic peptide ligands is a serious challenge. Here, we describe an approach for lead-optimization using the PD-1/PD-L1 system as a retrospective example of moving from initial lead compound to clinical candidate. We show how conformational restraints can be derived by exploiting NMR data to identify low-energy solution ensembles of a lead compound. Such restraints can be used to focus conformational search for analogs in order to accurately predict bound ligand poses through molecular docking and thereby estimate ligand strain and protein-ligand intermolecular binding energy. We also describe an analogous ligand-based approach that employs molecular similarity optimization to predict bound poses. Both approaches are shown to be effective for prioritizing lead-compound analogs. Surprisingly, relatively small ligand modifications, which may have minimal effects on predicted bound pose or intermolecular interactions, often lead to large changes in estimated strain that have dominating effects on overall binding energy estimates. Effective macrocyclic conformational search is crucial, whether in the context of NMR-based restraints, X-ray ligand refinement, partial torsional restraint for docking/ligand-similarity calculations or agnostic search for nominal global minima. Lead optimization for peptidic macrocycles can be made more productive using a multi-disciplinary approach that combines biophysical data with practical and efficient computational methods.
Collapse
Affiliation(s)
- Ajay N Jain
- Research and Development, BioPharmics LLC, Sonoma County, CA, USA.
| | | | | | - Ann E Cleves
- Research and Development, BioPharmics LLC, Sonoma County, CA, USA
| | | | | | | |
Collapse
|
9
|
Seidel T, Permann C, Wieder O, Kohlbacher SM, Langer T. High-Quality Conformer Generation with CONFORGE: Algorithm and Performance Assessment. J Chem Inf Model 2023; 63:5549-5570. [PMID: 37624145 PMCID: PMC10498443 DOI: 10.1021/acs.jcim.3c00563] [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] [Received: 04/13/2023] [Indexed: 08/26/2023]
Abstract
Knowledge of the putative bound-state conformation of a molecule is an essential prerequisite for the successful application of many computer-aided drug design methods that aim to assess or predict its capability to bind to a particular target receptor. An established approach to predict bioactive conformers in the absence of receptor structure information is to sample the low-energy conformational space of the investigated molecules and derive representative conformer ensembles that can be expected to comprise members closely resembling possible bound-state ligand conformations. The high relevance of such conformer generation functionality led to the development of a wide panel of dedicated commercial and open-source software tools throughout the last decades. Several published benchmarking studies have shown that open-source tools usually lag behind their commercial competitors in many key aspects. In this work, we introduce the open-source conformer ensemble generator CONFORGE, which aims at delivering state-of-the-art performance for all types of organic molecules in drug-like chemical space. The ability of CONFORGE and several well-known commercial and open-source conformer ensemble generators to reproduce experimental 3D structures as well as their computational efficiency and robustness has been assessed thoroughly for both typical drug-like molecules and macrocyclic structures. For small molecules, CONFORGE clearly outperformed all other tested open-source conformer generators and performed at least equally well as the evaluated commercial generators in terms of both processing speed and accuracy. In the case of macrocyclic structures, CONFORGE achieved the best average accuracy among all benchmarked generators, with RDKit's generator coming close in second place.
Collapse
Affiliation(s)
- Thomas Seidel
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Christian Permann
- NeGeMac
Research Platform, Department of Pharmaceutical Sciences, Division
of Pharmaceutical Chemistry, University
of Vienna, Josef-Holaubek-Platz
2, 1090 Vienna, Austria
| | - Oliver Wieder
- Christian
Doppler Laboratory for Molecular Informatics in the Biosciences, Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Stefan M. Kohlbacher
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Thierry Langer
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
- NeGeMac
Research Platform, Department of Pharmaceutical Sciences, Division
of Pharmaceutical Chemistry, University
of Vienna, Josef-Holaubek-Platz
2, 1090 Vienna, Austria
- Christian
Doppler Laboratory for Molecular Informatics in the Biosciences, Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| |
Collapse
|
10
|
Azad I, Khan T, Ahmad N, Khan AR, Akhter Y. Updates on drug designing approach through computational strategies: a review. Future Sci OA 2023; 9:FSO862. [PMID: 37180609 PMCID: PMC10167725 DOI: 10.2144/fsoa-2022-0085] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 04/12/2023] [Indexed: 05/16/2023] Open
Abstract
The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.
Collapse
Affiliation(s)
- Iqbal Azad
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Tahmeena Khan
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Naseem Ahmad
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Abdul Rahman Khan
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Yusuf Akhter
- Department of Biotechnology, Babasaheb Bhimrao Ambedkar University, Vidya Vihar, Raebareli Road, Lucknow, UP, 2260025, India
| |
Collapse
|
11
|
Kemisetti D, Amin R, Alam F, Gacem A, Emran TB, Alsufyani T, Alqahtani MS, Islam S, Matin MM, Jameel M. Novel Benzothiazole Derivatives Synthesis and its Analysis as Diuretic Agents. EVIDENCE-BASED COMPLEMENTARY AND ALTERNATIVE MEDICINE 2023; 2023:1-12. [DOI: 10.1155/2023/5460563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
Benzothiazoles, an anticonvulsant, antiviral, antihypertensive, and cancer-fighting medication of the heterocyclic scaffold family, also acts as antibacterial and antiviral agents. There is much interest in this chemical’s production because of the strong and vital biological action it possesses. Substituted aromatic aldehydes were combined with 2-amino-benzothiazole-6-sulfonic acid amides, or Schiff base derivatives, to create Schiff base derivatives. Recrystallized, characterized, and tested for diuretic efficacy in vivo using online tools, m.p. (melting point), Rf, FTIR (Fourier transform infrared), 1H-NMR (proton nuclear magnetic resonance) data The molecular characteristics of all the substances created were estimated using Lipinski’s rule of 5, OSIRIS (software) molecular property explorer, Molsoft, and Autodock 4.0 docking software. Male Wistar rats were used to make all the compounds traditionally in order to test for diuretic activity. Neither the elemental nor the spectral information for the synthesized compounds disagreed. There were five different methods used to evaluate these compounds: Lipinski rule of five, Molsoft to determine molecular characteristics, PASS (prediction of activity spectra for substances) values to determine the diuretic effect, and OSIRIS software to determine toxicology. In order to investigate the diuretic effects of the selected drugs, docking analysis was used. Acetazolamide was shown to have a diuretic effect that was superior to that of compounds IIIb and IIIe, whereas 2-{(E)-[(3-hydroxyphenyl)methylidene]amino}-1,3-benzothiazole-6-sulfonamide (IIIb) was found to be the most promising potential.
Collapse
Affiliation(s)
- Durgaprasad Kemisetti
- Faculty of Pharmaceutical Science, Assam Down Town University, Panikhaiti, Guwahati, Assam, India
| | - Ruhul Amin
- Faculty of Pharmaceutical Science, Assam Down Town University, Panikhaiti, Guwahati, Assam, India
| | - Faruk Alam
- Faculty of Pharmaceutical Science, Assam Down Town University, Panikhaiti, Guwahati, Assam, India
| | - Amel Gacem
- Department of Physics, Faculty of Sciences, University 20 Août 1955, Skikda, Algeria
| | - Talha Bin Emran
- Department of Pharmacy, BGC Trust University Bangladesh, Chittagong 4381, Bangladesh
- Department of Pharmacy, Faculty of Allied Health Sciences, Daffodil International University, Dhaka 1207, Bangladesh
| | - Taghreed Alsufyani
- Department of Chemistry, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
| | - Mohammed S. Alqahtani
- Radiological Sciences Department, College of Applied Medical Sciences, King Khalid University, Abha 61421, Saudi Arabia
- BioImaging Unit, Space Research Centre, Michael Atiyah Building, University of Leicester, Leicester, LE1 7RH, UK
- Research Center for Advanced Materials Science (RCAMS), King Khalid University, Postcode: 9004, Zip Code: 61413, Abha, Saudi Arabia
| | - Saiful Islam
- Civil Engineering Department, College of Engineering, King Khalid University, Abha 61421, Saudi Arabia
| | - Mohammed Mahbubul Matin
- Bioorganic and Medicinal Chemistry Laboratory, Faculty of Science, Department of Chemistry, University of Chittagong, Chittagong 4331, Bangladesh
| | - Mohammed Jameel
- Department of Civil Engineering, College of Engineering, King Khalid University, Abha, Saudi Arabia
| |
Collapse
|
12
|
Liebschner D, Moriarty NW, Poon BK, Adams PD. In situ ligand restraints from quantum-mechanical methods. Acta Crystallogr D Struct Biol 2023; 79:100-110. [PMID: 36762856 PMCID: PMC9912925 DOI: 10.1107/s2059798323000025] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2022] [Accepted: 01/02/2023] [Indexed: 01/21/2023] Open
Abstract
In macromolecular crystallographic structure refinement, ligands present challenges for the generation of geometric restraints due to their large chemical variability, their possible novel nature and their specific interaction with the binding pocket of the protein. Quantum-mechanical approaches are useful for providing accurate ligand geometries, but can be plagued by the number of minima in flexible molecules. In an effort to avoid these issues, the Quantum Mechanical Restraints (QMR) procedure optimizes the ligand geometry in situ, thus accounting for the influence of the macromolecule on the local energy minima of the ligand. The optimized ligand geometry is used to generate target values for geometric restraints during the crystallographic refinement. As demonstrated using a sample of >2330 ligand instances in >1700 protein-ligand models, QMR restraints generally result in lower deviations from the target stereochemistry compared with conventionally generated restraints. In particular, the QMR approach provides accurate torsion restraints for ligands and other entities.
Collapse
Affiliation(s)
- Dorothee Liebschner
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Nigel W. Moriarty
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Billy K. Poon
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| | - Paul D. Adams
- Molecular Biosciences and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
- Department of Bioengineering, University of California, Berkeley, CA 94720, USA
| |
Collapse
|
13
|
Jain AN, Brueckner AC, Cleves AE, Reibarkh M, Sherer EC. A Distributional Model of Bound Ligand Conformational Strain: From Small Molecules up to Large Peptidic Macrocycles. J Med Chem 2023; 66:1955-1971. [PMID: 36701387 PMCID: PMC9923749 DOI: 10.1021/acs.jmedchem.2c01744] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
The internal conformational strain incurred by ligands upon binding a target site has a critical impact on binding affinity, and expectations about the magnitude of ligand strain guide conformational search protocols. Estimates for bound ligand strain begin with modeled ligand atomic coordinates from X-ray co-crystal structures. By deriving low-energy conformational ensembles to fit X-ray diffraction data, calculated strain energies are substantially reduced compared with prior approaches. We show that the distribution of expected global strain energy values is dependent on molecular size in a superlinear manner. The distribution of strain energy follows a rectified normal distribution whose mean and variance are related to conformational complexity. The modeled strain distribution closely matches calculated strain values from experimental data comprising over 3000 protein-ligand complexes. The distributional model has direct implications for conformational search protocols as well as for directions in molecular design.
Collapse
Affiliation(s)
- Ajay N. Jain
- Research
& Development, BioPharmics LLC, Sonoma County, California95404, United States,
| | - Alexander C. Brueckner
- Molecular
Structure & Design, Bristol Myers Squibb, Princeton, New Jersey08543, United States
| | - Ann E. Cleves
- Research
& Development, BioPharmics LLC, Sonoma County, California95404, United States
| | - Mikhail Reibarkh
- Analytical
Research and Development, Merck & Co.
Inc., Rahway, New Jersey07065, United States
| | - Edward C. Sherer
- Analytical
Research and Development, Merck & Co.
Inc., Rahway, New Jersey07065, United States,
| |
Collapse
|
14
|
Breznik M, Ge Y, Bluck JP, Briem H, Hahn DF, Christ CD, Mortier J, Mobley DL, Meier K. Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods. ChemMedChem 2023; 18:e202200425. [PMID: 36240514 PMCID: PMC9868080 DOI: 10.1002/cmdc.202200425] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 10/10/2022] [Indexed: 01/26/2023]
Abstract
Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.
Collapse
Affiliation(s)
- Marko Breznik
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Yunhui Ge
- Department of Pharmaceutical Sciences, University of California, Irvine, CA 92697, USA
| | - Joseph P. Bluck
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Hans Briem
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - David F. Hahn
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | - Clara D. Christ
- Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Jérémie Mortier
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - David L. Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, CA 92697, USA,Department of Chemistry, University of California, Irvine, CA 92697, USA
| | - Katharina Meier
- Computational Life Science Technology Functions, Crop Science, R&D, Bayer AG, 40789 Monheim, Germany
| |
Collapse
|
15
|
Meixner M, Zachmann M, Metzler S, Scheerer J, Zacharias M, Antes I. Dynamic Docking of Macrocycles in Bound and Unbound Protein Structures with DynaDock. J Chem Inf Model 2022; 62:3426-3441. [PMID: 35796228 DOI: 10.1021/acs.jcim.2c00436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Macrocycles are interesting molecules with unique features due to their conformationally constrained yet flexible ring structure. This characteristic poses a difficult challenge for computational modeling studies since they rely on accurate structural descriptions. In particular, molecular docking calculations suffer from the lack of ring flexibility during pose generation, which is often compensated by using pregenerated ligand conformer ensembles. Moreover, receptor structures are mainly treated rigidly, which limits the use of many docking tools. In this study, we optimized our previous molecular dynamics-based sampling and docking pipeline specifically designed for the accurate prediction of macrocyclic compounds. We developed a dihedral classification procedure for in-depth conformational analysis of the macrocyclic rings and extracted structural ensembles that were subsequently docked in both bound and unbound protein structures employing a fully flexible approach. Our results suggest that including a ring conformer close to the bound state in the starting ensemble increases the chance of successful docking. The bioactive conformations of a diverse set of ligands could be predicted with high and decent accuracy in bound and unbound protein structures, respectively, due to the incorporation of full molecular flexibility in our approach. The remaining unsuccessful docking calculations were mainly caused by large flexible substituents that bind to surface-exposed binding sites, rather than the macrocyclic ring per se and could be further improved by explicit molecular dynamics simulations of the docked complex.
Collapse
Affiliation(s)
- Maximilian Meixner
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Martin Zachmann
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Sebastian Metzler
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Jonathan Scheerer
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Martin Zacharias
- Center of Functional Protein Assemblies, Technical University Munich, Ernst-Otto-Fischer-Straße 8, Garching bei München 85748, Germany
| | - Iris Antes
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| |
Collapse
|
16
|
Xie W, Wang F, Li Y, Lai L, Pei J. Advances and Challenges in De Novo Drug Design Using Three-Dimensional Deep Generative Models. J Chem Inf Model 2022; 62:2269-2279. [PMID: 35544331 DOI: 10.1021/acs.jcim.2c00042] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A persistent goal for de novo drug design is to generate novel chemical compounds with desirable properties in a labor-, time-, and cost-efficient manner. Deep generative models provide alternative routes to this goal. Numerous model architectures and optimization strategies have been explored in recent years, most of which have been developed to generate two-dimensional molecular structures. Some generative models aiming at three-dimensional (3D) molecule generation have also been proposed, gaining attention for their unique advantages and potential to directly design drug-like molecules in a target-conditioning manner. This review highlights current developments in 3D molecular generative models combined with deep learning and discusses future directions for de novo drug design.
Collapse
Affiliation(s)
- Weixin Xie
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Fanhao Wang
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Yibo Li
- Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Luhua Lai
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China.,Peking-Tsinghua Center for Life Science at BNLMS, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jianfeng Pei
- Center for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| |
Collapse
|
17
|
Kobauri P, Galenkamp NS, Schulte AM, de Vries J, Simeth NA, Maglia G, Thallmair S, Kolarski D, Szymanski W, Feringa BL. Hypothesis-Driven, Structure-Based Design in Photopharmacology: The Case of eDHFR Inhibitors. J Med Chem 2022; 65:4798-4817. [PMID: 35258959 PMCID: PMC8958501 DOI: 10.1021/acs.jmedchem.1c01962] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
Photopharmacology
uses light to regulate the biological activity
of drugs. This precise control is obtained through the incorporation
of molecular photoswitches into bioactive molecules. A major challenge
for photopharmacology is the rational design of photoswitchable drugs
that show light-induced activation. Computer-aided drug design is
an attractive approach toward more effective, targeted design. Herein,
we critically evaluated different structure-based approaches for photopharmacology
with Escherichia coli dihydrofolate reductase (eDHFR)
as a case study. Through the iterative examination of our hypotheses,
we progressively tuned the design of azobenzene-based, photoswitchable
eDHFR inhibitors in five design–make–switch–test–analyze
cycles. Targeting a hydrophobic subpocket of the enzyme and a specific
salt bridge only with the thermally metastable cis-isomer emerged as the most promising design strategy. We identified
three inhibitors that could be activated upon irradiation and reached
potencies in the low-nanomolar range. Above all, this systematic study
provided valuable insights for future endeavors toward rational photopharmacology.
Collapse
Affiliation(s)
- Piermichele Kobauri
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Nicole S Galenkamp
- Groningen Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Albert M Schulte
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Jisk de Vries
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Nadja A Simeth
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Institute for Organic and Biomolecular Chemistry, University of Goettingen, Tammannstr. 2, 37077 Göttingen, Germany
| | - Giovanni Maglia
- Groningen Biomolecular Science and Biotechnology Institute, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Sebastian Thallmair
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands.,Frankfurt Institute for Advanced Studies, Ruth-Moufang-Straße 1, 60438 Frankfurt am Main, Germany
| | - Dušan Kolarski
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,DWI-Leibniz Institut für interaktive Materialien e.V., RWTH Aachen University, Forckenbeckstraße 50, 52074 Aachen, Germany
| | - Wiktor Szymanski
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Department of Radiology, Medical Imaging Center, University of Groningen, University Medical Center Groningen, Hanzeplein 1, 9713 GZ Groningen, The Netherlands
| | - Ben L Feringa
- Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
18
|
Brown BP, Vu O, Geanes AR, Kothiwale S, Butkiewicz M, Lowe EW, Mueller R, Pape R, Mendenhall J, Meiler J. Introduction to the BioChemical Library (BCL): An Application-Based Open-Source Toolkit for Integrated Cheminformatics and Machine Learning in Computer-Aided Drug Discovery. Front Pharmacol 2022; 13:833099. [PMID: 35264967 PMCID: PMC8899505 DOI: 10.3389/fphar.2022.833099] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 01/24/2022] [Indexed: 01/31/2023] Open
Abstract
The BioChemical Library (BCL) cheminformatics toolkit is an application-based academic open-source software package designed to integrate traditional small molecule cheminformatics tools with machine learning-based quantitative structure-activity/property relationship (QSAR/QSPR) modeling. In this pedagogical article we provide a detailed introduction to core BCL cheminformatics functionality, showing how traditional tasks (e.g., computing chemical properties, estimating druglikeness) can be readily combined with machine learning. In addition, we have included multiple examples covering areas of advanced use, such as reaction-based library design. We anticipate that this manuscript will be a valuable resource for researchers in computer-aided drug discovery looking to integrate modular cheminformatics and machine learning tools into their pipelines.
Collapse
Affiliation(s)
- Benjamin P. Brown
- Chemical and Physical Biology Program, Medical Scientist Training Program, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
- *Correspondence: Jens Meiler, ; Jeffrey Mendenhall, ; Benjamin P. Brown,
| | - Oanh Vu
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Alexander R. Geanes
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Sandeepkumar Kothiwale
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Mariusz Butkiewicz
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Edward W. Lowe
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Ralf Mueller
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Richard Pape
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | - Jeffrey Mendenhall
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
- *Correspondence: Jens Meiler, ; Jeffrey Mendenhall, ; Benjamin P. Brown,
| | - Jens Meiler
- Department of Chemistry, Departments of Pharmacology and Biomedical Informatics, Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
- Institute for Drug Discovery, Leipzig University Medical School, Leipzig, Germany
- *Correspondence: Jens Meiler, ; Jeffrey Mendenhall, ; Benjamin P. Brown,
| |
Collapse
|
19
|
Trapl D, Krupička M, Višňovský V, Hozzová J, Ol'ha J, Křenek A, Spiwok V. Property Map Collective Variable as a Useful Tool for a Force Field Correction. J Chem Inf Model 2022; 62:567-576. [PMID: 35112877 DOI: 10.1021/acs.jcim.1c00651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accuracy of biomolecular simulations depends on the accuracy of an empirical molecular mechanics potential known as a force field: a set of parameters and expressions to estimate the potential from atomic coordinates. Accurate parametrization of force fields for small organic molecules is a challenge due to their high diversity. One of the possible approaches is to apply a correction to the existing force fields. Here, we propose an approach to estimate the density functional theory (DFT)-derived force field correction which is calculated during the run of molecular dynamics without significantly affecting its speed. Using the formula known as a property map collective variable, we approximate the force field correction by a weighted average of this force field correction calculated only for a small series of reference structures. To validate this method, we used seven AMBER force fields, and we show how it is possible to convert one force field to behave like the other one. We also present the force field correction for the important anticancer drug Imatinib as a use case example. Our method appears to be suitable for adjusting the force field for general drug-like molecules. We provide a pipeline that generates the correction; this pipeline is available at https://pmcvff-correction.cerit-sc.cz/.
Collapse
Affiliation(s)
- Dalibor Trapl
- Department of Biochemistry and Microbiology, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Martin Krupička
- Department of Organic Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Vladimír Višňovský
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Jana Hozzová
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Jaroslav Ol'ha
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Aleš Křenek
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Vojtěch Spiwok
- Department of Biochemistry and Microbiology, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| |
Collapse
|
20
|
Calculation of Crystal-Solution Dissociation Constants. Biomolecules 2022; 12:biom12020147. [PMID: 35204648 PMCID: PMC8961641 DOI: 10.3390/biom12020147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 01/08/2022] [Accepted: 01/10/2022] [Indexed: 11/17/2022] Open
Abstract
The calculation of dissociation constants is an important problem in molecular biophysics. For such a calculation, it is important to correctly calculate both terms of the binding free energy; that is, the enthalpy and entropy of binding. Both these terms can be computed using molecular dynamics simulations, but this approach is very computationally expensive, and entropy calculations are especially slow. We develop an alternative very fast method of calculating the binding entropy and dissociation constants. The main part of our approach is based on the evaluation of movement ranges of molecules in the bound state. Then, the range of molecular movements in the bound state (here, in molecular crystals) is used for the calculation of the binding entropies and, then (using, in addition, the experimentally measured sublimation enthalpies), the crystal-to-vapor dissociation constants. Previously, we considered the process of the reversible sublimation of small organic molecules from crystals to vapor. In this work, we extend our approach by considering the dissolution of molecules, in addition to their sublimation. Similar to the sublimation case, our method shows a good correlation with experimentally measured dissociation constants at the dissolution of crystals.
Collapse
|
21
|
Drug Discovery for Mycobacterium tuberculosis Using Structure-Based Computer-Aided Drug Design Approach. Int J Mol Sci 2021; 22:ijms222413259. [PMID: 34948055 PMCID: PMC8703488 DOI: 10.3390/ijms222413259] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Revised: 11/09/2021] [Accepted: 11/14/2021] [Indexed: 12/12/2022] Open
Abstract
Developing new, more effective antibiotics against resistant Mycobacterium tuberculosis that inhibit its essential proteins is an appealing strategy for combating the global tuberculosis (TB) epidemic. Finding a compound that can target a particular cavity in a protein and interrupt its enzymatic activity is the crucial objective of drug design and discovery. Such a compound is then subjected to different tests, including clinical trials, to study its effectiveness against the pathogen in the host. In recent times, new techniques, which involve computational and analytical methods, enhanced the chances of drug development, as opposed to traditional drug design methods, which are laborious and time-consuming. The computational techniques in drug design have been improved with a new generation of software used to develop and optimize active compounds that can be used in future chemotherapeutic development to combat global tuberculosis resistance. This review provides an overview of the evolution of tuberculosis resistance, existing drug management, and the design of new anti-tuberculosis drugs developed based on the contributions of computational techniques. Also, we show an appraisal of available software and databases on computational drug design with an insight into the application of this software and databases in the development of anti-tubercular drugs. The review features a perspective involving machine learning, artificial intelligence, quantum computing, and CRISPR combination with available computational techniques as a prospective pathway to design new anti-tubercular drugs to combat resistant tuberculosis.
Collapse
|
22
|
Foloppe N, Chen IJ. The reorganization energy of compounds upon binding to proteins, from dynamic and solvated bound and unbound states. Bioorg Med Chem 2021; 51:116464. [PMID: 34798378 DOI: 10.1016/j.bmc.2021.116464] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 09/26/2021] [Accepted: 10/01/2021] [Indexed: 11/30/2022]
Abstract
The intramolecular reorganization energy (ΔEReorg) of compounds upon binding to proteins is a component of the binding free energy, which has long received particular attention, for fundamental and practical reasons. Understanding ΔEReorg would benefit the science of molecular recognition and drug design. For instance, the tolerable strain energy of compounds upon binding has been elusive. Prior studies found some large ΔEReorg values (e.g. > 10 kcal/mol), received with skepticism since they imply excessive opposition to binding. Indeed, estimating ΔEReorg is technically difficult. Typically, ΔEReorg has been approached by taking two energy-minimized conformers representing the bound and unbound states, and subtracting their conformational energy. This is a drastic oversimplification, liable to conformational collapse of the unbound conformer. Instead, the present work applies extensive molecular dynamics (MD) and the modern OPLS3 force-field to simulate compounds bound and unbound states, in explicit solvent under physically relevant conditions. The thermalized unbound compounds populate multiple conformations, not reducible to one or a few energy-minimized conformers. The intramolecular energies in the bound and unbound states were averaged over pertinent conformational ensembles, and the reorganization enthalpy upon binding (ΔHReorg) deduced by subtraction. This was applied to 76 systems, including 43 approved drugs, carefully selected for i) the quality of the bioactive X-ray structures and ii) the diversity of the chemotypes, their properties and protein targets. It yielded comparatively low ΔHReorg values (median = 1.4 kcal/mol, mean = 3.0 kcal/mol). A new finding is the observation of negative ΔHReorg values. Indeed, reorganization energies do not have to oppose binding, e.g. when intramolecular interactions stabilize preferentially the bound state. Conversely, even with competing water molecules, intramolecular interactions can occur predominantly in the unbound compound, and be replaced by intermolecular counterparts upon protein binding. Such disruption of intramolecular interactions upon binding gives rise to occasional larger ΔHReorg values. Such counterintuitive larger ΔHReorg values may be rationalized as a redistribution of interactions upon binding, qualitatively compatible with binding.
Collapse
Affiliation(s)
- Nicolas Foloppe
- Vernalis (R&D) Ltd, Granta Park, Abington, Cambridge CB21 6GB, UK.
| | - I-Jen Chen
- Vernalis (R&D) Ltd, Granta Park, Abington, Cambridge CB21 6GB, UK.
| |
Collapse
|
23
|
Gu S, Smith MS, Yang Y, Irwin JJ, Shoichet BK. Ligand Strain Energy in Large Library Docking. J Chem Inf Model 2021; 61:4331-4341. [PMID: 34467754 DOI: 10.1021/acs.jcim.1c00368] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
While small molecule internal strain is crucial to molecular docking, using it in evaluating ligand scores has remained elusive. Here, we investigate a technique that calculates strain using relative torsional populations in the Cambridge Structural Database, enabling fast precalculation of these energies. In retrospective studies of large docking screens of the dopamine D4 receptor and of AmpC β-lactamase, where close to 600 docking hits were tested experimentally, including such strain energies improved hit rates by preferentially reducing the ranks of strained high-scoring decoy molecules. In a 40-target subset of the DUD-E benchmark, we found two thresholds that usefully distinguished between ligands and decoys: one based on the total strain energy of the small molecules and another based on the maximum strain allowed for any given torsion within them. Using these criteria, about 75% of the benchmark targets had improved enrichment after strain filtering. Relying on precalculated population distributions, this approach is rapid, taking less than 0.04 s to evaluate a conformation on a standard core, making it pragmatic for precalculating strain in even ultralarge libraries. Since it is scoring function agnostic, it may be useful to multiple docking approaches; it is openly available at http://tldr.docking.org.
Collapse
Affiliation(s)
- Shuo Gu
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States
| | - Matthew S Smith
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States.,Program of Biophysics, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States
| | - Ying Yang
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States
| | - John J Irwin
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States
| | - Brian K Shoichet
- Department of Pharmaceutical Chemistry, University of California, San Francisco, 1700 Fourth Street, San Francisco, California 94143-2550, United States
| |
Collapse
|
24
|
Castelli M, Serapian SA, Marchetti F, Triveri A, Pirota V, Torielli L, Collina S, Doria F, Freccero M, Colombo G. New perspectives in cancer drug development: computational advances with an eye to design. RSC Med Chem 2021; 12:1491-1502. [PMID: 34671733 PMCID: PMC8459323 DOI: 10.1039/d1md00192b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Accepted: 07/06/2021] [Indexed: 02/06/2023] Open
Abstract
Computational chemistry has come of age in drug discovery. Indeed, most pharmaceutical development programs rely on computer-based data and results at some point. Herein, we discuss recent applications of advanced simulation techniques to difficult challenges in drug discovery. These entail the characterization of allosteric mechanisms and the identification of allosteric sites or cryptic pockets determined by protein motions, which are not immediately evident in the experimental structure of the target; the study of ligand binding mechanisms and their kinetic profiles; and the evaluation of drug-target affinities. We analyze different approaches to tackle challenging and emerging biological targets. Finally, we discuss the possible perspectives of future application of computation in drug discovery.
Collapse
Affiliation(s)
- Matteo Castelli
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Stefano A Serapian
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Filippo Marchetti
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Alice Triveri
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Valentina Pirota
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Luca Torielli
- Department of Drug Sciences, Medicinal Chemistry and Pharmaceutical Technology Section, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Simona Collina
- Department of Drug Sciences, Medicinal Chemistry and Pharmaceutical Technology Section, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Filippo Doria
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Mauro Freccero
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| | - Giorgio Colombo
- Department of Chemistry, University of Pavia via Taramelli 12 27100 Pavia Italy
| |
Collapse
|
25
|
Mammoli A, Riccio A, Bianconi E, Coletti A, Camaioni E, Macchiarulo A. One Key and Multiple Locks: Substrate Binding in Structures of Tryptophan Dioxygenases and Hydroxylases. ChemMedChem 2021; 16:2732-2743. [PMID: 34137184 PMCID: PMC8518741 DOI: 10.1002/cmdc.202100312] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2021] [Revised: 06/14/2021] [Indexed: 12/18/2022]
Abstract
Since its discovery at the beginning of the past century, the essential nutrient l-Tryptophan (l-Trp) and its catabolic pathways have acquired an increasing interest in an ever wider scientific community for their pivotal roles in underlying many important physiological functions and associated pathological conditions. As a consequence, enzymes catalyzing rate limiting steps along l-Trp catabolic pathways - including IDO1, TDO, TPH1 and TPH2 - have turned to be interesting drug targets for the design and development of novel therapeutic agents for different disorders such as carcinoid syndrome, cancer and autoimmune diseases. This article provides a fresh comparative overview on the most recent advancements that crystallographic studies, biophysical and computational works have brought on structural aspects and molecular recognition patterns of these enzymes toward l-Trp. Finally, a conformational analysis of l-Trp is also discussed as part of the molecular recognition process governing the binding of a substrate to its cognate enzymes.
Collapse
Affiliation(s)
- Andrea Mammoli
- Department of Pharmaceutical SciencesUniversity of PerugiaVia del Liceo N. 106123PerugiaItaly
| | - Alessandra Riccio
- Department of Pharmaceutical SciencesUniversity of PerugiaVia del Liceo N. 106123PerugiaItaly
| | - Elisa Bianconi
- Department of Pharmaceutical SciencesUniversity of PerugiaVia del Liceo N. 106123PerugiaItaly
| | - Alice Coletti
- Department of Medicine and SurgeryUniversity of PerugiaP. le Gambuli06132PerugiaItaly
| | - Emidio Camaioni
- Department of Pharmaceutical SciencesUniversity of PerugiaVia del Liceo N. 106123PerugiaItaly
| | - Antonio Macchiarulo
- Department of Pharmaceutical SciencesUniversity of PerugiaVia del Liceo N. 106123PerugiaItaly
| |
Collapse
|
26
|
Seep L, Bonin A, Meier K, Diedam H, Göller AH. Ensemble completeness in conformer sampling: the case of small macrocycles. J Cheminform 2021; 13:55. [PMID: 34325738 PMCID: PMC8320181 DOI: 10.1186/s13321-021-00524-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 06/05/2021] [Indexed: 11/18/2022] Open
Abstract
In this study we compare the three algorithms for the generation of conformer ensembles Biovia BEST, Schrödinger Prime macrocycle sampling (PMM) and Conformator (CONF) form the University of Hamburg, with ensembles derived for exhaustive molecular dynamics simulations applied to a dataset of 7 small macrocycles in two charge states and three solvents. Ensemble completeness is a prerequisite to allow for the selection of relevant diverse conformers for many applications in computational chemistry. We apply conformation maps using principal component analysis based on ring torsions. Our major finding critical for all applications of conformer ensembles in any computational study is that maps derived from MD with explicit solvent are significantly distinct between macrocycles, charge states and solvents, whereas the maps for post-optimized conformers using implicit solvent models from all generator algorithms are very similar independent of the solvent. We apply three metrics for the quantification of the relative covered ensemble space, namely cluster overlap, variance statistics, and a novel metric, Mahalanobis distance, showing that post-optimized MD ensembles cover a significantly larger conformational space than the generator ensembles, with the ranking PMM > BEST >> CONF. Furthermore, we find that the distributions of 3D polar surface areas are very similar for all macrocycles independent of charge state and solvent, except for the smaller and more strained compound 7, and that there is also no obvious correlation between 3D PSA and intramolecular hydrogen bond count distributions.
Collapse
Affiliation(s)
- Lea Seep
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Anne Bonin
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Katharina Meier
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Holger Diedam
- Engineering & Technology, Applied Mathematics, Bayer AG, 51368, Leverkusen, Germany
| | - Andreas H Göller
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany.
| |
Collapse
|
27
|
Liebeschuetz JW. The Good, the Bad, and the Twisted Revisited: An Analysis of Ligand Geometry in Highly Resolved Protein-Ligand X-ray Structures. J Med Chem 2021; 64:7533-7543. [PMID: 34060310 DOI: 10.1021/acs.jmedchem.1c00228] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An analysis of the rotatable bond geometry of drug-like ligand models is reported for high-resolution (<1.1 Å) crystallographic protein-ligand complexes. In cases where the ligand fit to the electron density is very good, unusual torsional geometry is rare and, most often, though not exclusively, associated with strong polar, metal, or covalent ligand-protein interactions. It is rarely associated with a torsional strain of greater than 2 kcal mol-1 by calculation. An unusual torsional geometry is more prevalent where the fit to electron density is not perfect. Multiple low-strain conformer bindings were observed in 21% of the set and, it is suggested, may also lie behind many of the 35% of single-occupancy cases, where a poor fit to the e-density was found. It is concluded that multiple conformer ligand binding is an under-recognized phenomenon in structure-based drug design and that there is a need for more robust crystallographic refinement methods to better handle such cases.
Collapse
Affiliation(s)
- John W Liebeschuetz
- Skilos Chemoinformatics, 159 Water Street, Cambridge CB4 1PB, United Kingdom.,Astex Pharmaceuticals, 436 Cambridge Science Park, Milton Road, Milton, Cambridge CB4 0QA, United Kingdom
| |
Collapse
|
28
|
Roy S, Ali A, Bhattacharya S. Theoretical Insight into the Library Screening Approach for Binding of Intermolecular G-Quadruplex RNA and Small Molecules through Docking and Molecular Dynamics Simulation Studies. J Phys Chem B 2021; 125:5489-5501. [PMID: 34029082 DOI: 10.1021/acs.jpcb.0c10991] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The interactions of intermolecular G-quadruplex RNA and small molecules have been investigated by computational studies. Various anthraquinone, bisbenzimidazole, and carbazole-benzimidazole based ligands have shown a distinct preference to G-quadruplex structures as opposed to the corresponding duplex forms of DNA that were docked with telomeric G-quadruplex RNA. The comparative binding study of such ligands with G-quadruplex (G4) RNA showed higher binding affinities toward carbazole-benzimidazole ligands than those of the anthraquinone and bisbenzimidazole based ligands. A molecular dynamics simulation study was used to examine quadruplex-ligand interactions. Analysis of the binding free energy indicated the formation of the thermodynamically favorable RNA-ligand complex. The formation of several H-bonding interactions and the change of the solvent accessible surface area (SASA) also support the effective binding of the carbazole-benzimidazole ligands with G4 RNA structures. Thus, the library screening approach has assisted in getting a structure-activity relationship for the selected small molecules toward the G-quadruplex RNA binding, which can be applied in the targeting of G-quadruplex RNA medicated anticancer therapeutics.
Collapse
Affiliation(s)
- Soma Roy
- Department of Organic Chemistry, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | - Asfa Ali
- Department of Organic Chemistry, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | - Santanu Bhattacharya
- Department of Organic Chemistry, Indian Institute of Science, Bangalore, Karnataka 560012, India.,School of Applied & Interdisciplinary Sciences, Indian Association for the Cultivation of Science, Kolkata, West Bengal 700032, India
| |
Collapse
|
29
|
Nutschel C, Coscolín C, David B, Mulnaes D, Ferrer M, Jaeger KE, Gohlke H. Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones. J Chem Inf Model 2021; 61:2383-2395. [PMID: 33949194 DOI: 10.1021/acs.jcim.1c00152] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.
Collapse
Affiliation(s)
- Christina Nutschel
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Cristina Coscolín
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Benoit David
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Daniel Mulnaes
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Manuel Ferrer
- Institute of Catalysis, Consejo Superior de Investigaciones Científicas, 28049 Madrid, Spain
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, 52425 Jülich, Germany.,Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Holger Gohlke
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), Institute of Biological Information Processing (IBI-7: Structural Biochemistry), and Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.,Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
30
|
Heal JR, Sheridan JM, Kulkarni M, Liptrott N, Egan D, Khoo S, Matthews I, Hamilton WDO. AUTOSTERE: Systematic Search for Scaffold Replacement Opportunities within Structural Databases. J Chem Inf Model 2021; 61:1778-1788. [PMID: 33830752 DOI: 10.1021/acs.jcim.0c01261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Medicinal chemists often bias toward working with scaffolds with which previously they have had direct experience and successes. In this way, it is often the case that scaffolds which have proven tractable within a research group are "reused" across multiple and sometimes unrelated drug targets. With this concept in mind, we designed a new computer algorithm AUTOSTERE which could systematically assess the opportunities to replace any part of any molecule within an entire database of known ligand structures with a target scaffold and automatically evaluate the potential designs in the context of the original ligand's protein environment. As such, it performs scaffold replacement on an unprecedented scale and suggests new target opportunities for preferred chemistries rather than the conventional reverse situation. The results of this approach for one scaffold, a substituted triazolinone, applied to a set of 10 426 ligand conformations extracted from the PDB are described. This led to the identification of ∼600 novel ligands incorporating the triazolinone scaffolds in complex with their predicted drug targets. From these, design examples are provided for HSP-90, cathepsin K, and TIE-2 kinase. A further study involved the searching for possible drug targets for unusual pyridopyrimidine cores. This process resulted in the identification of potential novel HIV reverse transcriptase inhibitors which were synthesized and shown to exhibit similar in vitro potencies to marketed compounds. Overall, the methodology described provides a powerful new approach to identify new target opportunities for scaffolds of provenance.
Collapse
Affiliation(s)
- Jonathan R Heal
- RxCelerate Limited, Dorothy Hodgkin Building, Babraham Research Campus, Babraham, Cambridge CB22 3FH, United Kingdom
| | - Joseph M Sheridan
- RxCelerate Limited, Dorothy Hodgkin Building, Babraham Research Campus, Babraham, Cambridge CB22 3FH, United Kingdom
| | - Manisha Kulkarni
- ChemOvation Limited, 49 Arrivato Plaza, Hall Street, Saint Helens WA10 1GH, United Kingdom
| | - Neill Liptrott
- Pharmacology Research Laboratories, University of Liverpool, Block H, First Floor, 70 Pembroke Place, Liverpool L69 3GF, United Kingdom
| | - Deirdre Egan
- Pharmacology Research Laboratories, University of Liverpool, Block H, First Floor, 70 Pembroke Place, Liverpool L69 3GF, United Kingdom
| | - Saye Khoo
- Pharmacology Research Laboratories, University of Liverpool, Block H, First Floor, 70 Pembroke Place, Liverpool L69 3GF, United Kingdom
| | - Ian Matthews
- ChemOvation Limited, 49 Arrivato Plaza, Hall Street, Saint Helens WA10 1GH, United Kingdom
| | - William D O Hamilton
- RxCelerate Limited, Dorothy Hodgkin Building, Babraham Research Campus, Babraham, Cambridge CB22 3FH, United Kingdom
| |
Collapse
|
31
|
Borbulevych OY, Martin RI, Westerhoff LM. The critical role of QM/MM X-ray refinement and accurate tautomer/protomer determination in structure-based drug design. J Comput Aided Mol Des 2021; 35:433-451. [PMID: 33108589 PMCID: PMC8018927 DOI: 10.1007/s10822-020-00354-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Accepted: 10/12/2020] [Indexed: 12/29/2022]
Abstract
Conventional protein:ligand crystallographic refinement uses stereochemistry restraints coupled with a rudimentary energy functional to ensure the correct geometry of the model of the macromolecule-along with any bound ligand(s)-within the context of the experimental, X-ray density. These methods generally lack explicit terms for electrostatics, polarization, dispersion, hydrogen bonds, and other key interactions, and instead they use pre-determined parameters (e.g. bond lengths, angles, and torsions) to drive structural refinement. In order to address this deficiency and obtain a more complete and ultimately more accurate structure, we have developed an automated approach for macromolecular refinement based on a two layer, QM/MM (ONIOM) scheme as implemented within our DivCon Discovery Suite and "plugged in" to two mainstream crystallographic packages: PHENIX and BUSTER. This implementation is able to use one or more region layer(s), which is(are) characterized using linear-scaling, semi-empirical quantum mechanics, followed by a system layer which includes the balance of the model and which is described using a molecular mechanics functional. In this work, we applied our Phenix/DivCon refinement method-coupled with our XModeScore method for experimental tautomer/protomer state determination-to the characterization of structure sets relevant to structure-based drug design (SBDD). We then use these newly refined structures to show the impact of QM/MM X-ray refined structure on our understanding of function by exploring the influence of these improved structures on protein:ligand binding affinity prediction (and we likewise show how we use post-refinement scoring outliers to inform subsequent X-ray crystallographic efforts). Through this endeavor, we demonstrate a computational chemistry ↔ structural biology (X-ray crystallography) "feedback loop" which has utility in industrial and academic pharmaceutical research as well as other allied fields.
Collapse
Affiliation(s)
- Oleg Y Borbulevych
- QuantumBio Inc, 2790 West College Ave, Suite 900, State College, PA, 16801, USA
| | - Roger I Martin
- QuantumBio Inc, 2790 West College Ave, Suite 900, State College, PA, 16801, USA
| | - Lance M Westerhoff
- QuantumBio Inc, 2790 West College Ave, Suite 900, State College, PA, 16801, USA.
| |
Collapse
|
32
|
Gutten O, Jurečka P, Aliakbar Tehrani Z, Buděšínský M, Řezáč J, Rulíšek L. Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments. Phys Chem Chem Phys 2021; 23:7280-7294. [PMID: 33876088 DOI: 10.1039/d0cp05993e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.
Collapse
Affiliation(s)
- Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
33
|
Brueckner AC, Deng Q, Cleves AE, Lesburg CA, Alvarez JC, Reibarkh MY, Sherer EC, Jain AN. Conformational Strain of Macrocyclic Peptides in Ligand-Receptor Complexes Based on Advanced Refinement of Bound-State Conformers. J Med Chem 2021; 64:3282-3298. [PMID: 33724820 DOI: 10.1021/acs.jmedchem.0c02159] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Macrocyclic peptides are an important modality in drug discovery, but molecular design is limited due to the complexity of their conformational landscape. To better understand conformational propensities, global strain energies were estimated for 156 protein-macrocyclic peptide cocrystal structures. Unexpectedly large strain energies were observed when the bound-state conformations were modeled with positional restraints. Instead, low-energy conformer ensembles were generated using xGen that fit experimental X-ray electron density maps and gave reasonable strain energy estimates. The ensembles featured significant conformational adjustments while still fitting the electron density as well or better than the original coordinates. Strain estimates suggest the interaction energy in protein-ligand complexes can offset a greater amount of strain for macrocyclic peptides than for small molecules and non-peptidic macrocycles. Across all molecular classes, the approximate upper bound on global strain energies had the same relationship with molecular size, and bound-state ensembles from xGen yielded favorable binding energy estimates.
Collapse
Affiliation(s)
- Alexander C Brueckner
- Computational & Structural Chemistry, Merck & Co Inc, 2000 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Qiaolin Deng
- Computational & Structural Chemistry, Merck & Co Inc, 2000 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Ann E Cleves
- Bioengineering and Therapeutic Sciences, University of California San Francisco, Box 0128, San Francisco, California 94158, United States
| | - Charles A Lesburg
- Computational and Structural Chemistry, Merck and Co Inc, 33 Avenue Louis Pasteur, Boston, Massachusetts 02115, United States
| | - Juan C Alvarez
- Computational and Structural Chemistry, Merck and Co Inc, 33 Avenue Louis Pasteur, Boston, Massachusetts 02115, United States
| | - Mikhail Y Reibarkh
- Analytical Research and Development, Merck & Co Inc, 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Edward C Sherer
- Analytical Research and Development, Merck & Co Inc, 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Ajay N Jain
- Bioengineering and Therapeutic Sciences, University of California San Francisco, Box 0128, San Francisco, California 94158, United States
| |
Collapse
|
34
|
Tong J, Zhao S. Large-Scale Analysis of Bioactive Ligand Conformational Strain Energy by Ab Initio Calculation. J Chem Inf Model 2021; 61:1180-1192. [PMID: 33630603 DOI: 10.1021/acs.jcim.0c01197] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Ligand conformational strain energy (LCSE) plays an important role in virtual screening and lead optimization. While various studies have provided insights into LCSE for small-molecule ligands in the Protein Data Bank (PDB), conclusions are inconsistent mainly due to small datasets, poor quality control of crystal structures, and molecular mechanics (MM) or low-level quantum mechanics (QM) calculations. Here, we built a high-quality dataset (LigBoundConf) of 8145 ligand-bound conformations from PDB crystal structures and calculated LCSE at the M062X-D3/ma-TZVPP (SMD)//M062X-D3/def2-SVP(SMD) level for each case in the dataset. The mean/median LCSE is 4.6/3.7 kcal/mol for 6672 successfully calculated cases, which is significantly lower than the estimates based on molecular mechanics in many previous analyses. Especially, when removing ligands with nonaromatic ring(s) that are prone to have large LCSEs due to electron density overfitting, the mean/median LCSE was reduced to 3.3/2.5 kcal/mol. We further reveal that LCSE is correlated with several ligand properties, including formal atomic charge, molecular weight, number of rotatable bonds, and number of hydrogen-bond donors and acceptors. In addition, our results show that although summation of torsion strains is a good approximation of LCSE for most cases, for a small fraction (about 6%) of our dataset, it underestimates LCSEs if ligands could form nonlocal intramolecular interactions in the unbound state. Taken together, our work provides a comprehensive profile of LCSE for ligands in PDB, which could help ligand conformation generation, ligand docking pose evaluation, and lead optimization.
Collapse
Affiliation(s)
- Jiahui Tong
- iHuman Institute, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China.,School of Life Science and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China.,University of Chinese Academy of Sciences, No. 19A, Yuquan Road, Beijing 100049, China.,Shanghai Institute of Nutrition and Health, Chinese Academy of Sciences, 320 Yueyang Road, Shanghai 200031, China
| | - Suwen Zhao
- iHuman Institute, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China.,School of Life Science and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| |
Collapse
|
35
|
Abstract
Molecular descriptors encode a variety of molecular representations for computer-assisted drug discovery. Here, we focus on the Weighted Holistic Atom Localization and Entity Shape (WHALES) descriptors, which were originally designed for scaffold hopping from natural products to synthetic molecules. WHALES descriptors capture molecular shape and partial charges simultaneously. We introduce the key aspects of the WHALES concept and provide a step-by-step guide on how to use these descriptors for virtual compound screening and scaffold hopping. The results presented can be reproduced by using the code freely available from URL: github.com/ETHmodlab/scaffold_hopping_whales .
Collapse
Affiliation(s)
- Francesca Grisoni
- Department of Chemistry and Applied Biosciences, RETHINK, ETH Zurich, Zurich, Switzerland.
| | - Gisbert Schneider
- Department of Chemistry and Applied Biosciences, RETHINK, ETH Zurich, Zurich, Switzerland
| |
Collapse
|
36
|
Lewis-Atwell T, Townsend PA, Grayson MN. Comparisons of different force fields in conformational analysis and searching of organic molecules: A review. Tetrahedron 2021. [DOI: 10.1016/j.tet.2020.131865] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
37
|
Dubey A, Takeuchi K, Reibarkh M, Arthanari H. The role of NMR in leveraging dynamics and entropy in drug design. JOURNAL OF BIOMOLECULAR NMR 2020; 74:479-498. [PMID: 32720098 PMCID: PMC7686249 DOI: 10.1007/s10858-020-00335-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Accepted: 07/11/2020] [Indexed: 05/03/2023]
Abstract
Nuclear magnetic resonance (NMR) spectroscopy has contributed to structure-based drug development (SBDD) in a unique way compared to the other biophysical methods. The potency of a ligand binding to a protein is dictated by the binding free energy, which is an intricate interplay between entropy and enthalpy. In addition to providing the atomic resolution structural information, NMR can help to identify protein-ligand interactions that potentially contribute to the enthalpic component of the free energy. NMR can also illuminate dynamic aspects of the interaction, which correspond to the entropic term of the free energy. The ability of NMR to access both terms in the free energy equation stems from the suite of experiments developed to shed light on various aspects that contribute to both entropy and enthalpy, deepening our understanding of the biological function of macromolecules and assisting to target them in physiological conditions. Here we provide a brief account of the contribution of NMR to SBDD, highlighting hallmark examples and discussing the challenges that demand further method development. In the era of integrated biology, the unique ability of NMR to directly ascertain structural and dynamical aspects of macromolecule and monitor changes in these properties upon engaging a ligand can be combined with computational and other structural and biophysical methods to provide a more complete picture of the energetics of drug engagement with the target. Such efforts can be used to engineer better drugs.
Collapse
Affiliation(s)
- Abhinav Dubey
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA, 02115, USA
| | - Koh Takeuchi
- Cellular and Molecular Biotechnology Research Institute & Molecular Profiling Research Center for Drug Discovery (molprof), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, 135-0064, Japan.
| | - Mikhail Reibarkh
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ, 07065, USA
| | - Haribabu Arthanari
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA.
- Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA, 02115, USA.
| |
Collapse
|
38
|
A Potential Inhibition Process of Ricin Protein with the flavonoids Quercetin and Epigallocatechin Gallate. A Quantum-Chemical and Molecular Docking Study. Processes (Basel) 2020. [DOI: 10.3390/pr8111393] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Castor bean (Ricinus Communis) oil has been reported as one of the most important bio-based fuels; however, high amounts of toxic solid residue are generated in the production. This toxicity is due to several molecules, ricin protein being the most studied compound. The inhibition of the ricin protein is essential for eliminating its toxicity. The objective of this study is to predict the possible inhibition process via the interactions between the ricin protein and the flavonoids quercetin (Q) and epigallocatechin gallate (EGCG). The molecular structures of the complexes formed between the ricin protein and flavonoids were studied using quantum-chemical and molecular docking calculations to analyze the type of interaction, active site of the protein, binding energies, and different conformations in the inhibition process. Different methodologies were applied for the molecular structure determination; the best approximation was obtained with B3LYP/6-31G (d,p) theoretical methodology. Mappings of electrostatic potential (MEP) and frontier molecular orbitals were used for the identification of the probable sites of interaction, which were confirmed by molecular docking. The adjustment and alignment of flavonoid groups before and after the interaction, and charge transfer parameters, showed that Q and EGCG act as electron donors inside of the active site in ricin.
Collapse
|
39
|
Zheng Z, Borbulevych OY, Liu H, Deng J, Martin RI, Westerhoff LM. MovableType Software for Fast Free Energy-Based Virtual Screening: Protocol Development, Deployment, Validation, and Assessment. J Chem Inf Model 2020; 60:5437-5456. [PMID: 32791826 PMCID: PMC7781189 DOI: 10.1021/acs.jcim.0c00618] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
![]()
For decades, the
complicated energy surfaces found in macromolecular
protein:ligand structures, which require large amounts of computational
time and resources for energy state sampling, have been an inherent
obstacle to fast, routine free energy estimation in industrial drug
discovery efforts. Beginning in 2013, the Merz research group addressed
this cost with the introduction of a novel sampling methodology termed
“Movable Type” (MT). Using numerical integration methods,
the MT method reduces the computational expense for energy state sampling
by independently calculating each atomic partition function from an
initial molecular conformation in order to estimate the molecular
free energy using ensembles of the atomic partition functions. In
this work, we report a software package, the DivCon Discovery Suite
with the MovableType module from QuantumBio Inc., that performs this
MT free energy estimation protocol in a fast, fully encapsulated manner.
We discuss the computational procedures and improvements to the original
work, and we detail the corresponding settings for this software package.
Finally, we introduce two validation benchmarks to evaluate the overall
robustness of the method against a broad range of protein:ligand structural
cases. With these publicly available benchmarks, we show that the
method can use a variety of input types and parameters and exhibits
comparable predictability whether the method is presented with “expensive”
X-ray structures or “inexpensively docked” theoretical
models. We also explore some next steps for the method. The MovableType
software is available at http://www.quantumbioinc.com/
Collapse
Affiliation(s)
- Zheng Zheng
- QuantumBio Inc., 2790 West College Avenue, Suite 900, State College, Pennsylvania 16801, United States.,School of Chemistry, Chemical Engineering and Life Science, Wuhan University of Technology, 122 Luoshi Road, Wuhan 430070, P. R. China
| | - Oleg Y Borbulevych
- QuantumBio Inc., 2790 West College Avenue, Suite 900, State College, Pennsylvania 16801, United States
| | - Hao Liu
- School of Mechanical and Electronic Engineering, Wuhan University of Technology, 122 Luoshi Road, Wuhan 430070, P. R. China
| | - Jianpeng Deng
- School of Chemistry, Chemical Engineering and Life Science, Wuhan University of Technology, 122 Luoshi Road, Wuhan 430070, P. R. China
| | - Roger I Martin
- QuantumBio Inc., 2790 West College Avenue, Suite 900, State College, Pennsylvania 16801, United States
| | - Lance M Westerhoff
- QuantumBio Inc., 2790 West College Avenue, Suite 900, State College, Pennsylvania 16801, United States
| |
Collapse
|
40
|
Jain AN, Cleves AE, Brueckner AC, Lesburg CA, Deng Q, Sherer EC, Reibarkh MY. XGen: Real-Space Fitting of Complex Ligand Conformational Ensembles to X-ray Electron Density Maps. J Med Chem 2020; 63:10509-10528. [DOI: 10.1021/acs.jmedchem.0c01373] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ajay N. Jain
- Bioengineering and Therapeutic Sciences, University of California, San Francisco, California 94143 United States
| | - Ann E. Cleves
- BioPharmics LLC, Santa Rosa, California 95404 United States
| | | | | | - Qiaolin Deng
- Merck and Co., Inc., Kenilworth, New Jersey 07033 United States
| | | | | |
Collapse
|
41
|
Udvarhelyi A, Rodde S, Wilcken R. ReSCoSS: a flexible quantum chemistry workflow identifying relevant solution conformers of drug-like molecules. J Comput Aided Mol Des 2020; 35:399-415. [PMID: 32803515 DOI: 10.1007/s10822-020-00337-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 08/12/2020] [Indexed: 12/20/2022]
Abstract
Conformational equilibria are at the heart of drug design, yet their energetic description is often hampered by the insufficient accuracy of low-cost methods. Here we present a flexible and semi-automatic workflow based on quantum chemistry, ReSCoSS, designed to identify relevant conformers and predict their equilibria across different solvent environments in the Conductor-like Screening Model for Real Solvents (COSMO-RS) framework. We demonstrate the utility and accuracy of the workflow through conformational case studies on several drug-like molecules from literature where relevant conformations are known. We further show that including ReSCoSS conformers significantly improves COSMO-RS based predictions of physicochemical properties over single-conformation approaches. ReSCoSS has found broad adoption in the in-house drug discovery and development work streams and has contributed to establishing quantum-chemistry methods as a strategic pillar in ligand discovery.
Collapse
Affiliation(s)
- Anikó Udvarhelyi
- Technical Research and Development, Novartis Pharma AG, 4002, Basel, Switzerland.
| | - Stephane Rodde
- Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 4002, Basel, Switzerland
| | - Rainer Wilcken
- Global Discovery Chemistry, Novartis Institutes for BioMedical Research, 4002, Basel, Switzerland.
| |
Collapse
|
42
|
Zivanovic S, Colizzi F, Moreno D, Hospital A, Soliva R, Orozco M. Exploring the Conformational Landscape of Bioactive Small Molecules. J Chem Theory Comput 2020; 16:6575-6585. [DOI: 10.1021/acs.jctc.0c00304] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Sanja Zivanovic
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac, 10, 08028 Barcelona, Spain
| | - Francesco Colizzi
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac, 10, 08028 Barcelona, Spain
| | - David Moreno
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac, 10, 08028 Barcelona, Spain
| | - Adam Hospital
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac, 10, 08028 Barcelona, Spain
| | - Robert Soliva
- Nostrum Biodiscovery, Nexus II Building, 08034 Barcelona, Spain
| | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology (BIST), Baldiri Reixac, 10, 08028 Barcelona, Spain
- Departament de Bioquímica i Biomedicina, Facultat de Biologia, Universitat de Barcelona, E08028 Barcelona, Spain
| |
Collapse
|
43
|
Mandal M, Buevich A, Caldwell JP, Hyde L, Huang X, Liu X, McKittrick B, Mazzola RD, Pissarnitski D, Palani A, Zhang L, Parker E, Xiao L, Rindgen D, Zhu Z. Generation of Leads for γ-Secretase Modulation. J Med Chem 2020; 63:8216-8230. [DOI: 10.1021/acs.jmedchem.0c00446] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mihirbaran Mandal
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Alexei Buevich
- Department of NMR Structure Elucidation, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - John P. Caldwell
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Lynn Hyde
- Department of Neuroscience, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Xianhai Huang
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Xiaoxiang Liu
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Brian McKittrick
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Robert D. Mazzola
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Dmitri Pissarnitski
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Anandan Palani
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Lili Zhang
- Department of Neuroscience, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Eric Parker
- Department of Neuroscience, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Li Xiao
- Department of Modeling and Informatics, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Diane Rindgen
- Department of Pharmacokinetics, Pharmacodynamics, and Drug Metabolism, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Zhaoning Zhu
- Department of Medicinal Chemistry, Merck & Co., Inc., 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| |
Collapse
|
44
|
Senior T, Botha MJ, Kennedy AR, Calvo-Castro J. Understanding the Contribution of Individual Amino Acid Residues in the Binding of Psychoactive Substances to Monoamine Transporters. ACS OMEGA 2020; 5:17223-17231. [PMID: 32715208 PMCID: PMC7376891 DOI: 10.1021/acsomega.0c01370] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2020] [Accepted: 06/19/2020] [Indexed: 05/05/2023]
Abstract
The development of point-of-care detection methodologies for biologically relevant analytes that can facilitate rapid and appropriate treatment is at the forefront of current research efforts and interests. Among the various approaches, those exploiting host-guest chemistries where the optoelectronic signals of the chemical sensor can be modulated upon interaction with the target analyte are of particular interest. In aiding their rational development, judicious selection of peripheral functional groups anchored to core motifs with desired properties is critical. Herein, we report an in-depth investigation of the binding of three psychoactive substances, MDAI, mexedrone, and phenibut, to receptors of the monoamine transporters for dopamine, norepinephrine, and serotonin, particularly focusing on the role of individual amino acid residues. We first evaluated the conformational flexibility of the ligands by comparing their experimentally determined crystal structure geometries to those optimized by means of quantum as well as molecular mechanics, observing significant changes in the case of phenibut. Molecular docking studies were employed to identify preferential binding sites by means of calculated docking scores. In all cases, irrespective of the monoamine transporter, psychoactive substances exhibited preferred interaction with the S1 or central site of the proteins, in line with previous studies. However, we observed that experimental trends for their relative potency on the three transporters were only reproduced in the case of mexedrone. Subsequently, to further understand these findings and to pave the way for the rational development of superior chemical sensors for these substances, we computed the individual contributions of each nearest neighbor amino acid residue to the binding to the target analytes. Interestingly, these results are now in agreement with those experimental potency trends. In addition, these observations were in all cases associated with key intermolecular interactions with neighboring residues, such as tyrosine and aspartic acid, in the binding of the ligands to the monoamine transporter for dopamine. As a result, we believe this work will be of interest to those engaged in the rational development of chemical sensors for small molecule analytes as well as to those interested in the use of computational approaches to further understand protein-ligand interactions.
Collapse
Affiliation(s)
- Tamara Senior
- School
of Life and Medical Sciences, University
of Hertfordshire, Hatfield AL10 9AB, U.K.
| | - Michelle J. Botha
- School
of Life and Medical Sciences, University
of Hertfordshire, Hatfield AL10 9AB, U.K.
| | - Alan R. Kennedy
- Department
of Pure & Applied Chemistry, University
of Strathclyde, Glasgow G1 1XL, U.K.
| | - Jesus Calvo-Castro
- School
of Life and Medical Sciences, University
of Hertfordshire, Hatfield AL10 9AB, U.K.
| |
Collapse
|
45
|
Benchmarking Data Sets from PubChem BioAssay Data: Current Scenario and Room for Improvement. Int J Mol Sci 2020; 21:ijms21124380. [PMID: 32575564 PMCID: PMC7352161 DOI: 10.3390/ijms21124380] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 06/15/2020] [Accepted: 06/18/2020] [Indexed: 11/17/2022] Open
Abstract
Developing realistic data sets for evaluating virtual screening methods is a task that has been tackled by the cheminformatics community for many years. Numerous artificially constructed data collections were developed, such as DUD, DUD-E, or DEKOIS. However, they all suffer from multiple drawbacks, one of which is the absence of experimental results confirming the impotence of presumably inactive molecules, leading to possible false negatives in the ligand sets. In light of this problem, the PubChem BioAssay database, an open-access repository providing the bioactivity information of compounds that were already tested on a biological target, is now a recommended source for data set construction. Nevertheless, there exist several issues with the use of such data that need to be properly addressed. In this article, an overview of benchmarking data collections built upon experimental PubChem BioAssay input is provided, along with a thorough discussion of noteworthy issues that one must consider during the design of new ligand sets from this database. The points raised in this review are expected to guide future developments in this regard, in hopes of offering better evaluation tools for novel in silico screening procedures.
Collapse
|
46
|
Prosser K, Stokes RW, Cohen SM. Evaluation of 3-Dimensionality in Approved and Experimental Drug Space. ACS Med Chem Lett 2020; 11:1292-1298. [PMID: 32551014 PMCID: PMC7294711 DOI: 10.1021/acsmedchemlett.0c00121] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2020] [Accepted: 05/18/2020] [Indexed: 11/28/2022] Open
Abstract
The 3-dimensional (3D) structure of therapeutics and other bioactive molecules is an important factor in determining the strength and selectivity of their protein-ligand interactions. Previous efforts have considered the strain introduced and tolerated through conformational changes induced upon protein binding. Herein, we present an analysis of 3-dimentionality for energy-minimized structures from the DrugBank and ligands bound to proteins identified in the Protein Data Bank (PDB). This analysis reveals that the majority of molecules found in both the DrugBank and the PDB tend toward linearity and planarity, with few molecules having highly 3D conformations. Decidedly 3D geometries have been historically difficult to achieve, likely due to the synthetic challenge of making 3D organic molecules, and other considerations, such as adherence to the 'rule-of-five'. This has resulted in the dominance of planar and/or linear topologies of the molecules described here. Strategies to address the generally flat nature of these data sets are explored, including the use of 3D organic fragments and inorganic scaffolds as a means of accessing privileged 3D space. This work highlights the potential utility of libraries with greater 3D topological diversity so that the importance of molecular shape to biological behavior can be more fully understood in drug discovery campaigns.
Collapse
Affiliation(s)
- Kathleen
E. Prosser
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Ryjul W. Stokes
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Seth M. Cohen
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
47
|
Link BA, Sindt AJ, Shimizu LS, Do TD. Selective host-guest chemistry, self-assembly and conformational preferences of m-xylene macrocycles probed by ion-mobility spectrometry mass spectrometry. Phys Chem Chem Phys 2020; 22:9290-9300. [PMID: 32309846 DOI: 10.1039/c9cp06938k] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We demonstrated ion-mobility spectrometry mass spectrometry (IMS-MS) as a powerful tool for interrogating and preserving selective chemistry including non-covalent and host-guest complexes of m-xylene macrocycles formed in solution. The technique readily revealed the unique favorability of a thiourea-containing macrocycle MXT to Zn2+ to form a dimer complex with the cation in an off-axis sandwich structure having the Zn-S bonds in a tetrahedral coordination environment. Replacing thiourea with urea generates MXU which formed high-order oligomerization with weak binding interactions to neutral DMSO guests detected at every oligomer size. The self-assembly pathway observed for this macrocycle is consistent with the crystalline assembly. Further transformation of urea into squaramide produces MXS, a rare receptor for probing sulfate in solution. Tight complexes were observed for both monomeric and dimeric of MXS in which HSO4- bound stronger than SO42- to the host. The position of HSO4- at the binding cavity is a 180° inversion of the reported crystallographic SO42-. The MXS dimer formed a prism-like shape with HSO4- exhibiting strong contacts with the 8 amine protons of two MXS macrocycles. By eliminating intermolecular interferences, we detected the low energy structures of MXS with collisional cross section (CCS) matching cis-trans and cis-cis squaramides-amines, both were not observed in crystallization trials. The experiments collectively unravel multiple facets of macrocycle chemistry including conformational flexibility, self-assembly and ligand binding; all in one analysis. Our findings illustrate an inexpensive and widely applicable approach to investigate weak but important interactions that define the shape and binding of macrocycles.
Collapse
Affiliation(s)
- Benjamin A Link
- Department of Chemistry, University of Tennessee, Knoxville, TN 37996, USA.
| | | | | | | |
Collapse
|
48
|
Szűcs E, Marton J, Szabó Z, Hosztafi S, Kékesi G, Tuboly G, Bánki L, Horváth G, Szabó PT, Tömböly C, Varga ZK, Benyhe S, Ötvös F. Synthesis, biochemical, pharmacological characterization and in silico profile modelling of highly potent opioid orvinol and thevinol derivatives. Eur J Med Chem 2020; 191:112145. [PMID: 32092588 DOI: 10.1016/j.ejmech.2020.112145] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 01/22/2020] [Accepted: 02/12/2020] [Indexed: 02/02/2023]
Abstract
Morphine and its derivatives play inevitably important role in the μ-opioid receptor (MOR) targeted antinociception. A structure-activity relationship study is presented for novel and known orvinol and thevinol derivatives with varying 3-O, 6-O, 17-N and 20-alkyl substitutions starting from agonists, antagonists and partial agonists. In vitro competition binding experiments with [3H]DAMGO showed low subnanomolar affinity to MOR. Generally, 6-O-demethylation increased the affinity toward MOR and decreased the efficacy changing the pharmacological profile in some cases. In vivo tests in osteoarthritis inflammation model showed significant antiallodynic effects of thevinol derivatives while orvinol derivatives did not. The pharmacological character was modelled by computational docking to both active and inactive state models of MOR. Docking energy difference for the two states separates agonists and antagonists well while partial agonists overlapped with them. An interaction pattern of the ligands, involving the interacting receptor atoms, showed more efficient separation of the pharmacological profiles. In rats, thevinol derivatives showed antiallodynic effect in vivo. The orvinol derivatives, except for 6-O-desmethyl-dihydroetorfin (2c), did not show antiallodynic effect.
Collapse
Affiliation(s)
- Edina Szűcs
- Institute of Biochemistry, Biological Research Center, Temesvári krt. 62, H-6726, Szeged, Hungary; Doctoral School of Theoretical Medicine, Faculty of Medicine, University of Szeged, Dóm tér 10, H-6720, Szeged, Hungary
| | - János Marton
- ABX Advanced Biochemical Compounds, Biomedizinische Forschungsreagenzien GmbH, Heinrich-Glaeser-Strasse 10-14, D-01454, Radeberg, Germany
| | - Zoltán Szabó
- Royal Institute of Technology (KTH), School of Engineering Sciences in Chemistry, Biotechnology and Health, Department of Chemistry, Organic Chemistry, S-100 44, Stockholm, Sweden
| | - Sándor Hosztafi
- Institute of Pharmaceutical Chemistry, Semmelweis Medical University, Hőgyes Endre utca 9, H-1092, Budapest, Hungary
| | - Gabriella Kékesi
- Department of Physiology, Faculty of Medicine, University of Szeged, Dóm tér 10, H-6720, Szeged, Hungary
| | - Gábor Tuboly
- Department of Neurology, Faculty of Medicine, University of Szeged, Semmelweis u 6, H-6725, Szeged, Hungary
| | - László Bánki
- Department of Traumatology, Faculty of Medicine, University of Szeged, Semmelweis u 6, H-6725, Szeged, Hungary
| | - Gyöngyi Horváth
- Department of Physiology, Faculty of Medicine, University of Szeged, Dóm tér 10, H-6720, Szeged, Hungary
| | - Pál T Szabó
- Research Centre for Natural Sciences, MS Metabolomics Research Laboratory, H-1117, Budapest, Magyar tudósok krt. 2, Hungary
| | - Csaba Tömböly
- Institute of Biochemistry, Biological Research Center, Temesvári krt. 62, H-6726, Szeged, Hungary
| | - Zsuzsanna Katalin Varga
- Institute of Biochemistry, Biological Research Center, Temesvári krt. 62, H-6726, Szeged, Hungary; Doctoral School of Theoretical Medicine, Faculty of Medicine, University of Szeged, Dóm tér 10, H-6720, Szeged, Hungary
| | - Sándor Benyhe
- Institute of Biochemistry, Biological Research Center, Temesvári krt. 62, H-6726, Szeged, Hungary
| | - Ferenc Ötvös
- Institute of Biochemistry, Biological Research Center, Temesvári krt. 62, H-6726, Szeged, Hungary.
| |
Collapse
|
49
|
Macalino SJY, Billones JB, Organo VG, Carrillo MCO. In Silico Strategies in Tuberculosis Drug Discovery. Molecules 2020; 25:E665. [PMID: 32033144 PMCID: PMC7037728 DOI: 10.3390/molecules25030665] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 12/15/2019] [Accepted: 12/17/2019] [Indexed: 12/16/2022] Open
Abstract
Tuberculosis (TB) remains a serious threat to global public health, responsible for an estimated 1.5 million mortalities in 2018. While there are available therapeutics for this infection, slow-acting drugs, poor patient compliance, drug toxicity, and drug resistance require the discovery of novel TB drugs. Discovering new and more potent antibiotics that target novel TB protein targets is an attractive strategy towards controlling the global TB epidemic. In silico strategies can be applied at multiple stages of the drug discovery paradigm to expedite the identification of novel anti-TB therapeutics. In this paper, we discuss the current TB treatment, emergence of drug resistance, and the effective application of computational tools to the different stages of TB drug discovery when combined with traditional biochemical methods. We will also highlight the strengths and points of improvement in in silico TB drug discovery research, as well as possible future perspectives in this field.
Collapse
Affiliation(s)
- Stephani Joy Y. Macalino
- Chemistry Department, De La Salle University, 2401 Taft Avenue, Manila 0992, Philippines;
- OVPAA-EIDR Program, “Computer-Aided Discovery of Compounds for the Treatment of Tuberculosis in the Philippines”, Department of Physical Sciences and Mathematics, College of Arts and Sciences, University of the Philippines Manila, Manila 1000, Philippines; (V.G.O.); (M.C.O.C.)
| | - Junie B. Billones
- OVPAA-EIDR Program, “Computer-Aided Discovery of Compounds for the Treatment of Tuberculosis in the Philippines”, Department of Physical Sciences and Mathematics, College of Arts and Sciences, University of the Philippines Manila, Manila 1000, Philippines; (V.G.O.); (M.C.O.C.)
| | - Voltaire G. Organo
- OVPAA-EIDR Program, “Computer-Aided Discovery of Compounds for the Treatment of Tuberculosis in the Philippines”, Department of Physical Sciences and Mathematics, College of Arts and Sciences, University of the Philippines Manila, Manila 1000, Philippines; (V.G.O.); (M.C.O.C.)
| | - Maria Constancia O. Carrillo
- OVPAA-EIDR Program, “Computer-Aided Discovery of Compounds for the Treatment of Tuberculosis in the Philippines”, Department of Physical Sciences and Mathematics, College of Arts and Sciences, University of the Philippines Manila, Manila 1000, Philippines; (V.G.O.); (M.C.O.C.)
| |
Collapse
|
50
|
Conformational analysis of macrocycles: comparing general and specialized methods. J Comput Aided Mol Des 2020; 34:231-252. [PMID: 31965404 PMCID: PMC7036058 DOI: 10.1007/s10822-020-00277-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2019] [Accepted: 01/03/2020] [Indexed: 11/24/2022]
Abstract
Abstract Macrocycles represent an important class of medicinally relevant small molecules due to their interesting biological properties. Therefore, a firm understanding of their conformational preferences is important for drug design. Given the importance of macrocycle-protein modelling in drug discovery, we envisaged that a systematic study of both classical and recent specialized methods would provide guidance for other practitioners within the field. In this study we compare the performance of the general, well established conformational analysis methods Monte Carlo Multiple Minimum (MCMM) and Mixed Torsional/Low-Mode sampling (MTLMOD) with two more recent and specialized macrocycle sampling techniques: MacroModel macrocycle Baseline Search (MD/LLMOD) and Prime macrocycle conformational sampling (PRIME-MCS). Using macrocycles extracted from 44 macrocycle-protein X-ray crystallography complexes, we evaluated each method based on their ability to (i) generate unique conformers, (ii) generate unique macrocycle ring conformations, (iii) identify the global energy minimum, (iv) identify conformers similar to the X-ray ligand conformation after Protein Preparation Wizard treatment (X-rayppw), and (v) to the X-rayppw ring conformation. Computational speed was also considered. In addition, conformational coverage, as defined by the number of conformations identified, was studied. In order to study the relative energies of the bioactive conformations, the energy differences between the global energy minima and the energy minimized X-rayppw structures and, the global energy minima and the MCMM-Exhaustive (1,000,000 search steps) generated conformers closest to the X-rayppw structure, were calculated and analysed. All searches were performed using relatively short run times (10,000 steps for MCMM, MTLMOD and MD/LLMOD). To assess the performance of the methods, they were compared to an exhaustive MCMM search using 1,000,000 search steps for each of the 44 macrocycles (requiring ca 200 times more CPU time). Prior to our analysis, we also investigated if the general search methods MCMM and MTLMOD could also be optimized for macrocycle conformational sampling. Taken together, our work concludes that the more general methods can be optimized for macrocycle modelling by slightly adjusting the settings around the ring closure bond. In most cases, MCMM and MTLMOD with either standard or enhanced settings performed well in comparison to the more specialized macrocycle sampling methods MD/LLMOD and PRIME-MCS. When using enhanced settings for MCMM and MTLMOD, the X-rayppw conformation was regenerated with the greatest accuracy. The, MD/LLMOD emerged as the most efficient method for generating the global energy minima. Graphic abstract ![]()
Electronic supplementary material The online version of this article (10.1007/s10822-020-00277-2) contains supplementary material, which is available to authorized users.
Collapse
|