Design, synthesis and biological evaluation of 4-(4-aminophenoxy)picolinamide derivatives as potential antitumor agents

Cancer is a leading cause of death in humans. Molecular targeted therapy for cancer has become a research hotspot as it is associated with low toxicity and high efficiency. In this study, a total of 36 derivatives of 4-(4-aminophenoxy)pyridinamide were designed and synthesized, based on the analysis of the binding patterns of cabozantinib and BMS-777607 to MET protein. Most target compounds exhibited moderate to excellent antiproliferative activity against three different cell lines (A549, HeLa and MCF-7). A total of 7 compounds had stronger inhibitory activities than cabozantinib, and the IC₅₀ value of the most promising compound 46 was 0.26 μM against the A549 cells, which was 2.4 times more active than that of cabozantinib. The structure-activity relationship of the target compounds was analyzed and summarized, and the action mechanism was discussed. The acridine orange (AO) staining assay and cell cycle apoptosis revealed that compound 46 dose-dependently induced apoptosis of A549 cells, and blocked the cells mainly in G0/G1 phase. The IC₅₀ value of compound 46 on c-Met kinase was 46.5 nM. Further docking studies and molecular dynamics simulations signaled that compound 46 formed four key hydrogen bonds to c-Met kinase, and these key amino acids played a major role in binding free energy. In addition, compound 46 also showed good pharmacokinetic characteristics in rats. In conclusion, compound 46 is a promising antitumor agent.

The Conformational Transition Pathways and Hidden Intermediates in DFG-Flip Process of c-Met Kinase Revealed by Metadynamics Simulations

Protein kinases intrinsically translate their conformations between active and inactive states, which is key to their enzymatic activities. The conformational flipping of the three-residue conservative motif, Asp-Phe-Gly (DFG), is crucial for many kinases' biological functions. Obtaining a detailed demonstration of the DFG flipping process and its corresponding dynamical and thermodynamical features could broaden our understanding of kinases' conformation-activity relationship. In this study, we employed metadynamics simulation, a widely used enhanced sampling technique, to analyze the conformational transition pathways of the DFG flipping for the c-Met kinase. The corresponding free energy landscape suggested two distinct transition pathways between the "DFG-in" and "DFG-out" states of the DFG-flip from c-Met. On the basis of the orientation direction of the F1223 residue, we correspondingly named the two pathways the "DFG-up" path, featuring forming a commonly discovered "DFG-up" transition state, and the "DFG-down" path, a unique transition pathway with F1223 rotating along the opposite direction away from the hydrophobic cavity. The free energies along the two pathways were then calculated using the Path Collective Variable (PCV) metadynamics simulation. The simulation results showed that, though having similar free energy barriers, the free energy cuve for the DFG-down path suggested a two-step conformational transition mechanism, while that for the DFG-up path showed the one-step transition feature. The c-Met DFG flipping mechanism and the new intermediate state discovered in this work could provide a deeper understanding of the conformation-activity relationship for c-Met and, possibly, reveal a new conformational state as the drug target for c-Met and other similar kinases.

Lead generation of cysteine based mesenchymal epithelial transition (c-Met) kinase inhibitors: Using structure-based scaffold hopping, 3D-QSAR pharmacophore modeling, virtual screening, molecular docking, and molecular dynamics simulation

Cysteine-based mesenchymal-epithelial transition (c-Met) is a receptor tyrosine kinase that plays a definitive role during cancer progression and was identified as a possible target for anti-angiogenesis drugs. In the present study, different protocols of computer-based drug design were performed. Construction of predictive pharmacophore model using HypoGen algorithm resulted in a validated model of four features of positive ionizable, hydrogen bond acceptor, hydrophobic, and ring aromatic features with a correlation coefficient of 0.87, a configuration cost of 14.95, and a cost difference of 357.92. The model revealed a promising predictive power and had >90% probability of representing true correlation with the activity data. The model was established using Fisher's validation test at the 95% confidence level and test set prediction (r = 0.96), furthermore, the model was validated by mapping of set of compounds undergoing clinical trials as class Ⅱ c-met inhibitors. The generated valid pharmacophore model was then anticipated for virtual screening of three data bases. Moreover, scaffold hopping using replace fragments protocol was implemented. Hits generated were filtered according to Lipinski's rule; 510 selected hits were anatomized and subjected to molecular docking studies into the crystal structure of c-Met kinase. The good correlation between docking scores and ligand pharmacophore mapping fit values provided a reliable foundation for designing new potentially active candidates that may target c-Met kinase. Eventually, eight hits were selected as potential leads. Subsequently, seven (Hits) have displayed a higher dock score and demonstrated key residue interactions with stable molecular dynamics simulation. Therefore, these c-Met kinase inhibitors may further serve as new chemical spaces in designing new compounds.

TYROSINE KINASES: COMPLEX MOLECULAR SYSTEMS CHALLENGING COMPUTATIONAL METHODOLOGIES

Classical molecular dynamics (MD) simulations based on atomic models play an increasingly important role in a wide range of applications in physics, biology, and chemistry. Nonetheless, generating genuine knowledge about biological systems using MD simulations remains challenging. Protein tyrosine kinases are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Due to the large conformational changes and long timescales involved in their function, these kinases present particularly challenging problems to modern computational and theoretical frameworks aimed at elucidating the dynamics of complex biomolecular systems. Markov state models have achieved limited success in tackling the broader conformational ensemble and biased methods are often employed to examine specific long timescale events. Recent advances in machine learning continue to push the limitations of current methodologies and provide notable improvements when integrated with the existing frameworks. A broad perspective is drawn from a critical review of recent studies.

Comprehensive virtual screening of 4.8 k flavonoids reveals novel insights into allosteric inhibition of SARS-CoV-2 MPRO

SARS-CoV-2 main protease is a common target for inhibition assays due to its high conservation among coronaviruses. Since flavonoids show antiviral activity, several in silico works have proposed them as potential SARS-CoV-2 main protease inhibitors. Nonetheless, there is reason to doubt certain results given the lack of consideration for flavonoid promiscuity or main protease plasticity, usage of short library sizes, absence of control molecules and/or the limitation of the methodology to a single target site. Here, we report a virtual screening study where dorsilurin E, euchrenone a11, sanggenol O and CHEMBL2171598 are proposed to inhibit main protease through different pathways. Remarkably, novel structural mechanisms were observed after sanggenol O and CHEMBL2171598 bound to experimentally proven allosteric sites. The former drastically affected the active site, while the latter triggered a hinge movement which has been previously reported for an inactive SARS-CoV main protease mutant. The use of a curated database of 4.8 k flavonoids, combining two well-known docking software (AutoDock Vina and AutoDock4.2), molecular dynamics and MMPBSA, guaranteed an adequate analysis and robust interpretation. These criteria can be considered for future screening campaigns against SARS-CoV-2 main protease.

Research Progresses of Targeted Therapy and Immunotherapy for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide, with nearly one million new cases and deaths every year. Owing to the complex pathogenesis, hidden early symptoms, rapidly developing processes, and poor prognosis, the morbidity and mortality of HCC are increasing yearly. With the progress being made in modern medicine, the treatment of HCC is no longer limited to traditional methods. Targeted therapy and immunotherapy have emerged to treat advanced and metastatic HCC in recent years. Since Sorafenib is the first molecular targeting drug against angiogenesis, targeted drugs for HCC are continually emerging. Moreover, immunotherapy plays a vital role in clinical trials. In particular, the application of immune checkpoint inhibitors, which have received increasing attention in the field of cancer treatment, is a possible research path. Interestingly, these two therapies generally complement each other at some stages of HCC, bringing new hope for patients with advanced HCC. In this paper, we discuss the research progress of targeted therapy and immunotherapy for HCC in recent years, which will provide a reference for the further development of drugs for HCC.

Synthesis of triazolotriazine derivatives as c-Met inhibitors

Receptor tyrosine kinase c-Met is an important antitumor drug target. Triazolotriazine analogues 2-10 were prepared efficiently and evaluated the enzymatic and cellular c-Met activities. Brief structure-activity relationships of triazolotriazine core and CF₂-quinoline part were investigated, leading to the discovery of compound 8 with nanomolar enzymatic c-Met activity, and subnanomolar MKN45 and EBC-1 cellular potencies. The proposed binding model of 8 and c-Met unraveled that two canonical hydrogen bonds and a π-π stacking interaction formed between the inhibitor and the ATP binding site of c-Met kinase domain, which accounted for its potent c-Met activities.

Large-Scale Virtual Screening Against the MET Kinase Domain Identifies a New Putative Inhibitor Type

By using an ensemble-docking strategy, we undertook a large-scale virtual screening campaign in order to identify new putative hits against the MET kinase target. Following a large molecular dynamics sampling of its conformational space, a set of 45 conformers of the kinase was retained as docking targets to take into account the flexibility of the binding site moieties. Our screening funnel started from about 80,000 chemical compounds to be tested in silico for their potential affinities towards the kinase binding site. The top 100 molecules selected—thanks to the molecular docking results—were further analyzed for their interactions, and 25 of the most promising ligands were tested for their ability to inhibit MET activity in cells. F0514-4011 compound was the most efficient and impaired this scattering response to HGF (Hepatocyte Growth Factor) with an IC50 of 7.2 μM. Interestingly, careful docking analysis of this molecule with MET suggests a possible conformation halfway between classical type-I and type-II MET inhibitors, with an additional region of interaction. This compound could therefore be an innovative seed to be repositioned from its initial antiviral purpose towards the field of MET inhibitors. Altogether, these results validate our ensemble docking strategy as a cost-effective functional method for drug development.

Tyrosine Kinase Activation and Conformational Flexibility: Lessons from Src-Family Tyrosine Kinases

Protein kinases are enzymes that catalyze the covalent transfer of the γ-phosphate of an adenosine triphosphate (ATP) molecule onto a tyrosine, serine, threonine, or histidine residue in the substrate and thus send a chemical signal to networks of downstream proteins. They are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Unregulated protein kinase activity is often associated with a wide range of diseases, therefore making protein kinases major therapeutic targets. A prototypical system of central interest to understand the regulation of kinase activity is provided by tyrosine kinase c-Src, which belongs to the family of Src-related non-receptor tyrosine kinases (SFKs). Although the broad picture of autoinhibition via the regulatory domains and via the phosphorylation of the C-terminal tail is well characterized from a structural point of view, a detailed mechanistic understanding at the atomic-level is lacking. Advanced computational methods based on all-atom molecular dynamics (MD) simulations are employed to advance our understanding of tyrosine kinase activation. The computational studies suggest that the isolated kinase domain (KD) is energetically most favorable in the inactive conformation when the activation loop (A-loop) of the KD is not phosphorylated. The KD makes transient visits to a catalytically competent active-like conformation. The process of bimolecular trans-autophosphorylation of the A-loop eventually locks the KD in the active state. Activating point mutations may act by slightly increasing the population of the active-like conformation, enhancing the availability of the A-loop to be phosphorylated. The Src-homology 2 (SH2) and Src-homology 3 (SH3) regulatory domains, depending upon their configuration, either promote the inactive or the active state of the kinase domain. In addition to the roles played by the SH3, SH2, and KD, the Src-homology 4-Unique domain (SH4-U) region also serves as a key moderator of substrate specificity and kinase function. Thus, a fundamental understanding of the conformational propensity of the SH4-U region and how this affects the association to the membrane surface are likely to lead to the discovery of new intermediate states and alternate strategies for inhibition of kinase activity for drug discovery. The existence of a multitude of KD conformations poses a great challenge aimed at the design of specific inhibitors. One promising computational strategy to explore the conformational flexibility of the KD is to construct Markov state models from aggregated MD data.