Identification and comparison of plant-derived scaffolds as selective CDK5 inhibitors against standard molecules: Insights from umbrella sampling simulations

2D-QSAR, molecular docking and MD simulation based virtual screening of the herbal molecules against Alzheimer's disorder: an approach to predict CNS activity

Acetylcholinesterase (AChE) is one of the key enzyme targets that have been used clinically for the management of Alzheimer's Disorder (AD). Numerous reports in the literature predict and demonstrate in-vitro, and in-silico anticholinergic activity of herbal molecules, however, majority of them failed to find clinical application. To address these issues, we developed a 2D-QSAR model that could efficiently predict the AChE inhibitory activity of herbal molecules along with predicting their potential to cross the blood-brain-barrier (BBB) to exert their beneficial effects during AD. Virtual screening of the herbal molecules was performed and amentoflavone, asiaticoside, astaxanthin, bahouside, biapigenin, glycyrrhizin, hyperforin, hypericin, and tocopherol were predicted as the most promising herbal molecules for inhibiting AChE. Results were validated through molecular docking, atomistic molecular dynamics simulations and Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) studies against human AChE (PDB ID: 4EY7). To determine whether or not these molecules can cross BBB to inhibit AChE within the central nervous system (CNS) for being beneficial for the management of AD, we determined a CNS Multi-parameter Optimization (MPO) score, which was found in the range of 1 to 3.76. Overall, the best results were observed for amentoflavone and our results demonstrated a PIC50 value of 7.377 nM, molecular docking score of -11.5 kcal/mol, and CNS MPO score of 3.76. In conclusion, we successfully developed a reliable and efficient 2D-QSAR model and predicted amentoflavone to be the most promising molecule that could inhibit human AChE enzyme within the CNS and could prove beneficial for the management of AD.Communicated by Ramaswamy H. Sarma.

An insight from computational approach to explore novel, high-affinity phosphodiesterase 10A inhibitors for neurological disorders

The enzyme Phosphodiesterase 10A (PDE10A) plays a regulatory role in the cAMP/protein kinase A (PKA) signaling pathway by means of hydrolyzing cAMP and cGMP. PDE10A emerges as a relevant pharmacological drug target for neurological conditions such as psychosis, schizophrenia, Parkinson's, Huntington's disease, and other memory-related disorders. In the current study, we subjected a set of 1,2,3-triazoles to be explored as PDE10A inhibitors using diverse computational approaches, including molecular docking, classical molecular dynamics (MD) simulations, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations, steered MD, and umbrella sampling simulations. Molecular docking of cocrystallized ligands papaverine and PFJ, along with a set of in-house synthesized molecules, suggested that molecule 3i haded the highest binding affinity, followed by 3h and 3j. Furthermore, the structural stability studies using MD and MM-PBSA indicated that the 3h and 3j formed stable complexes with PDE10A. The binding free energy of -240.642 kJ/mol and -201.406 kJ/mol was observed for 3h and 3j, respectively. However, the cocrystallized ligands papaverine and PFJ exhibited comparitively higher binding free energy values of -202.030 kJ/mol and -138.764 kJ/mol, respectively. Additionally, steered MD and umbrella sampling simulations provided conclusive evidence that the molecules 3h and 3j could be exploited as promising candidates to target PDE10A.Communicated by Ramaswamy H. Sarma.

Identification of potential mutational hotspots in serratiopeptidase to address its poor pH tolerance issue

Serratiopeptidase is the multifunctionality metalloendopeptidase extensively employed in biopharmaceutical and industrial biotechnology. Despite its poor pH tolerance, serratiopeptidase must withstand the highly acidic environment of the gastrointestinal tract to be used as a potent anti-inflammatory and analgesic medication. In earlier studies, post-translational deamination related mutations showed alteration in the net charge of protein's surface. Therefore, the current study aimed to enhance the acid resistance of serratiopeptidase via implementing computational interventions to screen out the most stable mutational hotspot. The methodology used in this study is as follows: (a) Higher accessibility to surface (b) 4 Å away from active site region to avoid interference with its proteolytic activity, and (c) By converting non-conserved amide residues to acidic residues. A docking study has been conducted to establish the substrate specificity and binding affinity to native and mutant proteins. The docking outcomes were then validated using molecular dynamic simulations to clarify each mutant's molecular stability and conformation while preserving their activity. The results showed that N412D is the best-screened mutant with negative electrostatic potential that can alter the overall charge on the protein's surface with increased H+ ions. Alteration in overall charge leads to protein surface more acidic that causes a common ion effect in stomach pH and act as a buffer which could stabilize the serratiopeptidase amid extreme pH.Communicated by Ramaswamy H. Sarma.

Modulatory effects of phytol on the antiemetic property of domperidone, possibly through the D2 receptor interaction pathway: in vivo and in silico studies

The current study is designed to evaluate the antiemetic effect of the diterpenoid phytol (PHY) using in vivo and in silico studies. For this, emesis was induced in 4-day-old chicks by the oral administration of copper sulfate (CuSO4.5H2O) at 50 mg/kg. To see the possible antiemetic mechanism of PHY, we used a number of reference drugs such as domperidone (80 mg/kg), ondansetron (24 mg/kg) and hyoscine (100 mg/kg) as positive controls, while the vehicle served as a negative control group. PHY was administered orally at the doses of 50 and 75 mg/kg. Both PHY and reference drugs were given alone or in combined groups to evaluate their synergistic or antagonistic effects on the chicks. Molecular docking of PHY and reference drugs was carried out against 5HT3, D2, D3, H1, NK1, and mAChRs (M1-M5) receptors for estimating binding affinity to the receptors. Drug-receptor interactions and active sites of the receptors were observed with the aid of different computational tools. The drug-likeness and pharmacokinetics of all the drugs were predicted through the SwissADME online database. The results suggest that PHY reduces the mean number of retches and increases latency dose-dependently in the birds. In the combination groups, PHY75 showed better antiemetic effects with domperidone and ondansetron. In addition, PHY exhibited the highest binding affinity with the D2 receptor (6CM4) (- 7.3 kcal/mol). In conclusion, PHY showed an antiemetic activity in chicks, possibly through the D2 receptor interaction pathway.

Development of 9-(N-arylmethylamino) congeners of noscapine: the microtubule targeting drugs for the management of breast cancer

Noscapine is a natural lead molecule with anticancer activity at a higher concentrations. So, there is an urge for the development of more potent derivatives of noscapine. In this study, we have approached for development of 9-N-arylmethylamino derivatives of noscapine that kills cancer cells without affecting the normal cells. They were designed by substituting N-aryl methyl pharmacophore at the C-9 position and screened out top-ranked three derivatives 13a-c using molecular docking. Further, their theoretical free energy of binding with tubulin was calculated followed by chemical synthesis and experimental validation. In vitro antiproliferative activity of noscapine and its 9-N-arylmethylamino derivatives (13a-c) was carried out using MCF-7 (a triple receptors positive) and MDA-MB-231 (a triple receptor negative) breast cancer cell lines. Further, cytotoxicity to normal cells was examined using human embryonic kidney cells (HEK cells). Inhibition to cell cycle progression and induction of apoptosis was monitored using FACS. The binding of noscapine and 13a-c with tubulin was examined using fluorescence quenching assay. The 9-N-arylmethylamino derivatives of noscapine (13a-c) were found to inhibit the proliferation of cancer cells at a much lower concentration (IC50 values range between 9.1 to 47.3 µM) compared to noscapine (IC50 value is 45.8-59.3 µM). Surprisingly, the proliferation of HEK cells was not inhibited even at a concentration of 100 µM (cytotoxicity is < 5%). These derivatives induced apoptosis by arresting cells at G2/M-phase and also bind to tubulin. The 9-(N-arylmethylamino) noscapinoids have the potential to be a novel therapeutic agent for the treatment of breast cancer.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03445-3.

Network pharmacology-based analysis of signaling pathways of an anti-osteoporotic triterpenoid from Acyranthes bidentata Blume root

In Korea folk remedies, Acyranthes bidentata Blume is a functional food plant to treat bone diseases; especially, its roots have been used to alleviate osteoporosis (OP), but its key chemical compound(s) and mechanism of action against osteoporosis have not reported yet. This study suggests that Acyranthes bidentata Blume root (ABBR) has promising compound(s) against OP. We utilized network pharmacology to evaluate the therapeutic value. The chemical compounds from Acyranthes bidentata Blume root (ABBR) were identified by gas chromatography-mass spectrum (GC-MS); their physicochemical properties have been evaluated by SwissADME. Next, the target(s) related to a triterpenoid or OP-related targets were investigated by public databases. The signaling pathways from final targets were visualized, constructed, and analyzed by RPackage. Finally, we performed a molecular docking (MD) to explore key target(s) and compound(s) by employing AutoDockVina tools; the residues of amino acids interacted with ligands were identified by LigPlot + v.22. A total of 24 chemicals were accepted by the Lipinski's rules. We found a sole triterpenoid from ABBR via GC-MS, suggesting that might be a potent ligand to alleviate OP. Thereby, the 42 targets were associated with the triterpenoid; the 19 targets among them were connected to OP-targets (1426). The final 19 targets were related directly to 8 signaling pathways on STRING database. On Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and a key signaling pathway (PPAR signaling pathway), four key targets (PPARA, PPARD, FABP3, and FABP4) and a key compound (Methyl 3β-hydroxyolean-18-en-28-oate) were selected via MD. Collectively, the triterpenoid from ABBR might have potent anti-osteoporotic efficacy by activating PPARA, PPARD, FABP3, and FABP4 on PPAR signaling pathway.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03362-5.

The interaction analysis between human serum albumin with chlorpyrifos and its derivatives through sub-atomic docking and molecular dynamics simulation techniques

Chlorpyrifos (CPF) is an extensively used organophosphate pesticide for crop protection. However, there are concerns about it contaminating the environment and human health, with estimated three lakh deaths annually. The molecular modeling protocol was assisted in redesigning thirteen well-known CPF linkers and inserting them at five selectable CPF (R1-R5) positions of CPF to get 258 CPF derivatives. CPF and its derivatives were optimized using LigPrep and docked to a grid centralized on Trp214 using extra precision glide docking. The Binding free energy of complexes was calculated using molecular mechanics/generalized born surface area (MM-GBSA). CPF and CPFD-225 have glide scores of - 3.08 and - 6.152 kcal/mol, respectively, with human serum albumin and ΔG bind for CPF (- 33.041817 kcal/mol) (- 52.825 kcal/mol) for CPF-D225. The top ten CPF derivatives showed at least ninefold better binding free energy than the CPF proposed for polyclonal antibody production. Subsequently, molecular docking studies revealed that CPF and its derivatives could bind to human serum albumin (HSA). Furthermore, using the Desmond package, a 100-ns molecular dynamics (MD) simulation was performed on the potential binding site. The final systems of CPF-HSA and CPF-222D complexes consist of 76,014 and 76,026 atoms, respectively. The physical stability of both the systems (CPF-HSA and CPF-222D) was analyzed by considering the overall potential energy, RMSF, RMSD, Hydrophobic interactions, and water-mediated patterns, which showed total energy of - 141,610 kcal/mol and - 140,150 kcal/mol, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03344-7.

Computational targeting of allosteric site of MEK1 by quinoline-based molecules

MEK1 is an attractive target due to its role in selective extracellular-signal-regulated kinase phosphorylation, which plays a pivotal role in regulating cell proliferation. Another benefit of targeting the MEK protein is its unique hydrophobic pocket that can accommodate highly selective allosteric inhibitors. To date, various MEK1 inhibitors have reached clinical trials against several cancers, but they were discarded due to their severe toxicity and low efficacy. Thus, the development of allosteric inhibitors for MEK1 is the demand of the hour. In this in-silico study, molecular docking, long-term molecular dynamics (5 µs), and molecular mechanics Poisson-Boltzmann surface area analysis were undertaken to address the potential of quinolines as allosteric inhibitors. We selected four reference MEK1 inhibitors for the comparative analysis. The drug-likeness and toxicity of these molecules were also examined based on their ADMET and Toxicity Prediction by Komputer Assisted Technology profiles. The outcome of the analysis revealed that the quinolines (4m, 4o, 4s, and 4n) exhibited better stability and binding affinity while being nontoxic compared to reference inhibitors. We have reached the conclusion that these quinoline molecules could be checked by experimental studies to validate their use as allosteric inhibitors against MEK1.

Identification of 11β-HSD1 inhibitors through enhanced sampling methods

Aminoarylbenzosuberene (AAB) molecules were chosen for in silico analysis to develop effective and more competent 11β-hydroxysteroid dehydrogenase (11β-HSD1) protein inhibitors. The AAB4 molecule was shown to have stronger interactions and binding affinity than standard inhibitors (co-crystallized molecules). These results were based on conventional, steered and enhanced umbrella sampling simulations.