Molecular dynamics simulation and free energy calculation studies of kinase inhibitors binding to active and inactive conformations of VEGFR-2

Multiscale biomechanics and mechanotransduction from liver fibrosis to cancer

A growing body of multiscale biomechanical studies has been proposed to highlight the mechanical cues in the development of hepatic fibrosis and cancer. At the cellular level, changes in mechanical microenvironment induce phenotypic and functional alterations of hepatic cells, initiating a positive feedback loop that promotes liver fibrogenesis and hepatocarcinogenesis. Tumor mechanical microenvironment of hepatocellular carcinoma facilitates tumor cell growth and metastasis, and hinders the drug delivery and immunotherapy. At the molecular level, mechanical forces are sensed and transmitted into hepatic cells via allosteric activation of mechanoreceptors on the cell membrane, leading to the activation of various mechanotransduction pathways including integrin and YAP signaling and then regulating cell function. Thus, the application of mechanomedicine concept in the treatment of liver diseases is promising for rational design and cell-specific delivery of therapeutic drugs. This review mainly discusses the correlation between biomechanical cues and liver diseases from the viewpoint of mechanobiology.

Identification of novel allosteric binding sites and multi-targeted allosteric inhibitors of receptor and non-receptor tyrosine kinases using a computational approach

EGFR1, VEGFR2, Bcr-Abl and Src kinases are key drug targets in non-small cell lung cancer (NSCLC), bladder cancer, pancreatic cancer, CML, ALL, colorectal cancer, etc. The available drugs targeting these kinases have limited therapeutic efficacy due to novel mutations resulting in drug resistance and toxicity, as they target ATP binding site. Allosteric drugs have shown promising results in overcoming drug resistance, but the discovery of allosteric drugs is challenging. The allosteric binding pockets are difficult to predict, as they are generally associated with high energy conformations and regulate protein function in yet unknown mechanisms. In addition, the discovery of drugs using conventional methods takes long time and goes through several challenges, putting the lives of many cancer patients at risk. Therefore, the aim of the present work was to apply the most successful, drug repurposing approach in combination with computational methods to identify kinase inhibitors targeting novel allosteric sites on protein structure and assess their potential multi-kinase binding affinity. Multiple crystal structures belonging to EGFR1, VEGFR2, Bcr-Abl and Src tyrosine kinases were selected, including mutated, inhibitor bound and allosteric conformations to identify potential leads, close to physiological conditions. Interestingly the potential inhibitors identified were peptides. The drugs identified in this study could be used in therapy as a single multi-kinase inhibitor or in a combination of single kinase inhibitors after experimental validation. In addition, we have also identified new hot spots that are likely to be druggable allosteric sites for drug discovery of kinase-specific drugs in the future.Communicated by Ramaswamy H. Sarma.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild-type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the fluctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.

Gaining insight into crizotinib resistance mechanisms caused by L2026M and G2032R mutations in ROS1 via molecular dynamics simulations and free-energy calculations

ROS1 fusion kinase-highly expressed in a variety of human cancers-has emerged as an important and attractive target for anticancer drug design. Crizotinib, a well-known drug approved by the FDA as an ALK inhibitor to treat advanced NSCLC, also shows potent inhibitoy activity against ROS1. However, the development of serious resistance due to secondary mutations has been observed in clinical studies. To provide insight into the mechanisms of this drug resistance, molecular dynamics simulations and free-energy calculations were carried out for complexes of crizotinib with wild-type (WT) ROS1 as well as the mutated L2026M and G2032R forms. MD simulations indicated that the L2026M and G2032R systems are slightly less flexible than the WT system. Binding free energy calculations showed that the L2026M and G2032R mutations significantly reduce the binding affinity of crizotinib for ROS1, and that the resistance to crizotinib caused by the L2026M and G2032R mutations arises mostly from increases in entropic terms. Furthermore, calculations of per-residue binding free energies highlighted increased and decreased contributions of some residues in the L2026M and G2032R systems relative to those in the WT system. The present study therefore yielded detailed insight into the mechanisms of resistance to crizotinib caused by the L2026M and G2032R mutations, which should provide the basis for rational drug design to combat crizotinib resistance. Graphical Abstract Superposition of the average structures obtained from the last 10 ns of the molecular dynamics simulation trajectoriy for WT (green) and mutated ROS1 (cyan).

Bisarylureas Based on 1H-Pyrazolo[3,4-d]pyrimidine Scaffold as Novel Pan-RAF Inhibitors with Potent Anti-Proliferative Activities: Structure-Based Design, Synthesis, Biological Evaluation and Molecular Modelling Studies

RAF (Ras activating factor) kinases are important and attractive targets for cancer therapy. With the aim of discovering RAF inhibitors that bind to DFG-out inactive conformation created by the movement of Asp-Phe-Gly (DFG), we conducted structure-based drug design using the X-ray cocrystal structures of BRAF (v-raf murine sarcoma viral oncogene homolog B1), starting from bisarylurea derivative based on 1H-pyrazolo[3,4-d]pyrimidine scaffold 1a. Most of the synthesized compounds showed good to excellent inhibitory activities against BRAF^V600E kinase, possessed moderate to potent anti-proliferative activities against four tumor cell lines (A375, HT-29, PC-3 and A549) and good selectivity towards cancer cells rather normal cells (Madin-Darby canine kidney, MDCK). The most promising compound, 1v, exhibited potent inhibitory activity against not only BRAF^V600E (half maximal inhibitory concentration, IC₅₀ = 23.6 nM) but also wild-type BRAF (IC₅₀ = 51.5 nM) and C-RAF (IC₅₀ = 8.5 nM), and effective cellular anti-proliferative activities against A375, HT-29, PC-3 and A549 cell lines as well as a very good selectivity profile. Moreover, compound 1v mainly arrested the A375 cell line in the G0/G1 stage, and showed significant suppression of MEK (mitogen-activated protein kinase kinase) phosphorylation in A375 and HT-29 cell lines. Taken together, the optimal compound 1v showed excellent in vitro potency as a pan-RAF inhibitor. In addition, the promise of compound 1v was further confirmed by molecular dynamics simulation and binding free energy calculations.

Affinity enhancement of nanobody binding to EGFR: in silico site-directed mutagenesis and molecular dynamics simulation approaches

Epidermal growth factor receptor (EGFR), a transmembrane glycoprotein, is overexpressed in many cancers such as head-neck, breast, prostate, and skin cancers for this reason it is a good target in cancer therapy and diagnosis. In nanobody-based cancer diagnosis and treatment, nanobodies with high affinity toward receptor (e.g. EGFR) results in effective treatment or diagnosis of cancer. In this regard, the main aim of this study is to develop a method based on molecular dynamic (MD) simulations for designing of 7D12 based nanobody with high affinity compared with wild-type nanobody. By surveying electrostatic and desolvation interactions between different residues of 7D12 and EGFR, the critical residues of 7D12 that play the main role in the binding of 7D12 to EGFR were elucidated and based on these residues, five logical variants were designed. Following the 50 ns MD simulations, pull and umbrella sampling simulation were performed for 7D12 and all its variants in complex with EGFR. Binding free energy of 7D12 (and all its variants) with EGFR was obtained by weighted histogram analysis method. According to binding free energy results, GLY101 to GLU mutation showed the highest binding affinity but this variant is unstable after 50 ns MD simulations. ALA100 to GLU mutation shows suitable binding enhancement with acceptable structural stability. Suitable position and orientation of GLU in residue 100 of 7D12 against related amino acids of EGFR formed some extra hydrogen and electrostatic interactions which resulted in binding enhancement.

Inherent formulation issues of kinase inhibitors

The small molecular Kinase Inhibitor (smKI) drug class is very promising and rapidly expanding. All of these drugs are administered orally. The clear relationship between structure and function has led to drugs with a general low intrinsic solubility. The majority of the commercial pharmaceutical formulations of the smKIs are physical mixtures that are limited by the low drug solubility of a salt form. This class of drugs is therefore characterized by an impaired and variable bioavailability rendering them costly and their therapies suboptimal. New formulations are sparingly being reported in literature and patents. The presented data suggests that continued research into formulation design can help to develop more efficient and cost-effective smKI formulation. Moreover, it may also be of help in the future design of the formulations of new smKIs.

Advanced molecular dynamics simulation methods for kinase drug discovery

Interest in the application of molecular dynamics (MD) simulations has increased in the field of protein kinase (PK) drug discovery. PKs belong to an important drug target class because they are directly involved in a number of diseases, including cancer. MD methods simulate dynamic biological and chemical events at an atomic level. This information can be combined with other in silico and experimental methods to efficiently target selected receptors. In this review, we present common and advanced methods of MD simulations and we focus on the recent applications of MD-based methodologies that provided significant insights into the elucidation of biological mechanisms involving PKs and into the discovery of novel kinase inhibitors.

Discovery of vascular endothelial growth factor receptor tyrosine kinase inhibitors by quantitative structure–activity relationships, molecular dynamics simulation and free energy calculation

Employing the combined strategy to understand the features of KDR–ligands complexes, and provide a basis for rational design of inhibitors.