1
|
Fedotova MV, Chuev GN. The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands. Int J Mol Sci 2024; 25:3697. [PMID: 38612508 PMCID: PMC11011341 DOI: 10.3390/ijms25073697] [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: 03/09/2024] [Revised: 03/22/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.
Collapse
Affiliation(s)
- Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., 142290 Pushchino, Russia
| |
Collapse
|
2
|
Kriegel M, Muller YA. De novo prediction of explicit water molecule positions by a novel algorithm within the protein design software MUMBO. Sci Rep 2023; 13:16680. [PMID: 37794104 PMCID: PMC10550942 DOI: 10.1038/s41598-023-43659-w] [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: 08/01/2023] [Accepted: 09/26/2023] [Indexed: 10/06/2023] Open
Abstract
By mediating interatomic interactions, water molecules play a major role in protein-protein, protein-DNA and protein-ligand interfaces, significantly affecting affinity and specificity. This notwithstanding, explicit water molecules are usually not considered in protein design software because of high computational costs. To challenge this situation, we analyzed the binding characteristics of 60,000 waters from high resolution crystal structures and used the observed parameters to implement the prediction of water molecules in the protein design and side chain-packing software MUMBO. To reduce the complexity of the problem, we incorporated water molecules through the solvation of rotamer pairs instead of relying on solvated rotamer libraries. Our validation demonstrates the potential of our algorithm by achieving recovery rates of 67% for bridging water molecules and up to 86% for fully coordinated waters. The efficacy of our algorithm is highlighted further by the prediction of 3 different proteinligand complexes. Here, 91% of water-mediated interactions between protein and ligand are correctly predicted. These results suggest that the new algorithm could prove highly beneficial for structure-based protein design, particularly for the optimization of ligand-binding pockets or protein-protein interfaces.
Collapse
Affiliation(s)
- Mark Kriegel
- Division of Biotechnology, Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Yves A Muller
- Division of Biotechnology, Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany.
| |
Collapse
|
3
|
Yoon HR, Park GJ, Balupuri A, Kang NS. TWN-FS method: A novel fragment screening method for drug discovery. Comput Struct Biotechnol J 2023; 21:4683-4696. [PMID: 37841326 PMCID: PMC10568351 DOI: 10.1016/j.csbj.2023.09.037] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 09/26/2023] [Accepted: 09/27/2023] [Indexed: 10/17/2023] Open
Abstract
Fragment-based drug discovery (FBDD) is a well-established and effective method for generating diverse and novel hits in drug design. Kinases are suitable targets for FBDD due to their well-defined structure. Water molecules contribute to structure and function of proteins and also influence the environment within the binding pocket. Water molecules form a variety of hydrogen-bonded cyclic water-ring networks, collectively known as topological water networks (TWNs). Analyzing the TWNs in protein binding sites can provide valuable insights into potential locations and shapes for fragments within the binding site. Here, we introduce TWN-based fragment screening (TWN-FS) method, a novel screening method that suggests fragments through grouped TWN analysis within the protein binding site. We used this method to screen known CDK2, CHK1, IGF1R and ERBB4 inhibitors. Our findings suggest that TWN-FS method has the potential to effectively screen fragments. The TWN-FS method package is available on GitHub at https://github.com/pkj0421/TWN-FS.
Collapse
Affiliation(s)
- Hye Ree Yoon
- Graduate School of New Drug Discovery and Development, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, South Korea
| | - Gyoung Jin Park
- Graduate School of New Drug Discovery and Development, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, South Korea
| | - Anand Balupuri
- Graduate School of New Drug Discovery and Development, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, South Korea
| | - Nam Sook Kang
- Graduate School of New Drug Discovery and Development, Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon 34134, South Korea
| |
Collapse
|
4
|
Bryan DR, Kulp JL, Mahapatra MK, Bryan RL, Viswanathan U, Carlisle MN, Kim S, Schutte WD, Clarke KV, Doan TT, Kulp JL. BMaps: A Web Application for Fragment-Based Drug Design and Compound Binding Evaluation. J Chem Inf Model 2023; 63:4229-4236. [PMID: 37406353 DOI: 10.1021/acs.jcim.3c00209] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2023]
Abstract
Fragment-based drug design uses data about where, and how strongly, small chemical fragments bind to proteins, to assemble new drug molecules. Over the past decade, we have been successfully using fragment data, derived from thermodynamically rigorous Monte Carlo fragment-protein binding simulations, in dozens of preclinical drug programs. However, this approach has not been available to the broader research community because of the cost and complexity of doing simulations and using design tools. We have developed a web application, called BMaps, to make fragment-based drug design widely available with greatly simplified user interfaces. BMaps provides access to a large repository (>550) of proteins with 100s of precomputed fragment maps, druggable hot spots, and high-quality water maps. Users can also employ their own structures or those from the Protein Data Bank and AlphaFold DB. Multigigabyte data sets are searched to find fragments in bondable orientations, ranked by a binding-free energy metric. The designers use this to select modifications that improve affinity and other properties. BMaps is unique in combining conventional tools such as docking and energy minimization with fragment-based design, in a very easy to use and automated web application. The service is available at https://www.boltzmannmaps.com.
Collapse
Affiliation(s)
- Daniel R Bryan
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - John L Kulp
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
- Zymergen, Inc., 430 E. 29th Street, Suite 625, New York, New York 10016, United States
| | - Manoj K Mahapatra
- Kanak Manjari Institute of Pharmaceutical Sciences, Rourkela 769015, Odisha, India
| | - Richard L Bryan
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - Usha Viswanathan
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - Micah N Carlisle
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - Surim Kim
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
- Zymergen, Inc., 430 E. 29th Street, Suite 625, New York, New York 10016, United States
| | - William D Schutte
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - Kevaughn V Clarke
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - Tony T Doan
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| | - John L Kulp
- Conifer Point Pharmaceuticals, 3805 Old Easton Road, Doylestown, Pennsylvania 18902, United States
| |
Collapse
|
5
|
Zsidó BZ, Bayarsaikhan B, Börzsei R, Szél V, Mohos V, Hetényi C. The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering. Int J Mol Sci 2023; 24:11784. [PMID: 37511543 PMCID: PMC10381018 DOI: 10.3390/ijms241411784] [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: 06/20/2023] [Revised: 07/18/2023] [Accepted: 07/20/2023] [Indexed: 07/30/2023] Open
Abstract
Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.
Collapse
Affiliation(s)
- Balázs Zoltán Zsidó
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Bayartsetseg Bayarsaikhan
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Rita Börzsei
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Viktor Szél
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Violetta Mohos
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Csaba Hetényi
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| |
Collapse
|
6
|
Eberhardt J, Forli S. WaterKit: Thermodynamic Profiling of Protein Hydration Sites. J Chem Theory Comput 2023; 19:2535-2556. [PMID: 37094087 PMCID: PMC10732097 DOI: 10.1021/acs.jctc.2c01087] [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] [Indexed: 04/26/2023]
Abstract
Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.
Collapse
Affiliation(s)
- Jerome Eberhardt
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| |
Collapse
|
7
|
Leini R, Pantsar T. In Silico Evaluation of the Thr58-Associated Conserved Water with KRAS Switch-II Pocket Binders. J Chem Inf Model 2023; 63:1490-1505. [PMID: 36854010 PMCID: PMC10015465 DOI: 10.1021/acs.jcim.2c01479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
The KRAS switch-II pocket (SII-P) has proven to be one of the most successful tools for targeting KRAS with small molecules to date. This has been demonstrated with several KRAS(G12C)-targeting covalent inhibitors, already resulting in two FDA-approved drugs. Several earlier-stage compounds have also been reported to engage KRAS SII-P with other position 12 mutants, including G12D, G12S, and G12R. A highly conserved water molecule exists in the KRAS SII-P, linking Thr58 of switch-II and Gly10 of β1 sheet. This conserved water is also present in the cocrystal structures of most of the disclosed small-molecule inhibitors but is only displaced by a handful of SII-P binders. Here, we evaluated the conserved water molecule energetics by the WaterMap for the SII-P binders with publicly disclosed structures and studied the water behavior in the presence of selected inhibitors by microsecond timescale molecular dynamics (MD) simulations using two water models (total simulation time of 120 μs). Our data revealed the high-energy nature of this hydration site when coexisting with an SII-P binder and that there is a preference for a single isolated hydration site in this location within the most advanced compounds. Furthermore, water displacement was only achieved with a few disclosed compounds and was suboptimal, as for instance a cyanomethyl group as a water displacer appears to introduce repulsion with the native conformation of Thr58. These results suggested that this conserved water should be considered more central when designing new inhibitors, especially in the design of noncovalent inhibitors targeting the SII-P.
Collapse
Affiliation(s)
- Renne Leini
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, Yliopistonranta 1C, 70210 Kuopio, Finland
| | - Tatu Pantsar
- School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, Yliopistonranta 1C, 70210 Kuopio, Finland
| |
Collapse
|
8
|
Puch-Giner I, Molina A, Municoy M, Pérez C, Guallar V. Recent PELE Developments and Applications in Drug Discovery Campaigns. Int J Mol Sci 2022; 23:ijms232416090. [PMID: 36555731 PMCID: PMC9788188 DOI: 10.3390/ijms232416090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2022] [Revised: 12/12/2022] [Accepted: 12/13/2022] [Indexed: 12/23/2022] Open
Abstract
Computer simulation techniques are gaining a central role in molecular pharmacology. Due to several factors, including the significant improvements of traditional molecular modelling, the irruption of machine learning methods, the massive data generation, or the unlimited computational resources through cloud computing, the future of pharmacology seems to go hand in hand with in silico predictions. In this review, we summarize our recent efforts in such a direction, centered on the unconventional Monte Carlo PELE software and on its coupling with machine learning techniques. We also provide new data on combining two recent new techniques, aquaPELE capable of exhaustive water sampling and fragPELE, for fragment growing.
Collapse
Affiliation(s)
- Ignasi Puch-Giner
- Barcelona Supercomputing Center, Plaça d’Eusebi Güell, 1-3, 08034 Barcelona, Spain
| | - Alexis Molina
- Nostrum Biodiscovery S.L., Av. de Josep Tarradellas, 8-10, 3-2, 08029 Barcelona, Spain
| | - Martí Municoy
- Barcelona Supercomputing Center, Plaça d’Eusebi Güell, 1-3, 08034 Barcelona, Spain
- Nostrum Biodiscovery S.L., Av. de Josep Tarradellas, 8-10, 3-2, 08029 Barcelona, Spain
| | - Carles Pérez
- Nostrum Biodiscovery S.L., Av. de Josep Tarradellas, 8-10, 3-2, 08029 Barcelona, Spain
| | - Victor Guallar
- Barcelona Supercomputing Center, Plaça d’Eusebi Güell, 1-3, 08034 Barcelona, Spain
- Nostrum Biodiscovery S.L., Av. de Josep Tarradellas, 8-10, 3-2, 08029 Barcelona, Spain
- Correspondence:
| |
Collapse
|
9
|
Pierri M, Gazzillo E, Chini MG, Ferraro MG, Piccolo M, Maione F, Irace C, Bifulco G, Bruno I, Terracciano S, Lauro G. Introducing structure-based three-dimensional pharmacophore models for accelerating the discovery of selective BRD9 binders. Bioorg Chem 2021; 118:105480. [PMID: 34823196 DOI: 10.1016/j.bioorg.2021.105480] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 11/03/2021] [Accepted: 11/07/2021] [Indexed: 01/18/2023]
Abstract
A well-structured in silico workflow is here reported for disclosing structure-based pharmacophore models against bromodomain-containing protein 9 (BRD9), accelerating virtual screening campaigns and facilitating the identification of novel binders. Specifically, starting from 23 known ligands co-crystallized with BRD9, three-dimensional pharmacophore models, namely placed in a reference protein structure, were developed. Specifically, we here introduce a fragment-related pharmacophore model, useful for the identification of new promising small chemical probes targeting the protein region responsible of the acetyllysine recognition, and two further pharmacophore models useful for the selection of compounds featuring drug-like properties. A pharmacophore-driven virtual screening campaign was then performed to facilitate the selection of new selective BRD9 ligands, starting from a large library of commercially available molecules. The identification of a promising BRD9 binder (7) prompted us to re-iterate this computational workflow on a second focused in-house built library of synthesizable compounds and, eventually, three further novel BRD9 binders were disclosed (8-10). Moreover, all these compounds were tested among a panel comprising other nine bromodomains, showing a high selectivity for BRD9. Preclinical bioscreens for potential anticancer activity highlighted compound 7 as that showing the most promising biological effects, proving the reliability of this in silico pipeline and confirming the applicability of the here introduced structure-based three-dimensional (3D) pharmacophore models as straightforward tools for the selection of new BRD9 ligands.
Collapse
Affiliation(s)
- Martina Pierri
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy; PhD Program in Drug Discovery and Development, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy
| | - Erica Gazzillo
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy; PhD Program in Drug Discovery and Development, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy
| | - Maria Giovanna Chini
- Department of Biosciences and Territory, University of Molise, C.da Fonte Lappone, Pesche 86090, Italy
| | - Maria Grazia Ferraro
- Department of Pharmacy, School of Medicine and Surgery, University of Naples, Via Domenico Montesano 49, Naples 80131, Italy
| | - Marialuisa Piccolo
- Department of Pharmacy, School of Medicine and Surgery, University of Naples, Via Domenico Montesano 49, Naples 80131, Italy
| | - Francesco Maione
- Department of Pharmacy, School of Medicine and Surgery, University of Naples, Via Domenico Montesano 49, Naples 80131, Italy
| | - Carlo Irace
- Department of Pharmacy, School of Medicine and Surgery, University of Naples, Via Domenico Montesano 49, Naples 80131, Italy
| | - Giuseppe Bifulco
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy
| | - Ines Bruno
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy
| | - Stefania Terracciano
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy.
| | - Gianluigi Lauro
- Department of Pharmacy, University of Salerno, Via Giovanni Paolo II 132, Fisciano 84084, Italy.
| |
Collapse
|
10
|
Lukac I, Wyatt PG, Gilbert IH, Zuccotto F. Ligand binding: evaluating the contribution of the water molecules network using the Fragment Molecular Orbital method. J Comput Aided Mol Des 2021; 35:1025-1036. [PMID: 34458939 PMCID: PMC8523014 DOI: 10.1007/s10822-021-00416-3] [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: 02/21/2021] [Accepted: 08/09/2021] [Indexed: 11/15/2022]
Abstract
Water molecules play a crucial role in protein-ligand binding, and many tools exist that aim to predict the position and relative energies of these important, but challenging participants of biomolecular recognition. The available tools are, in general, capable of predicting the location of water molecules. However, predicting the effects of their displacement is still very challenging. In this work, a linear-scaling quantum mechanics-based approach was used to assess water network energetics and the changes in network stability upon ligand structural modifications. This approach offers a valuable way to improve understanding of SAR data and help guide compound design.
Collapse
Affiliation(s)
- Iva Lukac
- Drug Discovery Unit, Wellcome Centre for Anti-Infectives Research, Division of Biological Chemistry and Drug Discovery, University of Dundee, Dow Street, Dundee, DD1 5EH, UK
| | - Paul G Wyatt
- Drug Discovery Unit, Wellcome Centre for Anti-Infectives Research, Division of Biological Chemistry and Drug Discovery, University of Dundee, Dow Street, Dundee, DD1 5EH, UK.
| | - Ian H Gilbert
- Drug Discovery Unit, Wellcome Centre for Anti-Infectives Research, Division of Biological Chemistry and Drug Discovery, University of Dundee, Dow Street, Dundee, DD1 5EH, UK
| | - Fabio Zuccotto
- Drug Discovery Unit, Wellcome Centre for Anti-Infectives Research, Division of Biological Chemistry and Drug Discovery, University of Dundee, Dow Street, Dundee, DD1 5EH, UK.
| |
Collapse
|
11
|
Andreev S, Pantsar T, Tesch R, Kahlke N, El-Gokha A, Ansideri F, Grätz L, Romasco J, Sita G, Geibel C, Lämmerhofer M, Tarozzi A, Knapp S, Laufer SA, Koch P. Addressing a Trapped High-Energy Water: Design and Synthesis of Highly Potent Pyrimidoindole-Based Glycogen Synthase Kinase-3β Inhibitors. J Med Chem 2021; 65:1283-1301. [PMID: 34213342 DOI: 10.1021/acs.jmedchem.0c02146] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In small molecule binding, water is not a passive bystander but rather takes an active role in the binding site, which may be decisive for the potency of the inhibitor. Here, by addressing a high-energy water, we improved the IC50 value of our co-crystallized glycogen synthase kinase-3β (GSK-3β) inhibitor by nearly two orders of magnitude. Surprisingly, our results demonstrate that this high-energy water was not displaced by our potent inhibitor (S)-3-(3-((7-ethynyl-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propanenitrile ((S)-15, IC50 value of 6 nM). Instead, only a subtle shift in the location of this water molecule resulted in a dramatic decrease in the energy of this high-energy hydration site, as shown by the WaterMap analysis combined with microsecond timescale molecular dynamics simulations. (S)-15 demonstrated both a favorable kinome selectivity profile and target engagement in a cellular environment and reduced GSK-3 autophosphorylation in neuronal SH-SY5Y cells. Overall, our findings highlight that even a slight adjustment in the location of a high-energy water can be decisive for ligand binding.
Collapse
Affiliation(s)
- Stanislav Andreev
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Tatu Pantsar
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany.,School of Pharmacy, Faculty of Health Sciences, University of Eastern Finland, P.O. Box 1627, 70211 Kuopio, Finland
| | - Roberta Tesch
- Institute for Pharmaceutical Chemistry, Goethe University, Max-von-Laue-Straße 9, 60438 Frankfurt am Main, Germany.,Structural Genomics Consortium, Buchmann Institute for Life Sciences, Goethe University, Max-von-Laue-Straße 15, 60438 Frankfurt am Main, Germany
| | - Niclas Kahlke
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Ahmed El-Gokha
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany.,Chemistry Department, Faculty of Science, Menoufia University, Gamal Abdel-Nasser Street, 32511 Shebin El-Kom, Egypt
| | - Francesco Ansideri
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Lukas Grätz
- Department of Pharmaceutical/Medicinal Chemistry II, Institute of Pharmacy, University of Regensburg, Universitätsstraße 31, 93053 Regensburg, Germany
| | - Jenny Romasco
- Department for Life Quality Studies, Alma Mater Studiorum, University of Bologna, Corso D'Augusto 237, 47921 Rimini, Italy
| | - Giulia Sita
- Department of Pharmacy and Biotechnology, University of Bologna, Via Irnerio 48, 40126 Bologna, Italy
| | - Christian Geibel
- Institute of Pharmaceutical Sciences, Department of Pharmaceutical (Bio-)Analysis, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Michael Lämmerhofer
- Institute of Pharmaceutical Sciences, Department of Pharmaceutical (Bio-)Analysis, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Andrea Tarozzi
- Department for Life Quality Studies, Alma Mater Studiorum, University of Bologna, Corso D'Augusto 237, 47921 Rimini, Italy
| | - Stefan Knapp
- Institute for Pharmaceutical Chemistry, Goethe University, Max-von-Laue-Straße 9, 60438 Frankfurt am Main, Germany.,Structural Genomics Consortium, Buchmann Institute for Life Sciences, Goethe University, Max-von-Laue-Straße 15, 60438 Frankfurt am Main, Germany
| | - Stefan A Laufer
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany.,Tübingen Center for Academic Drug Discovery (TüCAD2), Auf der Morgenstelle 8, 72076 Tübingen, Germany
| | - Pierre Koch
- Institute of Pharmaceutical Sciences, Department of Medicinal and Pharmaceutical Chemistry, Eberhard Karls University Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany.,Department of Pharmaceutical/Medicinal Chemistry II, Institute of Pharmacy, University of Regensburg, Universitätsstraße 31, 93053 Regensburg, Germany
| |
Collapse
|
12
|
Starciuc T, Guinet Y, Hedoux A, Shalaev E. Water content thresholds in glycerol/water system: Low- and high-wavenumber Raman spectroscopy study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114678] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
13
|
Wang Y, Wu S, Wang L, Yang Z, Zhao J, Zhang L. Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations. RSC Adv 2020; 11:745-759. [PMID: 35423696 PMCID: PMC8693360 DOI: 10.1039/d0ra09469b] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2020] [Accepted: 12/17/2020] [Indexed: 12/12/2022] Open
Abstract
Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.
Collapse
Affiliation(s)
- Yan Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Shiliang Wu
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lifei Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University Nanchang 330045 China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lulu Zhang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| |
Collapse
|
14
|
Shanu-Wilson J, Evans L, Wrigley S, Steele J, Atherton J, Boer J. Biotransformation: Impact and Application of Metabolism in Drug Discovery. ACS Med Chem Lett 2020; 11:2087-2107. [PMID: 33214818 DOI: 10.1021/acsmedchemlett.0c00202] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Accepted: 08/13/2020] [Indexed: 02/07/2023] Open
Abstract
Biotransformation has a huge impact on the efficacy and safety of drugs. Ultimately the effects of metabolism can be the lynchpin in the discovery and development cycle of a new drug. This article discusses the impact and application of biotransformation of drugs by mammalian systems, microorganisms, and recombinant enzymes, covering active and reactive metabolites, the impact of the gut microbiome on metabolism, and how insights gained from biotransformation studies can influence drug design from the combined perspectives of a CRO specializing in a range of biotransformation techniques and pharma biotransformation scientists. We include a commentary on how biology-driven approaches can complement medicinal chemistry strategies in drug optimization and the in vitro and surrogate systems available to explore and exploit biotransformation.
Collapse
Affiliation(s)
- Julia Shanu-Wilson
- Hypha Discovery Ltd., 154B Brook Drive, Milton Park, Abingdon, Oxfordshire OX14 4SD, U.K
| | - Liam Evans
- Hypha Discovery Ltd., 154B Brook Drive, Milton Park, Abingdon, Oxfordshire OX14 4SD, U.K
| | - Stephen Wrigley
- Hypha Discovery Ltd., 154B Brook Drive, Milton Park, Abingdon, Oxfordshire OX14 4SD, U.K
| | - Jonathan Steele
- Hypha Discovery Ltd., 154B Brook Drive, Milton Park, Abingdon, Oxfordshire OX14 4SD, U.K
| | - James Atherton
- Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Jason Boer
- Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| |
Collapse
|
15
|
Shah B, Sindhikara D, Borrelli K, Leffler AE. Water Thermodynamics of Peptide Toxin Binding Sites on Ion Channels. Toxins (Basel) 2020; 12:toxins12100652. [PMID: 33053750 PMCID: PMC7599640 DOI: 10.3390/toxins12100652] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Revised: 09/21/2020] [Accepted: 10/09/2020] [Indexed: 12/20/2022] Open
Abstract
Peptide toxins isolated from venomous creatures, long prized as research tools due to their innate potency for ion channels, are emerging as drugs as well. However, it remains challenging to understand why peptide toxins bind with high potency to ion channels, to identify residues that are key for activity, and to improve their affinities via mutagenesis. We use WaterMap, a molecular dynamics simulation-based method, to gain computational insight into these three questions by calculating the locations and thermodynamic properties of water molecules in the peptide toxin binding sites of five ion channels. These include an acid-sensing ion channel, voltage-gated potassium channel, sodium channel in activated and deactivated states, transient-receptor potential channel, and a nicotinic receptor whose structures were recently determined by crystallography and cryo-electron microscopy (cryo-EM). All channels had water sites in the peptide toxin binding site, and an average of 75% of these sites were stable (low-energy), and 25% were unstable (medium or high energy). For the sodium channel, more unstable water sites were present in the deactivated state structure than the activated. Additionally, for each channel, unstable water sites coincided with the positions of peptide toxin residues that previous mutagenesis experiments had shown were important for activity. Finally, for the sodium channel in the deactivated state, unstable water sites were present in the peptide toxin binding pocket but did not overlap with the peptide toxin, suggesting that future experimental efforts could focus on targeting these sites to optimize potency.
Collapse
Affiliation(s)
- Binita Shah
- Schrödinger, Inc. 120 W. 45th St, New York, NY 10036, USA; (B.S.); (D.S.); (K.B.)
- PhD Program in Biological and Biomedical Sciences, Division of Medical Sciences, Harvard Medical School, Boston, MA 02115, USA
| | - Dan Sindhikara
- Schrödinger, Inc. 120 W. 45th St, New York, NY 10036, USA; (B.S.); (D.S.); (K.B.)
| | - Ken Borrelli
- Schrödinger, Inc. 120 W. 45th St, New York, NY 10036, USA; (B.S.); (D.S.); (K.B.)
| | - Abba E. Leffler
- Schrödinger, Inc. 120 W. 45th St, New York, NY 10036, USA; (B.S.); (D.S.); (K.B.)
- Correspondence:
| |
Collapse
|
16
|
Pavlovicz RE, Park H, DiMaio F. Efficient consideration of coordinated water molecules improves computational protein-protein and protein-ligand docking discrimination. PLoS Comput Biol 2020; 16:e1008103. [PMID: 32956350 PMCID: PMC7529342 DOI: 10.1371/journal.pcbi.1008103] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Revised: 10/01/2020] [Accepted: 06/29/2020] [Indexed: 12/25/2022] Open
Abstract
Highly coordinated water molecules are frequently an integral part of protein-protein and protein-ligand interfaces. We introduce an updated energy model that efficiently captures the energetic effects of these ordered water molecules on the surfaces of proteins. A two-stage method is developed in which polar groups arranged in geometries suitable for water placement are first identified, then a modified Monte Carlo simulation allows highly coordinated waters to be placed on the surface of a protein while simultaneously sampling amino acid side chain orientations. This “semi-explicit” water model is implemented in Rosetta and is suitable for both structure prediction and protein design. We show that our new approach and energy model yield significant improvements in native structure recovery of protein-protein and protein-ligand docking discrimination tests. Well-coordinated water molecules—those forming multiple hydrogen bonds with nearby polar groups—play an important role in the structure of biomolecular systems, yet the effect of these waters is often not considered in molecular energy computations. In this paper, we describe a method to efficiently consider these water molecules both implicitly and explicitly at the interfaces formed by two polar molecules. In computations related to determining how a protein interacts with binding partners, we show that the use of this new method significantly improves results. Future application of this approach may improve the design of new protein and small molecule drugs.
Collapse
Affiliation(s)
- Ryan E. Pavlovicz
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
| | - Hahnbeom Park
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
| | - Frank DiMaio
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- Institute for Protein Design, University of Washington, Seattle, Washington, United States of America
- * E-mail:
| |
Collapse
|
17
|
The role of water in ligand binding. Curr Opin Struct Biol 2020; 67:1-8. [PMID: 32942197 DOI: 10.1016/j.sbi.2020.08.002] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 08/11/2020] [Accepted: 08/13/2020] [Indexed: 02/06/2023]
|
18
|
Clegg MA, Bamborough P, Chung CW, Craggs PD, Gordon L, Grandi P, Leveridge M, Lindon M, Liwicki GM, Michon AM, Molnar J, Rioja I, Soden PE, Theodoulou NH, Werner T, Tomkinson NCO, Prinjha RK, Humphreys PG. Application of Atypical Acetyl-lysine Methyl Mimetics in the Development of Selective Inhibitors of the Bromodomain-Containing Protein 7 (BRD7)/Bromodomain-Containing Protein 9 (BRD9) Bromodomains. J Med Chem 2020; 63:5816-5840. [PMID: 32410449 DOI: 10.1021/acs.jmedchem.0c00075] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Non-BET bromodomain-containing proteins have become attractive targets for the development of novel therapeutics targeting epigenetic pathways. To help facilitate the target validation of this class of proteins, structurally diverse small-molecule ligands and methodologies to produce selective inhibitors in a predictable fashion are in high demand. Herein, we report the development and application of atypical acetyl-lysine (KAc) methyl mimetics to take advantage of the differential stability of conserved water molecules in the bromodomain binding site. Discovery of the n-butyl group as an atypical KAc methyl mimetic allowed generation of 31 (GSK6776) as a soluble, permeable, and selective BRD7/9 inhibitor from a pyridazinone template. The n-butyl group was then used to enhance the bromodomain selectivity of an existing BRD9 inhibitor and to transform pan-bromodomain inhibitors into BRD7/9 selective compounds. Finally, a solvent-exposed vector was defined from the pyridazinone template to enable bifunctional molecule synthesis, and affinity enrichment chemoproteomic experiments were used to confirm several of the endogenous protein partners of BRD7 and BRD9, which form part of the chromatin remodeling PBAF and BAF complexes, respectively.
Collapse
Affiliation(s)
- Michael A Clegg
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom.,WestCHEM, Department of Pure and Applied Chemistry, Thomas Graham Building, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| | - Paul Bamborough
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Chun-Wa Chung
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Peter D Craggs
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Laurie Gordon
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Paola Grandi
- Cellzome GmbH, R&D MST GlaxoSmithKline, Meyerhofstrasse 1 69117 Heidelberg, Germany
| | - Melanie Leveridge
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Matthew Lindon
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Gemma M Liwicki
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Anne-Marie Michon
- Cellzome GmbH, R&D MST GlaxoSmithKline, Meyerhofstrasse 1 69117 Heidelberg, Germany
| | - Judit Molnar
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Inmaculada Rioja
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Peter E Soden
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | - Natalie H Theodoulou
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom.,WestCHEM, Department of Pure and Applied Chemistry, Thomas Graham Building, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| | - Thilo Werner
- Cellzome GmbH, R&D MST GlaxoSmithKline, Meyerhofstrasse 1 69117 Heidelberg, Germany
| | - Nicholas C O Tomkinson
- WestCHEM, Department of Pure and Applied Chemistry, Thomas Graham Building, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, United Kingdom
| | - Rab K Prinjha
- GlaxoSmithKline R&D, Stevenage SG1 2NY, Hertfordshire, United Kingdom
| | | |
Collapse
|
19
|
Thapa B, Erickson J, Raghavachari K. Quantum Mechanical Investigation of Three-Dimensional Activity Cliffs Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Inf Model 2020; 60:2924-2938. [PMID: 32407081 DOI: 10.1021/acs.jcim.9b01123] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The concept of activity cliff (AC) (i.e., a small structural modification resulting in a substantial bioactivity change) is widely encountered in medicinal chemistry during compound design. Whereas the study of ACs is of high interest as it provides a wealth of opportunities for effective drug design, its practical application in the actual drug development process has been difficult because of significant computational challenges. To provide some understanding of the ACs, we have carried out a rigorous quantum-mechanical investigation of the electronic interactions of a wide range of ACs (205 cliffs formed by 261 protein-ligand complexes covering 37 different receptor types) using multilayer molecules-in-molecules (MIM) fragmentation-based methodology. The MIM methodology enables performing accurate high-level quantum mechanical (QM) calculations at a substantially lower computational cost, while allowing for a quantitative decomposition of the protein-ligand binding energy into the contributions from individual residues, solvation, and entropy. Our investigation in this study is mainly focused on whether the QM binding energy calculation can correctly identify the higher potency cliff partner for a given ligand pair having a sufficiently high activity difference. We have also analyzed the effect of including crystal water molecules as a part of the receptor as well as the impact of ligand desolvation energy on the correct identification of the more potent ligand in a cliff pair. Our analysis reveals that, in the majority of the cases, the AC prediction could be significantly improved by carefully identifying the critical crystal water molecules, whereas the contribution from the ligand desolvation also remains essential. Additionally, we have exploited the residue-specific interaction energies provided by MIM to identify the key residues and interaction hot-spots that are responsible for the experimentally observed drastic activity changes. The results show that our MIM fragmentation-based protocol provides comprehensive interaction energy profiles that can be employed to understand the distinctiveness of ligand modifications, for potential applications in structure-based drug design.
Collapse
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Jon Erickson
- Lilly Research Laboratories, Eli Lilly & Company, Indianapolis, Indiana 46285, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
20
|
Wang LF, Wang Y, Yang ZY, Zhao J, Sun HB, Wu SL. Revealing binding selectivity of inhibitors toward bromodomain-containing proteins 2 and 4 using multiple short molecular dynamics simulations and free energy analyses. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2020; 31:373-398. [PMID: 32496901 DOI: 10.1080/1062936x.2020.1748107] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Accepted: 03/24/2020] [Indexed: 06/11/2023]
Abstract
Emerging evidences indicate bromodomain-containing proteins 2 and 4 (BRD2 and BRD4) play critical roles in cancers, inflammations, cardiovascular diseases and other pathologies. Multiple short molecular dynamics (MSMD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were applied to investigate the binding selectivity of three inhibitors 87D, 88M and 89G towards BRD2 over BRD4. The root-mean-square fluctuation (RMSF) analysis indicates that the structural flexibility of BRD4 is stronger than that of BRD2. Moreover the calculated distances between the Cα atoms in the centres of the ZA_loop and BC_loop of BRD4 are also bigger than that of BRD2. The rank of binding free energies calculated using MM-GBSA method agrees well with that determined by experimental data. The results show that 87D can bind more favourably to BRD2 than BRD4, while 88M has better selectivity on BRD4 over BRD2. Residue-based free-energy decomposition method was utilized to estimate the inhibitor-residue interaction spectrum and the results not only identify the hot interaction spots of inhibitors with BRD2 and BRD4, but also demonstrate that several common residues, including (W370, W374), (P371, P375), (V376, V380) and (L381, L385) belonging to (BRD2, BRD4), generate significant binding difference of inhibitors to BRD2 and BRD4.
Collapse
Affiliation(s)
- L F Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - Y Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - Z Y Yang
- Department of Physics, Jiangxi Agricultural University , Nanchang, China
| | - J Zhao
- School of Science, Shandong Jiaotong University , Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University , Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University , Jinan, China
| |
Collapse
|
21
|
Bodnarchuk MS, Packer MJ, Haywood A. Utilizing Grand Canonical Monte Carlo Methods in Drug Discovery. ACS Med Chem Lett 2020; 11:77-82. [PMID: 31938467 DOI: 10.1021/acsmedchemlett.9b00499] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 12/11/2019] [Indexed: 12/17/2022] Open
Abstract
The concepts behind targeting waters for potency and selectivity gains have been well documented and explored, although maximizing such potential gains can prove to be challenging. This problem is exacerbated in cases where there are multiple interacting waters, wherein perturbation of one water can affect the free energy landscape of the remaining waters. Knowing the right modification a priori is challenging, and computational approaches are ideally suited to help answer the key question of which substitution is best to try. Here, we use Grand Canonical Monte Carlo and the recent Grand Canonical Alchemical Perturbation methods to both understand and predict the effect of ligand-mediated water displacement when more than one water molecule is involved, as well as to understand how exploiting water networks can help govern selectivity.
Collapse
Affiliation(s)
- Michael S. Bodnarchuk
- Computational Chemistry, R&D Oncology, AstraZeneca, Cambridge CB4 0WG, United Kingdom
| | - Martin J. Packer
- Computational Chemistry, R&D Oncology, AstraZeneca, Cambridge CB4 0WG, United Kingdom
| | - Alexe Haywood
- Computational Chemistry, R&D Oncology, AstraZeneca, Cambridge CB4 0WG, United Kingdom
| |
Collapse
|
22
|
A Molecular Dynamics Study of Vasoactive Intestinal Peptide Receptor 1 and the Basis of Its Therapeutic Antagonism. Int J Mol Sci 2019; 20:ijms20184348. [PMID: 31491880 PMCID: PMC6770453 DOI: 10.3390/ijms20184348] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2019] [Revised: 07/30/2019] [Accepted: 08/20/2019] [Indexed: 12/21/2022] Open
Abstract
Vasoactive intestinal peptide receptor 1 (VPAC1) is a member of a secretin-like subfamily of G protein-coupled receptors. Its endogenous neuropeptide (VIP), secreted by neurons and immune cells, modulates various physiological functions such as exocrine and endocrine secretions, immune response, smooth muscles relaxation, vasodilation, and fetal development. As a drug target, VPAC1 has been selected for therapy of inflammatory diseases but drug discovery is still hampered by lack of its crystal structure. In this study we presented the homology model of this receptor constructed with the well-known web service GPCRM. The VPAC1 model is composed of extracellular and transmembrane domains that form a complex with an endogenous hormone VIP. Using the homology model of VPAC1 the mechanism of action of potential drug candidates for VPAC1 was described. Only two series of small-molecule antagonists of confirmed biological activity for VPAC1 have been described thus far. Molecular docking and a series of molecular dynamics simulations were performed to elucidate their binding to VPAC1 and resulting antagonist effect. The presented work provides the basis for the possible binding mode of VPAC1 antagonists and determinants of their molecular recognition in the context of other class B GPCRs. Until the crystal structure of VPAC1 will be released, the presented homology model of VPAC1 can serve as a scaffold for drug discovery studies and is available from the author upon request.
Collapse
|
23
|
Nittinger E, Flachsenberg F, Bietz S, Lange G, Klein R, Rarey M. Placement of Water Molecules in Protein Structures: From Large-Scale Evaluations to Single-Case Examples. J Chem Inf Model 2018; 58:1625-1637. [DOI: 10.1021/acs.jcim.8b00271] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Eva Nittinger
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Florian Flachsenberg
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Stefan Bietz
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Gudrun Lange
- Bayer CropScience AG, Industriepark Hoechst G836, 65926 Frankfurt am Main, Germany
| | - Robert Klein
- Bayer CropScience AG, Industriepark Hoechst G836, 65926 Frankfurt am Main, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| |
Collapse
|