1
|
Kokkaliari S, Grauso L, Mangoni A, Seabra G, Paul VJ, Luesch H. Isolation, Structure Elucidation, and Biological Activity of the Selective TACR2 Antagonist Tumonolide and its Aldehyde from a Marine Cyanobacterium. Chemistry 2024:e202401393. [PMID: 39023398 DOI: 10.1002/chem.202401393] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
The macrocyclic tumonolide (1) with enamide functionality and the linear tumonolide aldehyde (2) are new interconverting natural products from a marine cyanobacterium with a peptide-polyketide skeleton, representing a hybrid of apratoxins and palmyrolides or laingolides. The planar structures were established by NMR and mass spectrometry. The relative configuration of the stereogenically-rich apratoxin-like polyketide portion was determined using J-based configuration analysis. The absolute configuration of tumonolide (1) was determined by chiral analysis of the amino acid units and computational methods, followed by NMR chemical shift and ECD spectrum prediction, indicating all-R configuration for the polyketide portion, as in palmyrolide A and contrary to the all-S configuration in apratoxins. Functional screening against a panel of 168 GPCR targets revealed tumonolide (1) as a selective antagonist of TACR2 with an IC50 of 7.0 μM, closely correlating with binding affinity. Molecular docking studies established the binding mode and rationalized the selectivity for TACR2 over TACR1 and TACR3. RNA sequencing upon treatment of HCT116 colorectal cancer cells demonstrated activation of the pulmonary fibrosis idiopathic signaling pathway and the insulin secretion signaling pathway at 20 μM, indicating its potential to modulate these pathways.
Collapse
Affiliation(s)
- Sofia Kokkaliari
- Department of Medicinal Chemistry and Center for Natural Products, Drug Discovery and Development (CNPD3), University of Florida, 1345 Center Drive, Gainesville, Florida 32610, United States
| | - Laura Grauso
- Dipartimento di Agraria, Università degli Studi di Napoli Federico II, 80055 Portici, Napoli, Italy
| | - Alfonso Mangoni
- Dipartimento di Farmacia, Università degli Studi di Napoli Federico II, 80131, Napoli, Italy
| | - Gustavo Seabra
- Department of Medicinal Chemistry and Center for Natural Products, Drug Discovery and Development (CNPD3), University of Florida, 1345 Center Drive, Gainesville, Florida 32610, United States
| | - Valerie J Paul
- Smithsonian Marine Station at Ft. Pierce, 701 Seaway Drive, Ft. Pierce, Florida 34949, United States
| | - Hendrik Luesch
- Department of Medicinal Chemistry and Center for Natural Products, Drug Discovery and Development (CNPD3), University of Florida, 1345 Center Drive, Gainesville, Florida 32610, United States
| |
Collapse
|
2
|
Cleves AE, Jain AN, Demeter DA, Buchan ZA, Wilmot J, Hancock EN. From UK-2A to florylpicoxamid: Active learning to identify a mimic of a macrocyclic natural product. J Comput Aided Mol Des 2024; 38:19. [PMID: 38630341 DOI: 10.1007/s10822-024-00555-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Accepted: 02/26/2024] [Indexed: 04/19/2024]
Abstract
Scaffold replacement as part of an optimization process that requires maintenance of potency, desirable biodistribution, metabolic stability, and considerations of synthesis at very large scale is a complex challenge. Here, we consider a set of over 1000 time-stamped compounds, beginning with a macrocyclic natural-product lead and ending with a broad-spectrum crop anti-fungal. We demonstrate the application of the QuanSA 3D-QSAR method employing an active learning procedure that combines two types of molecular selection. The first identifies compounds predicted to be most active of those most likely to be well-covered by the model. The second identifies compounds predicted to be most informative based on exhibiting low predicted activity but showing high 3D similarity to a highly active nearest-neighbor training molecule. Beginning with just 100 compounds, using a deterministic and automatic procedure, five rounds of 20-compound selection and model refinement identifies the binding metabolic form of florylpicoxamid. We show how iterative refinement broadens the domain of applicability of the successive models while also enhancing predictive accuracy. We also demonstrate how a simple method requiring very sparse data can be used to generate relevant ideas for synthetic candidates.
Collapse
Affiliation(s)
- Ann E Cleves
- BioPharmics Division, Optibrium Limited, Cambridge, CB25 9GL, UK.
| | - Ajay N Jain
- BioPharmics Division, Optibrium Limited, Cambridge, CB25 9GL, UK
| | | | | | | | | |
Collapse
|
3
|
Anderluzzi G, Ghitti M, Gasparri AM, Taiè G, Sacchi A, Gori A, Andolfo A, Pozzi F, Musco G, Curnis F, Corti A. A novel aminopeptidase N/CD13 inhibitor selectively targets an endothelial form of CD13 after coupling to proteins. Cell Mol Life Sci 2024; 81:68. [PMID: 38289472 PMCID: PMC10827914 DOI: 10.1007/s00018-023-05102-1] [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: 09/19/2023] [Revised: 11/28/2023] [Accepted: 12/20/2023] [Indexed: 02/01/2024]
Abstract
Aminopeptidase N/CD13, a membrane-bound enzyme upregulated in tumor vasculature and involved in angiogenesis, can be used as a receptor for the targeted delivery of drugs to tumors through ligand-directed targeting approaches. We describe a novel peptide ligand (VGCARRYCS, called "G4") that recognizes CD13 with high affinity and selectivity. Enzymological and computational studies showed that G4 is a competitive inhibitor that binds to the catalytic pocket of CD13 through its N-terminal region. Fusing the peptide C-terminus to tumor necrosis factor-alpha (TNF) or coupling it to a biotin/avidin complex causes loss of binding and inhibitory activity against different forms of CD13, including natural or recombinant ectoenzyme and a membrane form expressed by HL60 promyelocytic leukemia cells (likely due to steric hindrance), but not binding to a membrane form of CD13 expressed by endothelial cells (ECs). Furthermore, G4-TNF systemically administered to tumor-bearing mice exerted anticancer effects through a CD13-targeting mechanism, indicating the presence of a CD13 form in tumor vessels with an accessible binding site. Biochemical studies showed that most CD13 molecules expressed on the surface of ECs are catalytically inactive. Other functional assays showed that these molecules can promote endothelial cell adhesion to plates coated with G4-avidin complexes, suggesting that the endothelial form of CD13 can exert catalytically independent biological functions. In conclusion, ECs express a catalytically inactive form of CD13 characterized by an accessible conformation that can be selectively targeted by G4-protein conjugates. This form of CD13 may represent a specific target receptor for ligand-directed targeted delivery of therapeutics to tumors.
Collapse
Affiliation(s)
- Giulia Anderluzzi
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy
- Università Vita-Salute San Raffaele, Milan, Italy
| | - Michela Ghitti
- Biomolecular NMR Group, Division of Genetics and Cell Biology, IRCCS Ospedale San Raffaele Scientific Institute, Milan, Italy
| | - Anna Maria Gasparri
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy
| | - Giulia Taiè
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy
| | - Angelina Sacchi
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy
| | - Alessandro Gori
- Istituto di Scienze e Tecnologie Chimiche, C.N.R., Milan, Italy
| | - Annapaola Andolfo
- ProMeFa, Proteomics and Metabolomics Facility, Center for Omics Sciences, IRCCS Ospedale San Raffaele Scientific Institute, Milan, Italy
| | - Federica Pozzi
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy
| | - Giovanna Musco
- Biomolecular NMR Group, Division of Genetics and Cell Biology, IRCCS Ospedale San Raffaele Scientific Institute, Milan, Italy
| | - Flavio Curnis
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy.
| | - Angelo Corti
- Tumor Biology and Vascular Targeting Unit, Division of Experimental Oncology, IRCCS San Raffaele Scientific Institute, Via Olgettina 58, 20132, Milan, Italy.
- Università Vita-Salute San Raffaele, Milan, Italy.
| |
Collapse
|
4
|
Tang X, Kokot J, Waibl F, Fernández-Quintero ML, Kamenik AS, Liedl KR. Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity. J Chem Inf Model 2023; 63:7107-7123. [PMID: 37943023 PMCID: PMC10685455 DOI: 10.1021/acs.jcim.3c01123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2023] [Revised: 10/24/2023] [Accepted: 10/24/2023] [Indexed: 11/10/2023]
Abstract
We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.
Collapse
Affiliation(s)
- Xuechen Tang
- Department
of General, Inorganic and Theoretical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria
| | - Janik Kokot
- Department
of General, Inorganic and Theoretical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria
| | - Franz Waibl
- Department
of General, Inorganic and Theoretical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria
- Department
of Chemistry and Applied Biosciences, ETH
Zürich, 8093 Zürich, Switzerland
| | | | - Anna S. Kamenik
- Department
of General, Inorganic and Theoretical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria
| | - Klaus R. Liedl
- Department
of General, Inorganic and Theoretical Chemistry, University of Innsbruck, A-6020 Innsbruck, Austria
| |
Collapse
|
5
|
Yu L, Barros SA, Sun C, Somani S. Cyclic Peptide Linker Design and Optimization by Molecular Dynamics Simulations. J Chem Inf Model 2023; 63:6863-6876. [PMID: 37903231 DOI: 10.1021/acs.jcim.3c01359] [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: 11/01/2023]
Abstract
Cyclic peptides are an emerging therapeutic modality that can target protein-protein interaction sites with high affinity and selectivity. A common medicinal chemistry strategy for the optimization of peptide hits is conformational stabilization through macrocyclization. We present a method based on explicit solvent enhanced sampling molecular dynamics simulations for estimating the impact of varying linker lengths and chemistry on the conformational stability of a peptide. The method is demonstrated on three cyclic peptide series that bind to proteins PCSK9, trypsin, and MDM2 adopting loop, β-sheet, and helical secondary structures. In general, the simulations show greater solution stability of the receptor-bound conformation for the higher-affinity peptides, consistent with the idea that preorganizing a ligand for binding can enhance binding affinity. The impact of the force field and sampling is discussed for one series that does not follow this trend. We have successfully applied this method to internal discovery programs to design peptides with increased potency and chemical stability.
Collapse
Affiliation(s)
- Lei Yu
- Janssen Research & Development, LLC, Spring House, Pennsylvania 19477, United States
| | - Stephanie A Barros
- Janssen Research & Development, LLC, Spring House, Pennsylvania 19477, United States
| | - Chengzao Sun
- Janssen Research & Development, LLC, Spring House, Pennsylvania 19477, United States
| | - Sandeep Somani
- Janssen Research & Development, LLC, Spring House, Pennsylvania 19477, United States
| |
Collapse
|
6
|
Jain AN, Brueckner AC, Jorge C, Cleves AE, Khandelwal P, Cortes JC, Mueller L. Complex peptide macrocycle optimization: combining NMR restraints with conformational analysis to guide structure-based and ligand-based design. J Comput Aided Mol Des 2023; 37:519-535. [PMID: 37535171 PMCID: PMC10505130 DOI: 10.1007/s10822-023-00524-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/20/2023] [Indexed: 08/04/2023]
Abstract
Systematic optimization of large macrocyclic peptide ligands is a serious challenge. Here, we describe an approach for lead-optimization using the PD-1/PD-L1 system as a retrospective example of moving from initial lead compound to clinical candidate. We show how conformational restraints can be derived by exploiting NMR data to identify low-energy solution ensembles of a lead compound. Such restraints can be used to focus conformational search for analogs in order to accurately predict bound ligand poses through molecular docking and thereby estimate ligand strain and protein-ligand intermolecular binding energy. We also describe an analogous ligand-based approach that employs molecular similarity optimization to predict bound poses. Both approaches are shown to be effective for prioritizing lead-compound analogs. Surprisingly, relatively small ligand modifications, which may have minimal effects on predicted bound pose or intermolecular interactions, often lead to large changes in estimated strain that have dominating effects on overall binding energy estimates. Effective macrocyclic conformational search is crucial, whether in the context of NMR-based restraints, X-ray ligand refinement, partial torsional restraint for docking/ligand-similarity calculations or agnostic search for nominal global minima. Lead optimization for peptidic macrocycles can be made more productive using a multi-disciplinary approach that combines biophysical data with practical and efficient computational methods.
Collapse
Affiliation(s)
- Ajay N Jain
- Research and Development, BioPharmics LLC, Sonoma County, CA, USA.
| | | | | | - Ann E Cleves
- Research and Development, BioPharmics LLC, Sonoma County, CA, USA
| | | | | | | |
Collapse
|
7
|
Seidel T, Permann C, Wieder O, Kohlbacher SM, Langer T. High-Quality Conformer Generation with CONFORGE: Algorithm and Performance Assessment. J Chem Inf Model 2023; 63:5549-5570. [PMID: 37624145 PMCID: PMC10498443 DOI: 10.1021/acs.jcim.3c00563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Indexed: 08/26/2023]
Abstract
Knowledge of the putative bound-state conformation of a molecule is an essential prerequisite for the successful application of many computer-aided drug design methods that aim to assess or predict its capability to bind to a particular target receptor. An established approach to predict bioactive conformers in the absence of receptor structure information is to sample the low-energy conformational space of the investigated molecules and derive representative conformer ensembles that can be expected to comprise members closely resembling possible bound-state ligand conformations. The high relevance of such conformer generation functionality led to the development of a wide panel of dedicated commercial and open-source software tools throughout the last decades. Several published benchmarking studies have shown that open-source tools usually lag behind their commercial competitors in many key aspects. In this work, we introduce the open-source conformer ensemble generator CONFORGE, which aims at delivering state-of-the-art performance for all types of organic molecules in drug-like chemical space. The ability of CONFORGE and several well-known commercial and open-source conformer ensemble generators to reproduce experimental 3D structures as well as their computational efficiency and robustness has been assessed thoroughly for both typical drug-like molecules and macrocyclic structures. For small molecules, CONFORGE clearly outperformed all other tested open-source conformer generators and performed at least equally well as the evaluated commercial generators in terms of both processing speed and accuracy. In the case of macrocyclic structures, CONFORGE achieved the best average accuracy among all benchmarked generators, with RDKit's generator coming close in second place.
Collapse
Affiliation(s)
- Thomas Seidel
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Christian Permann
- NeGeMac
Research Platform, Department of Pharmaceutical Sciences, Division
of Pharmaceutical Chemistry, University
of Vienna, Josef-Holaubek-Platz
2, 1090 Vienna, Austria
| | - Oliver Wieder
- Christian
Doppler Laboratory for Molecular Informatics in the Biosciences, Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Stefan M. Kohlbacher
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| | - Thierry Langer
- Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
- NeGeMac
Research Platform, Department of Pharmaceutical Sciences, Division
of Pharmaceutical Chemistry, University
of Vienna, Josef-Holaubek-Platz
2, 1090 Vienna, Austria
- Christian
Doppler Laboratory for Molecular Informatics in the Biosciences, Department
of Pharmaceutical Sciences, Division of Pharmaceutical Chemistry, University of Vienna, Josef-Holaubek-Platz 2, 1090 Vienna, Austria
| |
Collapse
|
8
|
Franca TC, Goncalves ADS, Bérubé C, Voyer N, Aubry N, LaPlante SR. Determining the Predominant Conformations of Mortiamides A-D in Solution Using NMR Data and Molecular Modeling Tools. ACS OMEGA 2023; 8:25832-25838. [PMID: 37521620 PMCID: PMC10373451 DOI: 10.1021/acsomega.3c01206] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/22/2023] [Accepted: 06/07/2023] [Indexed: 08/01/2023]
Abstract
Macrocyclic peptidomimetics have been seriously contributing to our arsenal of drugs to combat diseases. The search for nature's discoveries led us to mortiamides A-D (found in a novel fungus from Northern Canada), which is a family of cyclic peptides that clearly have demonstrated impressive pharmaceutical potential. This prompted us to learn more about their solution-state properties as these are central for binding to target molecules. Here, we secured and isolated mortiamide D, and then acquired high-resolution nuclear magnetic resonance (NMR) data to learn more about its structure and dynamics attributes. Sets of two-dimensional NMR experiments provided atomic-level (through-bond and through-space) data to confirm the primary structure, and NMR-driven molecular dynamics (MD) simulations suggested that more than one predominant three-dimensional (3D) structure exist in solution. Further steps of MD simulations are consistent with the finding that the backbones of mortiamides A-C also have at least two prominent macrocyclic shapes, but the side-chain structures and dynamics differed significantly. Knowledge of these solution properties can be exploited for drug design and discovery.
Collapse
Affiliation(s)
- Tanos
C. C. Franca
- INRS
− Centre Armand-Frappier Santé Biotechnologie, Université de Québec, 531 Boulevard des Prairies, Laval, Quebec H7V 1B7, Canada
- Laboratory
of Molecular Modeling Applied to Chemical and Biological Defense, Military Institute of Engineering, 22290-270 Rio de Janeiro, Brazil
- Department
of Chemistry, Faculty of Science, University
of Hradec Králové, Rokitanskeho 62, 50003 Hradec Králové, Czech Republic
| | - Arlan da Silva Goncalves
- Department
of Chemistry, Federal Institute of Espírito
Santo − Unit Vila Velha, 29106-010 Vila Velha, ES, Brazil
- PPGQUI
(Graduate Program in Chemistry), Federal
University of Espírito Santo, Av. Fernando Ferrari, 514,, 29075-910 Vitória, ES, Brazil
| | - Christopher Bérubé
- Departement
de Chimie and PROTEO, Faculté des Sciences et de Génie, Université Laval, 1045 Avenue de la Médecine, Québec, Quebec G1V OA6, Canada
| | - Normand Voyer
- Departement
de Chimie and PROTEO, Faculté des Sciences et de Génie, Université Laval, 1045 Avenue de la Médecine, Québec, Quebec G1V OA6, Canada
| | - Norman Aubry
- NMR
consultant of Steven R. LaPlante’s Lab, INRS − Centre
Armand-Frappier Santé Biotechnologie, Université de Québec, 531 Boulevard des Prairies, Laval, Quebec H7V 1B7, Canada
| | - Steven R. LaPlante
- INRS
− Centre Armand-Frappier Santé Biotechnologie, Université de Québec, 531 Boulevard des Prairies, Laval, Quebec H7V 1B7, Canada
| |
Collapse
|
9
|
Hickey J, Sindhikara D, Zultanski SL, Schultz DM. Beyond 20 in the 21st Century: Prospects and Challenges of Non-canonical Amino Acids in Peptide Drug Discovery. ACS Med Chem Lett 2023; 14:557-565. [PMID: 37197469 PMCID: PMC10184154 DOI: 10.1021/acsmedchemlett.3c00037] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Accepted: 03/29/2023] [Indexed: 05/19/2023] Open
Abstract
Life is constructed primarily using a toolbox of 20 canonical amino acids-relying upon these building blocks for the assembly of proteins and peptides that regulate nearly every cellular task, including cell structure, function, and maintenance. While Nature continues to be a source of inspiration for drug discovery, medicinal chemists are not beholden to only 20 canonical amino acids and have begun to explore non-canonical amino acids (ncAAs) for the construction of designer peptides with improved drug-like properties. However, as our toolbox of ncAAs expands, drug hunters are encountering new challenges in approaching the iterative peptide design-make-test-analyze cycle with a seemingly boundless set of building blocks. This Microperspective focuses on new technologies that are accelerating ncAA interrogation in peptide drug discovery (including HELM notation, late-stage functionalization, and biocatalysis) while shedding light on areas where further investment could not only accelerate the discovery of new medicines but also improve downstream development.
Collapse
Affiliation(s)
- Jennifer
L. Hickey
- Department
of Medicinal Chemistry, Merck & Co.,
Inc., Kenilworth, New Jersey 07033, United States
| | - Dan Sindhikara
- Department
of Modeling and Informatics, Merck &
Co., Inc., Kenilworth, New Jersey 07033, United States
| | - Susan L. Zultanski
- Department
of Process Research & Development, Merck
& Co., Inc., Rahway, New Jersey 07065, United States
| | - Danielle M. Schultz
- Department
of Process Research & Development, Merck
& Co., Inc., Rahway, New Jersey 07065, United States
| |
Collapse
|
10
|
Thompson AA, Harbut MB, Kung PP, Karpowich NK, Branson JD, Grant JC, Hagan D, Pascual HA, Bai G, Zavareh RB, Coate HR, Collins BC, Côte M, Gelin CF, Damm-Ganamet KL, Gholami H, Huff AR, Limon L, Lumb KJ, Mak PA, Nakafuku KM, Price EV, Shih AY, Tootoonchi M, Vellore NA, Wang J, Wei N, Ziff J, Berger SB, Edwards JP, Gardet A, Sun S, Towne JE, Venable JD, Shi Z, Venkatesan H, Rives ML, Sharma S, Shireman BT, Allen SJ. Identification of small-molecule protein-protein interaction inhibitors for NKG2D. Proc Natl Acad Sci U S A 2023; 120:e2216342120. [PMID: 37098070 PMCID: PMC10160951 DOI: 10.1073/pnas.2216342120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2022] [Accepted: 03/10/2023] [Indexed: 04/26/2023] Open
Abstract
NKG2D (natural-killer group 2, member D) is a homodimeric transmembrane receptor that plays an important role in NK, γδ+, and CD8+ T cell-mediated immune responses to environmental stressors such as viral or bacterial infections and oxidative stress. However, aberrant NKG2D signaling has also been associated with chronic inflammatory and autoimmune diseases, and as such NKG2D is thought to be an attractive target for immune intervention. Here, we describe a comprehensive small-molecule hit identification strategy and two distinct series of protein-protein interaction inhibitors of NKG2D. Although the hits are chemically distinct, they share a unique allosteric mechanism of disrupting ligand binding by accessing a cryptic pocket and causing the two monomers of the NKG2D dimer to open apart and twist relative to one another. Leveraging a suite of biochemical and cell-based assays coupled with structure-based drug design, we established tractable structure-activity relationships with one of the chemical series and successfully improved both the potency and physicochemical properties. Together, we demonstrate that it is possible, albeit challenging, to disrupt the interaction between NKG2D and multiple protein ligands with a single molecule through allosteric modulation of the NKG2D receptor dimer/ligand interface.
Collapse
Affiliation(s)
- Aaron A. Thompson
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Michael B. Harbut
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Pei-Pei Kung
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Nathan K. Karpowich
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Jeffrey D. Branson
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Joanna C. Grant
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Deborah Hagan
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Heather A. Pascual
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Guoyun Bai
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | | | - Heather R. Coate
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Bernard C. Collins
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Marjorie Côte
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Christine F. Gelin
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | | | - Hadi Gholami
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Adam R. Huff
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Luis Limon
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Kevin J. Lumb
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Puiying A. Mak
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Kohki M. Nakafuku
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Edmund V. Price
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Amy Y. Shih
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Mandana Tootoonchi
- Discovery Immunology, Janssen Research & Development, San Diego, CA92121
| | - Nadeem A. Vellore
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Jocelyn Wang
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Na Wei
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Jeannie Ziff
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Scott B. Berger
- Discovery Immunology, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - James P. Edwards
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Agnès Gardet
- Discovery Immunology, Janssen Research & Development, San Diego, CA92121
| | - Siquan Sun
- Discovery Immunology, Janssen Research & Development, San Diego, CA92121
| | - Jennifer E. Towne
- Discovery Immunology, Janssen Research & Development, San Diego, CA92121
| | | | - Zhicai Shi
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | | | - Marie-Laure Rives
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Sujata Sharma
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| | - Brock T. Shireman
- Therapeutics Discovery, Janssen Research & Development, San Diego, CA92121
| | - Samantha J. Allen
- Therapeutics Discovery, Janssen Research & Development, Lower Gwynedd Township, PA19002
| |
Collapse
|
11
|
Xia Y, Jin R, Li M, Lan F, Zhu H, Yu Y, Miao D, Wang Q, Zhou Y, Selvaggi G, Ying S, Zhang J, Shen H, Le X, Li W. Potent antitumor activity of ensartinib in MET exon 14 skipping-mutated non-small cell lung cancer. Cancer Lett 2023; 561:216140. [PMID: 36948240 DOI: 10.1016/j.canlet.2023.216140] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Revised: 03/09/2023] [Accepted: 03/19/2023] [Indexed: 03/24/2023]
Abstract
Met proto-oncogene exon 14 skipping (METex14) mutations are targetable driver genes in approximately 3% of non-small-cell lung cancers (NSCLCs). Ensartinib, a type Ia MET inhibitor, is a multi-kinase inhibitor that has been approved for ALK-positive NSCLCs. Ensartinib was administered for compassionate use (cohort 1) and in a phase II clinical trial (cohort 2) to patients with METex14 mutant NSCLCs, with ORR as a primary endpoint. Molecular simulation was conducted to evaluate ensartinib c-MET interaction, and cell lines, patient-derived organoids (PDOs), and xenograft models were used to test the effectiveness of ensartinib. Among 29 evaluable patients, the ORR and DCR of ensartinib were 67% and 94% in cohort 1, and 73% and 91% in cohort 2. The median DoR was 6.8 months and median PFS was 6.1 months in the total population. Rash was the most common drug-related adverse event, and peripheral edema of any grade was reported in only 9% patients. Molecular simulations indicated favorable binding of ensartinib to c-MET. The kinase assay demonstrated an IC50 of 7.9 nM of ensartinib against METex14 protein. In vitro, Hs746T (METex14 mutation) and EBC-1 (MET amplification) cells were sensitive to ensartinib, with IC50 of 31 and 44 nM, respectively. Ensartinib exhibited comparable inhibitory effects on cell migration as crizotinib and tepotinib in both cell types. In vivo, ensartinib suppressed the growth of Hs746T cells. Ensartinib also potently inhibited the viability of PDOs. Overall, Ensartinib exhibited substantial antitumor effects against METex14 mutant NSCLCs in preclinical and clinical trials, with relatively low peripheral edema rates.
Collapse
Affiliation(s)
- Yang Xia
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Rui Jin
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Miao Li
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Fen Lan
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Hao Zhu
- Department of Respiratory and Critical Care Medicine, Wuyi First People's Hospital, Jinhua, Zhejiang, China
| | - Yinghui Yu
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Da Miao
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Qiyuan Wang
- Department of Radiology, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, China
| | - Yi Zhou
- Xcovery Holdings Inc, Palm Beach Gardens, FL, USA
| | | | - Songmin Ying
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Jianjun Zhang
- Department of Thoracic/Head and Neck Medical Oncology, MD Anderson Cancer Center, Houston, TX, USA
| | - Huahao Shen
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| | - Xiuning Le
- Department of Thoracic/Head and Neck Medical Oncology, MD Anderson Cancer Center, Houston, TX, USA.
| | - Wen Li
- Key Laboratory of Respiratory Disease of Zhejiang Province, Department of Respiratory and Critical Care Medicine, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, Zhejiang, China; Cancer Center, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| |
Collapse
|
12
|
Chen Q, Zhou X, Rehmel J, Steele JP, Svensson KA, Beck JP, Hembre EJ, Hao J. Ensemble Docking Approach to Mitigate Pregnane X Receptor-Mediated CYP3A4 Induction Risk. J Chem Inf Model 2023; 63:173-186. [PMID: 36473234 DOI: 10.1021/acs.jcim.2c01175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Three structurally closely related dopamine D1 receptor positive allosteric modulators (D1 PAMs) based on a tetrahydroisoquinoline (THIQ) scaffold were profiled for their CYP3A4 induction potentials. It was found that the length of the linker at the C5 position greatly affected the potentials of these D1 PAMs as CYP3A4 inducers, and the level of induction correlated well with the activation of the pregnane X receptor (PXR). Based on the published PXR X-ray crystal structures, we built a binding model specifically for these THIQ-scaffold-based D1 PAMs in the PXR ligand-binding pocket via an ensemble docking approach and found the model could explain the observed CYP induction disparity. Combined with our previously reported D1 receptor homology model, which identified the C5 position as pointing toward the solvent-exposed space, our PXR-binding model coincidentally suggested that structural modifications at the C5 position could productively modulate the CYP induction potential while maintaining the D1 PAM potency of these THIQ-based PAMs.
Collapse
Affiliation(s)
- Qi Chen
- Discovery Chemistry Research and Technologies, Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana46285, United States
| | - Xin Zhou
- Drug Disposition, Lilly Biotechnology Center, Eli Lilly and Company, 10290 Campus Point Drive, San Diego, California92121, United States
| | - Jessica Rehmel
- Drug Disposition, Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana46285, United States
| | - James P Steele
- Quantitative Biology, Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana46285, United States
| | - Kjell A Svensson
- Neuroscience Discovery, Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana46285, United States
| | - James P Beck
- Discovery Chemistry Research and Technologies, Lilly Biotechnology Center, Eli Lilly and Company, 10290 Campus Point Drive, San Diego, California92121, United States
| | - Erik J Hembre
- Discovery Chemistry Research and Technologies, Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana46285, United States
| | - Junliang Hao
- Discovery Chemistry Research and Technologies, Lilly Biotechnology Center, Eli Lilly and Company, 10290 Campus Point Drive, San Diego, California92121, United States
| |
Collapse
|
13
|
Hsueh SCC, Aina A, Plotkin SS. Ensemble Generation for Linear and Cyclic Peptides Using a Reservoir Replica Exchange Molecular Dynamics Implementation in GROMACS. J Phys Chem B 2022; 126:10384-10399. [PMID: 36410027 DOI: 10.1021/acs.jpcb.2c05470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The profile of shapes presented by a cyclic peptide modulates its therapeutic efficacy and is represented by the ensemble of its sampled conformations. Although some algorithms excel at creating a diverse ensemble of cyclic peptide conformations, they seldom address the entropic contribution of flexible conformations and often have significant practical difficulty producing an ensemble with converged and reliable thermodynamic properties. In this study, an accelerated molecular dynamics (MD) method, namely, reservoir replica exchange MD (R-REMD or Res-REMD), was implemented in GROMACS ver. 4.6.7 and benchmarked on two small cyclic peptide model systems: a cyclized furin cleavage site of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (cyclo-(CGPRRARSG)) and oxytocin (disulfide-bonded CYIQNCPLG). Additionally, we also benchmarked Res-REMD on alanine dipeptide and Trpzip2 to demonstrate its validity and efficiency over REMD. For Trpzip2, Res-REMD coupled with an umbrella-sampling-derived reservoir generated similar folded fractions as regular REMD but on a much faster time scale. For cyclic peptides, Res-REMD appeared to be marginally faster than REMD in ensemble generation. Finally, Res-REMD was more effective in sampling rare events such as trans to cis peptide bond isomerization. We provide a GitHub page with the modified GROMACS source code for running Res-REMD at https://github.com/PlotkinLab/Reservoir-REMD.
Collapse
Affiliation(s)
- Shawn C C Hsueh
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, BCV6T 1Z1, Canada
| | - Adekunle Aina
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, BCV6T 1Z1, Canada
| | - Steven S Plotkin
- Department of Physics and Astronomy, The University of British Columbia, Vancouver, BCV6T 1Z1, Canada.,Genome Science and Technology Program, The University of British Columbia, Vancouver, BCV6T 1Z1, Canada
| |
Collapse
|
14
|
Eastwood JRB, Jiang L, Bonneau R, Kirshenbaum K, Renfrew PD. Evaluating the Conformations and Dynamics of Peptoid Macrocycles. J Phys Chem B 2022; 126:5161-5174. [PMID: 35820178 DOI: 10.1021/acs.jpcb.2c01669] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Peptoid macrocycles are versatile and chemically diverse peptidomimetic oligomers. However, the conformations and dynamics of these macrocycles have not been evaluated comprehensively and require extensive further investigation. Recent studies indicate that two degrees of freedom, and four distinct conformations, adequately describe the behavior of each monomer backbone unit in most peptoid oligomers. On the basis of this insight, we conducted molecular dynamics simulations of model macrocycles using an exhaustive set of idealized possible starting conformations. Simulations of various sizes of peptoid macrocycles yielded a limited set of populated conformations. In addition to reproducing all relevant experimentally determined conformations, the simulations accurately predicted a cyclo-octamer conformation for which we now present the first experimental observation. Sets of three adjacent dihedral angles (ϕi, ψi, ωi+1) exhibited correlated crankshaft motions over the course of simulation for peptoid macrocycles of six residues and larger. These correlated motions may occur in the form of an inversion of one amide bond and the concerted rotation of the preceding ϕ and ψ angles to their mirror-image conformation, a variation on "crankshaft flip" motions studied in polymers and peptides. The energy landscape of these peptoid macrocycles can be described as a network of conformations interconnected by transformations of individual crankshaft flips. For macrocycles of up to eight residues, our mapping of the landscape is essentially complete.
Collapse
Affiliation(s)
- James R B Eastwood
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Linhai Jiang
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Richard Bonneau
- Center for Data Science, New York University, New York, New York 10011, United States.,Center for Computational Biology, Flatiron Institute, New York, New York 10010 United States
| | - Kent Kirshenbaum
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - P Douglas Renfrew
- Center for Computational Biology, Flatiron Institute, New York, New York 10010 United States
| |
Collapse
|
15
|
Meixner M, Zachmann M, Metzler S, Scheerer J, Zacharias M, Antes I. Dynamic Docking of Macrocycles in Bound and Unbound Protein Structures with DynaDock. J Chem Inf Model 2022; 62:3426-3441. [PMID: 35796228 DOI: 10.1021/acs.jcim.2c00436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Macrocycles are interesting molecules with unique features due to their conformationally constrained yet flexible ring structure. This characteristic poses a difficult challenge for computational modeling studies since they rely on accurate structural descriptions. In particular, molecular docking calculations suffer from the lack of ring flexibility during pose generation, which is often compensated by using pregenerated ligand conformer ensembles. Moreover, receptor structures are mainly treated rigidly, which limits the use of many docking tools. In this study, we optimized our previous molecular dynamics-based sampling and docking pipeline specifically designed for the accurate prediction of macrocyclic compounds. We developed a dihedral classification procedure for in-depth conformational analysis of the macrocyclic rings and extracted structural ensembles that were subsequently docked in both bound and unbound protein structures employing a fully flexible approach. Our results suggest that including a ring conformer close to the bound state in the starting ensemble increases the chance of successful docking. The bioactive conformations of a diverse set of ligands could be predicted with high and decent accuracy in bound and unbound protein structures, respectively, due to the incorporation of full molecular flexibility in our approach. The remaining unsuccessful docking calculations were mainly caused by large flexible substituents that bind to surface-exposed binding sites, rather than the macrocyclic ring per se and could be further improved by explicit molecular dynamics simulations of the docked complex.
Collapse
Affiliation(s)
- Maximilian Meixner
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Martin Zachmann
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Sebastian Metzler
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Jonathan Scheerer
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| | - Martin Zacharias
- Center of Functional Protein Assemblies, Technical University Munich, Ernst-Otto-Fischer-Straße 8, Garching bei München 85748, Germany
| | - Iris Antes
- TUM School of Life Sciences, Technical University Munich, Am Staudengarten 2, Freising 85354, Germany
| |
Collapse
|
16
|
Cansfield AD, Ator MA, Banerjee J, Bestwick M, Bortolato A, Brown GA, Brown J, Butkovic K, Cansfield JE, Christopher JA, Congreve M, Cseke G, Deflorian F, Dugan B, Hunjadi MP, Hutinec A, Inturi TK, Landek G, Mason J, O'Brien A, Ott GR, Rupcic R, Saxty G, Southall SM, Zadravec R, Watson SP. Novel Macrocyclic Antagonists of the Calcitonin Gene-Related Peptide Receptor: Design, Realization, and Structural Characterization of Protein-Ligand Complexes. ACS Chem Neurosci 2022; 13:751-765. [PMID: 35245037 DOI: 10.1021/acschemneuro.1c00696] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
A series of macrocyclic calcitonin gene-related peptide (CGRP) receptor antagonists identified using structure-based design principles, exemplified by HTL0028016 (1) and HTL0028125 (2), is described. Structural characterization by X-ray crystallography of the interaction of two of the macrocycle antagonists with the CGRP receptor ectodomain is described, along with structure-activity relationships associated with point changes to the macrocyclic antagonists. The identification of non-peptidic/natural product-derived, macrocyclic ligands for a G protein coupled receptor (GPCR) is noteworthy.
Collapse
Affiliation(s)
- Andrew D Cansfield
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Mark A Ator
- Teva Pharmaceuticals, 145 Brandywine Parkway, West Chester, Pennsylvania 19380, United States
| | - Joydeep Banerjee
- Syngene International, Biocon Park, Bommasandra, Bangalore 560099, India
| | - Michael Bestwick
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Andrea Bortolato
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Giles A Brown
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Jason Brown
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | | | - Julie E Cansfield
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - John A Christopher
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Miles Congreve
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Gabriella Cseke
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Francesca Deflorian
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Benjamin Dugan
- Teva Pharmaceuticals, 145 Brandywine Parkway, West Chester, Pennsylvania 19380, United States
| | | | - Antun Hutinec
- Fidelta d.o.o., Prilaz Baruna Filipovića 29, 10000 Zagreb, Croatia
| | | | - Goran Landek
- Fidelta d.o.o., Prilaz Baruna Filipovića 29, 10000 Zagreb, Croatia
| | - Jonathan Mason
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Alistair O'Brien
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Gregory R Ott
- Teva Pharmaceuticals, 145 Brandywine Parkway, West Chester, Pennsylvania 19380, United States
| | - Renata Rupcic
- Fidelta d.o.o., Prilaz Baruna Filipovića 29, 10000 Zagreb, Croatia
| | - Gordon Saxty
- Fidelta d.o.o., Prilaz Baruna Filipovića 29, 10000 Zagreb, Croatia
| | - Stacey M Southall
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| | - Rahela Zadravec
- Fidelta d.o.o., Prilaz Baruna Filipovića 29, 10000 Zagreb, Croatia
| | - Stephen P Watson
- Sosei Heptares, Steinmetz Building, Granta Park, Great Abington, Cambridge CB21 6DG, U.K
| |
Collapse
|
17
|
Robertson MJ, Meyerowitz JG, Panova O, Borrelli K, Skiniotis G. Plasticity in ligand recognition at somatostatin receptors. Nat Struct Mol Biol 2022; 29:210-217. [PMID: 35210615 PMCID: PMC11073612 DOI: 10.1038/s41594-022-00727-5] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 01/14/2022] [Indexed: 11/10/2022]
Abstract
Somatostatin is a signaling peptide that plays a pivotal role in physiologic processes relating to metabolism and growth through its actions at somatostatin receptors (SSTRs). Members of the SSTR subfamily, particularly SSTR2, are key drug targets for neuroendocrine neoplasms, with synthetic peptide agonists currently in clinical use. Here, we show the cryogenic-electron microscopy structures of active-state SSTR2 in complex with heterotrimeric Gi3 and either the endogenous ligand SST14 or the FDA-approved drug octreotide. Complemented by biochemical assays and molecular dynamics simulations, these structures reveal key details of ligand recognition and receptor activation at SSTRs. We find that SSTR ligand recognition is highly diverse, as demonstrated by ligand-induced conformational changes in ECL2 and substantial sequence divergence across subtypes in extracellular regions. Despite this complexity, we rationalize several known sources of SSTR subtype selectivity and identify an additional interaction for specific binding. These results provide valuable insights for structure-based drug discovery at SSTRs.
Collapse
Affiliation(s)
- Michael J Robertson
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Justin G Meyerowitz
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Anesthesiology, Perioperative and Pain Medicine, Stanford University School of Medicine, Stanford, CA, USA
| | - Ouliana Panova
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA
| | | | - Georgios Skiniotis
- Department of Molecular and Cellular Physiology, Stanford University School of Medicine, Stanford, CA, USA.
- Department of Structural Biology, Stanford University School of Medicine, Stanford, CA, USA.
| |
Collapse
|
18
|
Wang S, Krummenacher K, Landrum GA, Sellers BD, Di Lello P, Robinson SJ, Martin B, Holden JK, Tom JYK, Murthy AC, Popovych N, Riniker S. Incorporating NOE-Derived Distances in Conformer Generation of Cyclic Peptides with Distance Geometry. J Chem Inf Model 2022; 62:472-485. [PMID: 35029985 DOI: 10.1021/acs.jcim.1c01165] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nuclear magnetic resonance (NMR) data from NOESY (nuclear Overhauser enhancement spectroscopy) and ROESY (rotating frame Overhauser enhancement spectroscopy) experiments can easily be combined with distance geometry (DG) based conformer generators by modifying the molecular distance bounds matrix. In this work, we extend the modern DG based conformer generator ETKDG, which has been shown to reproduce experimental crystal structures from small molecules to large macrocycles well, to include NOE-derived interproton distances. In noeETKDG, the experimentally derived interproton distances are incorporated into the distance bounds matrix as loose upper (or lower) bounds to generate large conformer sets. Various subselection techniques can subsequently be applied to yield a conformer bundle that best reproduces the NOE data. The approach is benchmarked using a set of 24 (mostly) cyclic peptides for which NOE-derived distances as well as reference solution structures obtained by other software are available. With respect to other packages currently available, the advantages of noeETKDG are its speed and that no prior force-field parametrization is required, which is especially useful for peptides with unnatural amino acids. The resulting conformer bundles can be further processed with the use of structural refinement techniques to improve the modeling of the intramolecular nonbonded interactions. The noeETKDG code is released as a fully open-source software package available at www.github.com/rinikerlab/customETKDG.
Collapse
Affiliation(s)
- Shuzhe Wang
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Kajo Krummenacher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Gregory A Landrum
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Benjamin D Sellers
- Department of Discovery Chemistry, Genentech, Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Paola Di Lello
- Department of Structural Biology, Genentech, Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Sarah J Robinson
- Department of Discovery Chemistry, Genentech, Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Bryan Martin
- Department of Structural Biology, Genentech, Inc., 1 DNA Way, South San Francisco, California 94080, United States
| | - Jeffrey K Holden
- Department of Early Discovery Biochemistry, Genentech, Inc., South San Francisco, California 94080, United States
| | - Jeffrey Y K Tom
- Department of Early Discovery Biochemistry, Genentech, Inc., South San Francisco, California 94080, United States
| | - Anastasia C Murthy
- Department of Early Discovery Biochemistry, Genentech, Inc., South San Francisco, California 94080, United States
| | - Nataliya Popovych
- Department of Early Discovery Biochemistry, Genentech, Inc., South San Francisco, California 94080, United States
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
19
|
Novotná B, Holá L, Staś M, Gutten O, Smola M, Zavřel M, Vavřina Z, Buděšínský M, Liboska R, Chevrier F, Dobiaš J, Boura E, Rulíšek L, Birkuš G. Enzymatic Synthesis of 3'-5', 3'-5' Cyclic Dinucleotides, Their Binding Properties to the Stimulator of Interferon Genes Adaptor Protein, and Structure/Activity Correlations. Biochemistry 2021; 60:3714-3727. [PMID: 34788017 DOI: 10.1021/acs.biochem.1c00692] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The 3'-5', 3'-5' cyclic dinucleotides (3'3'CDNs) are bacterial second messengers that can also bind to the stimulator of interferon genes (STING) adaptor protein in vertebrates and activate the host innate immunity. Here, we profiled the substrate specificity of four bacterial dinucleotide synthases from Vibrio cholerae (DncV), Bacillus thuringiensis (btDisA), Escherichia coli (dgcZ), and Thermotoga maritima (tDGC) using a library of 33 nucleoside-5'-triphosphate analogues and then employed these enzymes to synthesize 24 3'3'CDNs. The STING affinity of CDNs was evaluated in cell-based and biochemical assays, and their ability to induce cytokines was determined by employing human peripheral blood mononuclear cells. Interestingly, the prepared heterodimeric 3'3'CDNs bound to the STING much better than their homodimeric counterparts and showed similar or better potency than bacterial 3'3'CDNs. We also rationalized the experimental findings by in-depth STING-CDN structure-activity correlations by dissecting computed interaction free energies into a set of well-defined and intuitive terms. To this aim, we employed state-of-the-art methods of computational chemistry, such as quantum mechanics/molecular mechanics (QM/MM) calculations, and complemented the computed results with the {STING:3'3'c-di-ara-AMP} X-ray crystallographic structure. QM/MM identified three outliers (mostly homodimers) for which we have no clear explanation of their impaired binding with respect to their heterodimeric counterparts, whereas the R2 = 0.7 correlation between the computed ΔG'int_rel and experimental ΔTm's for the remaining ligands has been very encouraging.
Collapse
Affiliation(s)
- Barbora Novotná
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic.,Faculty of Science, Charles University, Albertov 6, Prague 12800, Czech Republic
| | - Lucie Holá
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Monika Staś
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Miroslav Smola
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Martin Zavřel
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Zdeněk Vavřina
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic.,Faculty of Science, Charles University, Albertov 6, Prague 12800, Czech Republic
| | - Miloš Buděšínský
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Radek Liboska
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Florian Chevrier
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Juraj Dobiaš
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Evzen Boura
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, 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 and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| | - Gabriel Birkuš
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB, Flemingovo náměstí 2, Prague 16610, Czech Republic
| |
Collapse
|
20
|
Olanders G, Brandt P, Sköld C, Karlén A. Computational studies of molecular pre-organization through macrocyclization: Conformational distribution analysis of closely related non-macrocyclic and macrocyclic analogs. Bioorg Med Chem 2021; 49:116399. [PMID: 34601455 DOI: 10.1016/j.bmc.2021.116399] [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: 05/05/2021] [Revised: 08/24/2021] [Accepted: 09/01/2021] [Indexed: 11/29/2022]
Abstract
Macrocycles form an important compound class in medicinal chemistry due to their interesting structural and biological properties. To help design macrocycles, it is important to understand how the conformational preferences are affected upon macrocyclization of a lead compound. To address this, we collected a unique data set of protein-ligand complexes containing "non-macrocyclic" ("linear") ligands matched with macrocyclic analogs binding to the same protein in a similar pose. Out of the 39 co-crystallized ligands considered, 10 were linear and 29 were macrocyclic. To enable a more general analysis, 128 additional ligands from the publications associated with these protein data bank entries were added to the data set. Using in total 167 collected ligands, we investigated if the conformers in the macrocyclic conformational ensembles were more similar to the bioactive conformation in comparison to the conformers of their linear counterparts. Unexpectedly, in most cases the macrocycle conformational ensemble distributions were not very different from those of the linear compounds. Thus, care should be taken when designing macrocycles with the aim to focus their conformational preference towards the bioactive conformation. We also set out to investigate potential conformational flexibility differences between the two compound classes, computational energy window settings and evaluate a literature metric for approximating the conformational focusing on the bioactive conformation.
Collapse
Affiliation(s)
- Gustav Olanders
- Department of Medicinal Chemistry, Uppsala University, BMC, Box 574, SE-751 23 Uppsala, Sweden
| | - Peter Brandt
- Medicinal Chemistry, Research and Early Development Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Christian Sköld
- Department of Medicinal Chemistry, Uppsala University, BMC, Box 574, SE-751 23 Uppsala, Sweden
| | - Anders Karlén
- Department of Medicinal Chemistry, Uppsala University, BMC, Box 574, SE-751 23 Uppsala, Sweden.
| |
Collapse
|
21
|
Seep L, Bonin A, Meier K, Diedam H, Göller AH. Ensemble completeness in conformer sampling: the case of small macrocycles. J Cheminform 2021; 13:55. [PMID: 34325738 PMCID: PMC8320181 DOI: 10.1186/s13321-021-00524-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 06/05/2021] [Indexed: 11/18/2022] Open
Abstract
In this study we compare the three algorithms for the generation of conformer ensembles Biovia BEST, Schrödinger Prime macrocycle sampling (PMM) and Conformator (CONF) form the University of Hamburg, with ensembles derived for exhaustive molecular dynamics simulations applied to a dataset of 7 small macrocycles in two charge states and three solvents. Ensemble completeness is a prerequisite to allow for the selection of relevant diverse conformers for many applications in computational chemistry. We apply conformation maps using principal component analysis based on ring torsions. Our major finding critical for all applications of conformer ensembles in any computational study is that maps derived from MD with explicit solvent are significantly distinct between macrocycles, charge states and solvents, whereas the maps for post-optimized conformers using implicit solvent models from all generator algorithms are very similar independent of the solvent. We apply three metrics for the quantification of the relative covered ensemble space, namely cluster overlap, variance statistics, and a novel metric, Mahalanobis distance, showing that post-optimized MD ensembles cover a significantly larger conformational space than the generator ensembles, with the ranking PMM > BEST >> CONF. Furthermore, we find that the distributions of 3D polar surface areas are very similar for all macrocycles independent of charge state and solvent, except for the smaller and more strained compound 7, and that there is also no obvious correlation between 3D PSA and intramolecular hydrogen bond count distributions.
Collapse
Affiliation(s)
- Lea Seep
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Anne Bonin
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Katharina Meier
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany
| | - Holger Diedam
- Engineering & Technology, Applied Mathematics, Bayer AG, 51368, Leverkusen, Germany
| | - Andreas H Göller
- Pharmaceuticals R&D, Digital Technologies, Bayer AG, 42096, Wuppertal, Germany.
| |
Collapse
|
22
|
David L, Wenlock M, Barton P, Ritzén A. Prediction of Chameleonic Efficiency. ChemMedChem 2021; 16:2669-2685. [PMID: 34240561 DOI: 10.1002/cmdc.202100306] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 06/29/2021] [Indexed: 11/09/2022]
Abstract
Chameleonic properties, i. e., the capacity of a molecule to hide polarity in non-polar environments and expose it in water, help achieving sufficient permeability and solubility for drug molecules with high MW. We present models of experimental measures of polarity for a set of 24 FDA approved drugs (MW 405-1113) and one PROTAC (MW 1034). Conformational ensembles in aqueous and non-polar environments were generated using molecular dynamics. A linear regression model that predicts chromatographic apparent polarity (EPSA) with a mean unsigned error of 10 Å2 was derived based on separate terms for donor, acceptor, and total molecular SASA. A good correlation (R2 =0.92) with an experimental measure of hydrogen bond donor potential, Δlog Poct-tol , was found for the mean hydrogen bond donor SASA of the conformational ensemble scaled with Abraham's A hydrogen bond acidity. Two quantitative measures of chameleonic behaviour, the chameleonic efficiency indices, are introduced. We envision that the methods presented herein will be useful to triage designed molecules and prioritize those with the best chance of achieving acceptable permeability and solubility.
Collapse
Affiliation(s)
- Laurent David
- Computational Chemistry, H. Lundbeck A/S, Ottiliavej 9, 2300, Valby, Copenhagen, Denmark
| | - Mark Wenlock
- Physical Chemistry, Cyprotex Discovery Limited, Alderley Park, Nether Alderley, Cheshire, SK10 4TG, UK
| | - Patrick Barton
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK.,DMPK, UCB Celltech, Branch of UCB Pharma S.A., 208 Bath Road, Slough, Berkshire, SL1 3WE, UK
| | - Andreas Ritzén
- Drug Design, LEO Pharma A/S, Industriparken 55, 2550, Ballerup, Denmark.,Monte Rosa Therapeutics AG, Aeschenvorstadt 36, CH 4057, Basel, Switzerland
| |
Collapse
|
23
|
Pimková Polidarová M, Břehová P, Kaiser MM, Smola M, Dračínský M, Smith J, Marek A, Dejmek M, Šála M, Gutten O, Rulíšek L, Novotná B, Brázdová A, Janeba Z, Nencka R, Boura E, Páv O, Birkuš G. Synthesis and Biological Evaluation of Phosphoester and Phosphorothioate Prodrugs of STING Agonist 3',3'-c-Di(2'F,2'dAMP). J Med Chem 2021; 64:7596-7616. [PMID: 34019405 DOI: 10.1021/acs.jmedchem.1c00301] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Cyclic dinucleotides (CDNs) are second messengers that bind to the stimulator of interferon genes (STING) and trigger the expression of type I interferons and proinflammatory cytokines. Here we evaluate the activity of 3',3'-c-di(2'F,2'dAMP) and its phosphorothioate analogues against five STING allelic forms in reporter-cell-based assays and rationalize our findings with X-ray crystallography and quantum mechanics/molecular mechanics calculations. We show that the presence of fluorine in the 2' position of 3',3'-c-di(2'F,2'dAMP) improves its activity not only against the wild type (WT) but also against REF and Q STING. Additionally, we describe the synthesis of the acyloxymethyl and isopropyloxycarbonyl phosphoester prodrugs of CDNs. Masking the negative charges of the CDNs results in an up to a 1000-fold improvement of the activities of the prodrugs relative to those of their parent CDNs. Finally, the uptake and intracellular cleavage of pivaloyloxymethyl prodrugs to the parent CDN is rapid, reaching a peak intracellular concentration within 2 h.
Collapse
Affiliation(s)
- Markéta Pimková Polidarová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic.,Faculty of Science, Charles University, Albertov 6, Prague 128 00, Czech Republic
| | - Petra Břehová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Martin Maxmilian Kaiser
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Miroslav Smola
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Joshua Smith
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Aleš Marek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Milan Dejmek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Michal Šála
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Barbora Novotná
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic.,Faculty of Science, Charles University, Albertov 6, Prague 128 00, Czech Republic
| | - Andrea Brázdová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Zlatko Janeba
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Radim Nencka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Evzen Boura
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Ondřej Páv
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| | - Gabriel Birkuš
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo Náměstí 2, Prague 160 00, Czech Republic
| |
Collapse
|
24
|
Zhan W, Zhang H, Ginn J, Leung A, Liu YJ, Michino M, Toita A, Okamoto R, Wong TT, Imaeda T, Hara R, Yukawa T, Chelebieva S, Tumwebaze PK, Lafuente-Monasterio MJ, Martinez-Martinez MS, Vendome J, Beuming T, Sato K, Aso K, Rosenthal PJ, Cooper RA, Meinke PT, Nathan CF, Kirkman LA, Lin G. Development of a Highly Selective Plasmodium falciparum Proteasome Inhibitor with Anti-malaria Activity in Humanized Mice. Angew Chem Int Ed Engl 2021; 60:9279-9283. [PMID: 33433953 PMCID: PMC8087158 DOI: 10.1002/anie.202015845] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2020] [Revised: 12/29/2020] [Indexed: 01/01/2023]
Abstract
Plasmodium falciparum proteasome (Pf20S) inhibitors are active against Plasmodium at multiple stages-erythrocytic, gametocyte, liver, and gamete activation stages-indicating that selective Pf20S inhibitors possess the potential to be therapeutic, prophylactic, and transmission-blocking antimalarials. Starting from a reported compound, we developed a noncovalent, macrocyclic peptide inhibitor of the malarial proteasome with high species selectivity and improved pharmacokinetic properties. The compound demonstrates specific, time-dependent inhibition of the β5 subunit of the Pf20S, kills artemisinin-sensitive and artemisinin-resistant P. falciparum isolates in vitro and reduces parasitemia in humanized, P. falciparum-infected mice.
Collapse
Affiliation(s)
- Wenhu Zhan
- Department of Microbiology & Immunology, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - Hao Zhang
- Department of Microbiology & Immunology, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - John Ginn
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Annie Leung
- Department of Medicine, Division of Infectious Diseases, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - Yi J Liu
- Department of Medicine, Division of Infectious Diseases, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - Mayako Michino
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Akinori Toita
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Rei Okamoto
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Tzu-Tshin Wong
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Toshihiro Imaeda
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Ryoma Hara
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Takafumi Yukawa
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Sevil Chelebieva
- Department of Natural Sciences and Mathematics, Dominican University of California, San Rafael, CA, 94901, USA
| | | | - Maria Jose Lafuente-Monasterio
- Diseases of the Developing World (DDW), Tres Cantos Medicine Development Campus, GlaxoSmithKline, Severo Ochoa 2, 28760, Tres Cantos, Madrid, Spain
| | - Maria Santos Martinez-Martinez
- Diseases of the Developing World (DDW), Tres Cantos Medicine Development Campus, GlaxoSmithKline, Severo Ochoa 2, 28760, Tres Cantos, Madrid, Spain
| | | | | | - Kenjiro Sato
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Kazuyoshi Aso
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Philip J Rosenthal
- Department of Medicine, University of California, San Francisco, CA, 94143, USA
| | - Roland A Cooper
- Department of Natural Sciences and Mathematics, Dominican University of California, San Rafael, CA, 94901, USA
| | - Peter T Meinke
- Tri-Institutional Therapeutics Discovery Institute, 413 E 69th St, New York, NY, 10065, USA
| | - Carl F Nathan
- Department of Microbiology & Immunology, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - Laura A Kirkman
- Department of Medicine, Division of Infectious Diseases, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| | - Gang Lin
- Department of Microbiology & Immunology, Weill Cornell Medicine, 1300 York Ave, New York, NY, 10065, USA
| |
Collapse
|
25
|
Bepari AK, Reza HM. Identification of a novel inhibitor of SARS-CoV-2 3CL-PRO through virtual screening and molecular dynamics simulation. PeerJ 2021; 9:e11261. [PMID: 33954055 PMCID: PMC8051358 DOI: 10.7717/peerj.11261] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Accepted: 03/22/2021] [Indexed: 12/12/2022] Open
Abstract
Background The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has ravaged lives across the globe since December 2019, and new cases are still on the rise. Peoples’ ongoing sufferings trigger scientists to develop safe and effective remedies to treat this deadly viral disease. While repurposing the existing FDA-approved drugs remains in the front line, exploring drug candidates from synthetic and natural compounds is also a viable alternative. This study employed a comprehensive computational approach to screen inhibitors for SARS-CoV-2 3CL-PRO (also known as the main protease), a prime molecular target to treat coronavirus diseases. Methods We performed 100 ns GROMACS molecular dynamics simulations of three high-resolution X-ray crystallographic structures of 3CL-PRO. We extracted frames at 10 ns intervals to mimic conformational diversities of the target protein in biological environments. We then used AutoDock Vina molecular docking to virtual screen the Sigma–Aldrich MyriaScreen Diversity Library II, a rich collection of 10,000 druglike small molecules with diverse chemotypes. Subsequently, we adopted in silico computation of physicochemical properties, pharmacokinetic parameters, and toxicity profiles. Finally, we analyzed hydrogen bonding and other protein-ligand interactions for the short-listed compounds. Results Over the 100 ns molecular dynamics simulations of 3CL-PRO’s crystal structures, 6LZE, 6M0K, and 6YB7, showed overall integrity with mean Cα root-mean-square deviation (RMSD) of 1.96 (±0.35) Å, 1.98 (±0.21) Å, and 1.94 (±0.25) Å, respectively. Average root-mean-square fluctuation (RMSF) values were 1.21 ± 0.79 (6LZE), 1.12 ± 0.72 (6M0K), and 1.11 ± 0.60 (6YB7). After two phases of AutoDock Vina virtual screening of the MyriaScreen Diversity Library II, we prepared a list of the top 20 ligands. We selected four promising leads considering predicted oral bioavailability, druglikeness, and toxicity profiles. These compounds also demonstrated favorable protein-ligand interactions. We then employed 50-ns molecular dynamics simulations for the four selected molecules and the reference ligand 11a in the crystallographic structure 6LZE. Analysis of RMSF, RMSD, and hydrogen bonding along the simulation trajectories indicated that S51765 would form a more stable protein-ligand complexe with 3CL-PRO compared to other molecules. Insights into short-range Coulombic and Lennard-Jones potentials also revealed favorable binding of S51765 with 3CL-PRO. Conclusion We identified a potential lead for antiviral drug discovery against the SARS-CoV-2 main protease. Our results will aid global efforts to find safe and effective remedies for COVID-19.
Collapse
Affiliation(s)
- Asim Kumar Bepari
- Department of Pharmaceutical Sciences, North South University, Dhaka, Bangladesh
| | - Hasan Mahmud Reza
- Department of Pharmaceutical Sciences, North South University, Dhaka, Bangladesh
| |
Collapse
|
26
|
Moritsugu K, Takeuchi K, Kamiya N, Higo J, Yasumatsu I, Fukunishi Y, Fukuda I. Flexibility and Cell Permeability of Cyclic Ras-Inhibitor Peptides Revealed by the Coupled Nosé-Hoover Equation. J Chem Inf Model 2021; 61:1921-1930. [PMID: 33835817 DOI: 10.1021/acs.jcim.0c01427] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantifying the cell permeability of cyclic peptides is crucial for their rational drug design. However, the reasons remain unclear why a minor chemical modification, such as the difference between Ras inhibitors cyclorasin 9A5 and 9A54, can substantially change a peptide's permeability. To address this question, we performed enhanced sampling simulations of these two 11-mer peptides using the coupled Nosé-Hoover equation (cNH) we recently developed. The present cNH simulations realized temperature fluctuations over a wide range (240-600 K) in a dynamic manner, allowing structural samplings that were well validated by nuclear Overhauser effect measurements. The derived structural ensembles were comprehensively analyzed by all-atom structural clustering, mapping the derived clusters onto principal components (PCs) that characterize the cyclic structure, and calculating cluster-dependent geometric and chemical properties. The planar-open conformation was dominant in aqueous solvent, owing to inclusion of the Trp side chain in the main-chain ring, while the compact-closed conformation, which favors cell permeation due to its compactness and high polarity, was also accessible. Conformation-dependent cell permeability was observed in one of the derived PCs, demonstrating that decreased cell permeability in 9A54 is due to the high free energy barrier separating the two conformations. The origin of the change in free energy surface was determined to be loss of flexibility in the modified residues 2-3, resulting from the increased bulkiness of their side chains. The derived molecular mechanism of cell permeability highlights the significance of complete structural dynamics surveys for accelerating drug development with cyclic peptides.
Collapse
Affiliation(s)
- Kei Moritsugu
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehirocho, Tsurumi, Yokohama, Kanagawa 230-0045, Japan
| | - Koh Takeuchi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology, 2-3-26 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Isao Yasumatsu
- Structure-Based Drug Design Group, Organic Synthesis Department, Daiichi Sankyo RD Novare Co., Ltd., 1-2-58 Hiromachi, Shinagawa-ku, Tokyo 140-8710, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology, 2-3-26 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Ikuo Fukuda
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka 565-0871, Japan
| |
Collapse
|
27
|
Gutten O, Jurečka P, Aliakbar Tehrani Z, Buděšínský M, Řezáč J, Rulíšek L. Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments. Phys Chem Chem Phys 2021; 23:7280-7294. [PMID: 33876088 DOI: 10.1039/d0cp05993e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.
Collapse
Affiliation(s)
- Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
28
|
Golosov AA, Flyer AN, Amin J, Babu C, Gampe C, Li J, Liu E, Nakajima K, Nettleton D, Patel TJ, Reid PC, Yang L, Monovich LG. Design of Thioether Cyclic Peptide Scaffolds with Passive Permeability and Oral Exposure. J Med Chem 2021; 64:2622-2633. [PMID: 33629858 DOI: 10.1021/acs.jmedchem.0c01505] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Advances in the design of permeable peptides and in the synthesis of large arrays of macrocyclic peptides with diverse amino acids have evolved on parallel but independent tracks. Less precedent combines their respective attributes, thereby limiting the potential to identify permeable peptide ligands for key targets. Herein, we present novel 6-, 7-, and 8-mer cyclic peptides (MW 774-1076 g·mol-1) with passive permeability and oral exposure that feature the amino acids and thioether ring-closing common to large array formats, including DNA- and RNA-templated synthesis. Each oral peptide herein, selected from virtual libraries of partially N-methylated peptides using in silico methods, reflects the subset consistent with low energy conformations, low desolvation penalties, and passive permeability. We envision that, by retaining the backbone N-methylation pattern and consequent bias toward permeability, one can generate large peptide arrays with sufficient side chain diversity to identify permeability-biased ligands to a variety of protein targets.
Collapse
Affiliation(s)
- Andrei A Golosov
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Alec N Flyer
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Jakal Amin
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Charles Babu
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Christian Gampe
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Jingzhou Li
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Eugene Liu
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Katsumasa Nakajima
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - David Nettleton
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Tajesh J Patel
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Patrick C Reid
- PeptiDream, Inc., 3-25-23 Tonomachi, Kawasaki-Ku, Kawasaki, Kanagawa 210-0821, Japan
| | - Lihua Yang
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Lauren G Monovich
- Novartis Institutes for BioMedical Research, 181 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
29
|
Zhan W, Zhang H, Ginn J, Leung A, Liu YJ, Michino M, Toita A, Okamoto R, Wong T, Imaeda T, Hara R, Yukawa T, Chelebieva S, Tumwebaze PK, Lafuente‐Monasterio MJ, Martinez‐Martinez MS, Vendome J, Beuming T, Sato K, Aso K, Rosenthal PJ, Cooper RA, Meinke PT, Nathan CF, Kirkman LA, Lin G. Development of a Highly Selective
Plasmodium falciparum
Proteasome Inhibitor with Anti‐malaria Activity in Humanized Mice. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202015845] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Wenhu Zhan
- Department of Microbiology & Immunology Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - Hao Zhang
- Department of Microbiology & Immunology Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - John Ginn
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Annie Leung
- Department of Medicine Division of Infectious Diseases Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - Yi J. Liu
- Department of Medicine Division of Infectious Diseases Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - Mayako Michino
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Akinori Toita
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Rei Okamoto
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Tzu‐Tshin Wong
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Toshihiro Imaeda
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Ryoma Hara
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Takafumi Yukawa
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Sevil Chelebieva
- Department of Natural Sciences and Mathematics Dominican University of California San Rafael CA 94901 USA
| | | | - Maria Jose Lafuente‐Monasterio
- Diseases of the Developing World (DDW) Tres Cantos Medicine Development Campus GlaxoSmithKline Severo Ochoa 2 28760, Tres Cantos Madrid Spain
| | - Maria Santos Martinez‐Martinez
- Diseases of the Developing World (DDW) Tres Cantos Medicine Development Campus GlaxoSmithKline Severo Ochoa 2 28760, Tres Cantos Madrid Spain
| | | | | | - Kenjiro Sato
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Kazuyoshi Aso
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | | | - Roland A. Cooper
- Department of Natural Sciences and Mathematics Dominican University of California San Rafael CA 94901 USA
| | - Peter T. Meinke
- Tri-Institutional Therapeutics Discovery Institute 413 E 69th St New York NY 10065 USA
| | - Carl F. Nathan
- Department of Microbiology & Immunology Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - Laura A. Kirkman
- Department of Medicine Division of Infectious Diseases Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| | - Gang Lin
- Department of Microbiology & Immunology Weill Cornell Medicine 1300 York Ave New York NY 10065 USA
| |
Collapse
|
30
|
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
|
31
|
Álvarez RM, García AB, Riesco-Fagundo C, Martín JI, Varela C, Rodríguez Hergueta A, González Cantalapiedra E, Oyarzabal J, Di Geronimo B, Lorenzo M, Albarrán MI, Cebriá A, Cebrián D, Martínez-González S, Blanco-Aparicio C, Pastor J. Omipalisib inspired macrocycles as dual PI3K/mTOR inhibitors. Eur J Med Chem 2020; 211:113109. [PMID: 33360802 DOI: 10.1016/j.ejmech.2020.113109] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2020] [Revised: 11/25/2020] [Accepted: 12/13/2020] [Indexed: 02/08/2023]
Abstract
Activation of the phosphatidylinositol 3-kinase (PI3K)/mammalian target of rapamycin (mTOR) signaling pathway occurs frequently in a wide range of human cancers and is a main driver of cell growth, proliferation, survival, and chemoresistance of cancer cells. Compounds targeting this pathway are under active development as anticancer therapeutics and some of them have reached advanced clinical trials or been approved by the FDA. Dual PI3K/mTOR inhibitors combine multiple therapeutic efficacies in a single molecule by inhibiting the pathway both upstream and downstream of AKT. Herein, we report our efforts on the exploration of novel small molecule macrocycles (MCXs) as dual PI3K/mTOR inhibitors. Macrocyclization is an attractive approach used in drug discovery, as the semi-rigid character of these structures could provide improved potency, selectivity and favorable pharmacokinetic properties. Importantly, this strategy allows access to new chemical space thus obtaining a better intellectual property position. A series of MCXs based on GSK-2126458, a known clinical PI3K/mTOR inhibitor is described. These molecules showed potent biochemical and cellular dual PI3K/mTOR inhibition, demonstrated strong antitumoral effects in human cancer cell lines, and displayed good drug-like properties. Among them, MCX 83 presented remarkable selectivity against a panel of 468 kinases, high in vitro metabolic stability, and favorable pharmacokinetic parameters without significant CYP450 and h-ERG binding inhibition. This profile qualified this compound as a suitable candidate for future in vivo PK-PD and efficacy studies in mouse cancer models.
Collapse
Affiliation(s)
- Rosa M Álvarez
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Ana Belén García
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Concepción Riesco-Fagundo
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - José I Martín
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Carmen Varela
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Antonio Rodríguez Hergueta
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Esther González Cantalapiedra
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Julen Oyarzabal
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Bruno Di Geronimo
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Milagros Lorenzo
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - M Isabel Albarrán
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Antonio Cebriá
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - David Cebrián
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Sonia Martínez-González
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Carmen Blanco-Aparicio
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain
| | - Joaquín Pastor
- Experimental Therapeutics Programme, Spanish National Cancer Research Centre (CNIO), C/ Melchor Fernández Almagro 3, E-28029, Madrid, Spain.
| |
Collapse
|
32
|
Reyes Romero A, Ruiz-Moreno AJ, Groves MR, Velasco-Velázquez M, Dömling A. Benchmark of Generic Shapes for Macrocycles. J Chem Inf Model 2020; 60:6298-6313. [PMID: 33270455 PMCID: PMC7768607 DOI: 10.1021/acs.jcim.0c01038] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
![]()
Macrocycles
target proteins that are otherwise considered undruggable
because of a lack of hydrophobic cavities and the presence of extended
featureless surfaces. Increasing efforts by computational chemists
have developed effective software to overcome the restrictions of
torsional and conformational freedom that arise as a consequence of
macrocyclization. Moloc is an efficient algorithm, with an emphasis
on high interactivity, and has been constantly updated since 1986
by drug designers and crystallographers of the Roche biostructural
community. In this work, we have benchmarked the shape-guided algorithm
using a dataset of 208 macrocycles, carefully selected on the basis
of structural complexity. We have quantified the accuracy, diversity,
speed, exhaustiveness, and sampling efficiency in an automated fashion
and we compared them with four commercial (Prime, MacroModel, molecular
operating environment, and molecular dynamics) and four open-access
(experimental-torsion distance geometry with additional “basic
knowledge” alone and with Merck molecular force field minimization
or universal force field minimization, Cambridge Crystallographic
Data Centre conformer generator, and conformator) packages. With three-quarters
of the database processed below the threshold of high ring accuracy,
Moloc was identified as having the highest sampling efficiency and
exhaustiveness without producing thousands of conformations, random
ring splitting into two half-loops, and possibility to interactively
produce globular or flat conformations with diversity similar to Prime,
MacroModel, and molecular dynamics. The algorithm and the Python scripts
for full automatization of these parameters are freely available for
academic use.
Collapse
Affiliation(s)
- Atilio Reyes Romero
- Drug Design, Department of Pharmacy, University of Groningen, Antonius Deusinglaan 1, XB20, 9713 AV Groningen, The Netherlands
| | - Angel Jonathan Ruiz-Moreno
- Drug Design, Department of Pharmacy, University of Groningen, Antonius Deusinglaan 1, XB20, 9713 AV Groningen, The Netherlands.,Departamento de Farmacología y Unidad Periférica de Investigación en Biomedicina Trasnacional, Facultad de Medicina, Universidad Nacional Autónoma de México (UNAM), Av. Universidad 3000, Circuito Exterior S/N, Delegación Coyoacán, Ciudad Universitaria, 04510 Ciudad de México, Mexico.,Programa de Doctorado en Ciencias Biomédicas, UNAM, Av. Universidad 3000, Circuito Exterior S/N. Delegación Coyoacán, Ciudad Universitaria, 04510 Ciudad de México, Mexico
| | - Matthew R Groves
- Drug Design, Department of Pharmacy, University of Groningen, Antonius Deusinglaan 1, XB20, 9713 AV Groningen, The Netherlands
| | - Marco Velasco-Velázquez
- Departamento de Farmacología y Unidad Periférica de Investigación en Biomedicina Trasnacional, Facultad de Medicina, Universidad Nacional Autónoma de México (UNAM), Av. Universidad 3000, Circuito Exterior S/N, Delegación Coyoacán, Ciudad Universitaria, 04510 Ciudad de México, Mexico
| | - Alexander Dömling
- Drug Design, Department of Pharmacy, University of Groningen, Antonius Deusinglaan 1, XB20, 9713 AV Groningen, The Netherlands
| |
Collapse
|
33
|
Saha I, Dang EK, Svatunek D, Houk KN, Harran PG. Computational generation of an annotated gigalibrary of synthesizable, composite peptidic macrocycles. Proc Natl Acad Sci U S A 2020; 117:24679-24690. [PMID: 32948694 PMCID: PMC7547232 DOI: 10.1073/pnas.2007304117] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Peptidomimetic macrocycles have the potential to regulate challenging therapeutic targets. Structures of this type having precise shapes and drug-like character are particularly coveted, but are relatively difficult to synthesize. Our laboratory has developed robust methods that integrate small-peptide units into designed scaffolds. These methods create macrocycles and embed condensed heterocycles to diversify outcomes and improve pharmacological properties. The hypothetical scope of the methodology is vast and far outpaces the capacity of our experimental format. We now describe a computational rendering of our methodology that creates an in silico three-dimensional library of composite peptidic macrocycles. Our open-source platform, CPMG (Composite Peptide Macrocycle Generator), has algorithmically generated a library of 2,020,794,198 macrocycles that can result from the multistep reaction sequences we have developed. Structures are generated based on predicted site reactivity and filtered on the basis of physical and three-dimensional properties to identify maximally diverse compounds for prioritization. For conformational analyses, we also introduce ConfBuster++, an RDKit port of the open-source software ConfBuster, which allows facile integration with CPMG and ready parallelization for better scalability. Our approach deeply probes ligand space accessible via our synthetic methodology and provides a resource for large-scale virtual screening.
Collapse
Affiliation(s)
- Ishika Saha
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095
| | - Eric K Dang
- Department of Computer Science, University of California, Los Angeles, CA 90095
| | - Dennis Svatunek
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095
| | - Kendall N Houk
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095;
| | - Patrick G Harran
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095;
| |
Collapse
|
34
|
Sindhikara D, Wagner M, Gkeka P, Güssregen S, Tiwari G, Hessler G, Yapici E, Li Z, Evers A. Automated Design of Macrocycles for Therapeutic Applications: From Small Molecules to Peptides and Proteins. J Med Chem 2020; 63:12100-12115. [PMID: 33017535 DOI: 10.1021/acs.jmedchem.0c01500] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Macrocycles and cyclic peptides are increasingly attractive therapeutic modalities as they often have improved affinity, are able to bind to extended protein surfaces, and otherwise have favorable properties. Macrocyclization of a known binder may stabilize its bioactive conformation and improve its metabolic stability, cell permeability, and in certain cases oral bioavailability. Herein, we present implementation and application of an approach that automatically generates, evaluates, and proposes cyclizations utilizing a library of well-established chemical reactions and reagents. Using the three-dimensional (3D) conformation of the linear molecule in complex with a target protein as the starting point, this approach identifies attachment points, generates linkers, evaluates their geometric compatibility, and ranks the resulting molecules with respect to their predicted conformational stability and interactions with the target protein. As we show here with prospective and retrospective case studies, this procedure can be applied for the macrocyclization of small molecules and peptides and even PROteolysis TArgeting Chimeras (PROTACs) and proteins.
Collapse
Affiliation(s)
- Dan Sindhikara
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Michael Wagner
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| | - Paraskevi Gkeka
- Integrated Drug Discovery, Sanofi R&D, 1 Avenue Pierre Brossolette, 91385 Chilly-Mazarin, France
| | - Stefan Güssregen
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| | - Garima Tiwari
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| | - Gerhard Hessler
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| | - Engin Yapici
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Ziyu Li
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| | - Andreas Evers
- Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, Industriepark Hoechst, 65926 Frankfurt am Main, Germany
| |
Collapse
|
35
|
Caulfield TR, Hayes KE, Qiu Y, Coban M, Seok Oh J, Lane AL, Yoshimitsu T, Hazlehurst L, Copland JA, Tun HW. A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor. Biomolecules 2020; 10:E1407. [PMID: 33027969 PMCID: PMC7599554 DOI: 10.3390/biom10101407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 09/11/2020] [Accepted: 09/29/2020] [Indexed: 12/03/2022] Open
Abstract
Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.
Collapse
Affiliation(s)
- Thomas R. Caulfield
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
- Department of Neuroscience, Mayo Clinic, Jacksonville, FL 32224, USA
- Department of Neurosurgery, Mayo Clinic, Jacksonville, FL 32224, USA
- Department of Health Sciences Research, Mayo Clinic, Jacksonville, FL 32224, USA
| | - Karen E. Hayes
- Modulation Therapeutics, Inc., Morgantown, WV 26506, USA;
| | - Yushi Qiu
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
| | - Mathew Coban
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
| | - Joon Seok Oh
- Department of Chemistry, University of North Florida, Jacksonville, FL 32224, USA;
| | - Amy L. Lane
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
- Department of Chemistry, University of North Florida, Jacksonville, FL 32224, USA;
| | - Takehiko Yoshimitsu
- Division of Pharmaceutical Sciences, Graduate School of Medicine, Dentistry, and Pharmaceutical Sciences, Okayama University, 1-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan;
| | - Lori Hazlehurst
- Department of Pharmaceutical Sciences, West Virginia University, Morgantown, WV 26506, USA;
| | - John A. Copland
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
| | - Han W. Tun
- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL 32224, USA; (Y.Q.); (M.C.); (A.L.L.); (J.A.C.)
- Department of Hematology/Oncology, Mayo Clinic, Jacksonville, FL 32224, USA
| |
Collapse
|
36
|
Synthetic group A streptogramin antibiotics that overcome Vat resistance. Nature 2020; 586:145-150. [PMID: 32968273 PMCID: PMC7546582 DOI: 10.1038/s41586-020-2761-3] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 08/19/2020] [Indexed: 11/08/2022]
Abstract
Natural products serve as chemical blueprints for the majority of antibiotics in our clinical arsenal. The evolutionary process by which these molecules arise is inherently accompanied by the co-evolution of resistance mechanisms that shorten the clinical lifetime of any given class1. Virginiamycin acetyltransferases (Vats) are resistance proteins that provide protection against streptogramins2, potent Gram-positive antibiotics that inhibit the bacterial ribosome3. Due to the challenge of selectively modifying the chemically complex, 23-membered macrocyclic scaffold of group A streptogramins, analogs that overcome Vat resistance have not been previously accessed2. Here we report the design, synthesis, and antibacterial evaluation of group A streptogramin antibiotics with unprecedented structural variability. Using cryo-electron microscopy and forcefield-based refinement, we characterize the binding of eight analogs to the bacterial ribosome at high resolution, revealing new binding interactions that extend into the peptidyl tRNA binding site and towards synergistic binders that occupy the nascent peptide exit tunnel (NPET). One of these analogs has excellent activity against several streptogramin-resistant strains of S. aureus, exhibits decreased acetylation rates in vitro, and is effective at lowering bacterial load in a mouse model of infection. Our results demonstrate that the combination of rational design and modular chemical synthesis can revitalize classes of antibiotics that are limited by naturally arising resistance mechanisms.
Collapse
|
37
|
Hawkins PCD, Wlodek S. Decisions with Confidence: Application to the Conformation Sampling of Molecules in the Solid State. J Chem Inf Model 2020; 60:3518-3533. [DOI: 10.1021/acs.jcim.0c00358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Paul C. D. Hawkins
- OpenEye Scientific, 9 Bisbee Court, Suite D, Santa Fe, New Mexico 87508, United States
| | - Stanislaw Wlodek
- OpenEye Scientific, 9 Bisbee Court, Suite D, Santa Fe, New Mexico 87508, United States
| |
Collapse
|
38
|
Kamenik AS, Kraml J, Hofer F, Waibl F, Quoika PK, Kahler U, Schauperl M, Liedl KR. Macrocycle Cell Permeability Measured by Solvation Free Energies in Polar and Apolar Environments. J Chem Inf Model 2020; 60:3508-3517. [PMID: 32551643 PMCID: PMC7388155 DOI: 10.1021/acs.jcim.0c00280] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The relation of surface polarity and conformational preferences is decisive for cell permeability and thus bioavailability of macrocyclic drugs. Here, we employ grid inhomogeneous solvation theory (GIST) to calculate solvation free energies for a series of six macrocycles in water and chloroform as a measure of passive membrane permeability. We perform accelerated molecular dynamics simulations to capture a diverse structural ensemble in water and chloroform, allowing for a direct profiling of solvent-dependent conformational preferences. Subsequent GIST calculations facilitate a quantitative measure of solvent preference in the form of a transfer free energy, calculated from the ensemble-averaged solvation free energies in water and chloroform. Hence, the proposed method considers how the conformational diversity of macrocycles in polar and apolar solvents translates into transfer free energies. Following this strategy, we find a striking correlation of 0.92 between experimentally determined cell permeabilities and calculated transfer free energies. For the studied model systems, we find that the transfer free energy exceeds the purely water-based solvation free energies as a reliable estimate of cell permeability and that conformational sampling is imperative for a physically meaningful model. We thus recommend this purely physics-based approach as a computational tool to assess cell permeabilities of macrocyclic drug candidates.
Collapse
Affiliation(s)
- Anna S Kamenik
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Johannes Kraml
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Florian Hofer
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Patrick K Quoika
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Ursula Kahler
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Michael Schauperl
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| | - Klaus R Liedl
- Institute of General, Inorganic and Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80-82, Innsbruck A-6020 Austria
| |
Collapse
|
39
|
Paulsen JL, Yu HS, Sindhikara D, Wang L, Appleby T, Villaseñor AG, Schmitz U, Shivakumar D. Evaluation of Free Energy Calculations for the Prioritization of Macrocycle Synthesis. J Chem Inf Model 2020; 60:3489-3498. [DOI: 10.1021/acs.jcim.0c00132] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Janet L. Paulsen
- Schrödinger, Inc., 120 West 45th Street, 17th Floor, New York, New York 10036, United States
| | - Haoyu S. Yu
- Schrödinger, Inc., 120 West 45th Street, 17th Floor, New York, New York 10036, United States
| | - Dan Sindhikara
- Schrödinger, Inc., 120 West 45th Street, 17th Floor, New York, New York 10036, United States
| | - Lingle Wang
- Schrödinger, Inc., 120 West 45th Street, 17th Floor, New York, New York 10036, United States
| | - Todd Appleby
- Gilead, 333 Lakeside Drive, Foster City, California 94404, United States
| | | | - Uli Schmitz
- Gilead, 333 Lakeside Drive, Foster City, California 94404, United States
| | | |
Collapse
|
40
|
Le Roux A, Blaise É, Boudreault PL, Comeau C, Doucet A, Giarrusso M, Collin MP, Neubauer T, Kölling F, Göller AH, Seep L, Tshitenge DT, Wittwer M, Kullmann M, Hillisch A, Mittendorf J, Marsault E. Structure-Permeability Relationship of Semipeptidic Macrocycles-Understanding and Optimizing Passive Permeability and Efflux Ratio. J Med Chem 2020; 63:6774-6783. [PMID: 32453569 DOI: 10.1021/acs.jmedchem.0c00013] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We herein report the first thorough analysis of the structure-permeability relationship of semipeptidic macrocycles. In total, 47 macrocycles were synthesized using a hybrid solid-phase/solution strategy, and then their passive and cellular permeability was assessed using the parallel artificial membrane permeability assay (PAMPA) and Caco-2 assay, respectively. The results indicate that semipeptidic macrocycles generally possess high passive permeability based on the PAMPA, yet their cellular permeability is governed by efflux, as reported in the Caco-2 assay. Structural variations led to tractable structure-permeability and structure-efflux relationships, wherein the linker length, stereoinversion, N-methylation, and peptoids site-specifically impact the permeability and efflux. Extensive nuclear magnetic resonance, molecular dynamics, and ensemble-based three-dimensional polar surface area (3D-PSA) studies showed that ensemble-based 3D-PSA is a good predictor of passive permeability.
Collapse
Affiliation(s)
- Antoine Le Roux
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | - Émilie Blaise
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | - Pierre-Luc Boudreault
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | - Christian Comeau
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | - Annie Doucet
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | - Marilena Giarrusso
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| | | | - Thomas Neubauer
- Drug Discovery, Pharmaceuticals, Bayer AG, Wuppertal D-42096, Germany
| | - Florian Kölling
- Drug Discovery, Pharmaceuticals, Bayer AG, Wuppertal D-42096, Germany
| | - Andreas H Göller
- Drug Discovery, Pharmaceuticals, Bayer AG, Wuppertal D-42096, Germany
| | - Lea Seep
- Drug Discovery, Pharmaceuticals, Bayer AG, Wuppertal D-42096, Germany
| | | | - Matthias Wittwer
- Drug Discovery, Pharmaceuticals, Bayer AG, Wuppertal D-42096, Germany
| | | | | | | | - Eric Marsault
- Department of Pharmacology-Physiology, Institut de Pharmacologie de Sherbrooke, 3001, 12e av nord, Sherbrooke, Québec J1H 5N4, Canada
| |
Collapse
|
41
|
Mulligan VK. The emerging role of computational design in peptide macrocycle drug discovery. Expert Opin Drug Discov 2020; 15:833-852. [PMID: 32345066 DOI: 10.1080/17460441.2020.1751117] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Drug discovery is a laborious process with rising cost per new drug. Peptide macrocycles are promising therapeutics, though conformational flexibility can reduce target affinity and specificity. Recent computational advancements address this problem by enabling rational design of rigidly folded peptide macrocycles. AREAS COVERED This review summarizes currently approved peptide macrocycle therapeutics and discusses advantages of mesoscale drugs over small molecules or protein therapeutics. It describes the history, rationale, and state of the art of computational tools, such as Rosetta, that allow the design of rigidly structured peptide macrocycles. The emerging pipeline for designing peptide macrocycle drugs is described, including current challenges in designing permeable molecules that can emulate the chameleonic behavior of natural macrocycles. Prospects for reducing computational cost and improving accuracy with emerging computational technologies are also discussed. EXPERT OPINION To embrace computational design of peptide macrocycle drugs, we must shift current attitudes regarding the role of computation in drug discovery, and move beyond Lipinski's rules. This technology has the potential to shift failures to earlier in silico stages of the drug discovery process, improving success rates in costly clinical trials. Given the available tools, now is the time for drug developers to incorporate peptide macrocycle design into drug discovery pipelines.
Collapse
Affiliation(s)
- Vikram K Mulligan
- Systems Biology, Center for Computational Biology, Flatiron Institute , New York, NY, USA
| |
Collapse
|
42
|
Maryam A, Khalid RR, Siddiqi AR, Ece A. E-pharmacophore based virtual screening for identification of dual specific PDE5A and PDE3A inhibitors as potential leads against cardiovascular diseases. J Biomol Struct Dyn 2020; 39:2302-2317. [DOI: 10.1080/07391102.2020.1748718] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Arooma Maryam
- Department of Biosciences, COMSATS University, Islamabad, Pakistan
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Biruni University, Istanbul, Turkey
| | | | | | - Abdulilah Ece
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Biruni University, Istanbul, Turkey
| |
Collapse
|
43
|
Zin PPK, Williams G, Fourches D. SIME: synthetic insight-based macrolide enumerator to generate the V1B library of 1 billion macrolides. J Cheminform 2020; 12:23. [PMID: 33431002 PMCID: PMC7146965 DOI: 10.1186/s13321-020-00427-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 03/27/2020] [Indexed: 11/24/2022] Open
Abstract
We report on a new cheminformatics enumeration technology—SIME, synthetic insight-based macrolide enumerator—a new and improved software technology. SIME can enumerate fully assembled macrolides with synthetic feasibility by utilizing the constitutional and structural knowledge extracted from biosynthetic aspects of macrolides. Taken into account by the software are key information such as positions in macrolide structures at which chemical components can be inserted, and the types of structural motifs and sugars of interest that can be synthesized and incorporated at those positions. Additionally, we report on the chemical distribution analysis of the newly SIME-generated V1B (virtual 1 billion) library of macrolides. Those compounds were built based on the core of the Erythromycin structure, 13 structural motifs and a library of sugars derived from eighteen bioactive macrolides. This new enumeration technology can be coupled with cheminformatics approaches such as QSAR modeling and molecular docking to aid in drug discovery for rational designing of next generation macrolide therapeutics with desirable pharmacokinetic properties.![]()
Collapse
Affiliation(s)
- Phyo Phyo Kyaw Zin
- Department of Chemistry, North Carolina State University, Raleigh, NC, USA.,Bioinformatics Research Center, North Carolina State University, Raleigh, NC, USA
| | - Gavin Williams
- Department of Chemistry, North Carolina State University, Raleigh, NC, USA.,Comparative Medicine Institute, North Carolina State University, Raleigh, NC, USA
| | - Denis Fourches
- Department of Chemistry, North Carolina State University, Raleigh, NC, USA. .,Bioinformatics Research Center, North Carolina State University, Raleigh, NC, USA. .,Comparative Medicine Institute, North Carolina State University, Raleigh, NC, USA.
| |
Collapse
|
44
|
Wang S, Witek J, Landrum GA, Riniker S. Improving Conformer Generation for Small Rings and Macrocycles Based on Distance Geometry and Experimental Torsional-Angle Preferences. J Chem Inf Model 2020; 60:2044-2058. [PMID: 32155061 DOI: 10.1021/acs.jcim.0c00025] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The conformer generator ETKDG is a stochastic search method that utilizes distance geometry together with knowledge derived from experimental crystal structures. It has been shown to generate good conformers for acyclic, flexible molecules. This work builds on ETKDG to improve conformer generation of molecules containing small or large aliphatic (i.e., non-aromatic) rings. For one, we devise additional torsional-angle potentials to describe small aliphatic rings and adapt the previously developed potentials for acyclic bonds to facilitate the sampling of macrocycles. However, due to the larger number of degrees of freedom of macrocycles, the conformational space to sample is much broader than for small molecules, creating a challenge for conformer generators. We therefore introduce different heuristics to restrict the search space of macrocycles and bias the sampling toward more experimentally relevant structures. Specifically, we show the usage of elliptical geometry and customizable Coulombic interactions as heuristics. The performance of the improved ETKDG is demonstrated on test sets of diverse macrocycles and cyclic peptides. The code developed here will be incorporated into the 2020.03 release of the open-source cheminformatics library RDKit.
Collapse
Affiliation(s)
- Shuzhe Wang
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jagna Witek
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | | | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
45
|
Wallraven K, Holmelin FL, Glas A, Hennig S, Frolov AI, Grossmann TN. Adapting free energy perturbation simulations for large macrocyclic ligands: how to dissect contributions from direct binding and free ligand flexibility. Chem Sci 2020; 11:2269-2276. [PMID: 32180932 PMCID: PMC7057854 DOI: 10.1039/c9sc04705k] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 01/20/2020] [Indexed: 11/25/2022] Open
Abstract
Large and flexible ligands gain increasing interest in the development of bioactive agents. They challenge the applicability of computational ligand optimization strategies originally developed for small molecules. Free energy perturbation (FEP) is often used for predicting binding affinities of small molecule ligands, however, its use for more complex ligands remains limited. Herein, we report the structure-based design of peptide macrocycles targeting the protein binding site of human adaptor protein 14-3-3. We observe a surprisingly strong dependency of binding affinities on relatively small variations in substituent size. FEP was performed to rationalize observed trends. To account for insufficient convergence of FEP, restrained calculations were performed and complemented with extensive REST MD simulations of the free ligands. These calculations revealed that changes in affinity originate both from altered direct interactions and conformational changes of the free ligand. In addition, MD simulations provided the basis to rationalize unexpected trends in ligand lipophilicity. We also verified the anticipated interaction site and binding mode for one of the high affinity ligands by X-ray crystallography. The introduced fully-atomistic simulation protocol can be used to rationalize the development of structurally complex ligands which will support future ligand maturation efforts.
Collapse
Affiliation(s)
- Kerstin Wallraven
- Department of Chemistry & Pharmaceutical Sciences , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands .
| | - Fredrik L Holmelin
- Medicinal Chemistry, Research and Early Development Cardiovascular, Renal and Metabolism , BioPharmaceuticals R&D , AstraZeneca , Pepparedsleden 1, Mölndal , 431 83 , Sweden .
| | - Adrian Glas
- Department of Chemistry & Pharmaceutical Sciences , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands .
| | - Sven Hennig
- Department of Chemistry & Pharmaceutical Sciences , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands .
| | - Andrey I Frolov
- Medicinal Chemistry, Research and Early Development Cardiovascular, Renal and Metabolism , BioPharmaceuticals R&D , AstraZeneca , Pepparedsleden 1, Mölndal , 431 83 , Sweden .
| | - Tom N Grossmann
- Department of Chemistry & Pharmaceutical Sciences , VU University Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands .
| |
Collapse
|
46
|
Conformational analysis of macrocycles: comparing general and specialized methods. J Comput Aided Mol Des 2020; 34:231-252. [PMID: 31965404 PMCID: PMC7036058 DOI: 10.1007/s10822-020-00277-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2019] [Accepted: 01/03/2020] [Indexed: 11/24/2022]
Abstract
Abstract Macrocycles represent an important class of medicinally relevant small molecules due to their interesting biological properties. Therefore, a firm understanding of their conformational preferences is important for drug design. Given the importance of macrocycle-protein modelling in drug discovery, we envisaged that a systematic study of both classical and recent specialized methods would provide guidance for other practitioners within the field. In this study we compare the performance of the general, well established conformational analysis methods Monte Carlo Multiple Minimum (MCMM) and Mixed Torsional/Low-Mode sampling (MTLMOD) with two more recent and specialized macrocycle sampling techniques: MacroModel macrocycle Baseline Search (MD/LLMOD) and Prime macrocycle conformational sampling (PRIME-MCS). Using macrocycles extracted from 44 macrocycle-protein X-ray crystallography complexes, we evaluated each method based on their ability to (i) generate unique conformers, (ii) generate unique macrocycle ring conformations, (iii) identify the global energy minimum, (iv) identify conformers similar to the X-ray ligand conformation after Protein Preparation Wizard treatment (X-rayppw), and (v) to the X-rayppw ring conformation. Computational speed was also considered. In addition, conformational coverage, as defined by the number of conformations identified, was studied. In order to study the relative energies of the bioactive conformations, the energy differences between the global energy minima and the energy minimized X-rayppw structures and, the global energy minima and the MCMM-Exhaustive (1,000,000 search steps) generated conformers closest to the X-rayppw structure, were calculated and analysed. All searches were performed using relatively short run times (10,000 steps for MCMM, MTLMOD and MD/LLMOD). To assess the performance of the methods, they were compared to an exhaustive MCMM search using 1,000,000 search steps for each of the 44 macrocycles (requiring ca 200 times more CPU time). Prior to our analysis, we also investigated if the general search methods MCMM and MTLMOD could also be optimized for macrocycle conformational sampling. Taken together, our work concludes that the more general methods can be optimized for macrocycle modelling by slightly adjusting the settings around the ring closure bond. In most cases, MCMM and MTLMOD with either standard or enhanced settings performed well in comparison to the more specialized macrocycle sampling methods MD/LLMOD and PRIME-MCS. When using enhanced settings for MCMM and MTLMOD, the X-rayppw conformation was regenerated with the greatest accuracy. The, MD/LLMOD emerged as the most efficient method for generating the global energy minima. Graphic abstract ![]()
Electronic supplementary material The online version of this article (10.1007/s10822-020-00277-2) contains supplementary material, which is available to authorized users.
Collapse
|
47
|
Peng C, Atilaw Y, Wang J, Xu Z, Poongavanam V, Shi J, Kihlberg J, Zhu W, Erdélyi M. Conformation of the Macrocyclic Drug Lorlatinib in Polar and Nonpolar Environments: A MD Simulation and NMR Study. ACS OMEGA 2019; 4:22245-22250. [PMID: 31891108 PMCID: PMC6933765 DOI: 10.1021/acsomega.9b03797] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Accepted: 11/26/2019] [Indexed: 05/02/2023]
Abstract
The replica exchange molecular dynamics (REMD) simulation is demonstrated to readily predict the conformations of the macrocyclic drug lorlatinib, as validated by solution NMR studies. In aqueous solution, lorlatinib adopts a conformer identical to its target bound structure. This conformer is stabilized by an extensive hydrogen bond network to the solvents. In chloroform, lorlatinib populates two conformers with the second one being less polar, which may contribute to lorlatinib's ability to cross cell membranes.
Collapse
Affiliation(s)
- Cheng Peng
- Drug
Discovery and Design Center; CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy
of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
- University
of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, China
| | - Yoseph Atilaw
- Department
of Chemistry-BMC, Uppsala University, Box 576, SE-751 23 Uppsala, Sweden
| | - Jinan Wang
- Drug
Discovery and Design Center; CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy
of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Zhijian Xu
- Drug
Discovery and Design Center; CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy
of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
- University
of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, China
| | | | - Jiye Shi
- Drug
Discovery and Design Center; CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy
of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Jan Kihlberg
- Department
of Chemistry-BMC, Uppsala University, Box 576, SE-751 23 Uppsala, Sweden
| | - Weiliang Zhu
- Drug
Discovery and Design Center; CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy
of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
- University
of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, China
- E-mail: (W.Z.)
| | - Máté Erdélyi
- Department
of Chemistry-BMC, Uppsala University, Box 576, SE-751 23 Uppsala, Sweden
- E-mail: (M.E.)
| |
Collapse
|
48
|
Kadukova M, Chupin V, Grudinin S. Docking rigid macrocycles using Convex-PL, AutoDock Vina, and RDKit in the D3R Grand Challenge 4. J Comput Aided Mol Des 2019; 34:191-200. [PMID: 31784861 DOI: 10.1007/s10822-019-00263-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 11/22/2019] [Indexed: 12/15/2022]
Abstract
The D3R Grand Challenge 4 provided a brilliant opportunity to test macrocyclic docking protocols on a diverse high-quality experimental data. We participated in both pose and affinity prediction exercises. Overall, we aimed to use an automated structure-based docking pipeline built around a set of tools developed in our team. This exercise again demonstrated a crucial importance of the correct local ligand geometry for the overall success of docking. Starting from the second part of the pose prediction stage, we developed a stable pipeline for sampling macrocycle conformers. This resulted in the subangstrom average precision of our pose predictions. In the affinity prediction exercise we obtained average results. However, we could improve these when using docking poses submitted by the best predictors. Our docking tools including the Convex-PL scoring function are available at https://team.inria.fr/nano-d/software/.
Collapse
Affiliation(s)
- Maria Kadukova
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Vladimir Chupin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Sergei Grudinin
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France.
| |
Collapse
|
49
|
Novotná B, Vaneková L, Zavřel M, Buděšínský M, Dejmek M, Smola M, Gutten O, Tehrani ZA, Pimková Polidarová M, Brázdová A, Liboska R, Štěpánek I, Vavřina Z, Jandušík T, Nencka R, Rulíšek L, Bouřa E, Brynda J, Páv O, Birkuš G. Enzymatic Preparation of 2'-5',3'-5'-Cyclic Dinucleotides, Their Binding Properties to Stimulator of Interferon Genes Adaptor Protein, and Structure/Activity Correlations. J Med Chem 2019; 62:10676-10690. [PMID: 31715099 DOI: 10.1021/acs.jmedchem.9b01062] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Cyclic dinucleotides are second messengers in the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, which plays an important role in recognizing tumor cells and viral or bacterial infections. They bind to the STING adaptor protein and trigger expression of cytokines via TANK binding kinase 1 (TBK1)/interferon regulatory factor 3 (IRF3) and inhibitor of nuclear factor-κB (IκB) kinase (IKK)/nuclear factor-κB (NFκB) signaling cascades. In this work, we describe an enzymatic preparation of 2'-5',3'-5'-cyclic dinucleotides (2'3'CDNs) with use of cyclic GMP-AMP synthases (cGAS) from human, mouse, and chicken. We profile substrate specificity of these enzymes by employing a small library of nucleotide-5'-triphosphate (NTP) analogues and use them to prepare 33 2'3'CDNs. We also determine affinity of these CDNs to five different STING haplotypes in cell-based and biochemical assays and describe properties needed for their optimal activity toward all STING haplotypes. Next, we study their effect on cytokine and chemokine induction by human peripheral blood mononuclear cells (PBMCs) and evaluate their cytotoxic effect on monocytes. Additionally, we report X-ray crystal structures of two new CDNs bound to STING protein and discuss structure-activity relationship by using quantum and molecular mechanical (QM/MM) computational modeling.
Collapse
Affiliation(s)
- Barbora Novotná
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic.,Faculty of Science , Charles University , Prague 110 00 , Czech Republic
| | - Lenka Vaneková
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic.,Faculty of Science , Charles University , Prague 110 00 , Czech Republic
| | - Martin Zavřel
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Miloš Buděšínský
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Milan Dejmek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Miroslav Smola
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Zahra Aliakbar Tehrani
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Markéta Pimková Polidarová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic.,Faculty of Science , Charles University , Prague 110 00 , Czech Republic
| | - Andrea Brázdová
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Radek Liboska
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Ivan Štěpánek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Zdeněk Vavřina
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic.,Faculty of Science , Charles University , Prague 110 00 , Czech Republic
| | - Tomáš Jandušík
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic.,Faculty of Food and Biochemical Technology , University of Chemistry and Technology , Prague 166 28 , Czech Republic
| | - Radim Nencka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Evžen Bouřa
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Jiří Brynda
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Ondřej Páv
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| | - Gabriel Birkuš
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences and Gilead Sciences Research Centre at IOCB , Flemingovo nam. 2 , Prague 16610 , Czech Republic
| |
Collapse
|
50
|
Modelling the binding mode of macrocycles: Docking and conformational sampling. Bioorg Med Chem 2019; 28:115143. [PMID: 31771798 DOI: 10.1016/j.bmc.2019.115143] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Revised: 09/12/2019] [Accepted: 09/25/2019] [Indexed: 11/21/2022]
Abstract
Drug discovery is increasingly tackling challenging protein binding sites regarding molecular recognition and druggability, including shallow and solvent-exposed protein-protein interaction interfaces. Macrocycles are emerging as promising chemotypes to modulate such sites. Despite their chemical complexity, macrocycles comprise important drugs and offer advantages compared to non-cyclic analogs, hence the recent impetus in the medicinal chemistry of macrocycles. Elaboration of macrocycles, or constituent fragments, can strongly benefit from knowledge of their binding mode to a target. When such information from X-ray crystallography is elusive, computational docking can provide working models. However, few studies have explored docking protocols for macrocycles, since conventional docking methods struggle with the conformational complexity of macrocycles, and also potentially with the shallower topology of their binding sites. Indeed, macrocycle binding mode prediction with the mainstream docking software GOLD has hardly been explored. Here, we present an in-depth study of macrocycle docking with GOLD and the ChemPLP scores. First, we summarize the thorough curation of a test set of 41 protein-macrocycle X-ray structures, raising the issue of lattice contacts with such systems. Rigid docking of the known bioactive conformers was successful (three top ranked poses) for 92.7% of the systems, in absence of crystallographic waters. Thus, without conformational search issues, scoring performed well. However, docking success dropped to 29.3% with the GOLD built-in conformational search. Yet, the success rate doubled to 58.5% when GOLD was supplied with extensive conformer ensembles docked rigidly. The reasons for failure, sampling or scoring, were analyzed, exemplified with particular cases. Overall, binding mode prediction of macrocycles remains challenging, but can be much improved with tailored protocols. The analysis of the interplay between conformational sampling and docking will be relevant to the prospective modelling of macrocycles in general.
Collapse
|