1
|
Wang X, Chen X, Chen Z, Xu W, Lai R, Qiu X, Zeng Z, Wang C, Wang Z, Wang J. Integrated Anchored Stapling and Hierarchical Dynamics: MSICDA-Driven CREBBP Bromodomain Inhibition. J Chem Inf Model 2024; 64:4739-4758. [PMID: 38863138 DOI: 10.1021/acs.jcim.4c00381] [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: 06/13/2024]
Abstract
Despite recent success in the computational approaches of cyclic peptide design, current studies face challenges in modeling noncanonical amino acids and nonstandard cyclizations due to limited data. To address this challenge, we developed an integrated framework for the tailored design of stapled peptides (SPs) targeting the bromodomain of CREBBP (CREBBP-BrD). We introduce a powerful combination of anchored stapling and hierarchical molecular dynamics to design and optimize SPs by employing the MultiScale integrative conformational dynamics assessment (MSICDA) strategy, which involves an initial virtual screening of over 1.5 million SPs, followed by comprehensive simulations amounting to 154.54 μs across 5418 of instances. The MSICDA method provides a detailed and holistic stability view of peptide-protein interactions, systematically isolated optimized peptides and identified two leading candidates, DA#430 and DA#99409, characterized by their enhanced stability, optimized binding, and high affinity toward the CREBBP-BrD. In cell-free assays, DA#430 and DA#99409 exhibited 2- to 12-fold greater potency than inhibitor SGC-CBP30. Cell studies revealed higher peptide selectivity for cancerous versus normal cells over small molecules. DA#430 combined with (+)-JQ-1 showed promising synergistic effects. Our approach enables the identification of peptides with optimized binding, high affinity, and enhanced stability, leading to more precise and effective cyclic peptide design, thereby establishing MSICDA as a generalizable and transformative tool for uncovering novel targeted drug development in various therapeutic areas.
Collapse
Affiliation(s)
- Xinpei Wang
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| | - Xu Chen
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| | - Zhidong Chen
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| | - Wanting Xu
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| | - Ruizhi Lai
- Department of Pathology, The Eighth Affiliated Hospital, Sun Yat-sen University, Shenzhen 518033, China
| | - Xiaohui Qiu
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| | - Zekai Zeng
- Department of Pathology, The Eighth Affiliated Hospital, Sun Yat-sen University, Shenzhen 518033, China
| | - Chenglin Wang
- Shenzhen Qiyu Biotechnology Co., Ltd, Shenzhen 518107, China
| | - Zhe Wang
- Department of Pathology, The Eighth Affiliated Hospital, Sun Yat-sen University, Shenzhen 518033, China
| | - Junqing Wang
- School of Pharmaceutical Sciences, Shenzhen Campus of Sun Yat-sen University, Shenzhen 518107, China
| |
Collapse
|
2
|
Li L, Liu S, Wang B, Liu F, Xu S, Li P, Chen Y. An Updated Review on Developing Small Molecule Kinase Inhibitors Using Computer-Aided Drug Design Approaches. Int J Mol Sci 2023; 24:13953. [PMID: 37762253 PMCID: PMC10530957 DOI: 10.3390/ijms241813953] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 08/31/2023] [Accepted: 09/07/2023] [Indexed: 09/29/2023] Open
Abstract
Small molecule kinase inhibitors (SMKIs) are of heightened interest in the field of drug research and development. There are 79 (as of July 2023) small molecule kinase inhibitors that have been approved by the FDA and hundreds of kinase inhibitor candidates in clinical trials that have shed light on the treatment of some major diseases. As an important strategy in drug design, computer-aided drug design (CADD) plays an indispensable role in the discovery of SMKIs. CADD methods such as docking, molecular dynamic, quantum mechanics/molecular mechanics, pharmacophore, virtual screening, and quantitative structure-activity relationship have been applied to the design and optimization of small molecule kinase inhibitors. In this review, we provide an overview of recent advances in CADD and SMKIs and the application of CADD in the discovery of SMKIs.
Collapse
Affiliation(s)
- Linwei Li
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| | - Songtao Liu
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
- Key Laboratory of Pesticide, College of Plant Protection, Nanjing Agricultural University, Nanjing 210095, China
| | - Bi Wang
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| | - Fei Liu
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| | - Shu Xu
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| | - Pirui Li
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| | - Yu Chen
- Jiangsu Key Laboratory for the Research and Utilization of Plant Resources, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China; (L.L.); (S.L.); (B.W.); (F.L.); (S.X.)
- Jiangsu Province Engineering Research Center of Eco-Cultivation and High-Value Utilization of Chines Medicinal Materials, Institute of Botany, Jiangsu Province and Chinese Academy of Sciences, Nanjing 210014, China
| |
Collapse
|
3
|
Salutari I, Caflisch A. Dynamics of the Histone Acetyltransferase Lysine-Rich Loop in the Catalytic Core of the CREB-Binding Protein. J Chem Inf Model 2022; 62:1014-1024. [PMID: 35119862 DOI: 10.1021/acs.jcim.1c01423] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The tight control of transcriptional coactivators is a fundamental aspect of gene expression in cells. The regulation of the CREB-binding protein (CBP) and p300 coactivators, two paralog multidomain proteins, involves an autoinhibitory loop (AIL) of the histone acetyltransferase (HAT) domain. There is experimental evidence for the AIL engaging with the HAT binding site, thus interrupting the acetylation of histone tails or other proteins. Both CBP and p300 contain a domain of about 110 residues (called the bromodomain) that recognizes histone tails with one or more acetylated lysine side chains. Here, we investigate by molecular dynamics simulations whether the AIL of CBP (residues 1556-1618) acetylated at the side chain of Lys1595 can bind to the bromodomain. The structural instability and fast unbinding kinetics of the AIL from the bromodomain pocket suggest that the AIL is not a ligand of the bromodomain on the same protein chain. This is further supported by the absence of strong and persistent contacts at the binding interface. Furthermore, the simulations of unbinding show an initial fast detachment of the acetylated lysine and a slower phase necessary for complete AIL dissociation. We provide further evidence for the instability of the AIL intramolecular binding by comparison with a natural ligand, the histone peptide H3K56ac, which shows higher stability in the pocket.
Collapse
Affiliation(s)
- Ilaria Salutari
- Department of Biochemistry, University of Zürich, CH-8057 Zürich, Switzerland
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zürich, CH-8057 Zürich, Switzerland
| |
Collapse
|
4
|
Rodríguez Y, Gerona-Navarro G, Osman R, Zhou MM. In silico design and molecular basis for the selectivity of Olinone toward the first over the second bromodomain of BRD4. Proteins 2020; 88:414-430. [PMID: 31587361 PMCID: PMC6982606 DOI: 10.1002/prot.25818] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Revised: 09/12/2019] [Accepted: 09/17/2019] [Indexed: 01/11/2023]
Abstract
Bromodomains (BrDs), a conserved structural module in chromatin-associated proteins, are well known for recognizing ε-N-acetyl lysine residues on histones. One of the most relevant BrDs is BRD4, a tandem BrD containing protein (BrD1 and BrD2) that plays a critical role in numerous diseases including cancer. Growing evidence shows that the two BrDs of BRD4 have different biological functions; hence selective ligands that can be used to study their functions are of great interest. Here, as a follow-up of our previous work, we first provide a detailed characterization study of the in silico rational design of Olinone as part of a series of five tetrahydropyrido indole-based compounds as BRD4 BrD1 inhibitors. Additionally, we investigated the molecular basis for Olinone's selective recognition by BrD1 over BrD2. Molecular dynamics simulations, free energy calculations, and conformational analyses of the apo-BRD4-BrD1|2 and BRD4-BrD1|2/Olinone complexes showed that Olinone's selectivity is facilitated by five key residues: Leu92 in BrD1|385 in BrD2 of ZA loop, Asn140|433, Asp144|His437 and Asp145|Glu438 of BC loop, and Ile146|Val49 of helix C. Furthermore, the difference in hydrogen bonds number and in mobility of the ZA and BC loops of the acetyl-lysine binding site between BRD4 BrD1/Olinone and BrD2/Olinone complexes also contribute to the difference in Olinone's binding affinity and selectivity toward BrD1 over BrD2. Altogether, our computer-aided molecular design techniques can effectively guide the development of small-molecule BRD4 BrD1 inhibitors, explain their selectivity origin, and further open doors to the design of new therapeutically improved derivatives.
Collapse
Affiliation(s)
- Yoel Rodríguez
- Department of Natural Sciences, Hostos Community
College of CUNY, Bronx, NY 10451, USA,Department of Pharmacological Sciences, Icahn School
of Medicine at Mount Sinai, 1425 Madison Avenue, New York, NY 10029, USA,Corresponding Authors: Yoel Rodríguez.
Address: Department of Natural Sciences, Room A-507F, Hostos Community College
of CUNY, Bronx, NY 10451, USA. Phone: +1 (718) 518-4134, Fax: +1 (718) 518-1120.
- ; Ming-Ming Zhou. Address: Department
of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, 1425
Madison Avenue, Box 1677, New York, NY 10029, USA. Phone: +1 (212) 659-8652.
Fax: +1 (212) 849-2456.
| | - Guillermo Gerona-Navarro
- Department of Chemistry, Brooklyn College, 2900
Bedford Avenue, Room 351 NE, Brooklyn, NY 11210, USA,Ph.D. Program in Chemistry. The Graduate Center of
The City University of New York, NY 10016, USA
| | - Roman Osman
- Department of Pharmacological Sciences, Icahn School
of Medicine at Mount Sinai, 1425 Madison Avenue, New York, NY 10029, USA
| | - Ming-Ming Zhou
- Department of Pharmacological Sciences, Icahn School
of Medicine at Mount Sinai, 1425 Madison Avenue, New York, NY 10029, USA,Corresponding Authors: Yoel Rodríguez.
Address: Department of Natural Sciences, Room A-507F, Hostos Community College
of CUNY, Bronx, NY 10451, USA. Phone: +1 (718) 518-4134, Fax: +1 (718) 518-1120.
- ; Ming-Ming Zhou. Address: Department
of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, 1425
Madison Avenue, Box 1677, New York, NY 10029, USA. Phone: +1 (212) 659-8652.
Fax: +1 (212) 849-2456.
| |
Collapse
|
5
|
Cheminformatics Explorations of Natural Products. PROGRESS IN THE CHEMISTRY OF ORGANIC NATURAL PRODUCTS 2019; 110:1-35. [PMID: 31621009 DOI: 10.1007/978-3-030-14632-0_1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The chemistry of natural products is fascinating and has continuously attracted the attention of the scientific community for many reasons including, but not limited to, biosynthesis pathways, chemical diversity, the source of bioactive compounds and their marked impact on drug discovery. There is a broad range of experimental and computational techniques (molecular modeling and cheminformatics) that have evolved over the years and have assisted the investigation of natural products. Herein, we discuss cheminformatics strategies to explore the chemistry and applications of natural products. Since the potential synergisms between cheminformatics and natural products are vast, we will focus on three major aspects: (1) exploration of the chemical space of natural products to identify bioactive compounds, with emphasis on drug discovery; (2) assessment of the toxicity profile of natural products; and (3) diversity analysis of natural product collections and the design of chemical collections inspired by natural sources.
Collapse
|
6
|
Bacci M, Langini C, Vymětal J, Caflisch A, Vitalis A. Focused conformational sampling in proteins. J Chem Phys 2018; 147:195102. [PMID: 29166086 DOI: 10.1063/1.4996879] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.
Collapse
Affiliation(s)
- Marco Bacci
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Cassiano Langini
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Jiří Vymětal
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Andreas Vitalis
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
7
|
Polêto MD, Rusu VH, Grisci BI, Dorn M, Lins RD, Verli H. Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities. Front Pharmacol 2018; 9:395. [PMID: 29740321 PMCID: PMC5928326 DOI: 10.3389/fphar.2018.00395] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2017] [Accepted: 04/05/2018] [Indexed: 11/13/2022] Open
Abstract
The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.
Collapse
Affiliation(s)
- Marcelo D Polêto
- Grupo de Bioinformática Estrutural, Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Victor H Rusu
- Swiss National Supercomputing Centre, Lugano, Switzerland
| | - Bruno I Grisci
- Instituto de Informática, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Marcio Dorn
- Instituto de Informática, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Roberto D Lins
- Instituto Aggeu Magalhães, Fundação Oswaldo Cruz, Recife, Brazil
| | - Hugo Verli
- Grupo de Bioinformática Estrutural, Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| |
Collapse
|
8
|
Batiste L, Unzue A, Dolbois A, Hassler F, Wang X, Deerain N, Zhu J, Spiliotopoulos D, Nevado C, Caflisch A. Chemical Space Expansion of Bromodomain Ligands Guided by in Silico Virtual Couplings (AutoCouple). ACS CENTRAL SCIENCE 2018; 4:180-188. [PMID: 29532017 PMCID: PMC5833004 DOI: 10.1021/acscentsci.7b00401] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Indexed: 10/24/2023]
Abstract
Expanding the chemical space and simultaneously ensuring synthetic accessibility is of upmost importance, not only for the discovery of effective binders for novel protein classes but, more importantly, for the development of compounds against hard-to-drug proteins. Here, we present AutoCouple, a de novo approach to computational ligand design focused on the diversity-oriented generation of chemical entities via virtual couplings. In a benchmark application, chemically diverse compounds with low-nanomolar potency for the CBP bromodomain and high selectivity against the BRD4(1) bromodomain were achieved by the synthesis of about 50 derivatives of the original fragment. The binding mode was confirmed by X-ray crystallography, target engagement in cells was demonstrated, and antiproliferative activity was showcased in three cancer cell lines. These results reveal AutoCouple as a useful in silico coupling method to expand the chemical space in hit optimization campaigns resulting in potent, selective, and cell permeable bromodomain ligands.
Collapse
Affiliation(s)
- Laurent Batiste
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Andrea Unzue
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Aymeric Dolbois
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Fabrice Hassler
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Xuan Wang
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Nicholas Deerain
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Jian Zhu
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Dimitrios Spiliotopoulos
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Cristina Nevado
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Amedeo Caflisch
- Department
of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| |
Collapse
|
9
|
Fronik P, Gaiser BI, Sejer Pedersen D. Bitopic Ligands and Metastable Binding Sites: Opportunities for G Protein-Coupled Receptor (GPCR) Medicinal Chemistry. J Med Chem 2017; 60:4126-4134. [DOI: 10.1021/acs.jmedchem.6b01601] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Philipp Fronik
- Department of Drug Design
and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Birgit I. Gaiser
- Department of Drug Design
and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| | - Daniel Sejer Pedersen
- Department of Drug Design
and Pharmacology, Faculty of Health and Medical Sciences, University of Copenhagen, Jagtvej 162, 2100 Copenhagen, Denmark
| |
Collapse
|
10
|
Spiliotopoulos D, Kastritis PL, Melquiond ASJ, Bonvin AMJJ, Musco G, Rocchia W, Spitaleri A. dMM-PBSA: A New HADDOCK Scoring Function for Protein-Peptide Docking. Front Mol Biosci 2016; 3:46. [PMID: 27630991 PMCID: PMC5006095 DOI: 10.3389/fmolb.2016.00046] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2016] [Accepted: 08/16/2016] [Indexed: 11/13/2022] Open
Abstract
Molecular-docking programs coupled with suitable scoring functions are now established and very useful tools enabling computational chemists to rapidly screen large chemical databases and thereby to identify promising candidate compounds for further experimental processing. In a broader scenario, predicting binding affinity is one of the most critical and challenging components of computer-aided structure-based drug design. The development of a molecular docking scoring function which in principle could combine both features, namely ranking putative poses and predicting complex affinity, would be of paramount importance. Here, we systematically investigated the performance of the MM-PBSA approach, using two different Poisson-Boltzmann solvers (APBS and DelPhi), in the currently rising field of protein-peptide interactions (PPIs), identifying the correct binding conformations of 19 different protein-peptide complexes and predicting their binding free energies. First, we scored the decoy structures from HADDOCK calculation via the MM-PBSA approach in order to assess the capability of retrieving near-native poses in the best-scoring clusters and of evaluating the corresponding free energies of binding. MM-PBSA behaves well in finding the poses corresponding to the lowest binding free energy, however the built-in HADDOCK score shows a better performance. In order to improve the MM-PBSA-based scoring function, we dampened the MM-PBSA solvation and coulombic terms by 0.2, as proposed in the HADDOCK score and LIE approaches. The new dampened MM-PBSA (dMM-PBSA) outperforms the original MM-PBSA and ranks the decoys structures as the HADDOCK score does. Second, we found a good correlation between the dMM-PBSA and HADDOCK scores for the near-native clusters of each system and the experimental binding energies, respectively. Therefore, we propose a new scoring function, dMM-PBSA, to be used together with the built-in HADDOCK score in the context of protein-peptide docking simulations.
Collapse
Affiliation(s)
| | - Panagiotis L Kastritis
- Faculty of Science - Chemistry, Bijvoet Center, Utrecht UniversityUtrecht, Netherlands; European Molecular Biology Laboratory HeidelbergHeidelberg, Germany
| | - Adrien S J Melquiond
- Faculty of Science - Chemistry, Bijvoet Center, Utrecht University Utrecht, Netherlands
| | | | - Giovanna Musco
- Biomolecular Nuclear Magnetic Resonance Unit, Ospedale S. Raffaele Milan, Italy
| | - Walter Rocchia
- CONCEPT Lab, Istituto Italiano di Tecnologia Genoa, Italy
| | | |
Collapse
|
11
|
Mechanisms of histone lysine-modifying enzymes: A computational perspective on the role of the protein environment. J Mol Graph Model 2016; 67:69-84. [DOI: 10.1016/j.jmgm.2016.04.011] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2016] [Revised: 04/28/2016] [Accepted: 04/29/2016] [Indexed: 12/13/2022]
|
12
|
Unzue A, Xu M, Dong J, Wiedmer L, Spiliotopoulos D, Caflisch A, Nevado C. Fragment-Based Design of Selective Nanomolar Ligands of the CREBBP Bromodomain. J Med Chem 2015; 59:1350-6. [PMID: 26043365 DOI: 10.1021/acs.jmedchem.5b00172] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Novel ligands of the CREBBP bromodomain were identified by fragment-based docking. The in silico discovered hits have been optimized by chemical synthesis into selective nanomolar compounds, thereby preserving the ligand efficiency. The selectivity for the CREBBP bromodomain over other human bromodomain subfamilies has achieved by a benzoate moiety which was predicted by docking to be involved in favorable electrostatic interactions with the Arg1173 side chain, a prediction that could be verified a posteriori by the high-resolution crystal structure of the CREBBP bromodomain in complex with ligand 6 and also by MD simulations (see Xu, M.; Unzue, A.; Dong, J.; Spiliotopoulos, D.; Nevado, C.; Caflisch, A. Discovery of CREBBP bromodomain inhibitors by high-throughput docking and hit optimization guided by molecular dynamics. J. Med. Chem. 2015, DOI: 10.1021/acs.jmedchem.5b00171).
Collapse
Affiliation(s)
- Andrea Unzue
- Department of of Chemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Min Xu
- Department of of Biochemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Jing Dong
- Department of of Biochemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Lars Wiedmer
- Department of of Biochemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Dimitrios Spiliotopoulos
- Department of of Biochemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Amedeo Caflisch
- Department of of Biochemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| | - Cristina Nevado
- Department of of Chemistry, University of Zürich , Winterthurerstrasse 190, CH-8057, Zürich, Switzerland
| |
Collapse
|
13
|
Observed bromodomain flexibility reveals histone peptide- and small molecule ligand-compatible forms of ATAD2. Biochem J 2015; 466:337-46. [PMID: 25486442 DOI: 10.1042/bj20140933] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Preventing histone recognition by bromodomains emerges as an attractive therapeutic approach in cancer. Overexpression of ATAD2 (ATPase family AAA domain-containing 2 isoform A) in cancer cells is associated with poor prognosis making the bromodomain of ATAD2 a promising epigenetic therapeutic target. In the development of an in vitro assay and identification of small molecule ligands, we conducted structure-guided studies which revealed a conformationally flexible ATAD2 bromodomain. Structural studies on apo-, peptide-and small molecule-ATAD2 complexes (by co-crystallization) revealed that the bromodomain adopts a 'closed', histone-compatible conformation and a more 'open' ligand-compatible conformation of the binding site respectively. An unexpected conformational change of the conserved asparagine residue plays an important role in driving the peptide-binding conformation remodelling. We also identified dimethylisoxazole-containing ligands as ATAD2 binders which aided in the validation of the in vitro screen and in the analysis of these conformational studies.
Collapse
|
14
|
Ran T, Zhang Z, Liu K, Lu Y, Li H, Xu J, Xiong X, Zhang Y, Xu A, Lu S, Liu H, Lu T, Chen Y. Insight into the key interactions of bromodomain inhibitors based on molecular docking, interaction fingerprinting, molecular dynamics and binding free energy calculation. MOLECULAR BIOSYSTEMS 2015; 11:1295-304. [DOI: 10.1039/c4mb00723a] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The interaction mechanism of bromodomain inhibitors was investigated using interaction fingerprinting and binding free energy based methods.
Collapse
|
15
|
High-resolution visualisation of the states and pathways sampled in molecular dynamics simulations. Sci Rep 2014; 4:6264. [PMID: 25179558 PMCID: PMC4151098 DOI: 10.1038/srep06264] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2014] [Accepted: 08/15/2014] [Indexed: 11/20/2022] Open
Abstract
We have recently developed a scalable algorithm for ordering the instantaneous observations of a dynamical system evolving continuously in time. Here, we apply the method to long molecular dynamics trajectories. The procedure requires only a pairwise, geometrical distance as input. Suitable annotations of both structural and kinetic nature reveal the free energy basins visited by biomolecules. The profile is supplemented by a trace of the temporal evolution of the system highlighting the sequence of events. We demonstrate that the resultant SAPPHIRE (States And Pathways Projected with HIgh REsolution) plots provide a comprehensive picture of the thermodynamics and kinetics of complex, molecular systems exhibiting dynamics covering a range of time and length scales. Information on pathways connecting states and the level of recurrence are quickly inferred from the visualisation. The considerable advantages of our approach are speed and resolution: the SAPPHIRE plot is scalable to very large data sets and represents every single snapshot. This minimizes the risk of missing states because of overlap or prior coarse-graining of the data.
Collapse
|