1
|
Zhao M, Yu W, MacKerell AD. Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms. J Phys Chem B 2024; 128:7362-7375. [PMID: 39031121 PMCID: PMC11294009 DOI: 10.1021/acs.jpcb.4c03045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2024]
Abstract
In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.
Collapse
Affiliation(s)
- Mingtian Zhao
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| |
Collapse
|
2
|
Giugliano G, Gajo M, Marforio TD, Zerbetto F, Mattioli EJ, Calvaresi M. Identification of Potential Drug Targets of Calix[4]arene by Reverse Docking. Chemistry 2024; 30:e202400871. [PMID: 38777795 DOI: 10.1002/chem.202400871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 05/20/2024] [Accepted: 05/22/2024] [Indexed: 05/25/2024]
Abstract
Calixarenes are displaying great potential for the development of new drug delivery systems, diagnostic imaging, biosensing devices and inhibitors of biological processes. In particular, calixarene derivatives are able to interact with many different enzymes and function as inhibitors. By screening of the potential drug target database (PDTD) with a reverse docking procedure, we identify and discuss a selection of 100 proteins that interact strongly with calix[4]arene. We also discover that leucine (23.5 %), isoleucine (11.3 %), phenylalanines (11.3 %) and valine (9.5 %) are the most frequent binding residues followed by hydrophobic cysteines and methionines and aromatic histidines, tyrosines and tryptophanes. Top binders are peroxisome proliferator-activated receptors that already are targeted by commercial drugs, demonstrating the practical interest in calix[4]arene. Nuclear receptors, potassium channel, several carrier proteins, a variety of cancer-related proteins and viral proteins are prominent in the list. It is concluded that calix[4]arene, which is characterized by facile access, well-defined conformational characteristics, and ease of functionalization at both the lower and higher rims, could be a potential lead compound for the development of enzyme inhibitors and theranostic platforms.
Collapse
Affiliation(s)
- Giulia Giugliano
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| | - Margherita Gajo
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| | - Tainah Dorina Marforio
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| | - Francesco Zerbetto
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| | - Edoardo Jun Mattioli
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| | - Matteo Calvaresi
- Dipartimento di Chimica "Giacomo Ciamician", Alma Mater Studiorum - Università di Bologna, Via Francesco Selmi 2, 40126, Bologna, Italy E-Mail
| |
Collapse
|
3
|
Hu R, Zhang J, Kang Y, Wang Z, Pan P, Deng Y, Hsieh CY, Hou T. Comprehensive, Open-Source, and Automated Workflow for Multisite λ-Dynamics in Lead Optimization. J Chem Theory Comput 2024; 20:1465-1478. [PMID: 38300792 DOI: 10.1021/acs.jctc.3c01154] [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: 02/03/2024]
Abstract
Multisite λ-dynamics (MSLD) is a highly efficient binding free energy calculation method that samples multiple ligands in a single round by assigning different λ values to the alchemical part of each ligand. This method holds great promise for lead optimization (LO) in drug discovery. However, the complex data preparation and simulation process limits its widespread application in diverse protein-ligand systems. To address this challenge, we developed a comprehensive, open-source, and automated workflow for MSLD calculations based on the BLaDE dynamics engine. This workflow incorporates the Ligand Internal and Cartesian coordinate reconstruction-based alignment algorithm (LIC-align) and an optimized maximum common substructure (MCS) search algorithm to accurately generate MSLD multiple topologies with ideal perturbation patterns. Furthermore, our workflow is highly modularized, allowing straightforward integration and extension of various simulation techniques, and is highly accessible to nonexperts. This workflow was validated by calculating the relative binding free energies of large-scale congeneric ligands, many of which have large perturbing groups. The agreement between the calculations and experiments was excellent, with an average unsigned error of 1.08 ± 0.47 kcal/mol. More than 57.1% of the ligands had an error of less than 1.0 kcal/mol, and the perturbations of 6 targets were fully connected via the calculations, while those of 2 targets were connected via both calculations and experimental data. The Pearson correlation coefficient reached 0.88, indicating that the MSLD workflow provides accurate predictions that can guide lead optimization in drug discovery. We also examined the impact of single-site versus multisite perturbations, ligand grouping by perturbing group size, and the position of the anchor atom on the MSLD performance. By integrating our proposed LIC-align and optimized MCS search algorithm along with the coping strategies to handle challenging molecular substructures, our workflow can handle many realistic scenarios more reasonably than all previously published methods. Moreover, we observed that our MSLD workflow achieved similar accuracy to free energy perturbation (FEP) while improving computational efficiency by over 1 order of magnitude in speedup. These findings provide valuable insights and strategies for further MSLD development, making MSLD a competitive tool for lead optimization.
Collapse
Affiliation(s)
- Renling Hu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
- Polytechnic Institute, Zhejiang University, Hangzhou 310058, Zhejiang, China
- CarbonSilicon AI Technology Co., Ltd., Hangzhou 310018, Zhejiang, China
| | - Jintu Zhang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yu Kang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Zhe Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Peichen Pan
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yafeng Deng
- CarbonSilicon AI Technology Co., Ltd., Hangzhou 310018, Zhejiang, China
| | - Chang-Yu Hsieh
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| |
Collapse
|
4
|
Liu R, Li W, Yao Y, Wu Y, Luo HB, Li Z. Accelerating and Automating the Free Energy Perturbation Absolute Binding Free Energy Calculation with the RED-E Function. J Chem Inf Model 2023; 63:7755-7767. [PMID: 38048439 DOI: 10.1021/acs.jcim.3c01670] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/06/2023]
Abstract
The accurate prediction of the binding affinities between small molecules and biological macromolecules plays a fundamental role in structure-based drug design, which is still challenging. The free energy perturbation-based absolute binding free energy (FEP-ABFE) approach has shown potential in its reliability. To correctly calculate the energy related to the ligand being restrained by the receptor, additional restraints between the ligand and the receptor are needed. However, determining the restraint parameters for individual ligands empirically is too trivial to be automated, and usually gives rise to numerical instabilities, which set back the applications of FEP-ABFE. To address these issues, we derived the analytical expression for the probability distribution of energy differences, P(ΔU), during the process of restraint addition, which is called the RED-E (restraint energy distribution at equilibrium position) function. Simulations indicated that the RED-E function can accurately describe P(ΔU) when restraints are added at the equilibrium position. Based on the RED-E function, an automatic restraint selection method was proposed to select the best restraint. With this method, there is a high phase-space overlap between the free and restrained states, such that using a 2-λ perturbation can accurately calculate the free energy of the restraint addition, which is a nearly 6 times acceleration compared with current widely used 12-λ perturbation method. The RED-E function gives insight into the non-Gaussian behavior of the sampled P(ΔU) in certain FEP processes in an analytical way. The highly automated and accelerated restraint selection also makes it possible for the large-scale application of FEP-ABFE in real drug discovery practices.
Collapse
Affiliation(s)
- Runduo Liu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Wenchao Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Yufen Yao
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Yinuo Wu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Hai-Bin Luo
- Key Laboratory of Tropical Biological Resources of Ministry of Education, School of Pharmaceutical Sciences, Hainan University, Haikou, Hainan 570228, China
- Song Li' Academician Workstation of Hainan University (School of Pharmaceutical Sciences), Yazhou Bay, Sanya 572000, China
| | - Zhe Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| |
Collapse
|
5
|
Molani F, Webb S, Cho AE. Combining QM/MM Calculations with Classical Mining Minima to Predict Protein-Ligand Binding Free Energy. J Chem Inf Model 2023; 63:2728-2734. [PMID: 37079618 DOI: 10.1021/acs.jcim.2c01637] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/21/2023]
Abstract
We developed an effective binding free energy prediction protocol which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations to substitute the specified atomic charges of force fields with quantum-mechanically recalculated ones at a proposed pose using a mining minima approach with the VeraChem mining minima engine. We tested this protocol using seven well-known targets with 147 different ligands and compared it with classical mining minima and the most popular binding free energy (BFE) methods using different metrics. Our new protocol, dubbed Qcharge-VM2, yielded an overall Pearson correlation of 0.86, which was better than all the methods examined. Qcharge-VM2 performed significantly better than implicit solvent-based methods, such as MM-GBSA and MM-PBSA, but not as good as explicit water-based free energy perturbation methods, such as FEP+, in terms of root-mean-square error, RMSE (1.75 kcal/mol) and mean unsigned error, MUE (1.39 kcal/mol) on a limited set of targets. However, our protocol is substantially less computationally demanding compared with FEP+. The combined accuracy and efficiency of our method can be valuable in drug discovery campaigns.
Collapse
Affiliation(s)
- Farzad Molani
- Department of Bioinformatics, Korea University, 2511 Sejong-ro, Sejong 30119, Korea
| | - Simon Webb
- VeraChem LLC, 12850 Middlebrook Road STE 205, Germantown, Maryland 20874, United States
| | - Art E Cho
- Department of Bioinformatics, Korea University, 2511 Sejong-ro, Sejong 30119, Korea
| |
Collapse
|
6
|
Lee TS, Tsai HC, Ganguly A, York DM. ACES: Optimized Alchemically Enhanced Sampling. J Chem Theory Comput 2023; 19:10.1021/acs.jctc.2c00697. [PMID: 36630672 PMCID: PMC10333454 DOI: 10.1021/acs.jctc.2c00697] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We present an alchemical enhanced sampling (ACES) method implemented in the GPU-accelerated AMBER free energy MD engine. The methods hinges on the creation of an "enhanced sampling state" by reducing or eliminating selected potential energy terms and interactions that lead to kinetic traps and conformational barriers while maintaining those terms that curtail the need to otherwise sample large volumes of phase space. For example, the enhanced sampling state might involve transforming regions of a ligand and/or protein side chain into a noninteracting "dummy state" with internal electrostatics and torsion angle terms turned off. The enhanced sampling state is connected to a real-state end point through a Hamiltonian replica exchange (HREMD) framework that is facilitated by newly developed alchemical transformation pathways and smoothstep softcore potentials. This creates a counterdiffusion of real and enhanced-sampling states along the HREMD network. The effect of a differential response of the environment to the real and enhanced-sampling states is minimized by leveraging the dual topology framework in AMBER to construct a counterbalancing HREMD network in the opposite alchemical direction with the same (or similar) real and enhanced sampling states at inverted end points. The method has been demonstrated in a series of test cases of increasing complexity where traditional MD, and in several cases alternative REST2-like enhanced sampling methods, are shown to fail. The hydration free energy for acetic acid was shown to be independent of the starting conformation, and the values for four additional edge case molecules from the FreeSolv database were shown to have a significantly closer agreement with experiment using ACES. The method was further able to handle different rotamer states in a Cdk2 ligand identified as fractionally occupied in crystal structures. Finally, ACES was applied to T4-lysozyme and demonstrated that the side chain distribution of V111χ1 could be reliably reproduced for the apo state, bound to p-xylene, and in p-xylene→ benzene transformations. In these cases, the ACES method is shown to be highly robust and superior to a REST2-like enhanced sampling implementation alone.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Hsu-Chun Tsai
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Abir Ganguly
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
7
|
Ganguly A, Tsai HC, Fernández-Pendás M, Lee TS, Giese TJ, York DM. AMBER Drug Discovery Boost Tools: Automated Workflow for Production Free-Energy Simulation Setup and Analysis (ProFESSA). J Chem Inf Model 2022; 62:6069-6083. [PMID: 36450130 PMCID: PMC9881431 DOI: 10.1021/acs.jcim.2c00879] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
We report an automated workflow for production free-energy simulation setup and analysis (ProFESSA) using the GPU-accelerated AMBER free-energy engine with enhanced sampling features and analysis tools, part of the AMBER Drug Discovery Boost package that has been integrated into the AMBER22 release. The workflow establishes a flexible, end-to-end pipeline for performing alchemical free-energy simulations that brings to bear technologies, including new enhanced sampling features and analysis tools, to practical drug discovery problems. ProFESSA provides the user with top-level control of large sets of free-energy calculations and offers access to the following key functionalities: (1) automated setup of file infrastructure; (2) enhanced conformational and alchemical sampling with the ACES method; and (3) network-wide free-energy analysis with the optional imposition of cycle closure and experimental constraints. The workflow is applied to perform absolute and relative solvation free-energy and relative ligand-protein binding free-energy calculations using different atom-mapping procedures. Results demonstrate that the workflow is internally consistent and highly robust. Further, the application of a new network-wide Lagrange multiplier constraint analysis that imposes key experimental constraints substantially improves binding free-energy predictions.
Collapse
Affiliation(s)
- Abir Ganguly
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Hsu-Chun Tsai
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Mario Fernández-Pendás
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA,Donostia International Physics Center (DIPC), PK 1072, 20080 Donostia-San Sebastian, Spain
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA,
| |
Collapse
|
8
|
Wang S, Sun X, Cui W, Yuan S. MM/PB(GB)SA benchmarks on soluble proteins and membrane proteins. Front Pharmacol 2022; 13:1018351. [PMID: 36532746 PMCID: PMC9751045 DOI: 10.3389/fphar.2022.1018351] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Accepted: 11/17/2022] [Indexed: 08/27/2023] Open
Abstract
Predicting protein-ligand binding free energy rapidly and accurately remains a challenging question in modern drug discovery. Molecular mechanics/Poisson-Boltzmann (Generalized Born) surface area (MM/PB(GB)SA) has emerged as an essential tool for accelerating cost-efficient binding free energy calculation. This study presents benchmarks with three membrane-bound protein systems and six soluble protein systems. Different parameters were sampled for different benchmarks to explore the highest accuracy. These include ligand charges, protein force fields, extra points, GB models, nonpolar optimization methods, internal dielectric constants and membrane dielectric constants. Comparisons of accuracy were made between MM/PB(GB)SA, docking and free energy perturbation (FEP). The results reveal a competitive performance between MM/PB(GB)SA and FEP. In summary, MM/PB(GB)SA is a powerful approach to predict ligand binding free energy rapidly and accurately. Parameters of MM/PB(GB)SA calculations, such as the GB models and membrane dielectric constants, need to be optimized for different systems. This method can be served as a powerful tool for drug design.
Collapse
Affiliation(s)
- Shiyu Wang
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- College of Chemical Science, University of Chinese Academy of Sciences, Beijing, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Xiaolin Sun
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Wenqiang Cui
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Shuguang Yuan
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland
| |
Collapse
|
9
|
Applications and mechanisms of the cyclin-dependent kinase 4/6 inhibitor, PD-0332991, in solid tumors. Cell Oncol (Dordr) 2022; 45:1053-1071. [PMID: 36087253 DOI: 10.1007/s13402-022-00714-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/01/2022] [Indexed: 01/10/2023] Open
Abstract
Abnormal CDK4/6-Rb-E2F signal transduction is a common finding in tumors and is a driving factor for the excessive proliferation of various tumor cells. PD-0332991, a highly specific, small molecule inhibitor for CDK4 and 6, has been shown to inhibit tumor growth by abrogating the phosphorylating capacity of CDK4/6 and suppressing Rb phosphorylation. It has been promoted for the treatment of breast cancer and potentially for other tumor types such as liver cancers, lung cancers and sarcomas. Due to the risk of monotherapy resistance, PD-0332991 is commonly used in combination with other drugs. Such combination treatments have proved able to inhibit tumor proliferation more effectively, induce stronger senescence and apoptosis, and enhance the efficiency of immunotherapy. Therefore, tumor cells with senescence induced by PD-0332991 are now used as ideal screening tools of cytolytic drugs with more efficient and thorough anti-tumor properties. With more extensive understandings about the branching points between senescence and apoptosis, it is possible to refine the dosage of PD-0332991. Better characterization of resistant cells, of inhibitors and of adverse effects such as leukopenia are needed to overcome obstacles in the use of PD-0332991. In this review of PD-0332991 research, we hope to provide guidance of transitions from laboratory findings to clinical applications of PD-0332991 and to facilitate PD-0332991-based multi-inhibitor combination therapies for various tumors.
Collapse
|
10
|
Sokolsky A, Winterton S, Kennedy K, Drake K, Stump K, Huo L, Lo Y, Ye M, Covington M, Diamond S, Yang YO, Kim S, Yeleswaram S, Wu L, Yao W. Discovery of 5,7-Dihydro-6 H-pyrrolo[2,3- d]pyrimidin-6-ones as Highly Selective CDK2 Inhibitors. ACS Med Chem Lett 2022; 13:1797-1804. [PMID: 36385925 PMCID: PMC9661707 DOI: 10.1021/acsmedchemlett.2c00408] [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] [Received: 08/30/2022] [Accepted: 10/03/2022] [Indexed: 11/30/2022] Open
Abstract
A series of exceptionally selective CDK2 inhibitors are described. Starting from an HTS hit, we successfully scaffold hopped to a 5,7-dihydro-6H-pyrrolo[2,3-d]pyrimidin-6-one core structure, which imparted a promising initial selectivity within the CDK family. Extensive further SAR identified additional factors that drove selectivity to above 200× for CDKs 1/4/6/7/9. General kinome selectivity was also greatly improved. Finally, use of in vivo metabolite identification allowed us to pinpoint sulfonamide dealkylation as the primary metabolite, which was ameliorated through the deuterium effect.
Collapse
Affiliation(s)
- Alexander Sokolsky
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Sarah Winterton
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Keith Kennedy
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Katherine Drake
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Kristine Stump
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Lu Huo
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Yvonne Lo
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Min Ye
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Maryanne Covington
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Sharon Diamond
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Yan-ou Yang
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Sunkyu Kim
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Swamy Yeleswaram
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Liangxing Wu
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| | - Wenqing Yao
- Incyte
Research Institute, Incyte Corporation, 1801 Augustine Cut-off, Wilmington, Delaware 19803, United States
| |
Collapse
|
11
|
Hahn DF, Bayly CI, Boby ML, Macdonald HEB, Chodera JD, Gapsys V, Mey ASJS, Mobley DL, Benito LP, Schindler CEM, Tresadern G, Warren GL. Best practices for constructing, preparing, and evaluating protein-ligand binding affinity benchmarks [Article v0.1]. LIVING JOURNAL OF COMPUTATIONAL MOLECULAR SCIENCE 2022; 4:1497. [PMID: 36382113 PMCID: PMC9662604 DOI: 10.33011/livecoms.4.1.1497] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/01/2023]
Abstract
Free energy calculations are rapidly becoming indispensable in structure-enabled drug discovery programs. As new methods, force fields, and implementations are developed, assessing their expected accuracy on real-world systems (benchmarking) becomes critical to provide users with an assessment of the accuracy expected when these methods are applied within their domain of applicability, and developers with a way to assess the expected impact of new methodologies. These assessments require construction of a benchmark-a set of well-prepared, high quality systems with corresponding experimental measurements designed to ensure the resulting calculations provide a realistic assessment of expected performance when these methods are deployed within their domains of applicability. To date, the community has not yet adopted a common standardized benchmark, and existing benchmark reports suffer from a myriad of issues, including poor data quality, limited statistical power, and statistically deficient analyses, all of which can conspire to produce benchmarks that are poorly predictive of real-world performance. Here, we address these issues by presenting guidelines for (1) curating experimental data to develop meaningful benchmark sets, (2) preparing benchmark inputs according to best practices to facilitate widespread adoption, and (3) analysis of the resulting predictions to enable statistically meaningful comparisons among methods and force fields. We highlight challenges and open questions that remain to be solved in these areas, as well as recommendations for the collection of new datasets that might optimally serve to measure progress as methods become systematically more reliable. Finally, we provide a curated, versioned, open, standardized benchmark set adherent to these standards (PLBenchmarks) and an open source toolkit for implementing standardized best practices assessments (arsenic) for the community to use as a standardized assessment tool. While our main focus is free energy methods based on molecular simulations, these guidelines should prove useful for assessment of the rapidly growing field of machine learning methods for affinity prediction as well.
Collapse
Affiliation(s)
- David F. Hahn
- Computational Chemistry,Janssen Research & Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | | | - Melissa L. Boby
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065 USA
| | - Hannah E. Bruce Macdonald
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065 USA
- MSD R&D Innovation Centre, 120 Moorgate, London EC2M 6UR, United Kingdom
| | - John D. Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065 USA
| | - Vytautas Gapsys
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Antonia S. J. S. Mey
- EaStCHEM School of Chemistry, David Brewster Road, Joseph Black Building, The King’s Buildings, Edinburgh, EH9 3FJ, UK
| | - David L. Mobley
- Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine, CA USA
| | - Laura Perez Benito
- Computational Chemistry,Janssen Research & Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | | | - Gary Tresadern
- Computational Chemistry,Janssen Research & Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | | |
Collapse
|
12
|
Mader L, Hayward JJ, Porter LA, Trant JF. A revised synthesis of 6-alkoxy-2-aminopurines with late-stage convergence allowing for increased molecular complexity. NEW J CHEM 2022. [DOI: 10.1039/d2nj02204d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This streamlined synthesis allows the alcohol at the 6-position of 6-alkoxy-2-arylaminopurines to be used only in moderate excess, rather than as solvent, opening up the possibility of accessing more complicated molecules.
Collapse
Affiliation(s)
- Lavleen Mader
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, ON, N9B 3P4, Canada
| | - John J. Hayward
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, ON, N9B 3P4, Canada
| | - Lisa A. Porter
- Department of Biomedical Sciences, University of Windsor, 401 Sunset Avenue, Windsor, ON, N9B 3P4, Canada
| | - John F. Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, ON, N9B 3P4, Canada
| |
Collapse
|
13
|
Faber EB, Wang N, Georg GI. Review of rationale and progress toward targeting cyclin-dependent kinase 2 (CDK2) for male contraception†. Biol Reprod 2021; 103:357-367. [PMID: 32543655 PMCID: PMC7523694 DOI: 10.1093/biolre/ioaa107] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2020] [Revised: 06/03/2020] [Accepted: 04/19/2020] [Indexed: 12/30/2022] Open
Abstract
Cyclin-dependent kinase 2 (CDK2) is a member of the larger cell cycle regulating CDK family of kinases, activated by binding partner cyclins as its name suggests. Despite its canonical role in mitosis, CDK2 knockout mice are viable but sterile, suggesting compensatory mechanisms for loss of CDK2 in mitosis but not meiosis. Here, we review the literature surrounding the role of CDK2 in meiosis, particularly a cyclin-independent role in complex with another activator, Speedy 1 (SPY1). From this evidence, we suggest that CDK2 could be a viable nonhormonal male contraceptive target. Finally, we review the literature of pertinent CDK2 inhibitors from the preclinical to clinical stages, mostly developed to treat various cancers. To date, there is no potent yet selective CDK2 inhibitor that could be repurposed as a contraceptive without appreciable off-target toxicity. To achieve selectivity for CDK2 over closely related kinases, developing compounds that bind outside the conserved adenosine triphosphate-binding site may be necessary.
Collapse
Affiliation(s)
- Erik B Faber
- Department of Medicinal Chemistry, College of Pharmacy, University of Minnesota-Twin Cities, Minneapolis, MN, USA.,Medical-Scientist Training Program, University of Minnesota Medical School, University of Minnesota-Twin Cities, Minneapolis, MN, USA
| | - Nan Wang
- Department of Medicinal Chemistry, College of Pharmacy, University of Minnesota-Twin Cities, Minneapolis, MN, USA
| | - Gunda I Georg
- Department of Medicinal Chemistry, College of Pharmacy, University of Minnesota-Twin Cities, Minneapolis, MN, USA
| |
Collapse
|
14
|
Talapati SR, Goyal M, Nataraj V, Pothuganti M, R SM, Gore S, Ramachandra M, Antony T, More SS, Rao NK. Structural and binding studies of cyclin-dependent kinase 2 with NU6140 inhibitor. Chem Biol Drug Des 2021; 98:857-868. [PMID: 34423559 DOI: 10.1111/cbdd.13941] [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: 04/01/2021] [Revised: 07/21/2021] [Accepted: 08/16/2021] [Indexed: 11/30/2022]
Abstract
Cyclin-dependent kinase 2 (CDK2) is an established target protein for therapeutic intervention in various diseases, including cancer. Reported inhibitors of CDK2 target the ATP-binding pocket to inhibit the kinase activity. Many small molecule CDK2 inhibitors have been discovered, and their crystal structure with CDK2 or CDK2-cyclin A complex has been published. NU6140 is a CDK2 inhibitor with moderate potency and selectivity. Herein, we report the cocrystal structure determination of NU6140 in complex with CDK2 and confirmation of the binding using various biophysical methods. Our data show that NU6140 binds to CDK2 with a Kd of 800 nM as determined by SPR and stabilizes the protein against thermal denaturation (ΔTm -5°C). The cocrystal structure determined in our study shows that NU6140 binds in the ATP-binding pocket as expected for this class of compounds and interacts with Leu83 and Glu81 with regular hydrogen bonds and with Asp145 via water-mediated H-bond. Based on these data, we propose structural modifications of NU6140 to introduce new interactions with CDK2 that can improve its potency while retaining the selectivity.
Collapse
Affiliation(s)
- Sumalatha Rani Talapati
- Aurigene Discovery Technologies Ltd, Bangalore, India.,School of Basic and Applied Sciences, Dayananda Sagar University, Bangalore, India
| | - Megha Goyal
- Aurigene Discovery Technologies Ltd, Bangalore, India
| | | | | | - Sreevidya M R
- Aurigene Discovery Technologies Ltd, Bangalore, India
| | - Suraj Gore
- Aurigene Discovery Technologies Ltd, Bangalore, India
| | | | - Thomas Antony
- Aurigene Discovery Technologies Ltd, Bangalore, India
| | - Sunil S More
- School of Basic and Applied Sciences, Dayananda Sagar University, Bangalore, India
| | | |
Collapse
|
15
|
Goel H, Hazel A, Ustach VD, Jo S, Yu W, MacKerell AD. Rapid and accurate estimation of protein-ligand relative binding affinities using site-identification by ligand competitive saturation. Chem Sci 2021; 12:8844-8858. [PMID: 34257885 PMCID: PMC8246086 DOI: 10.1039/d1sc01781k] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 05/24/2021] [Indexed: 01/18/2023] Open
Abstract
Predicting relative protein-ligand binding affinities is a central pillar of lead optimization efforts in structure-based drug design. The site identification by ligand competitive saturation (SILCS) methodology is based on functional group affinity patterns in the form of free energy maps that may be used to compute protein-ligand binding poses and affinities. Presented are results obtained from the SILCS methodology for a set of eight target proteins as reported originally in Wang et al. (J. Am. Chem. Soc., 2015, 137, 2695-2703) using free energy perturbation (FEP) methods in conjunction with enhanced sampling and cycle closure corrections. These eight targets have been subsequently studied by many other authors to compare the efficacy of their method while comparing with the outcomes of Wang et al. In this work, we present results for a total of 407 ligands on the eight targets and include specific analysis on the subset of 199 ligands considered previously. Using the SILCS methodology we can achieve an average accuracy of up to 77% and 74% when considering the eight targets with their 199 and 407 ligands, respectively, for rank-ordering ligand affinities as calculated by the percent correct metric. This accuracy increases to 82% and 80%, respectively, when the SILCS atomic free energy contributions are optimized using a Bayesian Markov-chain Monte Carlo approach. We also report other metrics including Pearson's correlation coefficient, Pearlman's predictive index, mean unsigned error, and root mean square error for both sets of ligands. The results obtained for the 199 ligands are compared with the outcomes of Wang et al. and other published works. Overall, the SILCS methodology yields similar or better-quality predictions without a priori need for known ligand orientations in terms of the different metrics when compared to current FEP approaches with significant computational savings while additionally offering quantitative estimates of individual atomic contributions to binding free energies. These results further validate the SILCS methodology as an accurate, computationally efficient tool to support lead optimization and drug discovery.
Collapse
Affiliation(s)
- Himanshu Goel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy 20, Penn St. Baltimore Maryland 21201 USA +1-410-706-5017 +1-410-706-7442
| | - Anthony Hazel
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy 20, Penn St. Baltimore Maryland 21201 USA +1-410-706-5017 +1-410-706-7442
| | - Vincent D Ustach
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy 20, Penn St. Baltimore Maryland 21201 USA +1-410-706-5017 +1-410-706-7442
| | - Sunhwan Jo
- SilcsBio LLC 8 Market Place, Suite 300 Baltimore Maryland 21201 USA
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy 20, Penn St. Baltimore Maryland 21201 USA +1-410-706-5017 +1-410-706-7442
| | - Alexander D MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy 20, Penn St. Baltimore Maryland 21201 USA +1-410-706-5017 +1-410-706-7442
- SilcsBio LLC 8 Market Place, Suite 300 Baltimore Maryland 21201 USA
| |
Collapse
|
16
|
Min L, Wu Y, Cao G, Mi D, Chen C. A network pharmacology strategy to investigate the anti-osteoarthritis mechanism of main lignans components of Schisandrae Fructus. Int Immunopharmacol 2021; 98:107873. [PMID: 34182246 DOI: 10.1016/j.intimp.2021.107873] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2020] [Revised: 06/09/2021] [Accepted: 06/09/2021] [Indexed: 10/21/2022]
Abstract
Osteoarthritis (OA) is a chronic age-related progressive joint disorder. Degradation of the cartilage extracellular matrix (ECM) is considered a hallmark of OA and may be a target for new therapeutic methods. Schisandrae Fructus (SF) has been shown to be effective in treating OA. The major active components of SF are lignans. However, the targets of SF and the pharmacological mechanisms underlying the effects of SF lignans in the treatment of OA have not been elucidated. Therefore, based on network pharmacology, this research predicted the treatment targets of six lignans in SF, constructed a protein-protein interaction network and identified 15 hub genes in the OA-target protein-protein interaction network. Through Gene Ontology function and pathway analyses, the gene functions of lignans in the treatment of OA were determined. Finally, the anti-OA effects of lignans and underlying mechanisms identified in the network pharmacology analysis were verified by molecular docking, real-time PCR and western blotting in vitro. The biological processes of the genes and proteins targeted by lignans in the treatment of OA included the immune response, inflammatory response, cell signal transduction and phospholipid metabolism. Moreover, 20 metabolic pathways were enriched. Network pharmacology, molecular docking and in vitro and in vivo experimental results revealed that SF, schisanhenol and gamma-schisandrin inhibited EGFR and MAPK14 gene expression by inhibiting SRC gene expression and activity and then decreased MMP 13 and collagen II protein and gene expression. This research provides a basis for further study of the anti-OA effects and mechanisms of SF, schisanhenol and gamma-schisandrin.
Collapse
Affiliation(s)
- Lingtian Min
- Department of Orthopaedics, Nantong Hospital of Traditional Chinese Medicine, the Affiliated Hospital of Nanjing University of Traditional Chinese Medicine, Nantong 226000, China
| | - Yu Wu
- Department of Pharmacy, Nantong Hospital of Traditional Chinese Medicine, the Affiliated Hospital of Nanjing University of Traditional Chinese Medicine, Nantong 226000, China
| | - Gang Cao
- Department of Pharmacy, Nantong Hospital of Traditional Chinese Medicine, the Affiliated Hospital of Nanjing University of Traditional Chinese Medicine, Nantong 226000, China
| | - Daguo Mi
- Department of Orthopaedics, Nantong Hospital of Traditional Chinese Medicine, the Affiliated Hospital of Nanjing University of Traditional Chinese Medicine, Nantong 226000, China.
| | - Cheng Chen
- Department of Orthopaedics, Suqian First Hospital, Affiliated to Nanjing Medical University, Suqian 223800, China.
| |
Collapse
|
17
|
Mandal SK, Munshi P. Predicting Accurate Lead Structures for Screening Molecular Libraries: A Quantum Crystallographic Approach. Molecules 2021; 26:molecules26092605. [PMID: 33946965 PMCID: PMC8124947 DOI: 10.3390/molecules26092605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 04/22/2021] [Accepted: 04/24/2021] [Indexed: 11/17/2022] Open
Abstract
Optimization of lead structures is crucial for drug discovery. However, the accuracy of such a prediction using the traditional molecular docking approach remains a major concern. Our study demonstrates that the employment of quantum crystallographic approach-counterpoise corrected kernel energy method (KEM-CP) can improve the accuracy by and large. We select human aldose reductase at 0.66 Å, cyclin dependent kinase 2 at 2.0 Å and estrogen receptor β at 2.7 Å resolutions with active site environment ranging from highly hydrophilic to moderate to highly hydrophobic and several of their known ligands. Overall, the use of KEM-CP alongside the GoldScore resulted superior prediction than the GoldScore alone. Unlike GoldScore, the KEM-CP approach is neither environment-specific nor structural resolution dependent, which highlights its versatility. Further, the ranking of the ligands based on the KEM-CP results correlated well with that of the experimental IC50 values. This computationally inexpensive yet simple approach is expected to ease the process of virtual screening of potent ligands, and it would advance the drug discovery research.
Collapse
|
18
|
Giese TJ, York DM. Variational Method for Networkwide Analysis of Relative Ligand Binding Free Energies with Loop Closure and Experimental Constraints. J Chem Theory Comput 2021; 17:1326-1336. [PMID: 33528251 PMCID: PMC8011336 DOI: 10.1021/acs.jctc.0c01219] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
We describe an efficient method for the simultaneous solution of all free energies within a relative binding free-energy (RBFE) network with cycle closure and experimental/reference constraint conditions using Bennett Acceptance Ratio (BAR) and Multistate BAR (MBAR) analysis. Rather than solving the BAR or MBAR equations for each transformation independently, the simultaneous solution of all transformations are obtained by performing a constrained minimization of a global objective function. The nonlinear optimization of the objective function is subjected to affine linear constraints that couple the free energies between the network edges. The constraints are used to enforce the closure of thermodynamic cycles within the RBFE network, and to enforce an additional set of linear constraint conditions demonstrated here to be subsets of (1 or 2) experimental values. We describe details of the practical implementation of the network BAR/MBAR procedure, including use of generalized coordinates in the minimization of the free-energy objective function, propagation of bootstrap errors from those coordinates, and performance and memory optimization. In some cases it is found that use of restraints in the optimization is more practical than use of generalized coordinates for enforcing constraint conditions. The fast BARnet and MBARnet methods are used to analyze the RBFEs of six prototypical protein-ligand systems, and it is shown that enforcement of cycle closure conditions reduces the error in the predictions only modestly, and further reduction in errors can be achieved when one or two experimental RBFEs are included in the optimization procedure. These methods have been implemented into FE-ToolKit, a new free-energy analysis toolkit. The BARnet/MBARnet framework presented here opens the door to new, more efficient and robust free-energy analysis with enhanced predictive capability for drug discovery applications.
Collapse
Affiliation(s)
- Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854-8087 USA
| |
Collapse
|
19
|
Bieniek M, Bhati AP, Wan S, Coveney PV. TIES 20: Relative Binding Free Energy with a Flexible Superimposition Algorithm and Partial Ring Morphing. J Chem Theory Comput 2021; 17:1250-1265. [PMID: 33486956 PMCID: PMC7876800 DOI: 10.1021/acs.jctc.0c01179] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Indexed: 12/14/2022]
Abstract
The TIES (Thermodynamic Integration with Enhanced Sampling) protocol is a formally exact alchemical approach in computational chemistry to the calculation of relative binding free energies. The validity of TIES relies on the correctness of matching atoms across compared pairs of ligands, laying the foundation for the transformation along an alchemical pathway. We implement a flexible topology superimposition algorithm which uses an exhaustive joint-traversal for computing the largest common component(s). The algorithm is employed to enable matching and morphing of partial rings in the TIES protocol along with a validation study using 55 transformations and five different proteins from our previous work. We find that TIES 20 with the RESP charge system, using the new superimposition algorithm, reproduces the previous results with mean unsigned error of 0.75 kcal/mol with respect to the experimental data. Enabling the morphing of partial rings decreases the size of the alchemical region in the dual-topology transformations resulting in a significant improvement in the prediction precision. We find that increasing the ensemble size from 5 to 20 replicas per λ window only has a minimal impact on the accuracy. However, the non-normal nature of the relative free energy distributions underscores the importance of ensemble simulation. We further compare the results with the AM1-BCC charge system and show that it improves agreement with the experimental data by slightly over 10%. This improvement is partly due to AM1-BCC affecting only the charges of the atoms local to the mutation, which translates to even fewer morphed atoms, consequently reducing issues with sampling and therefore ensemble averaging. TIES 20, in conjunction with the enablement of ring morphing, reduces the size of the alchemical region and significantly improves the precision of the predicted free energies.
Collapse
Affiliation(s)
- Mateusz
K. Bieniek
- Centre for Computational Science, Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Agastya P. Bhati
- Centre for Computational Science, Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shunzhou Wan
- Centre for Computational Science, Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Peter V. Coveney
- Centre for Computational Science, Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
20
|
Hao D, He X, Ji B, Zhang S, Wang J. How Well Does the Extended Linear Interaction Energy Method Perform in Accurate Binding Free Energy Calculations? J Chem Inf Model 2020; 60:6624-6633. [PMID: 33213150 DOI: 10.1021/acs.jcim.0c00934] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
With continually increased computer power, molecular mechanics force field-based approaches, such as the endpoint methods of molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA), have been routinely applied in both drug lead identification and optimization. However, the MM-PB/GBSA method is not as accurate as the pathway-based alchemical free energy methods, such as thermodynamic integration (TI) or free energy perturbation (FEP). Although the pathway-based methods are more rigorous in theory, they suffer from slow convergence and computational cost. Moreover, choosing adequate perturbation routes is also crucial for the pathway-based methods. Recently, we proposed a new method, coined extended linear interaction energy (ELIE) method, to overcome some disadvantages of the MM-PB/GBSA method to improve the accuracy of binding free energy calculation. In this work, we have systematically assessed this approach using in total 229 protein-ligand complexes for eight protein targets. Our results showed that ELIE performed much better than the molecular docking and MM-PBSA method in terms of root-mean-square error (RMSE), correlation coefficient (R), predictive index (PI), and Kendall's τ. The mean values of PI, R, and τ are 0.62, 0.58, and 0.44 for ELIE calculations. We also explored the impact of the length of simulation, ranging from 1 to 100 ns, on the performance of binding free energy calculation. In general, extending simulation length up to 25 ns could significantly improve the performance of ELIE, while longer molecular dynamics (MD) simulation does not always perform better than short MD simulation. Considering both the computational efficiency and achieved accuracy, ELIE is adequate in filling the gap between the efficient docking methods and computationally demanding alchemical free energy methods. Therefore, ELIE provides a practical solution for the routine ranking of compounds in lead optimization.
Collapse
Affiliation(s)
- Dongxiao Hao
- School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.,Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Shengli Zhang
- School of Physics, Xi'an Jiaotong University, Xi'an 710049, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
21
|
Guglielmo S, Cortese D, Cano C, Fruttero R. Molecular dynamics simulations reveal the determinants of cyclin-dependent kinase 2 inhibition by 5-nitrosopyrimidine derivatives. J Biomol Struct Dyn 2020; 38:4016-4024. [PMID: 31498033 DOI: 10.1080/07391102.2019.1666032] [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/15/2019] [Revised: 09/02/2019] [Accepted: 09/03/2019] [Indexed: 10/26/2022]
Abstract
Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Stefano Guglielmo
- Dipartimento di Scienza e Tecnologia del Farmaco, Universita' degli Studi di Torino, Turin, Italy
- Scientific Computing Competence Centre (C3S), University of Turin, Turin, Italy
| | - Daniela Cortese
- Dipartimento di Scienza e Tecnologia del Farmaco, Universita' degli Studi di Torino, Turin, Italy
- Northern Institute for Cancer Research, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, UK
| | - Celine Cano
- Northern Institute for Cancer Research, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, UK
| | - Roberta Fruttero
- Dipartimento di Scienza e Tecnologia del Farmaco, Universita' degli Studi di Torino, Turin, Italy
| |
Collapse
|
22
|
Kuhn M, Firth-Clark S, Tosco P, Mey ASJS, Mackey M, Michel J. Assessment of Binding Affinity via Alchemical Free-Energy Calculations. J Chem Inf Model 2020; 60:3120-3130. [DOI: 10.1021/acs.jcim.0c00165] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Maximilian Kuhn
- Cresset, New Cambridge House, Bassingbourn Road, Litlington SG8 0SS, Cambridgeshire, U.K
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, U.K
| | - Stuart Firth-Clark
- Cresset, New Cambridge House, Bassingbourn Road, Litlington SG8 0SS, Cambridgeshire, U.K
| | - Paolo Tosco
- Cresset, New Cambridge House, Bassingbourn Road, Litlington SG8 0SS, Cambridgeshire, U.K
| | - Antonia S. J. S. Mey
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, U.K
| | - Mark Mackey
- Cresset, New Cambridge House, Bassingbourn Road, Litlington SG8 0SS, Cambridgeshire, U.K
| | - Julien Michel
- EaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, U.K
| |
Collapse
|
23
|
He X, Liu S, Lee TS, Ji B, Man VH, York DM, Wang J. Fast, Accurate, and Reliable Protocols for Routine Calculations of Protein-Ligand Binding Affinities in Drug Design Projects Using AMBER GPU-TI with ff14SB/GAFF. ACS OMEGA 2020; 5:4611-4619. [PMID: 32175507 PMCID: PMC7066661 DOI: 10.1021/acsomega.9b04233] [Citation(s) in RCA: 70] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 02/13/2020] [Indexed: 05/12/2023]
Abstract
Accurate prediction of the absolute or relative protein-ligand binding affinity is one of the major tasks in computer-aided drug design projects, especially in the stage of lead optimization. In principle, the alchemical free energy (AFE) methods such as thermodynamic integration (TI) or free-energy perturbation (FEP) can fulfill this task, but in practice, a lot of hurdles prevent them from being routinely applied in daily drug design projects, such as the demanding computing resources, slow computing processes, unavailable or inaccurate force field parameters, and difficult and unfriendly setting up and post-analysis procedures. In this study, we have exploited practical protocols of applying the CPU (central processing unit)-TI and newly developed GPU (graphic processing unit)-TI modules and other tools in the AMBER software package, combined with ff14SB/GAFF1.8 force fields, to conduct efficient and accurate AFE calculations on protein-ligand binding free energies. We have tested 134 protein-ligand complexes in total for four target proteins (BACE, CDK2, MCL1, and PTP1B) and obtained overall comparable performance with the commercial Schrodinger FEP+ program (WangJ. Am. Chem. Soc.2015, 137, 2695-2703). The achieved accuracy fits within the requirements for computations to generate effective guidance for experimental work in drug lead optimization, and the needed wall time is short enough for practical application. Our verified protocol provides a practical solution for routine AFE calculations in real drug design projects.
Collapse
Affiliation(s)
- Xibing He
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Shuhan Liu
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Tai-Sung Lee
- Laboratory
for Biomolecular Simulation Research, Center for Integrative Proteomics
Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Beihong Ji
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Viet H. Man
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Darrin M. York
- Laboratory
for Biomolecular Simulation Research, Center for Integrative Proteomics
Research, and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Junmei Wang
- Department
of Pharmaceutical Sciences and Computational Chemical Genomics Screening
Center, School of Pharmacy, University of
Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
- E-mail: . Phone: (412) 383-3268. Fax: (412) 383-7436
| |
Collapse
|
24
|
Jespers W, Esguerra M, Åqvist J, Gutiérrez-de-Terán H. QligFEP: an automated workflow for small molecule free energy calculations in Q. J Cheminform 2019; 11:26. [PMID: 30941533 PMCID: PMC6444553 DOI: 10.1186/s13321-019-0348-5] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 03/23/2019] [Indexed: 11/24/2022] Open
Abstract
The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.
Collapse
Affiliation(s)
- Willem Jespers
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Mauricio Esguerra
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Johan Åqvist
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
| | - Hugo Gutiérrez-de-Terán
- Department of Cell and Molecular Biology, Uppsala University, Uppsala, 75124 Sweden
- Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| |
Collapse
|
25
|
Roos K, Wu C, Damm W, Reboul M, Stevenson JM, Lu C, Dahlgren MK, Mondal S, Chen W, Wang L, Abel R, Friesner RA, Harder ED. OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules. J Chem Theory Comput 2019; 15:1863-1874. [PMID: 30768902 DOI: 10.1021/acs.jctc.8b01026] [Citation(s) in RCA: 633] [Impact Index Per Article: 126.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Katarina Roos
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
- Department of Cell and Molecular Biology, Uppsala University, Biomedical Centre, Box 596, SE-751 24 Uppsala, Sweden
| | - Chuanjie Wu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wolfgang Damm
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Mark Reboul
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - James M. Stevenson
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Chao Lu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Markus K. Dahlgren
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Sayan Mondal
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wei Chen
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Robert Abel
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Edward D. Harder
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| |
Collapse
|
26
|
Tadesse S, Caldon EC, Tilley W, Wang S. Cyclin-Dependent Kinase 2 Inhibitors in Cancer Therapy: An Update. J Med Chem 2018; 62:4233-4251. [PMID: 30543440 DOI: 10.1021/acs.jmedchem.8b01469] [Citation(s) in RCA: 139] [Impact Index Per Article: 23.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Cyclin-dependent kinase 2 (CDK2) drives the progression of cells into the S- and M-phases of the cell cycle. CDK2 activity is largely dispensable for normal development, but it is critically associated with tumor growth in multiple cancer types. Although the role of CDK2 in tumorigenesis has been controversial, emerging evidence proposes that selective CDK2 inhibition may provide a therapeutic benefit against certain tumors, and it continues to appeal as a strategy to exploit in anticancer drug development. Several small-molecule CDK2 inhibitors have progressed to the clinical trials. However, a CDK2-selective inhibitor is yet to be discovered. Here, we discuss the latest understandings of the role of CDK2 in normal and cancer cells, review the core pharmacophores used to target CDK2, and outline strategies for the rational design of CDK2 inhibitors. We attempt to provide an outlook on how CDK2-selective inhibitors may open new avenues for cancer therapy.
Collapse
Affiliation(s)
- Solomon Tadesse
- Centre for Drug Discovery and Development , University of South Australia Cancer Research Institute , Adelaide , SA 5000 , Australia
| | - Elizabeth C Caldon
- The Kinghorn Cancer Centre , Garvan Institute of Medical Research , Darlinghurst , NSW 2010 , Australia.,St Vincent's Clinical School, UNSW Medicine , UNSW Sydney , Darlinghurst , NSW 2010 , Australia
| | - Wayne Tilley
- Adelaide Medical School , University of Adelaide , Adelaide , SA 5000 , Australia
| | - Shudong Wang
- Centre for Drug Discovery and Development , University of South Australia Cancer Research Institute , Adelaide , SA 5000 , Australia
| |
Collapse
|
27
|
Chen W, Deng Y, Russell E, Wu Y, Abel R, Wang L. Accurate Calculation of Relative Binding Free Energies between Ligands with Different Net Charges. J Chem Theory Comput 2018; 14:6346-6358. [PMID: 30375870 DOI: 10.1021/acs.jctc.8b00825] [Citation(s) in RCA: 67] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Wei Chen
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Yuqing Deng
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Ellery Russell
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Yujie Wu
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Robert Abel
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| |
Collapse
|
28
|
Jorda R, Hendrychová D, Voller J, Řezníčková E, Gucký T, Kryštof V. How Selective Are Pharmacological Inhibitors of Cell-Cycle-Regulating Cyclin-Dependent Kinases? J Med Chem 2018; 61:9105-9120. [PMID: 30234987 DOI: 10.1021/acs.jmedchem.8b00049] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Cyclin-dependent kinases (CDKs) are an important and emerging class of drug targets for which many small-molecule inhibitors have been developed. However, there is often insufficient data available on the selectivity of CDK inhibitors (CDKi) to attribute the effects on the presumed target CDK to these inhibitors. Here, we highlight discrepancies between the kinase selectivity of CDKi and the phenotype exhibited; we evaluated 31 CDKi (claimed to target CDK1-4) for activity toward CDKs 1, 2, 4, 5, 7, 9 and for effects on the cell cycle. Our results suggest that most CDKi should be reclassified as pan-selective and should not be used as a tool. In addition, some compounds did not even inhibit CDKs as their primary cellular targets; for example, NU6140 showed potent inhibition of Aurora kinases. We also established an online database of commercially available CDKi for critical evaluation of their utility as molecular probes. Our results should help researchers select the most relevant chemical tools for their specific applications.
Collapse
Affiliation(s)
- Radek Jorda
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| | - Denisa Hendrychová
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| | - Jiří Voller
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| | - Eva Řezníčková
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| | - Tomáš Gucký
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| | - Vladimír Kryštof
- Laboratory of Growth Regulators, Centre of the Region Haná for Biotechnological and Agricultural Research , Palacký University and Institute of Experimental Botany ASCR , Šlechtitelů 27 , 78371 Olomouc , Czech Republic
| |
Collapse
|
29
|
Wu YZ, Ying HZ, Xu L, Cheng G, Chen J, Hu YZ, Liu T, Dong XW. Design, synthesis, and molecular docking study of 3H
-imidazole[4,5-c
]pyridine derivatives as CDK2 inhibitors. Arch Pharm (Weinheim) 2018; 351:e1700381. [DOI: 10.1002/ardp.201700381] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Revised: 04/04/2018] [Accepted: 04/05/2018] [Indexed: 11/08/2022]
Affiliation(s)
- Yi-Zhe Wu
- ZJU-ENS Joint Laboratory of Medicinal Chemistry; Zhejiang Province Key Laboratory of Anti-Cancer Drug Research; Hangzhou Institute of Innovative Medicine; College of Pharmaceutical Sciences; Zhejiang University; Hangzhou P. R. China
| | - Hua-Zhou Ying
- ZJU-ENS Joint Laboratory of Medicinal Chemistry; Zhejiang Province Key Laboratory of Anti-Cancer Drug Research; Hangzhou Institute of Innovative Medicine; College of Pharmaceutical Sciences; Zhejiang University; Hangzhou P. R. China
| | - Lei Xu
- School of Life Science and Technology; ShanghaiTech University; Shanghai P. R. China
- State Key Laboratory of Drug Research; Shanghai Institute of Materia Medica; Chinese Academy of Sciences; Shanghai P. R. China
| | - Gang Cheng
- College of Pharmaceutical Science; Zhejiang Chinese Medical University; Hangzhou P. R. China
| | - Jing Chen
- College of Pharmaceutical Science; Zhejiang Chinese Medical University; Hangzhou P. R. China
| | - Yong-Zhou Hu
- ZJU-ENS Joint Laboratory of Medicinal Chemistry; Zhejiang Province Key Laboratory of Anti-Cancer Drug Research; Hangzhou Institute of Innovative Medicine; College of Pharmaceutical Sciences; Zhejiang University; Hangzhou P. R. China
| | - Tao Liu
- ZJU-ENS Joint Laboratory of Medicinal Chemistry; Zhejiang Province Key Laboratory of Anti-Cancer Drug Research; Hangzhou Institute of Innovative Medicine; College of Pharmaceutical Sciences; Zhejiang University; Hangzhou P. R. China
| | - Xiao-Wu Dong
- ZJU-ENS Joint Laboratory of Medicinal Chemistry; Zhejiang Province Key Laboratory of Anti-Cancer Drug Research; Hangzhou Institute of Innovative Medicine; College of Pharmaceutical Sciences; Zhejiang University; Hangzhou P. R. China
| |
Collapse
|
30
|
Kamble NR, Sigurdsson ST. Purine-Derived Nitroxides for Noncovalent Spin-Labeling of Abasic Sites in Duplex Nucleic Acids. Chemistry 2018; 24:4157-4164. [DOI: 10.1002/chem.201705410] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2017] [Indexed: 12/27/2022]
Affiliation(s)
- Nilesh R. Kamble
- University of Iceland; Department of Chemistry; Science Institute; Dunhaga 3 107 Reykjavik Iceland
| | - Snorri Th. Sigurdsson
- University of Iceland; Department of Chemistry; Science Institute; Dunhaga 3 107 Reykjavik Iceland
| |
Collapse
|
31
|
Bogenberger J, Whatcott C, Hansen N, Delman D, Shi CX, Kim W, Haws H, Soh K, Lee YS, Peterson P, Siddiqui-Jain A, Weitman S, Stewart K, Bearss D, Mesa R, Warner S, Tibes R. Combined venetoclax and alvocidib in acute myeloid leukemia. Oncotarget 2017; 8:107206-107222. [PMID: 29291023 PMCID: PMC5739808 DOI: 10.18632/oncotarget.22284] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2017] [Accepted: 10/10/2017] [Indexed: 11/25/2022] Open
Abstract
More effective treatment options for elderly acute myeloid leukemia (AML) patients are needed as only 25-50% of patients respond to standard-of-care therapies, response duration is typically short, and disease progression is inevitable even with some novel therapies and ongoing clinical trials. Anti-apoptotic BCL-2 family inhibitors, such as venetoclax, are promising therapies for AML. Nonetheless, resistance is emerging. We demonstrate that venetoclax combined with cyclin-dependent kinase (CDK) inhibitor alvocidib is potently synergistic in venetoclax-sensitive and -resistant AML models in vitro, ex vivo and in vivo. Alvocidib decreased MCL-1, and/or increased pro-apoptotic proteins such as BIM or NOXA, often synergistically with venetoclax. Over-expression of BCL-XL diminished synergy, while knock-down of BIM almost entirely abrogated synergy, demonstrating that the synergistic interaction between alvocidib and venetoclax is primarily dependent on intrinsic apoptosis. CDK9 inhibition predominantly mediated venetoclax sensitization, while CDK4/6 inhibition with palbociclib did not potentiate venetoclax activity. Combined, venetoclax and alvocidib modulate the balance of BCL-2 family proteins through complementary, yet variable mechanisms favoring apoptosis, highlighting this combination as a promising therapy for AML or high-risk MDS with the capacity to overcome intrinsic apoptosis mechanisms of resistance. These results support clinical testing of combined venetoclax and alvocidib for the treatment of AML and advanced MDS.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - Raoul Tibes
- Mayo Clinic, Scottsdale, AZ, USA.,NYU School of Medicine, New York, NY, USA
| |
Collapse
|
32
|
Mahajan P, Chashoo G, Gupta M, Kumar A, Singh PP, Nargotra A. Fusion of Structure and Ligand Based Methods for Identification of Novel CDK2 Inhibitors. J Chem Inf Model 2017; 57:1957-1969. [PMID: 28723151 DOI: 10.1021/acs.jcim.7b00293] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Cyclin dependent kinases play a central role in cell cycle regulation which makes them a promising target with multifarious therapeutic potential. CDK2 regulates various events of the eukaryotic cell division cycle, and the pharmacological evidence indicates that overexpression of CDK2 causes abnormal cell-cycle regulation, which is directly associated with hyperproliferation of cancer cells. Therefore, CDK2 is regarded as a potential target molecule for anticancer medication. Thus, to decline CDK2 activity by potential lead compounds has proved to be an effective treatment for cancer. The availability of a large number of X-ray crystal structures and known inhibitors of CDK2 provides a gateway to perform efficient computational studies on this target. With the aim to identify new chemical entities from commercial libraries, with increased inhibitory potency for CDK2, ligand and structure based computational drug designing approaches were applied. A druglike library of 50,000 compounds from ChemDiv and ChemBridge databases was screened against CDK2, and 110 compounds were identified using the parallel application of these models. On in vitro evaluation of 40 compounds, seven compounds were found to have more than 50% inhibition at 10 μM. MD studies of the hits revealed the stability of these inhibitors and pivotal role of Glu81 and Leu83 for binding with CDK2. The overall study resulted in the identification of four new chemical entities possessing CDK2 inhibitory activity.
Collapse
Affiliation(s)
- Priya Mahajan
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| | - Gousia Chashoo
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| | - Monika Gupta
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| | - Amit Kumar
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| | - Parvinder Pal Singh
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| | - Amit Nargotra
- Discovery Informatics, ‡Cancer Pharmacology, §Medicinal Chemistry, and ∥Academy of Scientific and Innovative Research, CSIR-Indian Institute of Integrative Medicine , Canal Road, Jammu 180001, India
| |
Collapse
|
33
|
Coxon CR, Wong C, Bayliss R, Boxall K, Carr KH, Fry AM, Hardcastle IR, Matheson CJ, Newell DR, Sivaprakasam M, Thomas H, Turner D, Yeoh S, Wang LZ, Griffin RJ, Golding BT, Cano C. Structure-guided design of purine-based probes for selective Nek2 inhibition. Oncotarget 2017; 8:19089-19124. [PMID: 27833088 PMCID: PMC5386672 DOI: 10.18632/oncotarget.13249] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 10/17/2016] [Indexed: 01/23/2023] Open
Abstract
Nek2 (NIMA-related kinase 2) is a cell cycle-dependent serine/threonine protein kinase that regulates centrosome separation at the onset of mitosis. Overexpression of Nek2 is common in human cancers and suppression can restrict tumor cell growth and promote apoptosis. Nek2 inhibition with small molecules, therefore, offers the prospect of a new therapy for cancer. To achieve this goal, a better understanding of the requirements for selective-inhibition of Nek2 is required. 6-Alkoxypurines were identified as ATP-competitive inhibitors of Nek2 and CDK2. Comparison with CDK2-inhibitor structures indicated that judicious modification of the 6-alkoxy and 2-arylamino substituents could achieve discrimination between Nek2 and CDK2. In this study, a library of 6-cyclohexylmethoxy-2-arylaminopurines bearing carboxamide, sulfonamide and urea substituents on the 2-arylamino ring was synthesized. Few of these compounds were selective for Nek2 over CDK2, with the best result being obtained for 3-((6-(cyclohexylmethoxy)-9H-purin-2-yl)amino)-N,N-dimethylbenzamide (CDK2 IC50 = 7.0 μM; Nek2 IC50 = 0.62 μM) with >10-fold selectivity. Deletion of the 6-substituent abrogated activity against both Nek2 and CDK2. Nine compounds containing an (E)-dialkylaminovinyl substituent at C-6, all showed selectivity for Nek2, e.g. (E)-6-(2-(azepan-1-yl)vinyl)-N-phenyl-9H-purin-2-amine (CDK2 IC50 = 2.70 μM; Nek2 IC50 = 0.27 μM). Structural biology of selected compounds enabled a partial rationalization of the observed structure activity relationships and mechanism of Nek2 activation. This showed that carboxamide 11 is the first reported inhibitor of Nek2 in the DFG-in conformation.
Collapse
Affiliation(s)
- Christopher R. Coxon
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Christopher Wong
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Richard Bayliss
- Department of Molecular and Cell Biology, University of Leicester, Leicester, UK
| | - Kathy Boxall
- Cancer Research UK Cancer Therapeutics Unit, The Institute of Cancer Research, London, UK
| | - Katherine H. Carr
- Department of Molecular and Cell Biology, University of Leicester, Leicester, UK
| | - Andrew M. Fry
- Department of Molecular and Cell Biology, University of Leicester, Leicester, UK
| | - Ian R. Hardcastle
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Christopher J. Matheson
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - David R. Newell
- Northern Institute for Cancer Research, Newcastle University, Newcastle upon Tyne, UK
| | - Mangaleswaran Sivaprakasam
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Huw Thomas
- Northern Institute for Cancer Research, Newcastle University, Newcastle upon Tyne, UK
| | - David Turner
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Sharon Yeoh
- Department of Molecular and Cell Biology, University of Leicester, Leicester, UK
| | - Lan Z. Wang
- Northern Institute for Cancer Research, Newcastle University, Newcastle upon Tyne, UK
| | - Roger J. Griffin
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Bernard T. Golding
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| | - Céline Cano
- Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Newcastle upon Tyne, UK
| |
Collapse
|
34
|
Coxon C, Anscombe E, Harnor SJ, Martin MP, Carbain B, Golding BT, Hardcastle IR, Harlow LK, Korolchuk S, Matheson CJ, Newell DR, Noble MEM, Sivaprakasam M, Tudhope SJ, Turner DM, Wang LZ, Wedge SR, Wong C, Griffin RJ, Endicott JA, Cano C. Cyclin-Dependent Kinase (CDK) Inhibitors: Structure-Activity Relationships and Insights into the CDK-2 Selectivity of 6-Substituted 2-Arylaminopurines. J Med Chem 2017; 60:1746-1767. [PMID: 28005359 PMCID: PMC6111440 DOI: 10.1021/acs.jmedchem.6b01254] [Citation(s) in RCA: 69] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2016] [Indexed: 02/08/2023]
Abstract
Purines and related heterocycles substituted at C-2 with 4'-sulfamoylanilino and at C-6 with a variety of groups have been synthesized with the aim of achieving selectivity of binding to CDK2 over CDK1. 6-Substituents that favor competitive inhibition at the ATP binding site of CDK2 were identified and typically exhibited 10-80-fold greater inhibition of CDK2 compared to CDK1. Most impressive was 4-((6-([1,1'-biphenyl]-3-yl)-9H-purin-2-yl)amino) benzenesulfonamide (73) that exhibited high potency toward CDK2 (IC50 0.044 μM) but was ∼2000-fold less active toward CDK1 (IC50 86 μM). This compound is therefore a useful tool for studies of cell cycle regulation. Crystal structures of inhibitor-kinase complexes showed that the inhibitor stabilizes a glycine-rich loop conformation that shapes the ATP ribose binding pocket and that is preferred in CDK2 but has not been observed in CDK1. This aspect of the active site may be exploited for the design of inhibitors that distinguish between CDK1 and CDK2.
Collapse
Affiliation(s)
- Christopher
R. Coxon
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Elizabeth Anscombe
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Suzannah J. Harnor
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Mathew P. Martin
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - Benoit Carbain
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Bernard T. Golding
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Ian R. Hardcastle
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Lisa K. Harlow
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Svitlana Korolchuk
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - Christopher J. Matheson
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - David R. Newell
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - Martin E. M. Noble
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mangaleswaran Sivaprakasam
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Susan J. Tudhope
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - David M. Turner
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Lan Z. Wang
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - Stephen R. Wedge
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, Newcastle University Medical School, Paul O’Gorman Building, Framlington Place, Newcastle upon Tyne, NE2 4HH, U.K.
| | - Christopher Wong
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Roger J. Griffin
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| | - Jane A. Endicott
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Céline Cano
- Newcastle
Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Newcastle University, Bedson Building, Newcastle
upon Tyne NE1 7RU, U.K.
| |
Collapse
|
35
|
Bhati AP, Wan S, Wright DW, Coveney PV. Rapid, Accurate, Precise, and Reliable Relative Free Energy Prediction Using Ensemble Based Thermodynamic Integration. J Chem Theory Comput 2016; 13:210-222. [PMID: 27997169 DOI: 10.1021/acs.jctc.6b00979] [Citation(s) in RCA: 92] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The accurate prediction of the binding affinities of ligands to proteins is a major goal in drug discovery and personalized medicine. The time taken to make such predictions is of similar importance to their accuracy, precision, and reliability. In the past few years, an ensemble based molecular dynamics approach has been proposed that provides a route to reliable predictions of free energies based on the molecular mechanics Poisson-Boltzmann surface area method which meets the requirements of speed, accuracy, precision, and reliability. Here, we describe an equivalent methodology based on thermodynamic integration to substantially improve the speed, accuracy, precision, and reliability of calculated relative binding free energies. We report the performance of the method when applied to a diverse set of protein targets and ligands. The results are in very good agreement with experimental data (90% of calculations agree to within 1 kcal/mol), while the method is reproducible by construction. Statistical uncertainties of the order of 0.5 kcal/mol or less are achieved. We present a systematic account of how the uncertainty in the predictions may be estimated.
Collapse
Affiliation(s)
- Agastya P Bhati
- Centre for Computational Science, Department of Chemistry, University College London , 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Shunzhou Wan
- Centre for Computational Science, Department of Chemistry, University College London , 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - David W Wright
- Centre for Computational Science, Department of Chemistry, University College London , 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Peter V Coveney
- Centre for Computational Science, Department of Chemistry, University College London , 20 Gordon Street, London WC1H 0AJ, United Kingdom
| |
Collapse
|
36
|
Cortese D, Chegaev K, Guglielmo S, Wang LZ, Golding BT, Cano C, Fruttero R. Synthesis and Biological Evaluation of N(2) -Substituted 2,4-Diamino-6-cyclohexylmethoxy-5-nitrosopyrimidines and Related 5-Cyano-NNO-azoxy Derivatives as Cyclin-Dependent Kinase 2 (CDK2) Inhibitors. ChemMedChem 2016; 11:1705-8. [PMID: 27355194 DOI: 10.1002/cmdc.201600108] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Revised: 06/07/2016] [Indexed: 11/09/2022]
Abstract
The potent and selective cyclin-dependent kinase 2 (CDK2) inhibitor NU6027 (6-cyclohexylmethoxy-5-nitroso-2,4-diaminopyrimidine) was used as the lead for the synthesis of a series of analogues in order to provide further insight into the structure-activity relationships for 2,4-diaminopyrimidine CDK2 inhibitors. Aliphatic amino substituents were introduced at position 2. The use of linear or less sterically hindered amines gave rise to compounds endowed with slightly better activity than the lead; on the other hand, the compounds were less active if a bulkier amino substituent was used. Substitution of the 5-nitroso group with a 5-cyano-NNO-azoxy moiety afforded a new class of inhibitors, the activity of which against CDK2 was found to be similar to that of the nitroso series. The most active nitroso compound was 8 b ((2S)-2-[(4-amino-6-cyclohexylmethoxy-5-nitrosopyrimidin-2-yl)amino]propan-1-ol; IC50 =0.16 μm), while in the 5-cyano-NNO-azoxy series the most active compound was 9 b (4-amino-5-[(Z)-cyano-NNO-azoxy]-2-{[(2S)-1-hydroxypropan-2-yl]amino}-6-cyclohexylmethoxypyrimidine; IC50 =0.30 μm). Taken together, these new analogues of NU6027 enhance our understanding of the structure-activity relationships for 2,4-diaminopyrimidine CDK2 inhibitors.
Collapse
Affiliation(s)
- Daniela Cortese
- Department of Drug Science and Technology, Università degli Studi di Torino, Via P. Giuria 9, 10125, Turin, Italy
| | - Konstantin Chegaev
- Department of Drug Science and Technology, Università degli Studi di Torino, Via P. Giuria 9, 10125, Turin, Italy
| | - Stefano Guglielmo
- Department of Drug Science and Technology, Università degli Studi di Torino, Via P. Giuria 9, 10125, Turin, Italy
| | - Lan Z Wang
- Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle Upon Tyne, NE2 4HH, UK
| | - Bernard T Golding
- Northern Institute for Cancer Research, Bedson Building, School of Chemistry, Newcastle University, Newcastle Upon Tyne, NE1 7RU, UK
| | - Céline Cano
- Northern Institute for Cancer Research, Bedson Building, School of Chemistry, Newcastle University, Newcastle Upon Tyne, NE1 7RU, UK
| | - Roberta Fruttero
- Department of Drug Science and Technology, Università degli Studi di Torino, Via P. Giuria 9, 10125, Turin, Italy.
| |
Collapse
|
37
|
Wang DC, Xia R, Xie MS, Qu GR, Guo HM. Synthesis of cycloalkyl substituted purine nucleosides via a metal-free radical route. Org Biomol Chem 2016; 14:4189-93. [PMID: 27101306 DOI: 10.1039/c6ob00596a] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An efficient route to synthesize cycloalkyl substituted purine nucleosides was developed. This metal-free C-H activation was accomplished by a tBuOOtBu initiated radical reaction. By adjusting the amount of tBuOOtBu and reaction time, the selective synthesis of C6-monocycloalkyl or C6,C8-dicycloalkyl substituted purine nucleosides could be realized. Furthermore, uracil and related nucleosides were also suitable substrates, giving the C5-cyclohexyl substituted uracil derivatives in good yields with excellent regioselectivities.
Collapse
Affiliation(s)
- Dong-Chao Wang
- School of Environment, Henan Normal University, Xinxiang, Henan province 453007, P. R. China and School of Chemistry and Chemical Engineering, Key Laboratory of Green Chemical Media and Reactions of Ministry of Education, Henan Normal University, 46 Jianshe Road, Xinxiang, 453007, P. R. China.
| | - Ran Xia
- School of Environment, Henan Normal University, Xinxiang, Henan province 453007, P. R. China and School of Chemistry and Chemical Engineering, Xinxiang University, Xinxiang 453003, China
| | - Ming-Sheng Xie
- School of Chemistry and Chemical Engineering, Key Laboratory of Green Chemical Media and Reactions of Ministry of Education, Henan Normal University, 46 Jianshe Road, Xinxiang, 453007, P. R. China.
| | - Gui-Rong Qu
- School of Chemistry and Chemical Engineering, Key Laboratory of Green Chemical Media and Reactions of Ministry of Education, Henan Normal University, 46 Jianshe Road, Xinxiang, 453007, P. R. China.
| | - Hai-Ming Guo
- School of Environment, Henan Normal University, Xinxiang, Henan province 453007, P. R. China and School of Chemistry and Chemical Engineering, Key Laboratory of Green Chemical Media and Reactions of Ministry of Education, Henan Normal University, 46 Jianshe Road, Xinxiang, 453007, P. R. China.
| |
Collapse
|
38
|
Knight JL, Krilov G, Borrelli KW, Williams J, Gunn JR, Clowes A, Cheng L, Friesner RA, Abel R. Leveraging Data Fusion Strategies in Multireceptor Lead Optimization MM/GBSA End-Point Methods. J Chem Theory Comput 2015; 10:3207-20. [PMID: 26588291 DOI: 10.1021/ct500189s] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Accurate and efficient affinity calculations are critical to enhancing the contribution of in silico modeling during the lead optimization phase of a drug discovery campaign. Here, we present a large-scale study of the efficacy of data fusion strategies to leverage results from end-point MM/GBSA calculations in multiple receptors to identify potent inhibitors among an ensemble of congeneric ligands. The retrospective analysis of 13 congeneric ligand series curated from publicly available data across seven biological targets demonstrates that in 90% of the individual receptor structures MM/GBSA scores successfully identify subsets of inhibitors that are more potent than a random selection, and data fusion strategies that combine MM/GBSA scores from each of the receptors significantly increase the robustness of the predictions. Among nine different data fusion metrics based on consensus scores or receptor rankings, the SumZScore (i.e., converting MM/GBSA scores into standardized Z-Scores within a receptor and computing the sum of the Z-Scores for a given ligand across the ensemble of receptors) is found to be a robust and physically meaningful metric for combining results across multiple receptors. Perhaps most surprisingly, even with relatively low to modest overall correlations between SumZScore and experimental binding affinities, SumZScore tends to reliably prioritize subsets of inhibitors that are at least as potent as those that are prioritized from a "best" single receptor identified from known compounds within the congeneric series.
Collapse
Affiliation(s)
- Jennifer L Knight
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Goran Krilov
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Kenneth W Borrelli
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Joshua Williams
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - John R Gunn
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Alec Clowes
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Luciano Cheng
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Richard A Friesner
- Columbia University , Department of Chemistry, 3000 Broadway, MC 3110, New York, New York 10027, United States
| | - Robert Abel
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| |
Collapse
|
39
|
Harder E, Damm W, Maple J, Wu C, Reboul M, Xiang JY, Wang L, Lupyan D, Dahlgren MK, Knight JL, Kaus JW, Cerutti DS, Krilov G, Jorgensen WL, Abel R, Friesner RA. OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. J Chem Theory Comput 2015; 12:281-96. [PMID: 26584231 DOI: 10.1021/acs.jctc.5b00864] [Citation(s) in RCA: 2047] [Impact Index Per Article: 227.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.
Collapse
Affiliation(s)
- Edward Harder
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wolfgang Damm
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jon Maple
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Chuanjie Wu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Mark Reboul
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jin Yu Xiang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Dmitry Lupyan
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Markus K Dahlgren
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jennifer L Knight
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Joseph W Kaus
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - David S Cerutti
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Goran Krilov
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - William L Jorgensen
- Department of Chemistry, Yale University , New Haven, Connecticut 06520, United States
| | - Robert Abel
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
40
|
Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis. Future Med Chem 2015; 7:2317-31. [DOI: 10.4155/fmc.15.150] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]
Collapse
|
41
|
Anscombe E, Meschini E, Mora-Vidal R, Martin MP, Staunton D, Geitmann M, Danielson UH, Stanley WA, Wang LZ, Reuillon T, Golding BT, Cano C, Newell DR, Noble MEM, Wedge SR, Endicott JA, Griffin RJ. Identification and Characterization of an Irreversible Inhibitor of CDK2. CHEMISTRY & BIOLOGY 2015; 22:1159-64. [PMID: 26320860 PMCID: PMC4579270 DOI: 10.1016/j.chembiol.2015.07.018] [Citation(s) in RCA: 77] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 07/02/2015] [Accepted: 07/24/2015] [Indexed: 01/04/2023]
Abstract
Irreversible inhibitors that modify cysteine or lysine residues within a protein kinase ATP binding site offer, through their distinctive mode of action, an alternative to ATP-competitive agents. 4-((6-(Cyclohexylmethoxy)-9H-purin-2-yl)amino)benzenesulfonamide (NU6102) is a potent and selective ATP-competitive inhibitor of CDK2 in which the sulfonamide moiety is positioned close to a pair of lysine residues. Guided by the CDK2/NU6102 structure, we designed 6-(cyclohexylmethoxy)-N-(4-(vinylsulfonyl)phenyl)-9H-purin-2-amine (NU6300), which binds covalently to CDK2 as shown by a co-complex crystal structure. Acute incubation with NU6300 produced a durable inhibition of Rb phosphorylation in SKUT-1B cells, consistent with it acting as an irreversible CDK2 inhibitor. NU6300 is the first covalent CDK2 inhibitor to be described, and illustrates the potential of vinyl sulfones for the design of more potent and selective compounds.
Collapse
Affiliation(s)
- Elizabeth Anscombe
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | - Elisa Meschini
- Newcastle Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Bedson Building, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - Regina Mora-Vidal
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Mathew P Martin
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - David Staunton
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | | | - U Helena Danielson
- Beactica AB, Box 567, 751 22 Uppsala, Sweden; Department of Chemistry-BMC, Uppsala University, 751 23 Uppsala, Sweden
| | - Will A Stanley
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Lan Z Wang
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Tristan Reuillon
- Newcastle Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Bedson Building, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - Bernard T Golding
- Newcastle Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Bedson Building, Newcastle University, Newcastle upon Tyne NE1 7RU, UK.
| | - Celine Cano
- Newcastle Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Bedson Building, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | - David R Newell
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Martin E M Noble
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK
| | - Stephen R Wedge
- Newcastle Cancer Centre, Northern Institute for Cancer Research, Paul O'Gorman Building, Medical School, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
| | - Jane A Endicott
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, UK.
| | - Roger J Griffin
- Newcastle Cancer Centre, Northern Institute for Cancer Research, School of Chemistry, Bedson Building, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| |
Collapse
|
42
|
Wang L, Wu Y, Deng Y, Kim B, Pierce L, Krilov G, Lupyan D, Robinson S, Dahlgren MK, Greenwood J, Romero DL, Masse C, Knight JL, Steinbrecher T, Beuming T, Damm W, Harder E, Sherman W, Brewer M, Wester R, Murcko M, Frye L, Farid R, Lin T, Mobley DL, Jorgensen WL, Berne BJ, Friesner RA, Abel R. Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. J Am Chem Soc 2015; 137:2695-703. [PMID: 25625324 DOI: 10.1021/ja512751q] [Citation(s) in RCA: 801] [Impact Index Per Article: 89.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Designing tight-binding ligands is a primary objective of small-molecule drug discovery. Over the past few decades, free-energy calculations have benefited from improved force fields and sampling algorithms, as well as the advent of low-cost parallel computing. However, it has proven to be challenging to reliably achieve the level of accuracy that would be needed to guide lead optimization (∼5× in binding affinity) for a wide range of ligands and protein targets. Not surprisingly, widespread commercial application of free-energy simulations has been limited due to the lack of large-scale validation coupled with the technical challenges traditionally associated with running these types of calculations. Here, we report an approach that achieves an unprecedented level of accuracy across a broad range of target classes and ligands, with retrospective results encompassing 200 ligands and a wide variety of chemical perturbations, many of which involve significant changes in ligand chemical structures. In addition, we have applied the method in prospective drug discovery projects and found a significant improvement in the quality of the compounds synthesized that have been predicted to be potent. Compounds predicted to be potent by this approach have a substantial reduction in false positives relative to compounds synthesized on the basis of other computational or medicinal chemistry approaches. Furthermore, the results are consistent with those obtained from our retrospective studies, demonstrating the robustness and broad range of applicability of this approach, which can be used to drive decisions in lead optimization.
Collapse
Affiliation(s)
- Lingle Wang
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Xu W, Lucke AJ, Fairlie DP. Comparing sixteen scoring functions for predicting biological activities of ligands for protein targets. J Mol Graph Model 2015; 57:76-88. [PMID: 25682361 DOI: 10.1016/j.jmgm.2015.01.009] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2014] [Revised: 01/22/2015] [Accepted: 01/23/2015] [Indexed: 12/17/2022]
Abstract
Accurately predicting relative binding affinities and biological potencies for ligands that interact with proteins remains a significant challenge for computational chemists. Most evaluations of docking and scoring algorithms have focused on enhancing ligand affinity for a protein by optimizing docking poses and enrichment factors during virtual screening. However, there is still relatively limited information on the accuracy of commercially available docking and scoring software programs for correctly predicting binding affinities and biological activities of structurally related inhibitors of different enzyme classes. Presented here is a comparative evaluation of eight molecular docking programs (Autodock Vina, Fitted, FlexX, Fred, Glide, GOLD, LibDock, MolDock) using sixteen docking and scoring functions to predict the rank-order activity of different ligand series for six pharmacologically important protein and enzyme targets (Factor Xa, Cdk2 kinase, Aurora A kinase, COX-2, pla2g2a, β Estrogen receptor). Use of Fitted gave an excellent correlation (Pearson 0.86, Spearman 0.91) between predicted and experimental binding only for Cdk2 kinase inhibitors. FlexX and GOLDScore produced good correlations (Pearson>0.6) for hydrophilic targets such as Factor Xa, Cdk2 kinase and Aurora A kinase. By contrast, pla2g2a and COX-2 emerged as difficult targets for scoring functions to predict ligand activities. Although possessing a high hydrophobicity in its binding site, β Estrogen receptor produced reasonable correlations using LibDock (Pearson 0.75, Spearman 0.68). These findings can assist medicinal chemists to better match scoring functions with ligand-target systems for hit-to-lead optimization using computer-aided drug design approaches.
Collapse
Affiliation(s)
- Weijun Xu
- Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Andrew J Lucke
- Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia
| | - David P Fairlie
- Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia.
| |
Collapse
|
44
|
Exploring conformational search protocols for ligand-based virtual screening and 3-D QSAR modeling. J Comput Aided Mol Des 2014; 29:165-82. [DOI: 10.1007/s10822-014-9813-4] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2014] [Accepted: 11/06/2014] [Indexed: 10/24/2022]
|
45
|
Polyphony: superposition independent methods for ensemble-based drug discovery. BMC Bioinformatics 2014; 15:324. [PMID: 25265915 PMCID: PMC4261739 DOI: 10.1186/1471-2105-15-324] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2013] [Accepted: 09/17/2014] [Indexed: 12/04/2022] Open
Abstract
Background Structure-based drug design is an iterative process, following cycles of structural biology, computer-aided design, synthetic chemistry and bioassay. In favorable circumstances, this process can lead to the structures of hundreds of protein-ligand crystal structures. In addition, molecular dynamics simulations are increasingly being used to further explore the conformational landscape of these complexes. Currently, methods capable of the analysis of ensembles of crystal structures and MD trajectories are limited and usually rely upon least squares superposition of coordinates. Results Novel methodologies are described for the analysis of multiple structures of a protein. Statistical approaches that rely upon residue equivalence, but not superposition, are developed. Tasks that can be performed include the identification of hinge regions, allosteric conformational changes and transient binding sites. The approaches are tested on crystal structures of CDK2 and other CMGC protein kinases and a simulation of p38α. Known interaction - conformational change relationships are highlighted but also new ones are revealed. A transient but druggable allosteric pocket in CDK2 is predicted to occur under the CMGC insert. Furthermore, an evolutionarily-conserved conformational link from the location of this pocket, via the αEF-αF loop, to phosphorylation sites on the activation loop is discovered. Conclusions New methodologies are described and validated for the superimposition independent conformational analysis of large collections of structures or simulation snapshots of the same protein. The methodologies are encoded in a Python package called Polyphony, which is released as open source to accompany this paper [http://wrpitt.bitbucket.org/polyphony/].
Collapse
|
46
|
Hardy T, Lee M, Hames RS, Prosser SL, Cheary DM, Samant MD, Schultz F, Baxter JE, Rhee K, Fry AM. Multisite phosphorylation of C-Nap1 releases it from Cep135 to trigger centrosome disjunction. J Cell Sci 2014; 127:2493-506. [PMID: 24695856 PMCID: PMC4038944 DOI: 10.1242/jcs.142331] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2013] [Accepted: 03/12/2014] [Indexed: 01/17/2023] Open
Abstract
During mitotic entry, centrosomes separate to establish the bipolar spindle. Delays in centrosome separation can perturb chromosome segregation and promote genetic instability. However, interphase centrosomes are physically tethered by a proteinaceous linker composed of C-Nap1 (also known as CEP250) and the filamentous protein rootletin. Linker disassembly occurs at the onset of mitosis in a process known as centrosome disjunction and is triggered by the Nek2-dependent phosphorylation of C-Nap1. However, the mechanistic consequences of C-Nap1 phosphorylation are unknown. Here, we demonstrate that Nek2 phosphorylates multiple residues within the C-terminal domain of C-Nap1 and, collectively, these phosphorylation events lead to loss of oligomerization and centrosome association. Mutations in non-phosphorylatable residues that make the domain more acidic are sufficient to release C-Nap1 from the centrosome, suggesting that it is an increase in overall negative charge that is required for this process. Importantly, phosphorylation of C-Nap1 also perturbs interaction with the core centriolar protein, Cep135, and interaction of endogenous C-Nap1 and Cep135 proteins is specifically lost in mitosis. We therefore propose that multisite phosphorylation of C-Nap1 by Nek2 perturbs both oligomerization and Cep135 interaction, and this precipitates centrosome disjunction at the onset of mitosis.
Collapse
Affiliation(s)
- Tara Hardy
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Miseon Lee
- Department of Biological Sciences, Seoul National University, Seoul 151-747, Republic of Korea
| | - Rebecca S Hames
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Suzanna L Prosser
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Donna-Marie Cheary
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Mugdha D Samant
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Francisca Schultz
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Joanne E Baxter
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| | - Kunsoo Rhee
- Department of Biological Sciences, Seoul National University, Seoul 151-747, Republic of Korea
| | - Andrew M Fry
- Department of Biochemistry, University of Leicester, Lancaster Road, Leicester LE1 9HN, UK
| |
Collapse
|
47
|
Liu H, Li L, Qurat-ul-ain S, Jiang T. Synthesis and Inhibitory Activity Evaluation of 2,6-Disubstituted Purine Derivatives. J Heterocycl Chem 2014. [DOI: 10.1002/jhet.1959] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Hongxia Liu
- Key Laboratory of Marine Drugs, Chinese Ministry of Education, Shandong Provincial, School of Medicine and Pharmacy; Ocean University of China; Qingdao Shandong 266003 People's Republic of China
- College of Chemistry and Chemical Engineering; Qiqihar University; Qiqihar Heilongjiang 161006 China
| | - Libo Li
- Department of Pharmacology; Qiqihar Medical University; Qiqihar Heilongjiang 161006 People's Republic of China
| | - Shaikh Qurat-ul-ain
- Key Laboratory of Marine Drugs, Chinese Ministry of Education, Shandong Provincial, School of Medicine and Pharmacy; Ocean University of China; Qingdao Shandong 266003 People's Republic of China
| | - Tao Jiang
- Key Laboratory of Marine Drugs, Chinese Ministry of Education, Shandong Provincial, School of Medicine and Pharmacy; Ocean University of China; Qingdao Shandong 266003 People's Republic of China
| |
Collapse
|
48
|
Mojzych M, Šubertová V, Bielawska A, Bielawski K, Bazgier V, Berka K, Gucký T, Fornal E, Kryštof V. Synthesis and kinase inhibitory activity of new sulfonamide derivatives of pyrazolo[4,3-e][1,2,4]triazines. Eur J Med Chem 2014; 78:217-24. [PMID: 24681986 DOI: 10.1016/j.ejmech.2014.03.054] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2013] [Revised: 03/13/2014] [Accepted: 03/16/2014] [Indexed: 11/19/2022]
Abstract
A new series of sulfonamide derivatives of pyrazolo[4,3-e][1,2,4]triazine has been synthesized and characterized. Their anticancer activity was tested in vitro against multiple human cancer cell lines and were found to have dose-dependent antiproliferative effects. Furthermore, some of the new compounds inhibited the Abl protein kinase with low micromolar IC50 values and exhibited selective activity against the Bcr-Abl positive K562 and BV173 cell lines, providing starting points for the further development of this new kinase inhibitor scaffold.
Collapse
Affiliation(s)
- Mariusz Mojzych
- Department of Chemistry, Siedlce University of Natural Sciences and Humanities, ul. 3 Maja 54, Siedlce 08-110, Poland
| | - Veronika Šubertová
- Centre of the Region Haná for Biotechnological and Agricultural Research, Laboratory of Growth Regulators, Faculty of Science, Palacký University, Šlechtitelů 11, 783 71 Olomouc, Czech Republic
| | - Anna Bielawska
- Department of Medicinal Chemistry and Drug Technology, Medical University of Bialystok, Bialystok, Poland
| | - Krzysztof Bielawski
- Department of Medicinal Chemistry and Drug Technology, Medical University of Bialystok, Bialystok, Poland
| | - Václav Bazgier
- Centre of the Region Haná for Biotechnological and Agricultural Research, Laboratory of Growth Regulators, Faculty of Science, Palacký University, Šlechtitelů 11, 783 71 Olomouc, Czech Republic
| | - Karel Berka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Tomáš Gucký
- Centre of the Region Haná for Biotechnological and Agricultural Research, Laboratory of Growth Regulators, Faculty of Science, Palacký University, Šlechtitelů 11, 783 71 Olomouc, Czech Republic
| | - Emilia Fornal
- Department of Chemistry, Laboratory of Separation and Spectroscopic Method Applications, Center for Interdisciplinary Research, The John Paul II Catholic University of Lublin, al. Krasnicka 102, 20-718 Lublin, Poland
| | - Vladimír Kryštof
- Centre of the Region Haná for Biotechnological and Agricultural Research, Laboratory of Growth Regulators, Faculty of Science, Palacký University, Šlechtitelů 11, 783 71 Olomouc, Czech Republic.
| |
Collapse
|
49
|
Giese T, Chen H, Huang M, York DM. Parametrization of an Orbital-Based Linear-Scaling Quantum Force Field for Noncovalent Interactions. J Chem Theory Comput 2014; 10:1086-1098. [PMID: 24803856 PMCID: PMC3985928 DOI: 10.1021/ct401035t] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Indexed: 01/22/2023]
Abstract
We parametrize a linear-scaling quantum mechanical force field called mDC for the accurate reproduction of nonbonded interactions. We provide a new benchmark database of accurate ab initio interactions between sulfur-containing molecules. A variety of nonbond databases are used to compare the new mDC method with other semiempirical, molecular mechanical, ab initio, and combined semiempirical quantum mechanical/molecular mechanical methods. It is shown that the molecular mechanical force field significantly and consistently reproduces the benchmark results with greater accuracy than the semiempirical models and our mDC model produces errors twice as small as the molecular mechanical force field. The comparisons between the methods are extended to the docking of drug candidates to the Cyclin-Dependent Kinase 2 protein receptor. We correlate the protein-ligand binding energies to their experimental inhibition constants and find that the mDC produces the best correlation. Condensed phase simulation of mDC water is performed and shown to produce O-O radial distribution functions similar to TIP4P-EW.
Collapse
Affiliation(s)
- Timothy
J. Giese
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Haoyuan Chen
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | - Ming Huang
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
- Scientific
Computation, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455−0431, United States
| | - Darrin M. York
- BioMaPS
Institute and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| |
Collapse
|
50
|
Carbain B, Paterson DJ, Anscombe E, Campbell AJ, Cano C, Echalier A, Endicott JA, Golding BT, Haggerty K, Hardcastle IR, Jewsbury PJ, Newell DR, Noble MEM, Roche C, Wang LZ, Griffin RJ. 8-Substituted O(6)-cyclohexylmethylguanine CDK2 inhibitors: using structure-based inhibitor design to optimize an alternative binding mode. J Med Chem 2014; 57:56-70. [PMID: 24304238 DOI: 10.1021/jm401555v] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Evaluation of the effects of purine C-8 substitution within a series of CDK1/2-selective O(6)-cyclohexylmethylguanine derivatives revealed that potency decreases initially with increasing size of the alkyl substituent. Structural analysis showed that C-8 substitution is poorly tolerated, and to avoid unacceptable steric interactions, these compounds adopt novel binding modes. Thus, 2-amino-6-cyclohexylmethoxy-8-isopropyl-9H-purine adopts a "reverse" binding mode where the purine backbone has flipped 180°. This provided a novel lead chemotype from which we have designed more potent CDK2 inhibitors using, in the first instance, quantum mechanical energy calculations. Introduction of an ortho-tolyl or ortho-chlorophenyl group at the purine C-8 position restored the potency of these "reverse" binding mode inhibitors to that of the parent 2-amino-6-cyclohexylmethoxy-9H-purine. By contrast, the corresponding 8-(2-methyl-3-sulfamoylphenyl)-purine derivative exhibited submicromolar CDK2-inhibitory activity by virtue of engineered additional interactions with Asp86 and Lys89 in the reversed binding mode, as confirmed by X-ray crystallography.
Collapse
Affiliation(s)
- Benoit Carbain
- Department of Biochemistry, University of Oxford , South Parks Road, Oxford OX1 3QU, U.K
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|