1
|
Sun S, Fushimi M, Rossetti T, Kaur N, Ferreira J, Miller M, Quast J, van den Heuvel J, Steegborn C, Levin LR, Buck J, Myers RW, Kargman S, Liverton N, Meinke PT, Huggins DJ. Scaffold Hopping and Optimization of Small Molecule Soluble Adenyl Cyclase Inhibitors Led by Free Energy Perturbation. J Chem Inf Model 2023; 63:2828-2841. [PMID: 37060320 DOI: 10.1021/acs.jcim.2c01577] [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: 04/16/2023]
Abstract
Free energy perturbation is a computational technique that can be used to predict how small changes to an inhibitor structure will affect the binding free energy to its target. In this paper, we describe the utility of free energy perturbation with FEP+ in the hit-to-lead stage of a drug discovery project targeting soluble adenyl cyclase. The project was structurally enabled by X-ray crystallography throughout. We employed free energy perturbation to first scaffold hop to a preferable chemotype and then optimize the binding affinity to sub-nanomolar levels while retaining druglike properties. The results illustrate that effective use of free energy perturbation can enable a drug discovery campaign to progress rapidly from hit to lead, facilitating proof-of-concept studies that enable target validation.
Collapse
Affiliation(s)
- Shan Sun
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Makoto Fushimi
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Thomas Rossetti
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - Navpreet Kaur
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - Jacob Ferreira
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - Michael Miller
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Jonathan Quast
- Department of Biochemistry, University of Bayreuth, Bayreuth 95440, Germany
| | | | - Clemens Steegborn
- Department of Biochemistry, University of Bayreuth, Bayreuth 95440, Germany
| | - Lonny R Levin
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - Jochen Buck
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - Robert W Myers
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Stacia Kargman
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Nigel Liverton
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Peter T Meinke
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
- Department of Pharmacology, Weill Cornell Medicine, New York City, New York 10056, United States
| | - David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| |
Collapse
|
2
|
Losev TV, Gerasimov IS, Panova MV, Lisov AA, Abdyusheva YR, Rusina PV, Zaletskaya E, Stroganov OV, Medvedev MG, Novikov FN. Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design. J Chem Inf Model 2023; 63:1239-1248. [PMID: 36763797 DOI: 10.1021/acs.jcim.2c01212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Bioisosteres are molecules that differ in substituents but still have very similar shapes. Bioisosteric replacements are ubiquitous in modern drug design, where they are used to alter metabolism, change bioavailability, or modify activity of the lead compound. Prediction of relative affinities of bioisosteres with computational methods is a long-standing task; however, the very shape closeness makes bioisosteric substitutions almost intractable for computational methods, which use standard force fields. Here, we design a quantum mechanical (QM)-cluster approach based on the GFN2-xTB semi-empirical quantum-chemical method and apply it to a set of H → F bioisosteric replacements. The proposed methodology enables advanced prediction of biological activity change upon bioisosteric substitution of -H with -F, with the standard deviation of 0.60 kcal/mol, surpassing the ChemPLP scoring function (0.83 kcal/mol), and making QM-based ΔΔG estimation comparable to ∼0.42 kcal/mol standard deviation of in vitro experiment. The speed of the method and lack of tunable parameters makes it affordable in current drug research.
Collapse
Affiliation(s)
- Timofey V Losev
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, 119991 Moscow, Russian Federation.,A.N. Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Vavilov Str. 28, 119991 Moscow, Russian Federation
| | - Igor S Gerasimov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Maria V Panova
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Alexey A Lisov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Yana R Abdyusheva
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| | - Polina V Rusina
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Eugenia Zaletskaya
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| | - Oleg V Stroganov
- BioMolTech Corp., 226 York Mills Rd, Toronto, Ontario M2L 1L1, Canada
| | - Michael G Medvedev
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Fedor N Novikov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| |
Collapse
|
3
|
Bieniek M, Wade AD, Bhati AP, Wan S, Coveney PV. TIES 2.0: A Dual-Topology Open Source Relative Binding Free Energy Builder with Web Portal. J Chem Inf Model 2023; 63:718-724. [PMID: 36719676 PMCID: PMC9930115 DOI: 10.1021/acs.jcim.2c01596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Relative binding free energy (RBFE) calculations are widely used to aid the process of drug discovery. TIES, Thermodynamic Integration with Enhanced Sampling, is a dual-topology approach to RBFE calculations with support for NAMD and OpenMM molecular dynamics engines. The software has been thoroughly validated on publicly available datasets. Here we describe the open source software along with a web portal (https://ccs-ties.org) that enables users to perform such calculations correctly and rapidly.
Collapse
Affiliation(s)
- Mateusz
K. Bieniek
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom,School
of Natural and Environmental Sciences, Newcastle
University, Newcastle upon Tyne NE1 7RU, United
Kingdom
| | - Alexander D. Wade
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Agastya P. Bhati
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Shunzhou Wan
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Peter V. Coveney
- Centre
for Computational Science, Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom,Advanced
Research Computing Centre, University College
London, London WC1H 0AJ, United
Kingdom,Institute
for Informatics, Faculty of Science, University
of Amsterdam, 1098XH Amsterdam, The Netherlands,E-mail:
| |
Collapse
|
4
|
Miller M, Rossetti T, Ferreira J, Ghanem L, Balbach M, Kaur N, Levin LR, Buck J, Kehr M, Coquille S, van den Heuvel J, Steegborn C, Fushimi M, Finkin-Groner E, Myers RW, Kargman S, Liverton NJ, Huggins DJ, Meinke PT. Design, Synthesis, and Pharmacological Evaluation of Second-Generation Soluble Adenylyl Cyclase (sAC, ADCY10) Inhibitors with Slow Dissociation Rates. J Med Chem 2022; 65:15208-15226. [PMID: 36346696 PMCID: PMC9866367 DOI: 10.1021/acs.jmedchem.2c01133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Soluble adenylyl cyclase (sAC: ADCY10) is an enzyme involved in intracellular signaling. Inhibition of sAC has potential therapeutic utility in a number of areas. For example, sAC is integral to successful male fertility: sAC activation is required for sperm motility and ability to undergo the acrosome reaction, two processes central to oocyte fertilization. Pharmacologic evaluation of existing sAC inhibitors for utility as on-demand, nonhormonal male contraceptives suggested that both high intrinsic potency, fast on and slow dissociation rates are essential design elements for successful male contraceptive applications. During the course of the medicinal chemistry campaign described here, we identified sAC inhibitors that fulfill these criteria and are suitable for in vivo evaluation of diverse sAC pharmacology.
Collapse
Affiliation(s)
- Michael Miller
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Thomas Rossetti
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Jacob Ferreira
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Lubna Ghanem
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Melanie Balbach
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Navpreet Kaur
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Lonny R. Levin
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Jochen Buck
- Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| | - Maria Kehr
- Department of Biochemistry, University of Bayreuth, 95440 Bayreuth, Germany
| | - Sandrine Coquille
- Department of Biochemistry, University of Bayreuth, 95440 Bayreuth, Germany
| | - Joop van den Heuvel
- Helmholtz Centre for Infection Research, Recombinant Protein Expression, 38124 Braunschweig, Germany
| | - Clemens Steegborn
- Department of Biochemistry, University of Bayreuth, 95440 Bayreuth, Germany
| | - Makoto Fushimi
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Efrat Finkin-Groner
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Robert W. Myers
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Stacia Kargman
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - Nigel J. Liverton
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States
| | - David J. Huggins
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States; Department of Physiology and Biophysics, Weill Cornell Medicine, New York, New York 10021, United States
| | - Peter T. Meinke
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States; Department of Pharmacology, Weill Cornell Medicine, New York, New York 10021, United States
| |
Collapse
|
5
|
Abstract
Positional analogue scanning (PAS) is an accepted strategy for multiparameter lead optimization (MPO) in drug discovery. Small structural changes as introduced by PAS can lead to 10-fold changes in binding potency in ∼10-20% of cases, a significant parameter shift irrespective of other MPO objectives. Sometimes performing a complete PAS is challenging due to resource and time constraints, building block availability, or difficulty in synthesis. Calculating relative binding free energies (RBFEs) for all positions can contribute to prioritizing the most promising analogues for synthesis. We tested a well-established RBFE calculation method, Amber GPU-TI, for 20 positional analogue scans in 14 test systems (cyclin-dependent kinase 8 (CDK8), hepatitis C virus nonstructural protein 5B (HCV NS5B), tankyrase, RAC-α serine/threonine-protein kinase (Akt), phosphodiesterase 1B (PDE1B), orexin/hypocretin receptor type 1 (OX1R), orexin/hypocretin receptor type 2 (OX2R), histone acetyltransferase K (lysine) acetyltransferase 6A (KAT6A), peroxisome proliferator-activated receptor γ (PPARγ), extracellular signal-regulated kinases (ERK1/2), coactivator-associated arginine methyltransferase 1 (PRMT4), αvβ6, bromodomain 1 (BD1), human immunodeficiency virus-1 (HIV-1) entry) involving nitrogen, methyl, halogen, methoxy, and hydroxyl scans with at least four analogues per set. Among the 66 analogue positions explored, we found that in 18 cases Amber GPU-TI calculations predicted a more than 10-fold change in potency. In all of these cases, the experimentally observed direction of potency changes agreed with the predictions. In 16 cases, more than 10-fold changes in experimental potency were observed. Again, in all of these cases, Amber GPU-TI predicted the direction of the potency changes correctly. In none of these cases would a decision made for or against synthesis based on a 10-fold change in potency have resulted in missing an important analogue. Therefore, in silico RBFE calculations using Amber GPU-TI can meaningfully contribute to the prioritization of positional analogues before synthesis.
Collapse
Affiliation(s)
- Yuan Hu
- Alkermes, Inc., 852 Winter Street, Waltham, Massachusetts 02451-1420, United States
| | - Ingo Muegge
- Alkermes, Inc., 852 Winter Street, Waltham, Massachusetts 02451-1420, United States
| |
Collapse
|
6
|
Huggins DJ. Comparing the Performance of Different AMBER Protein Forcefields, Partial Charge Assignments, and Water Models for Absolute Binding Free Energy Calculations. J Chem Theory Comput 2022; 18:2616-2630. [PMID: 35266690 DOI: 10.1021/acs.jctc.1c01208] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Identifying chemical starting points is a vital first step in small molecule drug discovery and can take significant time and money. For this reason, computational approaches to virtual screening are of great interest as they can lower the cost and shorten timeframes. However, simple approaches such as molecular docking and pharmacophore screening are of limited accuracy and provide a low probability of success. Alchemical binding free energies represent a promising approach for virtual screening as they naturally incorporate the key effects of water molecules, protein flexibility, and binding entropy. However, the calculations are technically very challenging, with performance depending on the specific forcefield used. For this reason, it is important that the community has access to benchmark test sets to assess prediction accuracy. In this paper, we present an approach to alchemical binding free energies using OpenMM. We identify effective simulation parameters using an existing BRD4(1) test set and present two new benchmark sets (cMET and PDE2A) that can be used in the community for validation purposes. Our findings also highlight the effectiveness of some AMBER forcefields, in particular, AMBER ff15ipq.
Collapse
Affiliation(s)
- David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, New York, New York 10021, United States.,Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| |
Collapse
|
7
|
Titov IY, Stroylov VS, Rusina P, Svitanko IV. Preliminary modelling as the first stage of targeted organic synthesis. RUSSIAN CHEMICAL REVIEWS 2021. [DOI: 10.1070/rcr5012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The review aims to present a classification and applicability analysis of methods for preliminary molecular modelling for targeted organic, catalytic and biocatalytic synthesis. The following three main approaches are considered as a primary classification of the methods: modelling of the target – ligand coordination without structural information on both the target and the resulting complex; calculations based on experimentally obtained structural information about the target; and dynamic simulation of the target – ligand complex and the reaction mechanism with calculation of the free energy of the reaction. The review is meant for synthetic chemists to be used as a guide for building an algorithm for preliminary modelling and synthesis of structures with specified properties.
The bibliography includes 353 references.
Collapse
|
8
|
Abula A, Xu Z, Zhu Z, Peng C, Chen Z, Zhu W, Aisa HA. Substitution Effect of the Trifluoromethyl Group on the Bioactivity in Medicinal Chemistry: Statistical Analysis and Energy Calculations. J Chem Inf Model 2020; 60:6242-6250. [DOI: 10.1021/acs.jcim.0c00898] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Amina Abula
- The Key Laboratory of Plant Resources and Chemistry in Arid Regions and Key Laboratory of Xinjiang Indigenous Medicinal Plants Resource Utilization, Xinjiang Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, South Beijing Road 40-1, Urumqi, Xinjiang 830011, People’s Republic of China
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Zhijian Xu
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Zhengdan Zhu
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Cheng Peng
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Zhaoqiang Chen
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Weiliang Zhu
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Haji Akber Aisa
- The Key Laboratory of Plant Resources and Chemistry in Arid Regions and Key Laboratory of Xinjiang Indigenous Medicinal Plants Resource Utilization, Xinjiang Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, South Beijing Road 40-1, Urumqi, Xinjiang 830011, People’s Republic of China
- University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing 100049, People’s Republic of China
| |
Collapse
|
9
|
Wade AD, Huggins DJ. Identification of Optimal Ligand Growth Vectors Using an Alchemical Free-Energy Method. J Chem Inf Model 2020; 60:5580-5594. [PMID: 32810401 DOI: 10.1021/acs.jcim.0c00610] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, a novel method to rationally design inhibitors with improved steric contacts and enhanced binding free energies is presented. This new method uses alchemical single step perturbation calculations to rapidly optimize the van der Waals interactions of a small molecule in a protein-ligand complex in order to maximize its binding affinity. The results of the optimizer are used to predict beneficial growth vectors on the ligand, and good agreement is found between the predictions from the optimizer and a more rigorous free energy calculation, with a Spearman's rank order correlation of 0.59. The advantage of the method presented here is the significant speed up of over 10-fold compared to traditional free energy calculations and sublinear scaling with the number of growth vectors assessed. Where experimental data were available, mutations from hydrogen to a methyl group at sites highlighted by the optimizer were calculated with MBAR, and the mean unsigned error between experimental and calculated values of the binding free energy was 0.83 kcal/mol.
Collapse
Affiliation(s)
- Alexander D Wade
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, Belfer Research Building, 413 East 69th Street, 16th Floor, Box 300, New York, United States.,Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| |
Collapse
|
10
|
|
11
|
Mei H, Han J, White S, Graham DJ, Izawa K, Sato T, Fustero S, Meanwell NA, Soloshonok VA. Tailor-Made Amino Acids and Fluorinated Motifs as Prominent Traits in Modern Pharmaceuticals. Chemistry 2020; 26:11349-11390. [PMID: 32359086 DOI: 10.1002/chem.202000617] [Citation(s) in RCA: 69] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Revised: 04/23/2020] [Indexed: 12/11/2022]
Abstract
Structural analysis of modern pharmaceutical practices allows for the identification of two rapidly growing trends: the introduction of tailor-made amino acids and the exploitation of fluorinated motifs. Curiously, the former represents one of the most ubiquitous classes of naturally occurring compounds, whereas the latter is the most xenobiotic and comprised virtually entirely of man-made derivatives. Herein, 39 selected compounds, featuring both of these traits in the same molecule, are profiled. The total synthesis, source of the corresponding amino acids and fluorinated residues, and medicinal chemistry aspects and biological properties of the molecules are discussed.
Collapse
Affiliation(s)
- Haibo Mei
- Jiangsu Co-Innovation Center of Efficient Processing and Utilization of Forest Resources, College of Chemical Engineering, Nanjing Forestry University, Nanjing, 210037, P.R. China
| | - Jianlin Han
- Jiangsu Co-Innovation Center of Efficient Processing and Utilization of Forest Resources, College of Chemical Engineering, Nanjing Forestry University, Nanjing, 210037, P.R. China
| | - Sarah White
- Oakwood Chemical, Inc., 730 Columbia Hwy. N, Estill, SC, 29918, USA
| | - Daniel J Graham
- Oakwood Chemical, Inc., 730 Columbia Hwy. N, Estill, SC, 29918, USA
| | - Kunisuke Izawa
- Hamari Chemicals Ltd., 1-4-29 Kunijima, Higashi-Yodogawa-ku, Osaka, 533-0024, Japan
| | - Tatsunori Sato
- Hamari Chemicals Ltd., 1-4-29 Kunijima, Higashi-Yodogawa-ku, Osaka, 533-0024, Japan
| | - Santos Fustero
- Departamento de Química Orgánica, Universidad de Valencia, 46100, Burjassot, Valencia, Spain
| | - Nicholas A Meanwell
- Department of Small Molecule Drug Discovery, Bristol-Myers Squibb Research and Development, P.O. Box 4000, Princeton, NJ, 08543-4000, USA
| | - Vadim A Soloshonok
- Department of Organic Chemistry I, Faculty of Chemistry, University of the Basque Country UPV/EHU, Paseo Manuel Lardizábal 3, 20018, San Sebastián, Spain.,IKERBASQUE, Basque Foundation for Science, María Díaz de Haro 3, Plaza Bizkaia, 48013, Bilbao, Spain
| |
Collapse
|
12
|
de Ruiter A, Oostenbrink C. Advances in the calculation of binding free energies. Curr Opin Struct Biol 2020; 61:207-212. [PMID: 32088376 DOI: 10.1016/j.sbi.2020.01.016] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2019] [Revised: 01/21/2020] [Accepted: 01/24/2020] [Indexed: 01/19/2023]
Abstract
In recent years, calculations of binding affinities from molecular simulations seem to have matured significantly. While the number of applications of such methods in drug design and biotechnology increases, the number of truly new methodological developments decreases. This review provides an overview of the current status of the field as reflected in recent publications. The focus is on the challenges that remain when using endstate, alchemical and pathway methods. For endstate methods this is the calculation of entropic contributions. For alchemical methods there are unsolved problems associated with the solvation of the active site, sampling slow degrees of freedom and when modifying the net charge. For pathway methods achieving sufficient sampling remains challenging. New trends are also highlighted, including the use of pathway methods for the quantification of protein-protein interactions.
Collapse
Affiliation(s)
- Anita de Ruiter
- Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences (BOKU), Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences (BOKU), Vienna, Austria.
| |
Collapse
|
13
|
Wade AD, Huggins DJ. Optimization of Protein-Ligand Electrostatic Interactions Using an Alchemical Free-Energy Method. J Chem Theory Comput 2019; 15:6504-6512. [PMID: 31584802 PMCID: PMC7007198 DOI: 10.1021/acs.jctc.9b00976] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
![]()
We present an explicit solvent alchemical free-energy method for
optimizing the partial charges of a ligand to maximize the binding
affinity with a receptor. This methodology can be applied to known
ligand–protein complexes to determine an optimized set of ligand
partial atomic changes. Three protein–ligand complexes have
been optimized in this work: FXa, P38, and the androgen receptor.
The sets of optimized charges can be used to identify design principles
for chemical changes to the ligands which improve the binding affinity
for all three systems. In this work, beneficial chemical mutations
are generated from these principles and the resulting molecules tested
using free-energy perturbation calculations. We show that three quarters
of our chemical changes are predicted to improve the binding affinity,
with an average improvement for the beneficial mutations of approximately
1 kcal/mol. In the cases where experimental data are available, the
agreement between prediction and experiment is also good. The results
demonstrate that charge optimization in explicit solvent is a useful
tool for predicting beneficial chemical changes such as pyridinations,
fluorinations, and oxygen to sulfur mutations.
Collapse
Affiliation(s)
- Alexander D Wade
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , U.K
| | - David J Huggins
- TCM Group, Cavendish Laboratory , University of Cambridge , 19 J J Thomson Avenue , Cambridge CB3 0HE , U.K.,Tri-Institutional Therapeutics Discovery Institute , Belfer Research Building, 413 East 69th Street, 16th Floor, Box 300 , New York 10021 , United States.,Department of Physiology and Biophysics , Weill Cornell Medicine , 1300 York Avenue , New York , New York 10065 , United States
| |
Collapse
|