1
|
de Oliveira C, Leswing K, Feng S, Kanters R, Abel R, Bhat S. FEP Protocol Builder: Optimization of Free Energy Perturbation Protocols Using Active Learning. J Chem Inf Model 2023; 63:5592-5603. [PMID: 37594480 DOI: 10.1021/acs.jcim.3c00681] [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: 08/19/2023]
Abstract
Significant improvements have been made in the past decade to methods that rapidly and accurately predict binding affinity through free energy perturbation (FEP) calculations. This has been driven by recent advances in small-molecule force fields and sampling algorithms combined with the availability of low-cost parallel computing. Predictive accuracies of ∼1 kcal mol-1 have been regularly achieved, which are sufficient to drive potency optimization in modern drug discovery campaigns. Despite the robustness of these FEP approaches across multiple target classes, there are invariably target systems that do not display expected performance with default FEP settings. Traditionally, these systems required labor-intensive manual protocol development to arrive at parameter settings that produce a predictive FEP model. Due to the (a) relatively large parameter space to be explored, (b) significant compute requirements, and (c) limited understanding of how combinations of parameters can affect FEP performance, manual FEP protocol optimization can take weeks to months to complete, and often does not involve rigorous train-test set splits, resulting in potential overfitting. These manual FEP protocol development timelines do not coincide with tight drug discovery project timelines, essentially preventing the use of FEP calculations for these target systems. Here, we describe an automated workflow termed FEP Protocol Builder (FEP-PB) to rapidly generate accurate FEP protocols for systems that do not perform well with default settings. FEP-PB uses an active-learning workflow to iteratively search the protocol parameter space to develop accurate FEP protocols. To validate this approach, we applied it to pharmaceutically relevant systems where default FEP settings could not produce predictive models. We demonstrate that FEP-PB can rapidly generate accurate FEP protocols for the previously challenging MCL1 system with limited human intervention. We also apply FEP-PB in a real-world drug discovery setting to generate an accurate FEP protocol for the p97 system. FEP-PB is able to generate a more accurate protocol than the expert user, rapidly validating p97 as amenable to free energy calculations. Additionally, through the active-learning workflow, we are able to gain insight into which parameters are most important for a given system. These results suggest that FEP-PB is a robust tool that can aid in rapidly developing accurate FEP protocols and increasing the number of targets that are amenable to the technology.
Collapse
Affiliation(s)
- César de Oliveira
- Schrodinger, Inc., 9868 Scranton Road, Suite 3200, San Diego, California 92121, United States
| | - Karl Leswing
- Schrodinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Shulu Feng
- Schrodinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - René Kanters
- Schrodinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Robert Abel
- Schrodinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Sathesh Bhat
- Schrodinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
2
|
Dutkiewicz Z, Mikstacka R. Hydration and Structural Adaptations of the Human CYP1A1, CYP1A2, and CYP1B1 Active Sites by Molecular Dynamics Simulations. Int J Mol Sci 2023; 24:11481. [PMID: 37511239 PMCID: PMC10380238 DOI: 10.3390/ijms241411481] [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/29/2023] [Accepted: 07/11/2023] [Indexed: 07/30/2023] Open
Abstract
Cytochromes CYP1A1, CYP1A2, and CYP1B1, the members of the cytochrome P450 family 1, catalyze the metabolism of endogenous compounds, drugs, and non-drug xenobiotics which include substances involved in the process of carcinogenesis, cancer chemoprevention, and therapy. In the present study, the interactions of three selected polymethoxy-trans-stilbenes, analogs of a bioactive polyphenol trans-resveratrol (3,5,4'-trihydroxy-trans-stilbene) with the binding sites of CYP1 isozymes were investigated with molecular dynamics (MD) simulations. The most pronounced structural changes in the CYP1 binding sites were observed in two substrate recognition sites (SRS): SRS2 (helix F) and SRS3 (helix G). MD simulations show that the number and position of water molecules occurring in CYP1 APO and in the structures complexed with ligands are diverse. The presence of water in binding sites results in the formation of water-protein, water-ligand, and bridging ligand-water-protein hydrogen bonds. Analysis of the solvent and substrate channels opening during the MD simulation showed significant differences between cytochromes in relation to the solvent channel and the substrate channels 2c, 2ac, and 2f. The results of this investigation lead to a deeper understanding of the molecular processes that occur in the CYP1 binding sites and may be useful for further molecular studies of CYP1 functions.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs, Poznan University of Medical Sciences, Grunwaldzka 6, 60-780 Poznań, Poland
| | - Renata Mikstacka
- Department of Inorganic and Analytical Chemistry, Nicolaus Copernicus University, Collegium Medicum, Dr. A. Jurasza 2, 85-089 Bydgoszcz, Poland
| |
Collapse
|
3
|
Barros EP, Ries B, Champion C, Rieder SR, Riniker S. Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities. J Chem Inf Model 2023; 63:1794-1805. [PMID: 36917685 PMCID: PMC10052353 DOI: 10.1021/acs.jcim.2c01610] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/15/2023]
Abstract
Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.
Collapse
Affiliation(s)
- Emilia P Barros
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Benjamin Ries
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Candide Champion
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Salomé R Rieder
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
4
|
Melling O, Samways ML, Ge Y, Mobley DL, Essex JW. Enhanced Grand Canonical Sampling of Occluded Water Sites Using Nonequilibrium Candidate Monte Carlo. J Chem Theory Comput 2023; 19:1050-1062. [PMID: 36692215 PMCID: PMC9933432 DOI: 10.1021/acs.jctc.2c00823] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Indexed: 01/25/2023]
Abstract
Water molecules play a key role in many biomolecular systems, particularly when bound at protein-ligand interfaces. However, molecular simulation studies on such systems are hampered by the relatively long time scales over which water exchange between a protein and solvent takes place. Grand canonical Monte Carlo (GCMC) is a simulation technique that avoids this issue by attempting the insertion and deletion of water molecules within a given structure. The approach is constrained by low acceptance probabilities for insertions in congested systems, however. To address this issue, here, we combine GCMC with nonequilibium candidate Monte Carlo (NCMC) to yield a method that we refer to as grand canonical nonequilibrium candidate Monte Carlo (GCNCMC), in which the water insertions and deletions are carried out in a gradual, nonequilibrium fashion. We validate this new approach by comparing GCNCMC and GCMC simulations of bulk water and three protein binding sites. We find that not only is the efficiency of the water sampling improved by GCNCMC but that it also results in increased sampling of ligand conformations in a protein binding site, revealing new water-mediated ligand-binding geometries that are not observed using alternative enhanced sampling techniques.
Collapse
Affiliation(s)
- Oliver
J. Melling
- School
of Chemistry, University of Southampton, SouthamptonSO17 1BJ, U.K.
| | - Marley L. Samways
- School
of Chemistry, University of Southampton, SouthamptonSO17 1BJ, U.K.
| | - Yunhui Ge
- Department
of Pharmaceutical Sciences, University of
California, Irvine, California92697, United States
| | - David L. Mobley
- Department
of Pharmaceutical Sciences, University of
California, Irvine, California92697, United States
- Department
of Chemistry, University of California, Irvine, California92697, United States
| | - Jonathan W. Essex
- School
of Chemistry, University of Southampton, SouthamptonSO17 1BJ, U.K.
| |
Collapse
|
5
|
Chen W, He H, Wang J, Wang J, Chang CEA. Uncovering water effects in protein-ligand recognition: importance in the second hydration shell and binding kinetics. Phys Chem Chem Phys 2023; 25:2098-2109. [PMID: 36562309 PMCID: PMC9970846 DOI: 10.1039/d2cp04584b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Developing a ligand with high affinity for a specific protein target is essential for drug design, and water molecules are well known to play a key role in protein-drug recognition. However, predicting the role of particularly ordered water molecules in drug binding remains challenging. Furthermore, hydration free energy contributed from the water network, including the second shell of water molecules, is far from being well studied. In this research we focused on these aspects to accurately and efficiently evaluate water effects in protein-ligand binding affinity. We developed a new strategy using a free-energy calculation method, VM2. We successfully predicted the stable ordered water molecules in a number of protein systems: PDE 10a, HSP90, tryptophan synthase (TRPS), CDK2 and Factor Xa. In some of these, the second shell of water molecules appeared to be critical in protein-ligand binding. We also applied the strategy to largely improve binding free energy calculation using the MM/PBSA method. When applying MM/PBSA alone for two systems, CDK2 and Factor Xa, the computed binding free energy resulted in poor to moderate R2 values with experimental data. However, including water free energy correction greatly improved the free energy calculation. Furthermore, our work helped to explain how xk263 is a 1000 times faster binder to HIVp than ritonavir, a potentially useful tool for investigating binding kinetics. Our studies reveal the importance of fully considering water effects in therapeutic developments in pharmaceutical and biotechnology industries and for fundamental research in protein-ligand recognition.
Collapse
Affiliation(s)
- Wei Chen
- School of Pharmacy, Fuzhou Medical College of NanChang University, Fuzhou, JiangXi 344000, P. R. China.
| | - Huan He
- School of Pharmacy, Fuzhou Medical College of NanChang University, Fuzhou, JiangXi 344000, P. R. China.
| | - Jing Wang
- School of Pharmacy, Fuzhou Medical College of NanChang University, Fuzhou, JiangXi 344000, P. R. China.
| | - Jiahui Wang
- School of Pharmacy, Fuzhou Medical College of NanChang University, Fuzhou, JiangXi 344000, P. R. China.
| | - Chia-en A. Chang
- Department of Chemistry, University of California at Riverside, Riverside, CA 92521, USA,Corresponding Authors Wei Chen: School of Pharmacy, Fuzhou Medical College of NanChang University, Fuzhou, JiangXi 344000, China. , Chia-en A. Chang: Department of Chemistry, University of California at Riverside, Riverside, CA 92521, USA.
| |
Collapse
|
6
|
Samways M, Bruce Macdonald HE, Taylor RD, Essex JW. Water Networks in Complexes between Proteins and FDA-Approved Drugs. J Chem Inf Model 2023; 63:387-396. [PMID: 36469670 PMCID: PMC9832485 DOI: 10.1021/acs.jcim.2c01225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Water molecules at protein-ligand interfaces are often of significant pharmaceutical interest, owing in part to the entropy which can be released upon the displacement of an ordered water by a therapeutic compound. Protein structures may not, however, completely resolve all critical bound water molecules, or there may be no experimental data available. As such, predicting the location of water molecules in the absence of a crystal structure is important in the context of rational drug design. Grand canonical Monte Carlo (GCMC) is a computational technique that is gaining popularity for the simulation of buried water sites. In this work, we assess the ability of GCMC to accurately predict water binding locations, using a dataset that we have curated, containing 108 unique structures of complexes between proteins and Food and Drug Administration (FDA)-approved small-molecule drugs. We show that GCMC correctly predicts 81.4% of nonbulk crystallographic water sites to within 1.4 Å. However, our analysis demonstrates that the reported performance of water prediction methods is highly sensitive to the way in which the performance is measured. We also find that crystallographic water sites with more protein/ligand hydrogen bonds and stronger electron density are more reliably predicted by GCMC. An analysis of water networks revealed that more than half of the structures contain at least one ligand-contacting water network. In these cases, displacement of a water site by a ligand modification might yield unexpected results if the larger network is destabilized. Cooperative effects between waters should therefore be explicitly considered in structure-based drug design.
Collapse
Affiliation(s)
- Marley
L. Samways
- School
of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.
| | - Hannah E. Bruce Macdonald
- Computational
and Systems Biology Program, Memorial Sloan
Kettering Cancer Center, New York, New York 10065, United States
| | | | - Jonathan W. Essex
- School
of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.,
| |
Collapse
|
7
|
Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization. Life (Basel) 2022; 12:life12091438. [PMID: 36143474 PMCID: PMC9504344 DOI: 10.3390/life12091438] [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: 07/25/2022] [Revised: 09/02/2022] [Accepted: 09/03/2022] [Indexed: 11/27/2022] Open
Abstract
Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying “magic bullets” amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand–target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural–energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.
Collapse
|
8
|
Mitryukovskiy S, Vanpoucke DEP, Bai Y, Hannotte T, Lavancier M, Hourlier D, Roos G, Peretti R. On the influence of water on THz vibrational spectral features of molecular crystals. Phys Chem Chem Phys 2022; 24:6107-6125. [PMID: 35212691 DOI: 10.1039/d1cp03261e] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The nanoscale structure of molecular assemblies plays a major role in many (μ)-biological mechanisms. Molecular crystals are one of the most simple of these assemblies and are widely used in a variety of applications from pharmaceuticals and agrochemicals, to nutraceuticals and cosmetics. The collective vibrations in such molecular crystals can be probed using terahertz spectroscopy, providing unique characteristic spectral fingerprints. However, the association of the spectral features to the crystal conformation, crystal phase and its environment is a difficult task. We present a combined computational-experimental study on the incorporation of water in lactose molecular crystals, and show how simulations can be used to associate spectral features in the THz region to crystal conformations and phases. Using periodic DFT simulations of lactose molecular crystals, the role of water in the observed lactose THz spectrum is clarified, presenting both direct and indirect contributions. A specific experimental setup is built to allow the controlled heating and corresponding dehydration of the sample, providing the monitoring of the crystal phase transformation dynamics. Besides the observation that lactose phases and phase transformation appear to be more complex than previously thought - including several crystal forms in a single phase and a non-negligible water content in the so-called anhydrous phase - we draw two main conclusions from this study. Firstly, THz modes are spread over more than one molecule and require periodic computation rather than a gas-phase one. Secondly, hydration water does not only play a perturbative role but also participates in the facilitation of the THz vibrations.
Collapse
Affiliation(s)
- Sergey Mitryukovskiy
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| | - Danny E P Vanpoucke
- IMO, Hasselt University, 3590 Diepenbeek, Belgium./AMIBM, Maastricht University, 6167 Geleen, The Netherlands
| | - Yue Bai
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| | - Théo Hannotte
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| | - Mélanie Lavancier
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| | - Djamila Hourlier
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| | - Goedele Roos
- Univ. Lille, CNRS, UMR 8576 - UGSF - Unité de Glycobiologie Structurale et Fonctionnelle, F-59000 Lille, France
| | - Romain Peretti
- Institut d'Electronique de Microélectronique et de Nanotechnologie, Université Lille, CNRS, 59652 Villeneuve d'Ascq, France.
| |
Collapse
|
9
|
Kolos JM, Pomplun S, Jung S, Rieß B, Purder PL, Voll AM, Merz S, Gnatzy M, Geiger TM, Quist-Løkken I, Jatzlau J, Knaus P, Holien T, Bracher A, Meyners C, Czodrowski P, Krewald V, Hausch F. Picomolar FKBP inhibitors enabled by a single water-displacing methyl group in bicyclic [4.3.1] aza-amides. Chem Sci 2021; 12:14758-14765. [PMID: 34820091 PMCID: PMC8597852 DOI: 10.1039/d1sc04638a] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2021] [Accepted: 10/22/2021] [Indexed: 01/30/2023] Open
Abstract
Methyl groups can have profound effects in drug discovery but the underlying mechanisms are diverse and incompletely understood. Here we report the stereospecific effect of a single, solvent-exposed methyl group in bicyclic [4.3.1] aza-amides, robustly leading to a 2 to 10-fold increase in binding affinity for FK506-binding proteins (FKBPs). This resulted in the most potent and efficient FKBP ligands known to date. By a combination of co-crystal structures, isothermal titration calorimetry (ITC), density-functional theory (DFT), and 3D reference interaction site model (3D-RISM) calculations we elucidated the origin of the observed affinity boost, which was purely entropically driven and relied on the displacement of a water molecule at the protein-ligand-bulk solvent interface. The best compounds potently occupied FKBPs in cells and enhanced bone morphogenic protein (BMP) signaling. Our results show how subtle manipulation of the solvent network can be used to design atom-efficient ligands for difficult, solvent-exposed binding pockets.
Collapse
Affiliation(s)
- Jürgen M Kolos
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany .,Max Planck Institute of Psychiatry Kraepelinstr. 2-10 80804 München Germany
| | - Sebastian Pomplun
- Max Planck Institute of Psychiatry Kraepelinstr. 2-10 80804 München Germany
| | - Sascha Jung
- Technische Universität Dortmund, Fakultät für Chemie und Chemische Biologie Otto-Hahn-Straße 6 44227 Dortmund Germany
| | - Benedikt Rieß
- Max Planck Institute of Psychiatry Kraepelinstr. 2-10 80804 München Germany
| | - Patrick L Purder
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Andreas M Voll
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Stephanie Merz
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Monika Gnatzy
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Thomas M Geiger
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Ingrid Quist-Løkken
- Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology 7491 Trondheim Norway.,Department of Immunology and Transfusion Medicine, St. Olav's University Hospital 7030 Trondheim Norway.,Department of Hematology, St. Olav's University Hospital 7030 Trondheim Norway
| | - Jerome Jatzlau
- Institute for Chemistry and Biochemistry, Freie Universität Berlin 14195 Berlin Germany
| | - Petra Knaus
- Institute for Chemistry and Biochemistry, Freie Universität Berlin 14195 Berlin Germany
| | - Toril Holien
- Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology 7491 Trondheim Norway.,Department of Immunology and Transfusion Medicine, St. Olav's University Hospital 7030 Trondheim Norway.,Department of Hematology, St. Olav's University Hospital 7030 Trondheim Norway
| | - Andreas Bracher
- Research Department Cellular Biochemistry, Max Planck Institute of Biochemistry Am Klopferspitz 18, 82152 Planegg Germany
| | - Christian Meyners
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Paul Czodrowski
- Technische Universität Dortmund, Fakultät für Chemie und Chemische Biologie Otto-Hahn-Straße 6 44227 Dortmund Germany
| | - Vera Krewald
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| | - Felix Hausch
- Department of Chemistry, Technical University of Darmstadt Alarich-Weiss-Straße 4 64293 Darmstadt Germany
| |
Collapse
|
10
|
Biswal J, Jayaprakash P, Rayala SK, Venkatraman G, Rangaswamy R, Jeyaraman J. WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches. ACS OMEGA 2021; 6:26829-26845. [PMID: 34693105 PMCID: PMC8529594 DOI: 10.1021/acsomega.1c02032] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Indexed: 06/13/2023]
Abstract
p21-Activated kinase 1 (PAK1) is positioned at the nexus of several oncogenic signaling pathways. Currently, there are no approved inhibitors for disabling the transfer of phosphate in the active site directly, as they are limited by lower affinity, and poor kinase selectivity. In this work, a repurposing study utilizing FDA-approved drugs from the DrugBank database was pursued with an initial selection of 27 molecules out of ∼2162 drug molecules, based on their docking energies and molecular interaction patterns. From the molecules that were considered for WaterMap analysis, seven molecules, namely, Mitoxantrone, Labetalol, Acalabrutinib, Sacubitril, Flubendazole, Trazodone, and Niraparib, ascertained the ability to overlap with high-energy hydration sites. Considering many other displaced unfavorable water molecules, only Acalabrutinib, Flubendazole, and Trazodone molecules highlighted their prominence in terms of binding affinity gains through ΔΔG that ranges between 6.44 and 2.59 kcal/mol. Even if Mitoxantrone exhibited the highest docking score and greater interaction strength, it did not comply with the WaterMap and molecular dynamics simulation results. Moreover, detailed MD simulation trajectory analyses suggested that the drug molecules Flubendazole, Niraparib, and Acalabrutinib were highly stable, observed from their RMSD values and consistent interaction pattern with Glu315, Glu345, Leu347, and Asp407 including the hydrophobic interactions maintained in the three replicates. However, the drug molecule Trazodone displayed a loss of crucial interaction with Leu347, which was essential to inhibit the kinase activity of PAK1. The molecular orbital and electrostatic potential analyses elucidated the reactivity and strong complementarity potentials of the drug molecules in the binding pocket of PAK1. Therefore, the CADD-based reposition efforts, reported in this work, helped in the successful identification of new PAK1 inhibitors that requires further investigation by in vitro analysis.
Collapse
Affiliation(s)
- Jayashree Biswal
- Structural
Biology and Bio-Computing Laboratory, Department of Bioinformatics,
Science Block, Alagappa University, Karaikudi 630 004, Tamil Nadu, India
| | - Prajisha Jayaprakash
- Structural
Biology and Bio-Computing Laboratory, Department of Bioinformatics,
Science Block, Alagappa University, Karaikudi 630 004, Tamil Nadu, India
| | - Suresh Kumar Rayala
- Department
of Biotechnology, Indian Institute of Technology
Madras, Room No. BT 306, Chennai 600 036, Tamil Nadu, India
| | - Ganesh Venkatraman
- Department
of Human Genetics, College of Biomedical Sciences, Sri Ramachandra University, Porur, Chennai 600 116, Tamil Nadu, India
| | - Raghu Rangaswamy
- Structural
Biology and Bio-Computing Laboratory, Department of Bioinformatics,
Science Block, Alagappa University, Karaikudi 630 004, Tamil Nadu, India
| | - Jeyakanthan Jeyaraman
- Structural
Biology and Bio-Computing Laboratory, Department of Bioinformatics,
Science Block, Alagappa University, Karaikudi 630 004, Tamil Nadu, India
| |
Collapse
|
11
|
Xue H, Svatek H, Bertonha AF, Reisenauer K, Robinson J, Kim M, Ingros A, Ho M, Taube J, Romo D. Synthesis of agelastatin A and derivatives premised on a hidden symmetry element leading to analogs displaying anticancer activity. Tetrahedron 2021; 94. [DOI: 10.1016/j.tet.2021.132340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
12
|
Smola M, Gutten O, Dejmek M, Kožíšek M, Evangelidis T, Tehrani ZA, Novotná B, Nencka R, Birkuš G, Rulíšek L, Boura E. Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein. Angew Chem Int Ed Engl 2021; 60:10172-10178. [PMID: 33616279 PMCID: PMC8251555 DOI: 10.1002/anie.202016805] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Indexed: 12/19/2022]
Abstract
STING (stimulator of interferon genes) is a key regulator of innate immunity that has recently been recognized as a promising drug target. STING is activated by cyclic dinucleotides (CDNs) which eventually leads to expression of type I interferons and other cytokines. Factors underlying the affinity of various CDN analogues are poorly understood. Herein, we correlate structural biology, isothermal calorimetry (ITC) and computational modeling to elucidate factors contributing to binding of six CDNs-three pairs of natural (ribo) and fluorinated (2'-fluororibo) 3',3'-CDNs. X-ray structural analyses of six {STING:CDN} complexes did not offer any explanation for the different affinities of the studied ligands. ITC showed entropy/enthalpy compensation up to 25 kcal mol-1 for this set of similar ligands. The higher affinities of fluorinated analogues are explained with help of computational methods by smaller loss of entropy upon binding and by smaller strain (free) energy.
Collapse
Affiliation(s)
- Miroslav Smola
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Milan Dejmek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Milan Kožíšek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Thomas Evangelidis
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Zahra Aliakbar Tehrani
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Barbora Novotná
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Radim Nencka
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Gabriel Birkuš
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| | - Evzen Boura
- Gilead Sciences Research Centre at IOCBInstitute of Organic Chemistry and Biochemistry of the Czech Academy of SciencesFlemingovo náměstí 216610PragueCzech Republic
| |
Collapse
|
13
|
Smola M, Gutten O, Dejmek M, Kožíšek M, Evangelidis T, Tehrani ZA, Novotná B, Nencka R, Birkuš G, Rulíšek L, Boura E. Ligand Strain and Its Conformational Complexity Is a Major Factor in the Binding of Cyclic Dinucleotides to STING Protein. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202016805] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Miroslav Smola
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Milan Dejmek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Milan Kožíšek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Thomas Evangelidis
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Zahra Aliakbar Tehrani
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Barbora Novotná
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Radim Nencka
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Gabriel Birkuš
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Evzen Boura
- Gilead Sciences Research Centre at IOCB Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| |
Collapse
|
14
|
Vavřina Z, Gutten O, Smola M, Zavřel M, Aliakbar Tehrani Z, Charvát V, Kožíšek M, Boura E, Birkuš G, Rulíšek L. Protein-Ligand Interactions in the STING Binding Site Probed by Rationally Designed Single-Point Mutations: Experiment and Theory. Biochemistry 2021; 60:607-620. [PMID: 33586948 DOI: 10.1021/acs.biochem.0c00949] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
STING protein (stimulator of interferon genes) plays an important role in the innate immune system. A number of potent compounds regulating its activity have been reported, mostly derivatives of cyclic dinucleotides (CDNs), natural STING agonists. Here, we aim to provide complementary information to large-scale "ligand-profiling" studies by probing the importance of STING-CDN protein-ligand interactions on the protein side. We examined in detail six typical CDNs each in complex with 13 rationally devised mutations in STING: S162A, S162T, Y167F, G230A, R232K, R232H, A233L, A233I, R238K, T263A, T263S, R293Q, and G230A/R293Q. The mutations switch on and off various types of protein-ligand interactions: π-π stacking, hydrogen bonding, ionic pairing, and nonpolar contacts. We correlated experimental data obtained by differential scanning fluorimetry, X-ray crystallography, and isothermal titration calorimetry with theoretical calculations. This enabled us to provide a mechanistic interpretation of the differences in the binding of representative CDNs to STING. We observed that the G230A mutation increased the thermal stability of the protein-ligand complex, indicating an increased level of ligand binding, whereas R238K and Y167F led to a complete loss of stabilization (ligand binding). The effects of the other mutations depended on the type of ligand (CDN) and varied, to some extent. A very good correlation (R2 = 0.6) between the experimental binding affinities and interaction energies computed by quantum chemical methods enabled us to explain the effect of the studied mutations in detail and evaluate specific interactions quantitatively. Our work may inspire development of high-affinity ligands against the common STING haplotypes by targeting the key (sometimes non-intuitive) protein-ligand interactions.
Collapse
Affiliation(s)
- Zdeněk Vavřina
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic.,Department of Biochemistry, Faculty of Science, Charles University, Prague 110 00, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Miroslav Smola
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Martin Zavřel
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Zahra Aliakbar Tehrani
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic.,Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Průmyslová 595, Vestec 252 50, Czech Republic
| | - Vilém Charvát
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Milan Kožíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Evzen Boura
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Gabriel Birkuš
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, Prague 16610, Czech Republic
| |
Collapse
|
15
|
Fedorenko M, Alesio J, Fedorenko A, Slitt A, Bothun GD. Dominant entropic binding of perfluoroalkyl substances (PFASs) to albumin protein revealed by 19F NMR. CHEMOSPHERE 2021; 263:128083. [PMID: 33297081 PMCID: PMC8479757 DOI: 10.1016/j.chemosphere.2020.128083] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/29/2020] [Revised: 08/18/2020] [Accepted: 08/19/2020] [Indexed: 05/05/2023]
Abstract
Mechanistic insight into protein binding by poly- and perfluoroalkyl substances (PFASs) is critical to understanding how PFASs distribute and accumulate within the body and to developing predictive models within and across classes of PFASs. Fluorine nuclear magnetic resonance spectroscopy (19F NMR) has proven to be a powerful, yet underutilized tool to study PFAS binding; chemical shifts of each fluorine group reflect the local environment along the length of the PFAS molecule. Using bovine serum albumin (BSA), we report dissociation constants, Kd, for four common PFASs well below reported critical micelle concentrations (CMCs) - perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), perfluorohexanesulfonic acid (PFHxS), and perfluorooctanesulfonic acid (PFOS) - as a function of temperature in phosphate buffered saline. Kd values were determined based on the difluoroethyl group adjacent to the anionic headgroups and the terminal trifluoromethyl groups. Our results indicate that the hydrophobic tails exhibit greater binding affinity relative to the headgroup, and that the binding affinities are generally consistent with previous results showing that greater PFAS hydrophobicity leads to greater protein binding. However, the binding mechanism was dominated by entropic hydrophobic interactions attributed to desolvation of the PFAS tails within the hydrophobic cavities of the protein and on the surface of the protein. In addition, PFNA appears to form hemimicelles on the protein surfaces below reported CMC values. This work provides a renewed approach to utilizing 19F NMR for PFAS-protein binding studies and a new perspective on the role of solvent entropy.
Collapse
Affiliation(s)
- Michael Fedorenko
- Department of Chemical Engineering, University of Rhode Island, 2 East Alumni Ave, Kingston, RI, 02881, USA
| | - Jessica Alesio
- Department of Chemical Engineering, University of Rhode Island, 2 East Alumni Ave, Kingston, RI, 02881, USA
| | - Anatoliy Fedorenko
- Department of Chemical Engineering, University of Rhode Island, 2 East Alumni Ave, Kingston, RI, 02881, USA
| | - Angela Slitt
- Department of Biomedical & Pharmaceutical Sciences, University of Rhode Island, 7 Greenhouse Rd, Kingston, RI, 02881, USA
| | - Geoffrey D Bothun
- Department of Chemical Engineering, University of Rhode Island, 2 East Alumni Ave, Kingston, RI, 02881, USA.
| |
Collapse
|
16
|
Vázquez J, López M, Gibert E, Herrero E, Luque FJ. Merging Ligand-Based and Structure-Based Methods in Drug Discovery: An Overview of Combined Virtual Screening Approaches. Molecules 2020; 25:E4723. [PMID: 33076254 PMCID: PMC7587536 DOI: 10.3390/molecules25204723] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Revised: 10/06/2020] [Accepted: 10/11/2020] [Indexed: 12/20/2022] Open
Abstract
Virtual screening (VS) is an outstanding cornerstone in the drug discovery pipeline. A variety of computational approaches, which are generally classified as ligand-based (LB) and structure-based (SB) techniques, exploit key structural and physicochemical properties of ligands and targets to enable the screening of virtual libraries in the search of active compounds. Though LB and SB methods have found widespread application in the discovery of novel drug-like candidates, their complementary natures have stimulated continued efforts toward the development of hybrid strategies that combine LB and SB techniques, integrating them in a holistic computational framework that exploits the available information of both ligand and target to enhance the success of drug discovery projects. In this review, we analyze the main strategies and concepts that have emerged in the last years for defining hybrid LB + SB computational schemes in VS studies. Particularly, attention is focused on the combination of molecular similarity and docking, illustrating them with selected applications taken from the literature.
Collapse
Affiliation(s)
- Javier Vázquez
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona, Av. Prat de la Riba 171, E-08921 Santa Coloma de Gramanet, Spain
| | - Manel López
- AB Science, Parc Scientifique de Luminy, Zone Luminy Enterprise, Case 922, 163 Av. de Luminy, 13288 Marseille, France;
| | - Enric Gibert
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
| | - Enric Herrero
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
| | - F. Javier Luque
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona, Av. Prat de la Riba 171, E-08921 Santa Coloma de Gramanet, Spain
| |
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]
|