Elucidating the binding mechanisms of GABA and Muscimol as an avenue to discover novel GABA-mimetic small molecules

Gamma-aminobutyric acid (GABA) signaling is the principal inhibitory pathway in the central nervous system. It is critical in neuronal cell proliferation and fate determination. Any aberration in GABA inhibition results in psychiatric and neurological diseases. Thus, modulating GABAergic neurotransmission has become the basis of drug therapy for psychiatric and several neurological diseases. Though GABA and muscimol are classical inhibitors of GABA receptors, the search for novel inhibitors continues unabated. In this study, the binding mechanism of GABA and muscimol was elucidated and applied in the search for small molecule GABAergic inhibitors using comprehensive computational techniques. It was revealed that a high-affinity binding of GABA and muscimol was mediated by a water molecule involving α1Thr129 and then stabilized by strong interactions including salt bridges with β2Glu155 and α1Arg66 amidst hydrogen bonds, π-π stacking, and π -cation interactions with other residues. The binding of GABA and muscimol was also characterized by stability and deeper penetration into the hydrophobic core of the protein which resulted in conformational changes of the binding pocket and domain, by inducing correlated motions of the residues. Thermodynamics analysis showed GABA and muscimol exhibited total binding free energies of -19.85 ± 8.83 Kcal/mol and -26.55 ± 3.42 Kcal/mol, respectively. A pharmacophore model search, based on the energy contributions of implicating binding residues, resulted in the identification of ZINC68604167, ZINC19735138, ZINC04202466, ZINC00901626, and ZINC01532854 as potential GABA-mimetic compounds from metabolites and natural products libraries. This study has elucidated the binding mechanisms of GABA and muscimol and successfully applied in the identification of GABA-mimetic compounds.Communicated by Ramaswamy H. Sarma.

Iso-mukaadial acetate and ursolic acid acetate bind to Plasmodium Falciparum heat shock protein 70: towards targeting parasite protein folding pathway

Plasmodium falciparum is the most lethal malaria parasite. P. falciparum Hsp70 (PfHsp70) is an essential molecular chaperone (facilitates protein folding) and is deemed a prospective antimalarial drug target. The present study investigates the binding capabilities of select plant derivatives, iso-mukaadial acetate (IMA) and ursolic acid acetate (UAA), against P. falciparum using an in silico docking approach. The interaction between the ligands and PfHsp70 was evaluated using plasmon resonance (SPR) analysis. Molecular docking, binding free energy analysis and molecular dynamics simulations were conducted towards understanding the mechanisms by which the compounds bind to PfHsp70. The molecular docking results revealed ligand flexibilities, conformations and positions of key amino acid residues and protein-ligand interactions as crucial factors accounting for selective inhibition of Hsp70. The simulation results also suggest protein-ligand van der Waals forces as the driving force guiding the interaction of these compounds with PfHsp70. Of the two compounds, UAA and IMA bound to PfHsp70 within the micromolar range based on surface plasmon resonance (SPR) based binding assay. Our findings pave way for future rational design of new selective compounds targeting PfHsp70.

Insight into the corrosion inhibiting potential and anticancer activity of 1-(4-methoxyphenyl)-5-methyl-N'-(2-oxoindolin-3-ylidene)-1H-1,2,3-triazole-4-carbohydrazide via computational approaches

Cancer is a major health concern globally. Orthodox and traditional medicine have actively been explored to manage this disease. Also, corrosion is a natural catastrophe that weakens and deteriorates metallic structures and their alloys causing major structural failures and severe economic implications. Designing and exploring multi-functional materials are beneficial since they are adaptive to different fields including engineering and pharmaceutics. In this study, we examined the anti-corrosion and anti-cancer potentials of 1-(4-methoxyphenyl)-5-methyl-N'-(2-oxoindolin-3-ylidene)-1H-1,2,3-triazole-4-carbohydrazide (MAC) using computational approaches. The molecular reactivity descriptors and charge distribution parameters of MAC were studied in gas and water at density functional theory (DFT) at B3LYP/6-311++G(d,p) theory level. The binding and mechanism of interaction between MAC and iron surface was studied using Monte Carlo (MC) and molecular dynamics (MD) simulation in hydrochloric acid medium. From the DFT, MC, and MD simulations, it was observed that MAC interacted spontaneously with iron surface essentially via van der Waal and electrostatic interactions. The near-parallel alignment of the corrosion inhibitor on iron plane facilitates its adsorption and isolation of the metal surface from the acidic solution. Further, the compound was docked in the binding pocket of anaplastic lymphoma kinase (ALK: 4FNZ) protein to assess its anti-cancer potential. The binding score, pharmacokinetics, and drug-likeness of MAC were compared with the reference drug (Crizotinib). The MAC displayed binding scores of -5.729 kcal/mol while Crizotinib has -3.904 kcal/mol. MD simulation of the complexes revealed that MAC is more stable and exhibits more favourable hydrogen bonding with the ALK receptor's active site than Crizotinib.Communicated by Ramaswamy H. Sarma.

Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds

CONTEXT

[Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved.

METHODS

In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs.

RESULTS

Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

Fragment-based discovery of novel phenyltriazolyl derivatives as allosteric type-I1/2 ALK inhibitors with promising antitumor effects

Based on the high-throughput screening hit BY-1, a series of phenyltriazolyl derivatives were developed. Satisfyingly, most compounds were detected moderate to excellent antitumor effects against Karpas299 and H2228 cells. Among them, 12k bearing 4‑hydroxypiperidinyl group exhibited the optimal activities against tested cells with IC₅₀ values of 51 nM and 175 nM, as well as promising inhibitory effects on ALK^WT (3.7 nM) and ALK^L1196M (6.8 nM). Unlike the conventional type-I ALK inhibitors, molecular models identified 12k as an allosteric type-I^1/2 inhibitor by forming key interactions in both the ATP binding region and the hydrophobic back pocket of ALK. Intriguingly, 12k could dose-dependently induce apoptosis on H2228 cell and inhibit colony formation and tumor cell migration. Taken together, the rationalization of 12k may shed new light on the identification of novel allosteric type-I^1/2 ALK inhibitors.

EDock-ML: A web server for using ensemble docking with machine learning to aid drug discovery

EDock-ML is a web server that facilitates the use of ensemble docking with machine learning to help decide whether a compound is worthwhile to be considered further in a drug discovery process. Ensemble docking provides an economical way to account for receptor flexibility in molecular docking. Machine learning improves the use of the resulting docking scores to evaluate whether a compound is likely to be useful. EDock-ML takes a bottom-up approach in which machine-learning models are developed one protein at a time to improve predictions for the proteins included in its database. Because the machine-learning models are intended to be used without changing the docking and model parameters with which the models were trained, novice users can use it directly without worrying about what parameters to choose. A user simply submits a compound specified by an ID from the ZINC database (Sterling, T.; Irwin, J. J., J Chem Inf Model 2015, 55[11], 2,324-2,337.) or upload a file prepared by a chemical drawing program and receives an output helping the user decide the likelihood of the compound to be active or inactive for a drug target. EDock-ML can be accessed freely at edock-ml.umsl.edu.

Small-molecule activating SIRT6 elicits therapeutic effects and synergistically promotes anti-tumor activity of vitamin D3 in colorectal cancer

Colorectal cancer (CRC) is the leading cause of cancer death; however, targets with broad anti-CRC effects are limited. Sirtuin6 (SIRT6) is a conserved nicotinamide adenine dinucleotide (NAD⁺)-dependent deacetylase that is widely pathologically downregulated in CRC, but its pharmacological effect in CRC remains undefined due to the lack of small-molecule SIRT6 activators. We searched for a compound activating SIRT6 and investigated its anti-CRC effect in various models. Methods: We identified an allosteric SIRT6 activator, MDL-811. Its ability to enhance SIRT6 deacetylation at protein and cellular levels was evaluated by Fluor de Lys (FDL) and western blots. We assessed the proliferation of 26 CRC cell lines and patient-derived organoids (PDOs) treated with MDL-811. In vivo efficacy of MDL-811 was evaluated in HCT116 cell line- and patient-derived xenografts as well as a spontaneous CRC model. RNA sequencing and real-time quantitative PCR assays were performed to analyze gene expression changes in MDL-811-treated HCT116 cells. Along with controls in SIRT6-overexpressing HCT116 cells, the SIRT6-mediated histone H3 deacetylation at the Cytochrome P450 family 24 subfamily A member 1 (CYP24A1) gene locus was assessed by chromatin immunoprecipitation (ChIP) in MDL-811-treated HCT116 cells. A combination therapy against CRC based on the downstream gene of SIRT6 activation was evaluated in cells and mouse models. Results: MDL-811 significantly activated SIRT6 histone H3 deacetylation (H3K9Ac, H3K18Ac, and H3K56Ac) in vitro and had broad antiproliferative effects on diverse CRC cell lines and PDOs. More importantly, the in vivo anti-tumor efficacy of MDL-811 was demonstrated across cell line- and patient-derived xenografts and in the APC^min/+ spontaneous CRC model. Mechanically, we identified a new downstream target gene of SIRT6 in CRC, CYP24A1. Based on these findings, a combination drug strategy with MDL-811 to synergistically enhance the anti-CRC effect of vitamin D₃ was validated in vitro and in vivo. Conclusions: Our data provide proof of concept that targeting SIRT6 using a small-molecule activator is an attractive therapeutic strategy for CRC and that MDL-811 could be a promising lead compound for further preclinical and clinical studies of treatments for CRC.

A series of molecular modeling techniques to reveal selective mechanisms of inhibitors to β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) and β-site amyloid precursor protein cleaving enzyme 2 (BACE2)

Inhibition of β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) has been shown to be an effective treatment for Alzheimer's disease. A wealth of research has focused on finding highly selective small-molecule inhibitors targeting the BACE1 over its close homologue BACE2 to avoid potential side effects. However, given the highly structural similarities of BACE1 and BACE2, designing highly selective BACE1 inhibitors remains a huge challenge. Recently, it has been reported that a potential BACE1 inhibitor named C28 (∼52-fold selectivity) exhibited greater selectivity to BACE1 over BACE2 than the previously reported inhibitors AZD3293 and AZD3839 (∼1.5-fold and 14-fold selectivity). However, few computational studies have been performed to reveal its underlying mechanisms. In this study, a series of molecular modeling techniques were performed to reveal the selective mechanisms. Classical molecular dynamics (cMD) simulations indicated that the major variations appeared to be controlled by overall protein dynamics. Free energy calculations further suggested that the binding affinities of AZD3293 to BACE1 and BACE2 are similar, but the binding affinity of AZD3839 and C28 to BACE1 is much higher than to BACE2, and that the major variations are electrostatic interactions. The protein dynamics and energy differences were further observed in accelerated molecular dynamics (aMD) simulations. In addition, the umbrella sampling simulations revealed the inhibitors' different patterns of dissociation from the binding pockets of BACE1 and BACE2, and that different energy barriers were responsible for the selectivity. The physical principles revealed by this study may facilitate the rational design of more potent BACE1 selective inhibitors. Communicated by Ramaswamy H. Sarma.

Computational investigation of a ternary model of SnoN-SMAD3-SMAD4 complex

The transforming growth factor β (TGFβ) plays an essential role in the regulation of cellular processes such as cell proliferation, migration, differentiation, and apoptosis by association with SMAD transcriptional factors that are regulated by the transcriptional regulator SnoN. The crystal structure of SnoN-SMAD4 reveals that SnoN can adopt two binding modes, the open and closed forms, at the interfaces of SMAD4 subunits. Accumulating evidence indicates that SnoN can interact with both SMAD3 and SMAD4 to form a ternary SnoN-SMAD3-SMAD4 complex in the TGFβ signaling pathway. However, how the interaction of SnoN with the SMAD3 and SMAD4 remains unclear. Here, molecular dynamics (MD) simulations and molecular modeling methods were performed to figure out this issue. The simulations reveal that SnoN^open exists in two, open and semi-closed, conformations. Molecular modeling and MD simulation studies suggest that the SnoN^closed form interferes with the SMAD3-SMAD4 protein; in contract, the SnoN^open can form a stable SnoN-SMAD3-SMAD4 complex. These mechanistic mechanisms may help elucidate the detailed engagement of SnoN with two SMAD3 and SMAD4 transcriptional factors in the regulation of TGFβ signaling pathway.

Design, synthesis and biological evaluations of 2-amino-4-(1-piperidine) pyridine derivatives as novel anti crizotinib-resistant ALK/ROS1 dual inhibitors

ALK and ROS1 kinases have become promising therapeutic targets since Crizotinib was used to treat non-small-cell lung cancer clinically. Aiming to explore new potent inhibitors, a series of 2-amino-4-(1-piperidine) pyridine derivatives that stabilized a novel DFG-shifted conformation in the kinase domain of ALK were designed and synthesized on the base of lead compound A. Biological evaluation highlighted that most of these new compounds could also potently inhibit ROS1 kinase, leading to the promising inhibitors against both ROS1 and ALK. Among them, the representative compound 2e stood out potent anti-proliferative activity against ALK-addicted H3122 and ROS1-addicted HCC78 cell lines (IC₅₀ = 6.27 μM and 10.71 μM, respectively), which were comparable to that of Crizotinib. Moreover, 2e showed impressive enzyme activity against clinically Crizotinib-resistant ALK^L1196M with an IC₅₀ value of 41.3 nM, which was about 2-fold more potent than that of Crizotinib. 2e also showed potent inhibitory activity in about 6-fold superior to Crizotinib (IC₅₀: 104.7 nM vs. 643.5 nM) in Ba/F3 cell line harboring ROS1^G2032R. Furthermore, molecular modeling disclosed that all the representative inhibitors could dock into the active site of ALK and ROS1, which gave a probable explanation of anti Crizotinib-resistant mutants. These results indicated that our work has established a path forward for the generation of anti Crizotinib-resistant ALK/ROS1 dual inhibitors.

Discovery of novel natural compound inhibitors targeting estrogen receptor α by an integrated virtual screening strategy

Estrogen receptor (ER) is a nuclear hormone receptor and plays an important role in mediating the cellular effects of estrogen. ER can be classified into two receptors: estrogen receptor alpha (ERα) and beta (ERβ), and the former is expressed in 50~80% of breast tumors and has been extensively investigated in breast cancer for decades. Excessive exposure to estrogen can obviously stimulate the growth of breast cancers primarily mediated by ERα, and thus anti-estrogen therapies by small molecules are of concern to clinicians and pharmaceutical industry in the treatment of ERα-positive breast cancers. Although a series of estrogen receptor modulators have been developed, these drugs can lead to resistance and side effects. Therefore, the development of small molecule inhibitors with high target specificity has been intensified. In this pursuit, an integrated computer-aided virtual screening technique, including molecular docking and pharmacophore model screening, was used to screen traditional Chinese medicine (TCM) databases. The compounds with high docking scores and fit values were subjected to ADME (adsorption, distribution, metabolism, excretion) and toxicity prediction, and ten hits were identified as potential inhibitors targeting ERα. Molecular docking was used to investigate the binding modes between ERα and three most potent hits, and molecular dynamic simulations were chosen to explore the stability of these complexes. The rank of the predicted binding free energies evaluated by MM/GBSA is consistent with the docking score. These novel scaffolds discovered in the present study can be used as critical starting point in the drug discovery process for treating ERα-positive breast cancer. Graphical abstract .

Drug Discovery Targeting Anaplastic Lymphoma Kinase (ALK)

As a receptor tyrosine kinase of insulin receptor (IR) subfamily, anaplastic lymphoma kinase (ALK) has been validated to play important roles in various cancers, especially anaplastic large cell lymphoma (ALCL), nonsmall cell lung cancer (NSCLC), and neuroblastomas. Currently, five small-molecule inhibitors of ALK, including Crizotinib, Ceritinib, Alectinib, Brigatinib, and Lorlatinib, have been approved by the U.S. Food and Drug Administration (FDA) against ALK-positive NSCLCs. Novel type-I¹/₂ and type-II ALK inhibitors with improved kinase selectivity and enhanced capability to combat drug resistance have also been reported. Moreover, the "proteolysis targeting chimera" (PROTAC) technique has been successfully applied in developing ALK degraders, which opened a new avenue for targeted ALK therapies. This review provides an overview of the physiological and biological functions of ALK, the discovery and development of drugs targeting ALK by focusing on their chemotypes, activity, selectivity, and resistance as well as potential therapeutic strategies to overcome drug resistance.

The structural basis of the autoinhibition mechanism of glycogen synthase kinase 3β (GSK3β): molecular modeling and molecular dynamics simulation studies

The autoinhibition phenomenon has been frequently observed in enzymes and represents an important regulatory strategy to fine-tune enzyme activity. Evolution has exploited this mechanism to reduce enzymatic activity. Glycogen synthase kinase 3β (GSK3β) undergoes autoinhibition via the phosphorylation of Ser9 at the N-terminus of the kinase, which, acting as a pseudosubstrate, occupies the catalytic domain of GSK3β and subsequently blocks primed substrates from having access to the catalytic domain. The detailed structural basis of the autoinhibition mechanism of GSK3β by the pseudosubstrate, however, has not yet been fully resolved. Here, a three-dimensional model of the binary GSK3β-pseudosubstrate complex was built via the molecular modeling method. Based on the constructed model, extensive molecular dynamics (MD) simulations and subsequent molecular mechanics generalized Born/surface area (MM_GBSA) calculations were performed on the wild-type GSK3β-pseudosubstrate complex and three mutated systems (R4A, R6A, and S9A). Analyses of MD simulations and binding free energies revealed that the phosphorylation of Ser9 is the prerequisite for the autoinhibition of GSK3β, and both mutations of Arg4 and Arg6 to alanine markedly reduced the binding affinities of the mutated pseudosubstrate to the GSK3β catalytic domain, thereby disrupting the autoinhibition of the kinase. This study highlights the importance of Ser9, Arg6, and Arg4 in modulating the autoinhibition mechanism of GSK3β, contributing to a deeper understanding of GSK3β biology.Communicated by Ramaswamy H. Sarma.

Importance of protein flexibility in ranking ERK2 Type I1/2 inhibitor affinities: a computational study

Extracellular-regulated kinase (ERK2) has been regarded as an essential target for various cancers, especially melanoma. Recently, pyrrolidine piperidine derivatives were reported as Type I1/2 inhibitors of ERK2, which occupy both the ATP binding pocket and the allosteric pocket. Due to the dynamic behavior of ERK2 upon the binding of Type I1/2 inhibitors, it is difficult to predict the binding structures and relative binding potencies of these inhibitors with ERK2 accurately. In this work, the binding mechanism of pyrrolidine piperidines was discussed by using different simulation techniques, including molecular docking, ensemble docking based on multiple receptor conformation, molecular dynamics simulations and free energy calculations. Our computational results show that the traditional docking method cannot predict the relative binding ability of the studied inhibitors with high accuracy, but incorporating ERK2 protein flexibility into docking is an effective method to improve the prediction accuracy. It is worth noting that the binding free energies predicted by MM/GBSA or MM/PBSA based on the MD simulations for the docked poses have the highest correlation with the experimental data, which highlights the importance of protein flexibility for accurately predicting the binding ability of Type I1/2 inhibitors of ERK2. In addition, the comprehensive analysis of several representative inhibitors indicates that hydrogen bonds and hydrophobic interactions are of significance for improving the binding affinities of the inhibitors. We hope this work will provide valuable information for further design of novel and efficient Type I1/2 ERK2 inhibitors.

Importance of Incorporating Protein Flexibility in Molecule Modeling: A Theoretical Study on Type I1/2 NIK Inhibitors

NF-κB inducing kinase (NIK), which is considered as the central component of the non-canonical NF-κB pathway, has been proved to be an important target for the regulation of the immune system. In the past few years, NIK inhibitors with various scaffolds have been successively reported, among which type I^1/2 inhibitors that can not only bind in the ATP-binding pocket at the DFG-in state but also extend into an additional back pocket, make up the largest proportion of the NIK inhibitors, and are worthy of more attention. In this study, an integration protocol that combines molecule docking, MD simulations, ensemble docking, MM/GB(PB)SA binding free energy calculations, and decomposition was employed to understand the binding mechanism of 21 tricyclic type I^1/2 NIK inhibitors. It is found that the docking accuracy is largely dependent on the selection of docking protocols as well as the crystal structures. The predictions given by the ensemble docking based on multiple receptor conformations (MRCs) and the MM/GB(PB)SA calculations based on MD simulations showed higher linear correlations with the experimental data than those given by conventional rigid receptor docking (RRD) methods (Glide, GOLD, and Autodock Vina), highlighting the importance of incorporating protein flexibility in predicting protein–ligand interactions. Further analysis based on MM/GBSA demonstrates that the hydrophobic interactions play the most essential role in the ligand binding to NIK, and the polar interactions also make an important contribution to the NIK-ligand recognition. A deeper comparison of several pairs of representative derivatives reveals that the hydrophobic interactions are vitally important in the structural optimization of analogs as well. Besides, the H-bond interactions with some key residues and the large desolvation effect in the back pocket devote to the affinity distinction. It is expected that our study could provide valuable insights into the design of novel and potent type I^1/2 NIK inhibitors.

Insights into the structure-affinity relationships and solvation effects between OfHex1 and inhibitors using molecular dynamics simulations

OfHex1 is a potential target for the rational design of pesticides. TMG-chitotriomycin is one of the most highly specific known inhibitors of chitinolytic β-GlcNAcases from bacteria, fungi and insects. TMG-chitotriomycin and its analogues show different activities to OfHex1, dependent on the number of GlcNAc units. Subsequently, it is essential to explore how these GlcNAc unit number changes cause alterations in activity. In this study, we examined the free energy patterns and per residue decomposition of binding within the complexes of OfHex1 and a series of inhibitors, utilizing restricted molecular dynamics (MD) and water-mediated MM/GBSA calculations. The results indicated Glu328 could form a stronger polar interaction with OfHex1 inhibitors, while Trp448 and Trp490 had important non-polar contributions. Interestingly, the conformation of Trp448 was different in the open or closed state, when OfHex1 bound different inhibitors. Moreover, the water molecule that mediates the GlcNAc Ⅱ and Trp490 may be critical to stabilizing the hydrophobic interaction. Further study showed that isomerization of TMG-chitotriomycin analogs did not decrease binding affinity, however, there was a highly positive correlation between the calculated binding affinities and the experimental activity data (r² = 0.92) when water molecules were explicitly taken into account. Moreover, the water molecules that mediated GlcNAc II and Trp490 might be critical to the stabilization of the hydrophobic interaction and cause the activity difference between TMG-(GlcNAc)₂ and TMG-(GlcNAc).

Insight on mutation-induced resistance to anaplastic lymphoma kinase inhibitor ceritinib from molecular dynamics simulations

Ceritinib, an advanced anaplastic lymphoma kinase (ALK) next-generation inhibitor, has been proved excellent antitumor activity in the treatment of ALK-associated cancers. However, the accumulation of acquired resistance mutations compromise the therapeutic efficacy of ceritinib. Despite abundant mutagenesis data, the structural determinants for reduced ceritinib binding in mutants remains elusive. Focusing on the G1123S and F1174C mutations, we applied molecular dynamics (MD) simulations to study possible reasons for drug resistance caused by these mutations. The MD simulations predict that the studied mutations allosterically impact the configurations of the ATP-binding pocket. An important hydrophobic cluster is identified that connects P-loop and the αC-helix, which has effects on stabilizing the conformation of ATP-binding pocket. It is suggested, in this study, that the G1123S and F1174C mutations can induce the conformational change of P-loop thereby causing the reduced ceritinib affinity and causing drug resistance.

Reliability of Docking-Based Virtual Screening for GPCR Ligands with Homology Modeled Structures: A Case Study of the Angiotensin II Type I Receptor

The number of solved G-protein-coupled receptor (GPCR) crystal structures has expanded rapidly, but most GPCR structures remain unsolved. Therefore, computational techniques, such as homology modeling, have been widely used to produce the theoretical structures of various GPCRs for structure-based drug design (SBDD). Due to the low sequence similarity shared by the transmembrane domains of GPCRs, accurate prediction of GPCR structures by homology modeling is quite challenging. In this study, angiotensin II type I receptor (AT1R) was taken as a typical case to assess the reliability of class A GPCR homology models for SBDD. Four homology models of angiotensin II type I receptor (AT1R) at the inactive state were built based on the crystal structures of CXCR4 chemokine receptor, CCR5 chemokine receptor, and δ-opioid receptor, and refined through molecular dynamics (MD) simulations and induced-fit docking, to allow for backbone and side-chain flexibility. Then, the quality of the homology models was assessed relative to the crystal structures in terms of two criteria commonly used in SBDD: prediction accuracy of ligand-binding poses and screening power of docking-based virtual screening. It was found that the crystal structures outperformed the homology models prior to any refinement in both assessments. MD simulations could generally improve the docking results for both the crystal structures and homology models. Moreover, the optimized homology model refined by MD simulations and induced-fit docking even shows a similar performance of the docking assessment to the crystal structures. Our results indicate that it is possible to establish a reliable class A GPCR homology model for SBDD through the refinement by integrating multiple molecular modeling techniques.

Comprehensive assessment of nine docking programs on type II kinase inhibitors: prediction accuracy of sampling power, scoring power and screening power

Protein kinases have been regarded as important therapeutic targets for many diseases. Currently, a total of 41 kinase inhibitors have been approved by the Food and Drug Administration, along with a large number of kinase inhibitors being evaluated in clinical and preclinical trials. Among all, allosteric inhibitors, such as type II kinase inhibitors, have attracted extensive attention owing to their potential high selectivity. Nowadays, molecular docking has become a powerful tool to search for novel kinase inhibitors. However, as for type II kinase inhibitors, their allosteric characteristics may exert a deep influence on docking accuracy. In this study, a comprehensive assessment was conducted to evaluate the effectiveness of nine docking algorithms towards type II kinase inhibitors. The calculation results showed that most tested docking programs, especially Glide with XP scoring, LeDock and Surflex-Dock, succeeded in the accurate identification of near-native binding poses, with the success rates ranging from 0.80 to 0.90, and the scoring functions in GOLD and LeDock outperformed the others in the prediction of relative binding affinities. In terms of the P-values, areas under the curve and enrichment factors, Glide with XP scoring, Surflex-Dock, GOLD with Astex Statistical Potential scoring and LeDock had better screening power to discriminate between active compounds and decoys. However, the screening power is sensitive to different initial conformations of the same target. It is expected that our study can provide some guidance for docking-based virtual screening to discover novel type II kinase inhibitors, as well as other allosteric inhibitors.

Resistance mechanism of the oncogenic β3-αC deletion mutation in BRAF kinase to dabrafenib and vemurafenib revealed by molecular dynamics simulations and binding free energy calculations

BRAF kinase is an essential target for anti-cancer drug development. Emergence of the β3-αC loop deletion mutation (ΔNVTAP) in BRAF kinase frequently occurred in human cancers seriously compromises the therapeutic efficacy of some BRAF kinase inhibitors, such as dabrafenib and vemurafenib. However, the mechanism of this resistance is still not well understood. In this study, the influence of the β3-αC deletion mutation on the binding profiles of three BRAF kinase inhibitors (AZ628, dabrafenib, and vemurafenib) with BRAF^V600E or BRAF^ΔNVTAP was explored by conventional molecular dynamics (MD) simulations and binding free energy calculations. The simulation results indicated that the β3-αC deletion mutation enhances the flexibility of the αC helix and alters their conformations, which amplify the conformational entropy change (-TΔS) and weaken the interactions between the inhibitors and BRAF. The further per-residue binding free energy decomposition analysis revealed that the ΔNVTAP mutation changed the contributions of a few key residues to the bindings of dabrafenib or vemurafenib, such as L57, L66, W83, C84, F135, G145, and F147, but did not have obvious impact on the contributions of these residues to AZ628. Our results provide valuable clues to understand the mechanisms of drug resistance conferred by the β3-αC deletion mutation.

Exploring the binding mechanisms of PDE5 with chromeno[2,3-c]pyrrol-9(2H)-one by theoretical approaches

Cyclic nucleotide phosphodiesterase type 5 (PDE5), exclusively specific for the cyclic guanosine monophosphate (cGMP), is an important drug target for the treatment of erectile dysfunction and pulmonary arterial hypertension (PAH). Although many PDE5 inhibitors have been approved, such as sildenafil, vardenafil, tadalafil and so on, extensive studies have reported some side effects, such as vision disturbance and hearing loss as a result of the amino acid sequence and the secondary structural similarity of other PDEs to the catalytic domain of PDE5. In this study, multiple docking strategies, molecular dynamics (MD) simulations, free energy calculations and decomposition were employed to explore the structural determinants of PDE5 with a series of chromeno[2,3-c]pyrrol-9(2H)-one derivatives. First, reliable docking results were obtained using quantum mechanics/molecular mechanics (QM/MM) docking. Then, MD simulations and MM/GBSA free energy calculations were used to explore the dynamic binding process and characterize the binding modes of the inhibitors with different activities. The predicted binding free energies are in good agreement with the experimental data, and the MM/GBSA free energy decomposition analysis sheds light on the importance of hydrogen bonds with Gln817, π-π stacks against Phe820 and hydrophobic residues for the PDE5 binding of the studied inhibitors. The structural and energetic insights obtained here are useful for understanding the molecular mechanism of ligand binding and designing novel potent and selective PDE5 inhibitors with new scaffolds.

Conformational Transition of Key Structural Features Involved in Activation of ALK Induced by Two Neuroblastoma Mutations and ATP Binding: Insight from Accelerated Molecular Dynamics Simulations

Deregulated kinase activity of anaplastic lymphoma kinase (ALK) has been observed to be implicated in the development of tumor progression. The activation mechanism of ALK is proposed to be similar to other receptor tyrosine kinases (RTKs), but the distinct static X-ray crystal conformation of ALK suggests its unique conformational transition. Herein, we have illustrated the dynamic conformational property of wild-type ALK as well as the kinase activation equilibrium variation induced by two neuroblastoma mutations (R1275Q and Y1278S) and ATP binding by performing enhanced sampling accelerated Molecular Dynamics (aMD) simulations. The results suggest that the wild-type ALK is mostly favored in the inactive state, whereas the mutations and ATP binding promote a clear shift toward the active-like conformation. The R1275Q mutant stabilizes the active conformation by rigidifying the αC-in conformation. The Y1278S mutant promotes activation at the expense of a π-stacking hydrophobic cluster, which plays a critical role in the stabilization of the inactive conformation of native ALK. ATP produces a more compact active site and thereby facilitates the activation of ALK. Taken together, these findings not only elucidate the diverse conformations in different ALKs but can also shed light on new strategies for protein engineering and structural-based drug design for ALK.

Insights into Resistance Mechanisms of Inhibitors to Mps1 C604Y Mutation via a Comprehensive Molecular Modeling Study

Mono-polar spindle 1 (Mps1/TTK) represents a protein kinase reported to be vital for cell division processes and is generally regarded as an attractive target for the treatment of hepatocellular carcinoma, breast carcinoma, and colon cancer. However, the C604Y mutation has been linked to acquired resistance. Recently, three potential small-molecule inhibitors of Mps1 (i.e., reversine, NMS-P715, and its derivative Cpd-5) were reported for the C604Y mutation that exhibit significant resistance to NMS-P715 and Cpd-5, but retain affinity for reversine. In this study, classical molecular dynamic (MD) simulations, accelerated MD (aMD) simulations, and umbrella sampling (US) simulations were performed to illustrate the resistance mechanisms of inhibitors to Mps1. The classical MD simulations combined with free energy calculations revealed that reversine features similar binding affinity characteristics to both Mps1^WT and Mps1^C604Y, but both NMS-P715 and Cpd-5 feature much higher binding affinities to Mps1^WT than to Mps1^C604Y. The major variations were shown to be controlled by electrostatic energy and the conformational change of A-loop-induced entropy increased. The large conformational changes of Mps1^C604Y bound to NMS-P715 and Cpd-5 were also observed in aMD simulations. The US simulation results further suggest that reversine and Cpd-5 both exhibit similar dissociation processes from both Mps1^WT and Mps1^C604Y, but Cpd-5 and NMS-P715 were found to dissociate more easily from Mps1^C604Y than from Mps1^WT, thus a reduced residence time was responsible for the inhibitors resistance to the C604Y mutation. The physical principles provided by the present study may provide important clues for the discovery and rational design of novel inhibitors to combat the C604Y mutation of Mps1.

Insight into tetrodotoxin blockade and resistance mechanisms of Nav 1.2 sodium channel by theoretical approaches

Na_v 1.2, a member of voltage-gated sodium channels (Na_v s) that are responsible for the generation and propagation of action potentials along the cell membrane, and play a vital role in the process of information transmission within the nervous system and muscle contraction, is preferentially expressed in the central nervous system. As a potent and selective blocker of Na_v s, tetrodotoxin (TTX) has been extensively studied in biological and chemical sciences, whereas the detailed mechanism by which it blocks nine Na_v 1 channel subtypes remain elusive. Despite the high structural similarity, the TTX metabolite 4,9-anhydro-TTX is 161 times less effective toward the mammalian Na_v 1.2, which puzzled us to ask a question why such a subtle structural variation results in the largely binding affinity difference. In the current work, an integrated computational strategy, including homology modeling, induced fit docking, explicit-solvent MD simulations, and free energy calculations, was employed to investigate the binding mechanism and conformational determinants of TTX analogs. Based on the computational results, the H-bond interactions between C4-OH and C9-OH of TTX and the outer ring carboxylates of the selectivity-filter residues, and the cation-π interaction between the primary amine of guanidinium of TTX and Phe385 determine the difference of their binding affinities. Moreover, the computationally simulations were carried out for the D384N and E945K mutants of hNa_v 1.2-TTX, and the rank of the predicted binding free energies is in accordance with the experimental data. These observations provide a valuable model to design potent and selective neurotoxins of Na_v 1.2 and shed light on the blocking mechanism of TTX to sodium channels.

Importance of protein flexibility on molecular recognition: modeling binding mechanisms of aminopyrazine inhibitors to Nek2

NIMA-related kinase 2 (Nek2) plays a significant role in cell cycle regulation, and overexpression of Nek2 has been observed in several types of carcinoma, suggesting it is a potential target for cancer therapy. In this study, we attempted to gain more insight into the binding mechanisms of a series of aminopyrazine inhibitors of Nek2 through multiple molecular modeling techniques, including molecular docking, molecular dynamics (MD) simulations and free energy calculations. The simulation results showed that the induced fit docking and ensemble docking based on multiple protein structures yield better predictions than conventional rigid receptor docking, highlighting the importance of incorporating receptor flexibility into the accurate predictions of the binding poses and binding affinities of Nek2 inhibitors. Additionally, we observed that the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculations did not show better performance than the docking scoring to rank the binding affinities of the studied inhibitors, suggesting that MM/GBSA is system-dependent and may not be the best choice for the Nek2 systems. Moreover, the detailed information on protein-ligand binding was characterized by the MM/GBSA free energy decomposition, and a number of derivatives with improved docking scores were designed. It is expected that our study can provide valuable information for the future rational design of novel and potent inhibitors of Nek2.

Understanding the resistance mechanism of penicillin binding protein 1a mutant against cefotaxime using molecular dynamic simulation

Antibiotic resistance is a threatening challenge for global health, as the expansion of resistance to current antibiotics has made serious therapeutic problems. Genome mutations are key evolutionary mechanisms conferring antibiotic resistance in bacterial pathogens. For example, penicillin and cephalosporins resistance is mostly mediated by mutations in penicillin binding proteins to change the affinity of the drug. Accordingly, threonine point mutations were reported to develop antibiotic resistance in various bacterial infections including pneumococcal infections. In this study, conventional molecular dynamics simulations, umbrella sampling simulations and MM/GBSA free energy calculations were applied to figure out how the Threonine to Alanine mutation (T to A) at STMK motif affects the binding of cefotaxime to Penicillin Binding Protein 1a and to reveal the resistance mechanism induced by the T to A mutation. The results obtained from the computational methods demonstrate that the T to A mutation increases the flexibility of the binding pocket and changes its conformation, which leads to increased conformational entropy change (-TΔS) and attenuates the bonds between the ligand and the receptor. In brief, our findings indicate that both of the alterations of the conformational enthalpy and entropy contribute to the T to A-induced resistance in the binding of cefotaxime into penicillin binding protein 1a.

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.

Combating Drug-Resistant Mutants of Anaplastic Lymphoma Kinase with Potent and Selective Type-I1/2 Inhibitors by Stabilizing Unique DFG-Shifted Loop Conformation

Targeted inhibition of anaplastic lymphoma kinase (ALK) dramatically improved therapeutic outcomes in the treatment of ALK-positive cancers, but unfortunately patients invariably progressed due to acquired resistance mutations in ALK. Currently available drugs are all type-I inhibitors bound to the ATP-binding pocket and are most likely to be resistant in patients harboring genetic mutations surrounding the ATP pocket. To overcome drug resistance, we rationally designed a novel kind of "bridge" inhibitor, which specially bind into an extended hydrophobic back pocket adjacent to the ATP-binding site of ALK. The novel type-I1/2 inhibitors display excellent antiproliferation activity against ALK-positive cancer cells and appear superior to two clinically used drugs, crizotinib and ceritinib. Structural and molecular modeling analyses indicate that the inhibitor induces dramatic conformational transition and stabilizes unique DFG-shifted loop conformation, enabling persistent sensitivity to different genetic mutations in ALK. These data highlight a rationale for further development of next-generation ALK inhibitors to combat drug resistance.

How Does the L884P Mutation Confer Resistance to Type-II Inhibitors of JAK2 Kinase: A Comprehensive Molecular Modeling Study

Janus kinase 2 (JAK2) has been regarded as an essential target for the treatment of myeloproliferative neoplasms (MPNs). BBT594 and CHZ868, Type-II inhibitors of JAK2, illustrate satisfactory efficacy in preclinical MPNs and acute lymphoblastic leukemia (ALL) models. However, the L884P mutation of JAK2 abrogates the suppressive effects of BBT594 and CHZ868. In this study, conventional molecular dynamics (MD) simulations, umbrella sampling (US) simulations and MM/GBSA free energy calculations were employed to explore how the L884P mutation affects the binding of BBT594 and CHZ868 to JAK2 and uncover the resistance mechanism induced by the L884P mutation. The results provided by the US and MD simulations illustrate that the L884P mutation enhances the flexibility of the allosteric pocket and alters their conformations, which amplify the conformational entropy change (−TΔS) and weaken the interactions between the inhibitors and target. Additionally, the structural analyses of BBT594 and CHZ868 in complex with the WT JAK2 illustrate that the drug tail with strong electronegativity and small size located in the allosteric pocket of JAK2 may enhance anti-resistance capability. In summary, our results highlight that both of the changes of the conformational entropies and enthalpies contribute to the L884P-induced resistance in the binding of two Type-II inhibitors into JAK2 kinase.

Understanding the microscopic binding mechanism of hydroxylated and sulfated polybrominated diphenyl ethers with transthyretin by molecular docking, molecular dynamics simulations and binding free energy calculations

Polybrominated diphenyl ethers (PBDEs), one typical type of persistent environmental contaminant, have toxicological effects such as disrupting thyroid homeostasis in the human body. The high binding affinities of hydroxylated metabolites of PBDEs (OH-PBDEs) with transthyretin (TTR) were considered to be one major reason for their extraordinary capacity of passing through the blood-brain barrier via competitive thyroid hormone (T4) transport protein binding. Recent findings showed that sulfated PBDEs can be formed in human liver cytosol as phase-II metabolites. However, experimentally determined data for the TTR binding potential of the sulfated PBDEs are still not available. Therefore, molecular docking and molecular dynamics (MD) simulations were employed in the present study to probe the molecular basis of TTR interacting with hydroxylated and sulfated PBDEs at the atomic level. The docking scores of LeDock were used to construct the structure-based predictive model. The calculated results showed that the sulfated PBDEs have stronger affinity for TTR than the corresponding OH-PBDEs. Further analysis of structural characteristics based on MD simulations indicated that upon the binding of PBDE metabolites, the stability of TTR was enhanced and the dissociation rate of the tetrameric protein structure was potentially decreased. Subsequent binding free energy calculations implied that van der Waals interactions are the dominant forces for the binding of these metabolites of PBDEs at the T4 site of TTR. The residues Ser117/Ser117' and Lys15/Lys15' were identified, by both residue energy decomposition and computational alanine-scanning mutagenesis methods, as key residues which play an important role in determining the binding orientations of the -OSO₃^- group of sulfated PBDEs by formation of either hydrogen bonds or electrostatic interactions, respectively. In general, the combination of docking calculations with MD simulations provided a theoretically toxicological assessment for the metabolites of PBDEs, deep insight into the recognition mechanism of TTR for these compounds, and thus more comprehensive understanding of the thyroid-related toxic effects of PBDEs as well.

Theoretical studies on the selective mechanisms of GSK3β and CDK2 by molecular dynamics simulations and free energy calculations

Glycogen synthase kinase 3 (GSK3) is a serine/threonine protein kinase which is widely involved in cell signaling and controls a broad number of cellular functions. GSK3 contains α and β isoforms, and GSK3β has received more attention and becomes an attractive drug target for the treatment of several diseases. The binding pocket of cyclin-dependent kinase 2 (CDK2) shares high sequence identity to that of GSK3β, and therefore, the design of highly selective inhibitors toward GSK3β remains a big challenge. In this study, a computational strategy, which combines molecular docking, molecular dynamics simulations, free energy calculations, and umbrella sampling simulations, was employed to explore the binding mechanisms of two selective inhibitors to GSK3β and CDK2. The simulation results highlighted the key residues critical for GSK3β selectivity. It was observed that although GSK3β and CDK2 share the conserved ATP-binding pockets, some different residues have significant contributions to protein selectivity. This study provides valuable information for understanding the GSK3β-selective binding mechanisms and the rational design of selective GSK3β inhibitors.

Constructing and Validating High-Performance MIEC-SVM Models in Virtual Screening for Kinases: A Better Way for Actives Discovery

The MIEC-SVM approach, which combines molecular interaction energy components (MIEC) derived from free energy decomposition and support vector machine (SVM), has been found effective in capturing the energetic patterns of protein-peptide recognition. However, the performance of this approach in identifying small molecule inhibitors of drug targets has not been well assessed and validated by experiments. Thereafter, by combining different model construction protocols, the issues related to developing best MIEC-SVM models were firstly discussed upon three kinase targets (ABL, ALK, and BRAF). As for the investigated targets, the optimized MIEC-SVM models performed much better than the models based on the default SVM parameters and Autodock for the tested datasets. Then, the proposed strategy was utilized to screen the Specs database for discovering potential inhibitors of the ALK kinase. The experimental results showed that the optimized MIEC-SVM model, which identified 7 actives with IC50 < 10 μM from 50 purchased compounds (namely hit rate of 14%, and 4 in nM level) and performed much better than Autodock (3 actives with IC50 < 10 μM from 50 purchased compounds, namely hit rate of 6%, and 2 in nM level), suggesting that the proposed strategy is a powerful tool in structure-based virtual screening.

Flexible receptor docking for drug discovery

INTRODUCTION

Molecular docking has become a popular method for virtual screening. Docking small molecules to a rigid biological receptor is fast but could produce many false negatives and identify less diverse compounds. Flexible receptor docking has alleviated this problem.

AREAS COVERED

This article focuses on reviewing ensemble docking as an approximate but inexpensive method to incorporate receptor flexibility in molecular docking. It outlines key features and recent advances of this method and points out problem areas that need to be addressed to make it even more useful in drug discovery.

EXPERT OPINION

Among the different methods introduced for flexible receptor docking, ensemble docking represents one of the most popular approaches, especially for high-throughput virtual screening. One can generate structural ensembles by using experimental structures, by structural modeling and by various types of molecular simulations. In building a structural ensemble, a judicious choice of the structures to be included can improve performance. Furthermore, reducing the size of the structural ensemble can cut computational costs, and removing the structures that can bind few ligands well could enrich the number of true actives identified by ensemble docking. The ability of ensemble docking to identify more true positives at the top of a rank-ordered list also depends on the choice of the methods to score and rank compounds, an area that needs further research.

The origin of specificity and insight into recognition between an aminoacyl carrier protein and its partner ligase

Acyl carrier proteins (ACPs) are among the most promiscuous proteins in terms of protein-protein interactions and it is quite puzzling how ACPs select the correct partner between many possible upstream and downstream binding proteins. To address this question, we performed molecular dynamics simulations on dimeric Bradyrhizobium japonicum Gly:CP ligase 1 to inspect the origin of its selectivity towards the three types of carrier proteins, namely holoCP, apoCP, and holoCP-Gly, which only differ in the attached prosthetic group. In line with experiments, MM-GBSA analysis revealed that the ligase preferentially binds the holoCP form to both subunits with the binding free energies of -20.7 and -19.1 kcal mol(-1), while the apoCP form, without the prosthetic group, is also recognized, but the binding values of -9.2 and -3.6 kcal mol(-1) suggest that there is no competition for the ligase binding as long as the holoCP is present. After the prosthetic group becomes glycylated, the holoCP-Gly dissociates from the ligase, as supported by its endergonic binding free energies of 2.9 and 20.9 kcal mol(-1). Our results indicate that these affinity differences are influenced by three aspects: the form of the prosthetic group and the specific non-polar hydrophobic interactions, as well as charge complementarity dominantly manifested through Arg220-Glu53 ion pair within the binding region among proteins. A careful examination of the bonding patterns within the ligase active site elucidated the interactions with Arg258, Asp215 and Tyr132 as being predominant in stabilizing the prosthetic group, which are significantly diminished upon glycation, thus promoting complex dissociation.

Graphene Oxide Selectively Enhances Thermostability of Trypsin

In the past few years, graphene and its derivative, graphene oxide (GO), have been extensively studied for their applications in biotechnology. In our previous work, we reported certain PEGylated GOs (GO-PEGs) can selectively promote trypsin activity and enhance its thermostability. To further explore this, here we synthesized a series of GO-PEGs with varying PEGylation degrees. Enzymatic activity assay shows that both GO and GO-PEGs can protect trypsin, but not chymotrypsin, from thermal denaturation at high temperature. Surprisingly, the lower the PEGylation degree, the better the protection, and GO as well as the GO-PEG with the lowest PEGylation degree show the highest protection efficiency (∼70% retained activity at 70 °C). Fluorescence spectroscopy analysis shows that GO/GO-PEGs have strong interactions with trypsin. Molecular Dynamics (MD) simulation results reveal that trypsin is adsorbed onto the surface of GO through its cationic residues and hydrophilic residues. Different from chymotrypsin adsorbed on GO, the active site of trypsin is covered by GO. MD simulation at high temperature shows that, through such interaction with GO, trypsin's active site is therefore stabilized and protected by GO. Our work not only illustrates the promising potential of GO/GO-PEGs as efficient, selective modulators for trypsin, but also provides the interaction mechanism of GO with specific proteins at the nano-bio interface.