Application of Monte Carlo-based receptor ensemble docking to virtual screening for GPCR ligands

Identification of abemaciclib derivatives targeting cyclin-dependent kinase 4 and 6 using molecular dynamics, binding free energy calculation, synthesis, and pharmacological evaluation

CDK4/6 plays a crucial role in various cancers and is an effective anticancer drug target. However, the gap between clinical requirements and approved CDK4/6 drugs is unresolved. Thus, there is an urgent need to develop selective and oral CDK4/6 inhibitors, particularly for monotherapy. Here, we studied the interaction between abemaciclib and human CDK6 using molecular dynamics simulations, binding free energy calculations, and energy decomposition. V101 and H100 formed stable hydrogen bonds with the amine-pyrimidine group, and K43 interacted with the imidazole ring via an unstable hydrogen bond. Meanwhile, I19, V27, A41, and L152 interacted with abemaciclib through π-alkyl interactions. Based on the binding model, abemaciclib was divided into four regions. With one region modification, 43 compounds were designed and evaluated using molecular docking. From each region, three favorable groups were selected and combined with each other to obtain 81 compounds. Among them, C2231-A, which was obtained by removing the methylene group from C2231, showed better inhibition than C2231. Kinase profiling revealed that C2231-A showed inhibitory activity similar to that of abemaciclib; additionally, C2231-A inhibited the growth of MDA-MB-231 cells to a greater extent than did abemaciclib. Based on molecular dynamics simulation, C2231-A was identified as a promising candidate compound with considerable inhibitory effects on human breast cancer cell lines.

Molecular Docking, Molecular Dynamics Simulations, and Free Energy Calculation Insights into the Binding Mechanism between VS-4718 and Focal Adhesion Kinase

Focal adhesion kinase (FAK) is a 125 kDa nonreceptor tyrosine kinase that plays an important role in many carcinomas. Thus, the targeting of FAK by small molecules is considered to be promising for cancer therapy. Some FAK inhibitors have been reported as potential anticancer drugs and have entered into clinical development; for example, VS-4718 is currently undergoing clinical trials. However, the lack of crystal structural data for the binding of VS-4718 with FAK has hindered the optimization of this anticancer agent. In this work, the VS-4718/FAK interaction model was obtained by molecular docking and molecular dynamics simulations. The binding free energies of VS-4718/FAK were also calculated using the molecular mechanics generalized Born surface area method. It was found that the aminopyrimidine group formed hydrogen bonds with the C502 residue of the hinge loop, while the D564 residue of the T-loop interacted with the amide group. In addition, I428, A452, V484, M499, G505, and L553 residues formed hydrophobic interactions with VS-4718. The obtained results therefore provide an improved understanding of the interaction between human FAK and VS-4718. Based on the obtained binding mechanism, 47 novel compounds were designed to target the adenosine 5'-triphosphate-binding pocket of human FAK, and ensemble docking was performed to assess the effects of these modifications on the inhibitor binding affinity. This work is also expected to provide additional insights into potential future target design strategies based on VS-4718.

G protein-coupled receptors: structure- and function-based drug discovery

As one of the most successful therapeutic target families, G protein-coupled receptors (GPCRs) have experienced a transformation from random ligand screening to knowledge-driven drug design. We are eye-witnessing tremendous progresses made recently in the understanding of their structure-function relationships that facilitated drug development at an unprecedented pace. This article intends to provide a comprehensive overview of this important field to a broader readership that shares some common interests in drug discovery.

Computational methods in drug discovery

The process for drug discovery and development is challenging, time consuming and expensive. Computer-aided drug discovery (CADD) tools can act as a virtual shortcut, assisting in the expedition of this long process and potentially reducing the cost of research and development. Today CADD has become an effective and indispensable tool in therapeutic development. The human genome project has made available a substantial amount of sequence data that can be used in various drug discovery projects. Additionally, increasing knowledge of biological structures, as well as increasing computer power have made it possible to use computational methods effectively in various phases of the drug discovery and development pipeline. The importance of in silico tools is greater than ever before and has advanced pharmaceutical research. Here we present an overview of computational methods used in different facets of drug discovery and highlight some of the recent successes. In this review, both structure-based and ligand-based drug discovery methods are discussed. Advances in virtual high-throughput screening, protein structure prediction methods, protein-ligand docking, pharmacophore modeling and QSAR techniques are reviewed.

Identifying ligands at orphan GPCRs: current status using structure-based approaches

GPCRs are the most successful pharmaceutical targets in history. Nevertheless, the pharmacology of many GPCRs remains inaccessible as their endogenous or exogenous modulators have not been discovered. Tools that explore the physiological functions and pharmacological potential of these 'orphan' GPCRs, whether they are endogenous and/or surrogate ligands, are therefore of paramount importance. Rates of receptor deorphanization determined by traditional reverse pharmacology methods have slowed, indicating a need for the development of more sophisticated and efficient ligand screening approaches. Here, we discuss the use of structure-based ligand discovery approaches to identify small molecule modulators for exploring the function of orphan GPCRs. These studies have been buoyed by the growing number of GPCR crystal structures solved in the past decade, providing a broad range of template structures for homology modelling of orphans. This review discusses the methods used to establish the appropriate signalling assays to test orphan receptor activity and provides current examples of structure-based methods used to identify ligands of orphan GPCRs. Linked Articles This article is part of a themed section on Molecular Pharmacology of G Protein-Coupled Receptors. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v173.20/issuetoc.

Improvements, trends, and new ideas in molecular docking: 2012-2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.

Structure-based virtual screening of the nociceptin receptor: hybrid docking and shape-based approaches for improved hit identification

The antagonist-bound crystal structure of the nociceptin receptor (NOP), from the opioid receptor family, was recently reported along with those of the other opioid receptors bound to opioid antagonists. We recently reported the first homology model of the ‘active-state’ of the NOP receptor, which when docked with ‘agonist’ ligands showed differences in the TM helices and residues, consistent with GPCR activation after agonist binding. In this study, we explored the use of the active-state NOP homology model for structure-based virtual screening to discover NOP ligands containing new chemical scaffolds. Several NOP agonist and antagonist ligands previously reported are based on a common piperidine scaffold. Given the structure–activity relationships for known NOP ligands, we developed a hybrid method that combines a structure-based and ligand-based approach, utilizing the active-state NOP receptor as well as the pharmacophoric features of known NOP ligands, to identify novel NOP binding scaffolds by virtual screening. Multiple conformations of the NOP active site including the flexible second extracellular loop (EL2) loop were generated by simulated annealing and ranked using enrichment factor (EF) analysis and a ligand–decoy dataset containing known NOP agonist ligands. The enrichment factors were further improved by combining shape-based screening of this ligand–decoy dataset and calculation of consensus scores. This combined structure-based and ligand-based EF analysis yielded higher enrichment factors than the individual methods, suggesting the effectiveness of the hybrid approach. Virtual screening of the CNS Permeable subset of the ZINC database was carried out using the above-mentioned hybrid approach in a tiered fashion utilizing a ligand pharmacophore-based filtering step, followed by structure-based virtual screening using the refined NOP active-state models from the enrichment analysis. Determination of the NOP receptor binding affinity of a selected set of top-scoring hits resulted in identification of several compounds with measurable binding affinity at the NOP receptor, one of which had a new chemotype for NOP receptor binding. The hybrid ligand-based and structure-based methodology demonstrates an effective approach for virtual screening that leverages existing SAR and receptor structure information for identifying novel hits for NOP receptor binding. The refined active-state NOP homology models obtained from the enrichment studies can be further used for structure-based optimization of these new chemotypes to obtain potent and selective NOP receptor ligands for therapeutic development.