Molecular Dynamics Simulations of Wild Type and Mutants of SAPAP in Complexed with Shank3

Screening of Phytocompounds for Identification of Prospective Histone Deacetylase 1 (HDAC1) Inhibitor: An In Silico Molecular Docking, Molecular Dynamics Simulation, and MM-GBSA Approach

The upregulation of HDAC1 facilitate the induction of epigenetic repression of genes responsible for suppressing tumourigenesis, thereby triggering the development of cancer. HDAC1 inhibitors have thus emerged as possible therapeutic approaches against a variety of human malignancies, as they can inhibit the activity of certain HDACs, repair the overexpression of tumour suppressor genes, and induce cell differentiation, cell cycle arrest, and apoptosis. In this study, among 810 virtually screened compounds, Pinocembrin (PHUB000396) had a significant binding affinity (-7.99 kcal/mol). In molecular dynamics simulation (MD) studies for 200 ns time scale, the compound Pinocembrin effectively undergoes conformational optimization, thereby enabling its accommodation within the active site of the receptor. This outcome serves as a rational for the observed binding affinity. The optimal binding free energy calculations using the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) (-35.86 ± 7.52 kcal/mol) showed the significant role of van der Waals forces and Coulomb interactions in the stability of the respective complex. The pharmacokinetic study showed its potential as a lead compound. The in-silico cytotoxicity prediction also confirmed its potential as an active anticancer phytocompound in lung and brain cancer. Therefore, it can be predicted that Pinocembrin could be a useful bioactive compound as an HDAC1 inhibitor and could be used in developing epigenetic therapy in cancer such as brain cancer and lung cancer to regulate gene expression.

The road to evolution of ProTx2: how to be a subtype-specific inhibition of human Nav1.7

The human voltage-gated sodium channel Nav1.7 is a widely proven target for analgesic drug studies. ProTx2, a 30-residue polypeptide from Peruvian green tarantula venom, shows high specificity to activity against human Nav1.7, suggesting its potential to become a non-addictive analgesic. However, its high sensitivity to human Nav1.4 raises concerns about muscle side effects. Here, we engineered three mutants (R13A, R13D, and K27Y) of ProTx2 to evaluate their pharmacological activities toward Nav1.7 and Nav1.4. It is demonstrated that the mutant R13D maintained the analgesic effect in mice while dramatically reducing its muscle toxicity compared with ProTx2. The main reason is the formation of a strong electrostatic interaction between R13D and the negatively charged amino acid residues in DII/S3-S4 of Nav1.7, which is absent in Nav1.4. This study advances our understanding and insights on peptide toxins, paving the way for safer, effective non-addictive analgesic development.

Identification of potential inhibitor against Leishmania donovani mitochondrial DNA primase through in-silico and in vitro drug repurposing approaches

Leishmania donovani is the causal organism of leishmaniasis with critical health implications affecting about 12 million people around the globe. Due to less efficacy, adverse side effects, and resistance, the available therapeutic molecules fail to control leishmaniasis. The mitochondrial primase of Leishmania donovani (LdmtPRI1) is a vital cog in the DNA replication mechanism, as the enzyme initiates the replication of the mitochondrial genome of Leishmania donovani. Hence, we target this protein as a probable drug target against leishmaniasis. The de-novo approach enabled computational prediction of the three-dimensional structure of LdmtPRI1, and its active sites were identified. Ligands from commercially available drug compounds were selected and docked against LdmtPRI1. The compounds were chosen for pharmacokinetic study and molecular dynamics simulation based on their binding energies and protein interactions. The LdmtPRI1 gene was cloned, overexpressed, and purified, and a primase activity assay was performed. The selected compounds were verified experimentally by the parasite and primase inhibition assay. Capecitabine was observed to be effective against the promastigote form of Leishmania donovani, as well as inhibiting primase activity. This study's findings suggest capecitabine might be a potential anti-leishmanial drug candidate after adequate further studies.

Netting into the Sophoretin pool: An approach to trace GSTP1 inhibitors for reversing chemoresistance

Chemoresistance, a significant challenge in cancer treatment, is often associated with the cellular glutathione-related detoxification system. The GSTP1 isoenzyme (glutathione S-transferases) plays a critical role in the cytoplasmic inactivation of anticancer drugs. This suggests the identification of GSTP1 inhibitors to combat chemoresistance. We screened Sophoretin (also called quercetin) derivatives for molecular properties, pharmacokinetics, and toxicity profiles. Following that, we conducted molecular docking and simulations between selected derivatives and GSTP1. The best-docked complex, GSTP1-quercetin 7-O-β-D-glucoside, exhibited a binding affinity of -8.1 kcal/mol, with no predicted toxicity and good pharmacokinetic properties. Molecular dynamics simulations confirmed the stability of this complex. Quercetin 7-O-β-D-glucoside shows promise as a lead candidate for addressing chemoresistance in cancer patients, although further experimental studies are needed to validate its efficacy and therapeutic potential.

Integration of network pharmacology, molecular docking, and simulations to evaluate phytochemicals from Drymaria cordata against cervical cancer

Introduction: Cervical cancer is prevalent among women worldwide. It is a type of cancer that occurs in the cells of the cervix, the lower part of the uterus. Mostly, it is observed in developing nations due to limited access to screening tools. Natural products with anticancer properties and fewer side effects have gained attention. Therefore, this study evaluates the potential of Drymaria cordata as a natural source for treating cervical cancer. Methodology: Phytocompounds present in Drymaria cordata were screened for their molecular properties and drug-likeness. The selected compounds were studied using systems biology tools such as network pharmacology, molecular docking, and molecular dynamics simulations, including MMGBSA studies. Results: Through network pharmacology, molecular docking, and molecular dynamics simulations, quercetin 3-O-β-d-glucopyranosyl-(1→2)-rhamnopyranoside was identified as a hit compound targeting HRAS and VEGFA proteins. These proteins were found to be responsible for the maximum number of pathway modulations in cervical cancer. Conclusion: Drymaria cordata exhibits potential for treating cervical cancer due to the presence of quercetin 3-O-β-d-glucopyranosyl-(1→2)-rhamnopyranoside. Further validation of these findings through in vitro and in vivo studies is required.

Assessing the potential of Psidium guajava derived phytoconstituents as anticholinesterase inhibitor to combat Alzheimer's disease: an in-silico and in-vitro approach

Acetylcholinesterase (AChE) inhibitors play a crucial role in the treatment of Alzheimer's disease. These drugs increase acetylcholine levels by inhibiting the enzyme responsible for its degradation, which is a vital neurotransmitter involved in memory and cognition. This intervention intermittently improves cognitive symptoms and augments neurotransmission. This study investigates the potential of Psidium guajava fruit extract as an acetylcholinesterase (AChE) inhibitor for Alzheimer's disease treatment. Molecular characteristics and drug-likeness were analyzed after HR-LCMS revealed phytocompounds in an ethanolic extract of Psidium guajava fruit. Selected phytocompounds were subjected to molecular docking against AChE, with the best-docked compound then undergoing MD simulation, MMGBSA, DCCM, FEL, and PCA investigations to evaluate the complex stability. The hit compound's potential toxicity and further pharmacokinetic features were also predicted. Anticholinesterase activity was also studied using in vitro assay. The HR-LCMS uncovered 68 compounds. Based on computational analysis, Fluspirilene was determined to have the highest potential to inhibit AChE. It was discovered that the Fluspirilene-AChE complex is stable and that Fluspirilene has a high binding affinity for AChE. Extract of Psidium guajava fruit significantly inhibits AChE (88.37% at 200 μg/ml). It is comparable to the standard AChE inhibitor Galantamine. Fluspirilene exhibited remarkable binding to AChE. Psidium guajava fruit extract demonstrated substantial AChE inhibitory activity, indicating its potential for Alzheimer's treatment. The study underscores natural sources' significance in drug discovery.Communicated by Ramaswamy H. Sarma.

Rhamnetin, a nutraceutical flavonoid arrests cell cycle progression of human ovarian cancer (SKOV3) cells by inhibiting the histone deacetylase 2 protein

Overexpression of HDAC 2 promotes cell proliferation in ovarian cancer. HDAC 2 is involved in chromatin remodeling, transcriptional repression, and the formation of condensed chromatin structures. Targeting HDAC 2 presents a promising therapeutic approach for correcting cancer-associated epigenetic abnormalities. Consequently, HDAC 2 inhibitors have evolved as an attractive class of anti-cancer agents. This work intended to investigate the anti-cancer abilities and underlying molecular mechanisms of Rhamnetin in human epithelial ovarian carcinoma cells (SKOV3), which remain largely unexplored. We employed various in vitro methods, including MTT, apoptosis study, cell cycle analysis, fluorescence microscopy imaging, and in vitro enzymatic HDAC 2 protein inhibition, to examine the chemotherapeutic sensitivity of Rhamnetin in SKOV3 cells. Additionally, we conducted in silico studies using molecular docking, MD simulation, MM-GBSA, DFT, and pharmacokinetic analysis to investigate the binding interaction mechanism within Rhamnetin and HDAC 2, alongside the compound's prospective as a lead candidate. The in vitro assay confirmed the cytotoxic effects of Rhamnetin on SKOV3 cells, through its inhibition of HDAC 2 activity. Rhamnetin, a nutraceutical flavonoid, halted at the G1 phase of the cell cycle and triggered apoptosis in SKOV3 cells. Furthermore, computational studies provided additional evidence of its stable binding to the HDAC 2 protein's binding site cavity. Based on our findings, we conclude that Rhamnetin effectively promotes apoptosis and mitigates the proliferation of SKOV3 cells through HDAC 2 inhibition. These results highlight Rhamnetin as a potential lead compound, opening a new therapeutic strategy for human epithelial ovarian cancer.Communicated by Ramaswamy H. Sarma.

New lead compounds identification against KRas mediated cancers through pharmacophore-based virtual screening and in vitro assays

Cancer remains the leading cause of mortality and morbidity in the world, with 19.3 million new diagnoses and 10.1 million deaths in 2020. Cancer is caused due to mutations in proto-oncogenes and tumor-suppressor genes. Genetic analyses found that Ras (Rat sarcoma) is one of the most deregulated oncogenes in human cancers. The Ras oncogene family members including NRas (Neuroblastoma ras viral oncogene homolog), HRas (Harvey rat sarcoma) and KRas are involved in different types of human cancers. The mutant KRas is considered as the most frequent oncogene implicated in the development of lung, pancreatic and colon cancers. However, there is no efficient clinical drug even though it has been identified as an oncogene for 30 years. Therefore there is an emerging need to develop potent, new anticancer drugs. In this study, computer-aided drug designing approaches as well as experimental methods were employed to find new and potential anti-cancer drugs. The pharmacophore model was developed from an already known FDA approved anti-cancer drug Bortezomib using the software MOE. The validated pharmacophore model was then used to screen the in-house and commercially available databases. The pharmacophore-based virtual screening resulted in 26 and 86 hits from in-house and commercial databases respectively. Finally, 6/13 (in-house database) and 24/64 hits (commercial databases) were selected with different scaffolds having good interactions with the significant active residues of KRasG12D protein that were predicted as potent lead compounds. Finally, the results of pharmacophore-based virtual screening were further validated by molecular dynamics simulation analysis. The 6 hits of the in-house database were further evaluated experimentally. The experimental results showed that these compounds have good anti-cancer activity which validate the protocol of our in silico studies. KRasG12D protein is a very important anti-cancer target and potent inhibitors for this target are still not available, so small lead compound inhibitors were identified to inhibit the activity of this protein by blocking the GTP-binding pocket.Communicated by Ramaswamy H. Sarma.

Mutational analysis of catalytic site domain of CCHFV L RNA segment

INTRODUCTION

Crimean-Congo haemorrhagic fever virus (CCHFV) has tripartite RNA genome and is endemic in various countries of Asia, Africa and Europe.

METHOD

The present study is focused on mutation profiling of CCHFV L segment and phylogenetic clustering of protein dataset into six CCHFV genotypes.

RESULTS

Phylogenetic tree rooted with NCBI reference sequence (YP_325663.1) indicated less divergence from genotype III and the sequences belonging to same genotypes have shown less divergence among each other. Mutation frequency at 729 mutated positions was calculated and 563, 49, 33, 46 and 38 amino acid positions were found to be mutated at mutation frequency intervals of 0-0.2, 0.21-0.4, 0.41-0.6, 0.61-0.8 and 0.81-1.0 respectively. Thirty-eight highly frequent mutations (0.81-1.0 interval) were found in all genotypes and mapping in L segment (encoded for RdRp) revealed four mutations (V2074I, I2134T/A, V2148A and Q2695H/R) in catalytic site domain and no mutation in OTU domain. Molecular dynamic simulation and in silico analysis showed that catalytic site domain displayed large deviation and fluctuation upon introduction of these point mutations.

CONCLUSION

Overall study provides strong evidence that OTU domain is highly conserved and less prone to mutation whereas point mutations recorded in catalytic domain have affected the stability of protein and were found to be persistent in the large population.

Unraveling the Structural Changes in the DNA-Binding Region of Tumor Protein p53 (TP53) upon Hotspot Mutation p53 Arg248 by Comparative Computational Approach

The vital tissue homeostasis regulator p53 forms a tetramer when it binds to DNA and regulates the genes that mediate essential biological processes such as cell-cycle arrest, senescence, DNA repair, and apoptosis. Missense mutations in the core DNA-binding domain (109-292) simultaneously cause the loss of p53 tumor suppressor function and accumulation of the mutant p53 proteins that are carcinogenic. The most common p53 hotspot mutation at codon 248 in the DNA-binding region, where arginine (R) is substituted by tryptophan (W), glycine (G), leucine (L), proline (P), and glutamine (Q), is reported in various cancers. However, it is unclear how the p53 Arg248 mutation with distinct amino acid substitution affects the structure, function, and DNA binding affinity. Here, we characterized the pathogenicity and protein stability of p53 hotspot mutations at codon 248 using computational tools PredictSNP, Align GVGD, HOPE, ConSurf, and iStable. We found R248W, R248G, and R248P mutations highly deleterious and destabilizing. Further, we subjected all five R248 mutant-p53-DNA and wt-p53-DNA complexes to molecular dynamics simulation to investigate the structural stability and DNA binding affinity. From the MD simulation analysis, we observed increased RMSD, RMSF, and Rg values and decreased protein-DNA intermolecular hydrogen bonds in the R248-p53-DNA than the wt-p53-DNA complexes. Likewise, due to high SASA values, we observed the shrinkage of proteins in R248W, R248G, and R248P mutant-p53-DNA complexes. Compared to other mutant p53-DNA complexes, the R248W, R248G, and R248P mutant-p53-DNA complexes showed more structural alteration. MM-PBSA analysis showed decreased binding energies with DNA in all five R248-p53-DNA mutants than the wt-p53-DNA complexes. Henceforth, we conclude that the amino acid substitution of Arginine with the other five amino acids at codon 248 reduces the p53 protein's affinity for DNA and may disrupt cell division, resulting in a gain of p53 function. The proposed study influences the development of rationally designed molecular-targeted treatments that improve p53-based therapeutic outcomes in cancer.

Molecular Dynamic Simulation of Neurexin1α Mutations Associated with Mental Disorder

Neurexin1 gene is essential for formulating synaptic cell adhesion to establish synapses. In a previous work, 38 SNPs in Neurexin1 recoded in mental disorder patients have been collected. Five computational prediction tools have been used to predict the effect of SNPs on protein function and stability. Only four SNPs in Neurexin1α have deleterious prediction results from at least four tools. The current work aims to use molecular dynamic simulation (MD) to study the effects of the four mutations on Neurexin1α both on the whole protein as well as identifying affected domains by mutations. A protein model that consists of five domains out of six domains in the real protein was used; missing residues were added, and model was tested for quality. The MD experiment has last for 1.5 μs where four parameters have been used for studying the whole protein in addition to three more parameters for the domain analysis. The whole protein study has shown that two mutations E427I for Autism and R525C for non-syndromic intellectual disability (NSID) have distinctive behavior across the four used parameters. Domain study has confirmed the previous results where the five domains of R525C have acted differently from wild type (WT), while E427I has acted differently for four domains from wild type. The other two mutations D104H and G379E have three domains that only acted differently from wild type. The fourth domain of all mutations has an obvious distinctive behavior from wild type. Further study of E427I and R525C mutations can lead to better understanding of autism and NSID.

Abetting host immune response by inhibiting rhipicephalus sanguineus Evasin-1: An in silico approach

The brown dog tick (Rhipicephalus sanguineus) is the most prevalent tick in the world and a well-recognized vector of many pathogens affecting dogs and occasionally humans. Pathogens exploit tick salivary molecules for their survival and multiplication in the vector and transmission to and establishment in the hosts. Tick saliva contains various non-proteinaceous substances and secreted proteins that are differentially produced during feeding and comprise of inhibitors of blood congealing and platelet aggregation, vasodilatory and immunomodulatory substances, and compounds preventing itch and pain. One of these proteins is Evasin-1, which has a high binding affinity to certain types of chemokines. The binding of Evasin-1 to chemokines prevents the detection and immune response of the host to R. sanguineus, which may result in the successful transmission of pathogens. In this study, we screened potential Evasin-1 inhibitor based on the pharmacophore model derived from the binding site residues. Hit ligands were further screened via molecular docking and virtual ADMET prediction, which resulted in ZINC8856727 as the top ligand (binding affinity: -9.1 kcal/mol). Molecular dynamics simulation studies, coupled with MM-GBSA calculations and principal component analysis revealed that ZINC8856727 plays a vital role in the stability of Evasin-1. We recommend continuing the study by developing a formulation that serves as a potential medicine aid immune response during R. sanguineus infestation.

Repurposing food molecules as a potential BACE1 inhibitor for Alzheimer's disease

Alzheimer's disease (AD) is a severe neurodegenerative disorder of the brain that manifests as dementia, disorientation, difficulty in speech, and progressive cognitive and behavioral impairment. The emerging therapeutic approach to AD management is the inhibition of β-site APP cleaving enzyme-1 (BACE1), known to be one of the two aspartyl proteases that cleave β-amyloid precursor protein (APP). Studies confirmed the association of high BACE1 activity with the proficiency in the formation of β-amyloid-containing neurotic plaques, the characteristics of AD. Only a few FDA-approved BACE1 inhibitors are available in the market, but their adverse off-target effects limit their usage. In this paper, we have used both ligand-based and target-based approaches for drug design. The QSAR study entails creating a multivariate GA-MLR (Genetic Algorithm-Multilinear Regression) model using 552 molecules with acceptable statistical performance (R 2 = 0.82, Q 2 loo = 0.81). According to the QSAR study, the activity has a strong link with various atoms such as aromatic carbons and ring Sulfur, acceptor atoms, sp2-hybridized oxygen, etc. Following that, a database of 26,467 food compounds was primarily used for QSAR-based virtual screening accompanied by the application of the Lipinski rule of five; the elimination of duplicates, salts, and metal derivatives resulted in a truncated dataset of 8,453 molecules. The molecular descriptor was calculated and a well-validated 6-parametric version of the QSAR model was used to predict the bioactivity of the 8,453 food compounds. Following this, the food compounds whose predicted activity (pKi) was observed above 7.0 M were further docked into the BACE1 receptor which gave rise to the Identification of 4-(3,4-Dihydroxyphenyl)-2-hydroxy-1H-phenalen-1-one (PubChem I.D: 4468; Food I.D: FDB017657) as a hit molecule (Binding Affinity = -8.9 kcal/mol, pKi = 7.97 nM, Ki = 10.715 M). Furthermore, molecular dynamics simulation for 150 ns and molecular mechanics generalized born and surface area (MMGBSA) study aided in identifying structural motifs involved in interactions with the BACE1 enzyme. Molecular docking and QSAR yielded complementary and congruent results. The validated analyses can be used to improve a drug/lead candidate's inhibitory efficacy against the BACE1. Thus, our approach is expected to widen the field of study of repurposing nutraceuticals into neuroprotective as well as anti-cancer and anti-viral therapeutic interventions.

Identification of phytochemicals from Eclipta alba and assess their potentiality against Hepatitis C virus envelope glycoprotein: virtual screening, docking, and molecular dynamics simulation study

Hepatitis C virus has a major role in spreading chronic liver disease and hepatocellular carcinoma. Factors such as high costs, pharmacological side effects, and the development of drug resistance strains require the development of new and potentially effective antiviral to treat the various stages of Hepatitis C. Bioactive chemicals have been extracted from medicinal plants and are utilized by humans for the goal of maintaining a healthy lifestyle. The goal of this work is to recognize phytochemicals from Eclipta alba and assess their potentiality activity against the hepatitis C virus envelope glycoprotein using in silico approaches. Phytochemicals from Eclipta alba were virtually screened by Auto dock raccoon and 12 compounds were selected for molecular docking to probe the active binding site. The top two compounds based on the binding score like ecliptalbine and oleanolic acid with HCV E2 glycoprotein exhibit binding energy -8.88 and -8.02 kcal/mol, respectively. The chemicals' usefulness was reinforced by positive pharmacokinetic data. The phytocompounds were identified as potent HCV inhibitors based on the drug likeness and ADMET properties. Both ecliptalbine and oleanolic acid underwent molecular dynamics simulations to determine features such as RMSD, RMSF, SASA, hydrogen-bond number, and MM-PBSA-based binding free energy. From the molecular docking and molecular dynamics simulation study revealed that oleanolic acid obtained from Eclipta alba can be used as inhibitors against Hepatitis C. The identified inhibitor from our study will be study in vitro and in vivo studies to check their efficacy against Hepatitis C.Communicated by Ramaswamy H. Sarma.

A cholesterol analog stabilizes the human β2-adrenergic receptor nonlinearly with temperature

In cell membranes, G protein-coupled receptors (GPCRs) interact with cholesterol, which modulates their assembly, stability, and conformation. Previous studies have shown how cholesterol modulates the structural properties of GPCRs at ambient temperature. Here, we characterized the mechanical, kinetic, and energetic properties of the human β2-adrenergic receptor (β2AR) in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) at room temperature (25°C), at physiological temperature (37°C), and at high temperature (42°C). We found that CHS stabilized various structural regions of β2AR differentially, which changed nonlinearly with temperature. Thereby, the strongest effects were observed for structural regions that are important for receptor signaling. Moreover, at 37°C, but not at 25° or 42°C, CHS caused β2AR to increase and stabilize conformational substates to adopt to basal activity. These findings indicate that the nonlinear, temperature-dependent action of CHS in modulating the structural and functional properties of this GPCR is optimized for 37°C.

Potent Bioactive Compounds From Seaweed Waste to Combat Cancer Through Bioinformatics Investigation

The seaweed industries generate considerable amounts of waste that must be appropriately managed. This biomass from marine waste is a rich source of high-value bioactive compounds. Thus, this waste can be adequately utilized by recovering the compounds for therapeutic purposes. Histone deacetylases (HDACs) are key epigenetic regulators established as one of the most promising targets for cancer chemotherapy. In the present study, our objective is to find the HDAC 2 inhibitor. We performed top-down in silico methodologies to identify potential HDAC 2 inhibitors by screening compounds from edible seaweed waste. One hundred ninety-three (n = 193) compounds from edible seaweeds were initially screened and filtered with drug-likeness properties using SwissADME. After that, the filtered compounds were followed to further evaluate their binding potential with HDAC 2 protein by using Glide high throughput virtual screening (HTVS), standard precision (SP), extra precision (XP), and quantum polarized ligand docking (QPLD). One compound with higher negative binding energy was selected, and to validate the binding mode and stability of the complex, molecular dynamics (MD) simulations using Desmond were performed. The complex-binding free energy calculation was performed using molecular mechanics-generalized born surface area (MM-GBSA) calculation. Post-MD simulation analyses such as PCA, DCCM, and free energy landscape were also evaluated. The quantum mechanical and electronic properties of the potential bioactive compounds were assessed using the density functional theory (DFT) study. These findings support the use of marine resources like edible seaweed waste for cancer drug development by using its bioactive compounds. The obtained results encourage further in vitro and in vivo research. Our in silico findings show that the compound has a high binding affinity for the catalytic site of the HDAC 2 protein and has drug-likeness properties, and can be utilized in drug development against cancer.

Nonlinear molecular dynamics of quercetin in Gynocardia odorata and Diospyros malabarica fruits: Its mechanistic role in hepatoprotection

Liver performs number of critical physiological functions in human system. Intoxication of liver leads to accumulation of free radicals that eventually cause damage, fibrosis, cirrhosis and cancer. Carbon tetrachloride (CCl₄) belongs to hepatotoxin is converted to a highly reactive free radical by cytochrome P450 enzymes that causes liver damage. Plant extracts derived quercetin has substantial role in hepatoprotection. This study highlights the possible mechanism by which quercetin plays significant role in hepatoprotection. HPLC analysis revealed the abundance of quercetin in the fruit extracts of Gynocardia odorata and Diospyros malabarica, were isolated, purified and subjected to liver function analysis on Wistar rats. Post quercetin treatment improved liver function parameters in the hepatotoxic Wistar rats by augmenting bilirubin content, SGOT and SGPT activity. Gene expression profile of quercetin treated rats revealed down regulation of HGF, TIMP1 and MMP2 expressed during CCl₄ induction. In silico molecular mechanism prediction suggested that quercetin has a high affinity for cell signaling pathway proteins BCL-2, JAK2 and Cytochrome P450 Cyp2E1, which all play a significant role in CCl₄ induced hepatotoxicity. In silico molecular docking and molecular dynamics simulation have shown that quercetin has a plausible affinity for major signaling proteins in liver. MMGBSA studies have revealed high binding of quercetin (ΔG) -41.48±11.02, -43.53±6.55 and -39.89±5.78 kcal/mol, with BCL-2, JAK2 and Cyp2E1, respectively which led to better stability of the quercetin bound protein complexes. Therefore, quercetin can act as potent inhibitor against CCl₄ induced hepatic injury by regulating BCL-2, JAK2 and Cyp2E1.

Identification of Potent Aldose Reductase Inhibitors as Antidiabetic (Anti-hyperglycemic) agents using QSAR Based Virtual Screening, Molecular Docking, MD Simulation and MMGBSA Approaches

The aldose reductase (AR) enzyme is an important target enzyme in the development of therapeutics against hyperglycaemia induced health complications such as retinopathy, etc. In the present study, a quantitative structure activity relationship (QSAR) evaluation of a dataset of 226 reported AR inhibitor (ARi) molecules is performed using a genetic algorithm – multi linear regression (GA-MLR) technique. Multi-criteria decision making (MCDM) analysis furnished two five variables based QSAR models with acceptably high performance reflected in various statistical parameters such as, R² = 0.79–0.80, Q²_LOO = 0.78–0.79, Q²_LMO = 0.78–0.79. The QSAR model analysis revealed some of the molecular features that play crucial role in deciding inhibitory potency of the molecule against AR such as; hydrophobic Nitrogen within 2 Å of the center of mass of the molecule, non-ring Carbon separated by three and four bonds from hydrogen bond donor atoms, number of sp2 hybridized Oxygen separated by four bonds from sp2 hybridized Carbon atoms, etc. 14 in silico generated hits, using a compound 18 (a most potent ARi from present dataset with pIC₅₀ = 8.04 M) as a template, on QSAR based virtual screening (QSAR-VS) furnished a scaffold 5 with better ARi activity (pIC₅₀ = 8.05 M) than template compound 18. Furthermore, molecular docking of compound 18 (Docking Score = –7.91 kcal/mol) and scaffold 5 (Docking Score = –8.08 kcal/mol) against AR, divulged that they both occupy the specific pocket(s) in AR receptor binding sites through hydrogen bonding and hydrophobic interactions. Molecular dynamic simulation (MDS) and MMGBSA studies right back the docking results by revealing the fact that binding site residues interact with scaffold 5 and compound 18 to produce a stable complex similar to co-crystallized ligand's conformation. The QSAR analysis, molecular docking, and MDS results are all in agreement and complementary. QSAR-VS successfully identified a more potent novel ARi and can be used in the development of therapeutic agents to treat diabetes.

Target Specific Inhibition of Protein Tyrosine Kinase in Conjunction With Cancer and SARS-COV-2 by Olive Nutraceuticals

The fact that viruses cause human cancer dates back to the early 1980s. By reprogramming cellular signaling pathways, viruses encoded protein that can regulate altered control of cell cycle events. Viruses can interact with a superfamily of membrane bound protein, receptor tyrosine kinase to modulate their activity in order to increase virus entrance into cells and promotion of viral replication within the host. Therefore, our study aimed at screening of inhibitors of tyrosine kinase using natural compounds from olive. Protein tyrosine kinase (PTK) is an important factor for cancer progression and can be linked to coronavirus. It is evident that over expression of Protein tyrosine kinase (PTK) enhance viral endocytosis and proliferation and the use of tyrosine kinase inhibitors reduced the period of infection period. Functional network studies were carried out using two major PTKs viz. Anaplastic lymphoma kinase (ALK) and B-lymphocytic kinase (BTK). They are associated with coronavirus in regulation of cell signaling proteins for cellular processes. We virtually screened for 161 library of natural compounds from olive found overexpressed in ALK and BTK in metastatic as well as virus host cells. We have employed both ligand and target-based approach for drug designing by high throughput screening using Multilinear regression model based QSAR and docking. The QSAR based virtual screening of 161 olive nutraceutical compounds has successfully identified certain new hit; Wedelosin, in which, the descriptor rsa (ratio of molecular surface area to the solvent accessible surface area) plays crucial role in deciding Wedelosin’s inhibitory potency. The best-docked olive nutraceuticals further investigated for the stability and effectivity of the BTK and ALK during in 150 ns molecular dynamics and simulation. Post simulation analysis and binding energy estimation in MMGBSA further revealed the intensive potential of the olive nutraceuticals in PTK inhibition. This study is therefore expected to widen the use of nutraceuticals from olive in cancer as well as SARS-CoV2 alternative therapy.

Dynamics of natural product Lupenone as a potential fusion inhibitor against the spike complex of novel Semliki Forest Virus

The Semliki Forest Virus (SFV) is an RNA virus with a positive-strand that belongs to the Togaviridae family’s Alphavirus genus. An epidemic was observed among French troops stationed in the Central African Republic, most likely caused by the SFV virus. The two transmembrane proteins El and E2 and the peripheral protein E3 make up the viral spike protein. The virus binds to the host cell and is internalized via endocytosis; endosome acidification causes the E1/E2 heterodimer to dissociate and the E1 subunits to trimerize. Lupenone was evaluated against the E1 spike protein of SFV in this study based on state-of-the-art cheminformatics approaches, including molecular docking, molecular dynamics simulation, and binding free energy calculation. The molecular docking study envisaged major interactions of Lupenone with binding cavity residues involved non-bonded van der Waal’s and Pi-alkyl interactions. Molecular dynamic simulation of a time scale 200 ns corroborated interaction pattern with molecular docking studies between Lupenone and E1 spike protein. Nevertheless, Lupenone intearcation with the E1 spike protein conforming into a stable complex substantiated by free energy landscape (FEL), PCA analysis. Free energy decomposition of the binding cavity resdiues of E1 spike protein also ensured the efficient non-bonded van der Waal’s interaction contributing most energy to interact with the Lupenone. Therefore, Lupenone interacted strongly at the active site conforming into higher structural stability throughout the dynamic evolution of the complex. Thus, this study perhaps comprehend the efficiency of Lupenone as lead molecule against SFV E1 spike protein for future therapeutic purpose.

Effect of Early Swimming on the Behavior and Striatal Transcriptome of the Shank3 Knockout Rat Model of Autism

BACKGROUND

Autism spectrum disorder (ASD) is a developmental disorder characterized by social behavior deficits and stereotyped behaviors in childhood that lacks satisfactory medical intervention. Early swimming intervention is a noninvasive method combining enriched environment and exercise, which has been proven to improve brain development in young children and to treat neurodevelopmental diseases.

METHODS

In this study, we tested the autism-like behavior of rats with deletions in exons 11-21 of the Shank3 gene and evaluated the effect of early swimming intervention (from postnatal day 8 to 60) on the behavior of this animal model of autism. In addition, the transcriptomes of the striatal tissues of wild-type, Shank3 knockout and Shank3 knockout swimming groups rats were analyzed.

RESULTS

Shank3 knockout rats exhibit core symptoms of autism, and early swimming improved the social and stereotyped behaviors in this autism rat model. Transcriptomics results revealed that compared to the wild-type group, 291 differentially expressed genes (DEGs) were identified in the striatum of the Shank3 knockout group. Compared to Shank3 knockout group, 534 DEGs were identified in the striatum of Shank3 knockout swimming group. The DEGs annotated by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway shows that the impacts of Shank3 deletion were primarily reflected in synaptic structure, development, morphology, receptor function and signaling, and swimming primarily changed the terms related to the synapses in the striatum of Shank3 knockout rats, including the morphology, structure, composition, development and regulation of synapses.

CONCLUSION

Early swimming intervention can ameliorate behavioral abnormalities caused by Shank3 knockout, by a mechanism that may involve the process of striatal synaptic development and should be further investigated.

A Conservative Replacement in the Transmembrane Domain of SARS-CoV-2 ORF7a as a Putative Risk Factor in COVID-19

Simple Summary

The pathogenicity and transmissibility of the COVID-19 pandemic causative agent, the SARS-CoV-2 virus, is related to the functions of the proteins synthesized intracellularly, as guided by viral RNA. It is vitally important to accurately pinpoint novel variants of concern of the SARS-CoV-2 virus, in order to understand the molecular features of novel mutations and manage the on-going battle against the COVID-19 pandemic. We focused on A105V mutation in the ORF7a accessory protein. Sequencing and clinical data showed that this mutation is associated with increased severity and lethality in a group of Romanian patients, despite a lower viral copy number and a lower number of associated comorbidities. This effect is primarily due to increased protein stability through allosteric effects as shown by molecular dynamics analyses. This behavior manifests especially among residues 39–56, and the ones adjacent to 26–30 loop, placed in direct contact with potential interaction partners. Together, the results provide novel insights into the role of ORF7a in the pathogenicity of SARS-CoV-2.

Abstract

The ongoing COVID-19 pandemic follows an unpredictable evolution, driven by both host-related factors such as mobility, vaccination status, and comorbidities and by pathogen-related ones. The pathogenicity of its causative agent, SARS-CoV-2 virus, relates to the functions of the proteins synthesized intracellularly, as guided by viral RNA. These functions are constantly altered through mutations resulting in increased virulence, infectivity, and antibody-evasion abilities. Well-characterized mutations in the spike protein, such as D614G, N439K, Δ69–70, E484K, or N501Y, are currently defining specific variants; however, some less studied mutations outside the spike region, such as p. 3691 in NSP6, p. 9659 in ORF-10, 8782C > T in ORF-1ab, or 28144T > C in ORF-8, have been proposed for altering SARS-CoV-2 virulence and pathogenicity. Therefore, in this study, we focused on A105V mutation of SARS-CoV-2 ORF7a accessory protein, which has been associated with severe COVID-19 clinical manifestation. Molecular dynamics and computational structural analyses revealed that this mutation differentially alters ORF7a dynamics, suggesting a gain-of-function role that may explain its role in the severe form of COVID-19 disease.

Type 2 Diabetes Mellitus Mediation by the Disruptive Activity of Environmental Toxicants on Sex Hormone Receptors: In Silico Evaluation

This study investigates the disruptive activity of environmental toxicants on sex hormone receptors mediating type 2 diabetes mellitus (T2DM). Toxicokinetics, gene target prediction, molecular docking, molecular dynamics, and gene network analysis were applied in silico techniques. From the results, permethrin, perfluorooctanoic acid, dichlorodiphenyltrichloroethane, O-phenylphenol, bisphenol A, and diethylstilbestrol were the active toxic compounds that could modulate androgen (AR) and estrogen-α and -β receptors (ER) to induce T2DM. Early growth response 1 (EGR1), estrogen receptor 1 (ESR1), and tumour protein 63 (TP63) were the major transcription factors, while mitogen-activated protein kinases (MAPK) and cyclin-dependent kinases (CDK) were the major kinases upregulated by these toxicants via interactions with intermediary proteins such as PTEN, AKT1, NfKβ1, SMAD3 and others in the gene network analysis to mediate T2DM. These toxicants pose a major challenge to public health; hence, monitoring their manufacture, use, and disposal should be enforced. This would ensure reduced interaction between people and these toxic chemicals, thereby reducing the incidence and prevalence of T2DM.

Mutations of SARS-CoV-2 RBD May Alter Its Molecular Structure to Improve Its Infection Efficiency

The receptor-binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mediates the viral-host interaction and is a target for most neutralizing antibodies. Nevertheless, SARS-CoV-2 RBD mutations pose a threat due to their role in host cell entry via the human angiotensin-converting enzyme 2 receptor that might strengthen SARS-CoV-2 infectivity, viral load, or resistance against neutralizing antibodies. To understand the molecular structural link between RBD mutations and infectivity, the top five mutant RBDs (i.e., N501Y, E484K L452R, S477N, and N439K) were selected based on their recorded case numbers. These mutants along with wild-type (WT) RBD were studied through all-atom molecular dynamics (MD) simulations of 100 ns. The principal component analysis and the free energy landscape were used too. Interestingly, N501Y, N439K, and E484K mutations were observed to increase the rigidity in some RBD regions while increasing the flexibility of the receptor-binding motif (RBM) region, suggesting a compensation of the entropy penalty. However, S477N and L452R RBDs were observed to increase the flexibility of the RBM region while maintaining similar flexibility in other RBD regions in comparison to WT RBD. Therefore, both mutations (especially S477N) might destabilize the RBD structure, as loose conformation compactness was observed. The destabilizing effect of S477N RBD was consistent with previous work on S477N mutation. Finally, the free energy landscape results showed that mutations changed WT RBD conformation while local minima were maintained for all mutant RBDs. In conclusion, RBD mutations definitely impact the WT RBD structure and conformation as well as increase the binding affinity to angiotensin-converting enzyme receptor.

Macrolactin A as a Novel Inhibitory Agent for SARS-CoV-2 Mpro: Bioinformatics Approach

COVID-19 is a disease that puts most of the world on lockdown and the search for therapeutic drugs is still ongoing. Therefore, this study used in silico screening to identify natural bioactive compounds from fruits, herbaceous plants, and marine invertebrates that are able to inhibit protease activity in SARS-CoV-2 (PDB: 6LU7). We have used extensive screening strategies such as drug likeliness, antiviral activity value prediction, molecular docking, ADME, molecular dynamics (MD) simulation, and MM/GBSA. A total of 17 compounds were shortlisted using Lipinski’s rule in which 5 compounds showed significant predicted antiviral activity values. Among these 5, only 2 compounds, Macrolactin A and Stachyflin, showed good binding energy of −9.22 and −8.00 kcal/mol, respectively, within the binding pocket of the M^pro catalytic residues (HIS 41 and CYS 145). These two compounds were further analyzed to determine their ADME properties. The ADME evaluation of these 2 compounds suggested that they could be effective in developing therapeutic drugs to be used in clinical trials. MD simulations showed that protein–ligand complexes of Macrolactin A and Stachyflin with the target receptor (6LU7) were stable for 100 nanoseconds. The MM/GBSA calculations of M^pro–Macrolactin A complex indicated higher binding free energy (−42.58 ± 6.35 kcal/mol). Dynamic cross-correlation matrix (DCCM) and principal component analysis (PCA) on the residual movement in the MD trajectories further confirmed the stability of Macrolactin A bound state with 6LU7. In conclusion, this study showed that marine natural compound Macrolactin A could be an effective therapeutic inhibitor against SARS-CoV-2 protease (6LU7). Additional in vitro and in vivo validations are strongly needed to determine the efficacy and therapeutic dose of Macrolactin A in biological systems.

Sesamin and sesamolin rescues Caenorhabditis elegans from Pseudomonas aeruginosa infection through the attenuation of quorum sensing regulated virulence factors

Pseudomonas aeruginosa is an opportunistic pathogen emerging as a public health threat owing to their multidrug resistance profiles. The quorum sensing systems of P. aeruginosa play a pivotal role in the regulation of virulence and act as the target for the development of alternative therapeutics. The study discussed about anti-quorum sensing and antibiofilm properties of lignans (sesamin and sesamolin) found in Sesamum indicum (L.) against P. aeruginosa. The effect of lignans, sesamin and sesamolin on LasR/RhlR mediated virulence factor production, biofilm formation and bacterial motility were studied. To elucidate the mechanism of action of lignans on QS pathways, QS gene expression and in depth in silico analysis were performed. Both the lignans exerted anti-quorum sensing activity at 75 μg/ml without affecting the growth of bacteria. SA and SO exhibited decreased production of virulence factors such as pyocyanin, proteases, elastase and chitinase. The important biofilm constituents of P. aeruginosa including alginate, exopolysaccharides and rhamnolipids were strongly affected by the lignans. Likewise, plausible mechanism of action of lignans were determined through the down regulation of QS regulated gene expression, molecular docking and molecular simulation studies. The in vitro analysis was supported by C. elegans infection model. SA and SO rescued pre-infected worms within 8 days of post infection and reduced the colonization of bacteria inside the intestine due to the anti-infective properties of lignans. The lignans exhibited profound action on Las pathway rather than Rhl which was elucidated through in vitro and in silico assays. In silico pharmacokinetic analysis portrayed the opportunities to employ ligands as potential therapeutics for human use. The deep insights into the anti-QS, anti-biofilm and mechanism of action of lignans can contribute to the development of novel anti-infectives against pseuodmonal infections.

In silico analysis of Phyllanthus amarus phytochemicals as potent drugs against SARS-CoV-2 main protease

Phyllanthus amarus, also known as Bhui Korma in India, is well known for its medicinal properties and is used to treat several diseases worldwide. This study aims to identify phytochemicals from P. amarus and assess their anti-viral activity through in silico methods against the main protease (3CL^Pro/M^Pro) enzyme of the novel coronavirus. 190 compounds were obtained from literature and docked against 3CLPro and 16 compounds showed higher binding affinity with 3CL^Pro with their values lying between -8.9 ​kcal/mol to -9.6 ​kcal/mol. The top two compounds, Myricitrin (CID: 5352000) and Quercetin-3-O-glucuronide (CID: 12004528) gave high binding affinity values of -9.6 ​kcal/mol and -9.4 ​kcal/mol respectively and also display favourable binding interactions with the 3CL^Pro. Both the compounds were further subjected to molecular dynamics simulation and MM-PBSA based binding free energy calculations. ADMET and drug-likeness properties were studied to assess the pharmacokinetic properties of the compounds. Favourable pharmacokinetic results reinforced the applicability of the compounds assessed. Along with continuous studies being carried out with chemical compounds, research needs to expand into all areas, including the use of natural compounds as drug compounds. The identified hits from this study can be taken further for in vitro and in vivo studies to examine their efficacy against COVID-19.

Understand the Functions of Scaffold Proteins in Cell Signaling by a Mesoscopic Simulation Method

Scaffold proteins are central players in regulating the spatial-temporal organization of many important signaling pathways in cells. They offer physical platforms to downstream signaling proteins so that their transient interactions in a crowded and heterogeneous environment of cytosol can be greatly facilitated. However, most scaffold proteins tend to simultaneously bind more than one signaling molecule, which leads to the spatial assembly of multimeric protein complexes. The kinetics of these protein oligomerizations are difficult to quantify by traditional experimental approaches. To understand the functions of scaffold proteins in cell signaling, we developed a, to our knowledge, new hybrid simulation algorithm in which both spatial organization and binding kinetics of proteins were implemented. We applied this new technique to a simple network system that contains three molecules. One molecule in the network is a scaffold protein, whereas the other two are its binding targets in the downstream signaling pathway. Each of the three molecules in the system contains two binding motifs that can interact with each other and are connected by a flexible linker. By applying the new simulation method to the model, we show that the scaffold proteins will promote not only thermodynamics but also kinetics of cell signaling given the premise that the interaction between the two signaling molecules is transient. Moreover, by changing the flexibility of the linker between two binding motifs, our results suggest that the conformational fluctuations in a scaffold protein play a positive role in recruiting downstream signaling molecules. In summary, this study showcases the capability of computational simulation in understanding the general principles of scaffold protein functions.

Using structural analysis to explore the role of hepatitis B virus mutations in immune escape from liver cancer in Chinese, European and American populations

Hepatitis B virus (HBV) infection is an important problem threatening human health. After HBV virus invades human body, it may assemble a complete virus particle in the cytoplasm to trigger the immune reaction, especially the interaction between the HBV virus and the host that mediated by CD8⁺ T cell. We collected the sequences of HBV from the HBVdb database, then screened candidate mutation sites in Chinese, European and American populations based on conservation and physicochemical properties. After that we constructed the three-dimensional structure of Major histocompatibility complex class I (MHC I) -peptide complexes, performed molecular docking, run molecular dynamics to compare the binding free energy, stability, and affinity of MHC I-peptide complexes with the aim to estimate the effect of peptide mutation. The specific HBV virus subtypes of the Chinese, European and American population were studied and the candidate mutation sites were used to predict the mutant peptide antigen. Finally, based on physical and chemical properties and peptide antigen prediction scores, 21 HBV mutation sites were selected. Then combined with specific Human lymphocyte antigen (HLA) subtypes, 11 mutations were found to have a significant negative impact on affinity, stability and binding free energy. Overall, our work found important potential mutations, which provide an evaluation of HBV mutations and a clue of it in immunotherapy.Communicated by Ramaswamy H. Sarma.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Microsecond Simulation of the Proteoglycan-like Region of Carbonic Anhydrase IX and Design of Chemical Inhibitors Targeting pH Homeostasis in Cancer Cells

Carbonic anhydrase IX (CAIX) is a membrane-bound enzyme associated with tumor hypoxia and found to be over expressed in various tumor conditions. Targeting CAIX catalytic activity is proven to be efficient modality in modulating pH homeostasis in cancer cells. Proteoglycan-like (PG) region is unique to CAIX and is proposed to serve as an antenna enhancing the export of protons in conjunction with facilitated efflux of lactate ions via monocarboxylate transporters. Moreover, the PG region is also reported to contribute to the assembly and maturation of focal adhesion links during cellular attachment and dispersion on solid supports. Thus, drug targeting of this region shall efficiently modulate pH homeostasis and cell adhesion in cancer cells. As the PG region is intrinsically disordered, the complete crystal structure is not elucidated. Hence, in this study, we intend to sample the conformational landscape of the PG region at microsecond scale simulation in order to sample the most probable conformations that shall be utilized for structure-based drug design. In addition, the sampled conformations were subjected to high-throughput virtual screening against NCI and Maybridge datasets to identify potential hits based on consensus scoring and validation by molecular dynamics simulation. Further, the identified hits were experimentally validated for efficacy by in vitro and direct enzymatic assays. The results reveal 5-(2-aminoethyl)-1,2,3-benzenetriol to be the most promising hit as it showed significant CAIX inhibition at all levels of in silico and experimental validation.

Exploration of the selective binding mechanism of protein kinase Aurora A selectivity via a comprehensive molecular modeling study

Background

The kinase of Aurora A has been regarded as a promising therapeutic target due to its altered expression in various human cancers. However, given the high similarity of the active binding site of Aurora A to other kinases, designing highly selective inhibitors towards Aurora A remains a challenge. Recently, two potential small-molecule inhibitors named AT9283 and Danusertib were reported to exhibit significant selectivity to Aurora A, but not to Gleevec. It was argued that protein dynamics is crucial for drug selectivity to Aurora A. However, little computational research has been conducted to shed light on the underlying mechanisms.

Methods

In this study, MM/GBSA calculations based on conventional molecular dynamics (cMD) simulations and enhanced sampling simulations including Gaussian accelerated MD (GaMD) simulations and umbrella sampling were carried out to illustrate the selectivity of inhibitors to Aurora A.

Results

The calculation results from cMD simulation showed that the binding specificity is primarily controlled by conformational change of the kinase hinge. The protein dynamics and energetic differences were further supported by the GaMD simulations. Umbrella sampling further proved that AT9283 and Danusertib have similar potential of mean force (PMF) profiles toward Aurora A in terms of PMF depth. Compared with AT9283 and Danusertib, Gleevec has much lower PMF depth, indicating that Gleevec is more easily dissociated from Aurora A than AT9283 and Danusertib. These results not only show the selective determinants of Aurora A, but also provide valuable clues for the further development of novel potent Aurora A selective inhibitors.