Molecular docking of substituted pteridinones and pyrimidines to the ATP-binding site of the N-terminal domain of RSK2 and associated MM/GBSA and molecular field datasets

Trypsin is inhibited by phytocompounds liquiritin and terpinen-4-ol from the herb Glycyrrhiza glabra: in vitro and in silico studies

Serine proteases are a class of hydrolytic enzymes involved in various physiological functions like digestion, coagulation, fibrinolysis and immunity. The present study evaluates the serine protease inhibitory potential of phytochemicals liquiritin and terpinen-4-ol present in the herb Glycyrrhiza glabra L. using trypsin as the model enzyme. In silico studies showed that both the compounds have a significant binding affinity towards trypsin with a binding energy of -26.66 kcal/mol and -19.79 kcal/mol for liquiritin and terpinen-4-ol, respectively. Their binding affinity was confirmed through in vitro enzyme inhibition assays. The mode of inhibition was found to be uncompetitive. In order to explain the mode of inhibition, docking of the ligands to the enzyme-substrate complex was also done and binding energy was calculated after MD simulation. The energy values showed that the binding affinities of these compounds towards the enzyme substrate complex are more than that towards the enzyme alone. This explains the uncompetitive mode of inhibition.Communicated by Ramaswamy H. Sarma.

Discovery of small molecule entry inhibitors targeting the linoleic acid binding pocket of SARS-CoV-2 spike protein

Spike glycoprotein has a significant role in the entry of SARS-CoV-2 to host cells, which makes it a potential drug target. Continued accumulation of non-synonymous mutations in the receptor binding domain of spike protein poses great challenges in identifying antiviral drugs targeting this protein. This study aims to identify potential entry inhibitors of SARS-CoV-2 using virtual screening and molecular dynamics (MD) simulations from three distinct chemical libraries including Pandemic Response Box, Drugbank and DrugCentral, comprising 6971 small molecules. The molecules were screened against a binding pocket identified in the receptor-binding domain (RBD) region of the spike protein which is known as the linoleic acid binding pocket, a highly conserved motif among several SARS-CoV-2 variants. Through virtual screening and binding free energy calculations, we identified four top-scoring compounds, MMV1579787 ([2-Oxo-2-[2-(3-phenoxyphenyl)ethylamino]ethyl]phosphonic acid), Tretinoin, MMV1633963 ((2E,4E)-5-[3-(3,5-dichlorophenoxy)phenyl]penta-2,4-dienoic acid) and Polydatin, which were previously reported to have antibacterial, antifungal or antiviral properties. These molecules showed stable binding on MD simulations over 100 ns and maintained stable interactions with TYR365, PHE338, PHE342, PHE377, TYR369, PHE374 and LEU368 of the spike protein RBD that are found to be conserved among SARS-CoV-2 variants. Our findings were further validated with free energy landscape, principal component analysis and dynamic cross-correlation analysis. Our in silico analysis of binding mode and MD simulation analyses suggest that the identified compounds may impede viral entrance by interacting with the linoleic acid binding site of the spike protein of SARS-CoV-2 regardless of its variants, and they thus demand for further in vitro and in vivo research.Communicated by Ramaswamy H. Sarma.

Screening of phytochemicals from Clerodendrum inerme (L.) Gaertn as potential anti-breast cancer compounds targeting EGFR: an in-silico approach

Breast cancer (BC) is the most prevalent malignancy among women around the world. The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor (RTK) of the ErbB/HER family. It is essential for triggering the cellular signaling cascades that control cell growth and survival. However, perturbations in EGFR signaling lead to cancer development and progression. Hence, EGFR is regarded as a prominent therapeutic target for breast cancer. Therefore, in the current investigation, EGFR was targeted with phytochemicals from Clerodendrum inerme (L.) Gaertn (C. inerme). A total of 121 phytochemicals identified by gas chromatography-mass spectrometry (GC-MS) analysis were screened against EGFR through molecular docking, ADMET analyses (Absorption, Distribution, Metabolism, Excretion, and Toxicity), PASS predictions, and molecular dynamics simulation, which revealed three potential hit compounds with CIDs 10586 [i.e. alpha-bisabolol (-6.4 kcal/mol)], 550281 [i.e. 2,(4,4-Trimethyl-3-hydroxymethyl-5a-(3-methyl-but-2-enyl)-cyclohexene) (-6.5 kcal/mol)], and 161271 [i.e. salvigenin (-7.4 kcal/mol)]. The FDA-approved drug gefitinib was used to compare the inhibitory effects of the phytochemicals. The top selected compounds exhibited good ADMET properties and obeyed Lipinski's rule of five (ROF). The molecular docking analysis showed that salvigenin was the best among the three compounds and formed bonds with the key residue Met 793. Furthermore, the molecular mechanics generalized born surface area (MMGBSA) calculations, molecular dynamics simulation, and normal mode analysis validated the binding affinity of the compounds and also revealed the strong stability and compactness of phytochemicals at the docked site. Additionally, DFT and DOS analyses were done to study the reactivity of the compounds and to further validate the selected phytochemicals. These results suggest that the identified phytochemicals possess high inhibitory potential against the target EGFR and can treat breast cancer. However, further in vitro and in vivo investigations are warranted towards the development of these constituents into novel anti-cancer drugs.Communicated by Ramaswamy H. Sarma.

Development of novel thiazolidine-2,4-dione derivatives as PPAR-γ agonists through design, synthesis, computational docking, MD simulation, and comprehensive in vitro and in vivo evaluation

The present study was conducted to develop new novel 2,4-thiazolidinedione derivatives (3h-3j) as peroxisome proliferator-activated receptor-γ (PPAR-γ) modulators for antidiabetic activity. The objective was to overcome the adverse effects of existing thiazolidinediones while maintaining their pharmacological benefits. The synthesized compounds were elucidated based on FT-IR, 1H-NMR, 13C-NMR, and MS techniques. Molecular docking was utilized to investigate the interaction binding modes, binding free energy, and amino acids engaged in the compounds' interactions with the target protein. Subsequently, molecular dynamics modelling was used to assess the stability of the top-docked complexes and an assay was utilized to assess the cytotoxicity of the compounds to C2C12 myoblasts. Compounds 3h-3j exhibited PPAR-γ modulatory activity and demonstrated significant hypoglycaemic effects when compared to the reference drug pioglitazone. The new compounds were evaluated for their in vivo blood glucose-lowering potential by using a dexamethasone-induced diabetic rat model. All the compounds showed a hypoglycaemic effect of 108.04 ± 4.39, 112.55 ± 6.10, and 117.48 ± 43.93, respectively, along with pioglitazone (153.93 ± 4.61) compared to the diabetic control. Additionally, all the compounds significantly reduced AST and ALT levels and did not cause liver damage.

Identification of coastal pesticide pollutants as potent inhibitors of Bacillus pasteurii urease mediated calcium carbonate precipitation: a computational approach

Microbially induced calcite precipitation (MICP) through urease enzyme has attained a lot of recognition in various fields of civil engineering and geotechnology for stabilizing the strength of soil and various concrete materials. The activity of urease has been found to be affected by various factors like temperature, substrate concentrations, pH of the medium, presence of inhibitors, etc. Through this study, the outcome of the interaction of pesticides (commonly found in Indian coastal regions) on Bacillus pasteurii urease, a major organism reported for MICP studies has been investigated in silico. The results from the study revealed that the enzyme has higher interactions of -4.1, -3.2, and -3.4 kJ/mol with common pesticides like dichloro diphenyl dichloro ethane(DDD), dichloro diphenyl trichloroe thane (DDT), and methyl parathion of organochlorides and organophosphates class. From the molecular dynamics simulation analysis, complex 1 (DDD -receptor) has been found to have the highest and more compact structure followed by methyl parathion -receptor. Prime MM-GBSA analysis also revealed the highest binding energy of -27.8 kcal/mol with the protein and DDD. Thus, it can be inferred from the current study that pesticides, particularly, DDD, DDT, and methyl parathion present in the coastal areas may have an impact on urease. This interaction can result in the inhibition of the urease activity of B. pasteurii, thus preventing the biomineralization process. This study would be the first report on the computational approach to understanding the interaction of prominent pesticides on the coastal region and B. pasteurii urease.Communicated by Ramaswamy H. Sarma.

Comprehensive molecular docking and dynamic simulations for drug repurposing of clinical drugs against multiple cancer kinase targets

Drug repurposing is a method to identify novel therapeutic agents from the existing drugs and clinical compounds. In the present comprehensive work, molecular docking, virtual screening and dynamics simulations were carried out for ten cancer types viz breast, colon, central nervous system, leukaemia, melanoma, ovarian, prostate, renal and lung (non-small and small cell) against validated eighteen kinase targets. The study aims to understand the action of chemotherapy drugs mechanism through binding interactions against selected targets via comparative docking simulations with the state-art molecular modelling suits such as MOE, Cresset-Flare, AutoDock Vina, GOLD and GLIDE. Chemotherapeutic drugs (n = 112) were shortlisted from standard drug databases with appropriate chemoinformatic filters. Based on docking studies it was revealed that leucovorin, nilotinib, ellence, thalomid and carfilzomib drugs possessed potential against other cancer targets. A library was built to enumerate novel molecules based on the scaffold and functional groups extracted from known drugs and clinical compounds. Twenty novel molecules were prioritised further based on drug-like attributes. These were cross docked against 1MQ4 Aurora-A Protein Kinase for prostate cancer and 4UYA Mitogen-activated protein kinase for renal cancer. All docking programs yielded similar results but interestingly AutoDock Vina yielded the lowest RMSD with the native ligand. To further validate the final docking results at atomistic level, molecular dynamics simulations were performed to ascertain the stability of the protein-ligand complex. The study enables repurposing of drugs and lead identification by employing a host of structure and ligand based virtual screening tools and techniques.Communicated by Ramaswamy H. Sarma.

Role of indole curcumin in the epigenetic activation of apoptosis and cell cycle regulating genes

Background

Head-and-neck squamous cell carcinoma is associated with the epigenetic silencing of various genes such as DAPK, ataxia telangiectasia mutated (ATM), BRCA1, p16INK4a, pVHL, p16, and RASSF1A. The most common epigenetic change observed in these genes is DNA methylation that directs the studies toward finding inhibitors for DNA methyltransferases (DNMTs), the protagonist in the action. The present study focuses on analyzing the possibility whether indole curcumin can reverse epigenetic changes of the various tumor suppressor genes, characteristically silenced by methylation, by inhibiting the major methylation enzyme DNA methyltransferase 1 or DNMT1.

Materials and Methods

The cytotoxic effects of indole curcumin were studied through the MTT and lactate dehydrogenase assays. To determine the apoptosis-mediated death of HEp-2 cells, fluorescence imaging using different stains was done. Gene or mRNA expression analysis was done for p53, ATM, and DAPK genes.

Results

The results obtained from this study clearly indicate that the indole analog of curcumin plays a remarkable role in activating genes involved in cell cycle regulation and apoptosis induction through epigenetic regulation. The influence that the drug has on the methylation status of gene promoter sequence of the ATM gene is also very significant.

Conclusion

Indole curcumin, being an analog of curcumin, promises to be a very useful drug molecule having various potential targets. The target selected for this study was DNMT1 enzyme and the drug seems to actually show the effects; it was predicted to be having on the target molecule.

Computational Target-Based Screening of Anti-MRSA Natural Products Reveals Potential Multitarget Mechanisms of Action through Peptidoglycan Synthesis Proteins

Methicillin-resistant Staphylococcus aureus (MRSA) is one of the leading causes of bacterial infections in both healthcare and community settings. MRSA can acquire resistance to any current antibiotic, which has major implications for its current and future treatment options. As such, it is globally a major focus for infection control efforts. The mechanical rigidity provided by peptidoglycans in the bacteria cell walls makes it a promising target for broad-spectrum antibacterial drug discovery. The development of drugs that can target different stages of the synthesis of peptidoglycan in MRSA may compromise the integrity of its cell wall and consequently result in the rapid decline of diseases associated with this drug-resistant bacteria. The present study is aimed at screening natural products with known in vitro activities against MRSA to identify their potential to inhibit the proteins involved in the biosynthesis of the peptidoglycan cell wall. A total of 262 compounds were obtained when a literature survey was conducted on anti-MRSA natural products (AMNPs). Virtual screening of the AMNPs was performed against various proteins (targets) that are involved in the biosynthesis of the peptidoglycan (PPC) cell wall using Schrödinger software (release 2020-3) to determine their binding affinities. Nine AMNPs were identified as potential multitarget inhibitors against peptidoglycan biosynthesis proteins. Among these compounds, DB211 showed the strongest binding affinity and interactions with six protein targets, representing three stages of peptidoglycan biosynthesis, and thus was selected as the most promising compound. The MD simulation results for DB211 and its proteins indicated that the protein-ligand complexes were relatively stable over the simulation period of 100 ns. In conclusion, DB211 showed the potential to inhibit six proteins involved in the biosynthesis of the peptidoglycan cell wall in MRSA, thus reducing the chance of MRSA developing resistance to this compound. Therefore, DB211 provided a starting point for the design of new compounds that can inhibit multiple targets in the biosynthesis of the peptidoglycan layer in MRSA.

Cu(II)-Catalysed Hydrocarboxylation of Imines Utilizing CO2 to Synthesize α-Unsaturated Aminocarboxylic Acids

Here, we report the Cu(II)-photocatalysed hydrocarboxylation of imines (C=N) from a series of synthesized Schiff Base derivatives, namely (E)-1-(4-((4-methylbenzylidene)amino)phenyl)ethanone, (E)-1-(3-((5-bromo-2-hydroxybenzylidene)amino)phenyl)ethanone, (E)-4-((5-bromo-2-hydroxybenzylidene)amino)-1,5-dimethyl-2-phenyl-1H-pyrazol-3(2H)-one, and (E)-1,5-dimethyl-4-((4-methylbenzylidene)amino)-2-phenyl-1H-pyrazol-3(2H)-one, with carbon dioxide (CO₂) to generate disubstituted amino acids. Under mild conditions (atmospheric pressure of CO₂, room temperature, and 30 W Blue LED light), good to excellent yields confirming the formation of substituted amino acid unsaturated acid derivatives were obtained. Single crystal X-ray diffraction (SC-XRD) and UV-Vis diffuse reflectance spectroscopy (UV-Vis-DRS) confirmed the square pyramidal geometry of the Cu(II) photocatalyst. Docking and DFT calculations of the substituted amino acid unsaturated acid derivatives showed their potential as antimicrobial molecules.

Optimizing the Sunitinib for cardio-toxicity and thyro-toxicity by scaffold hopping approach

Sunitinib is a potent anti-cancer scaffold that acts as a VEGFR-2 inhibitor. Although the scaffold exhibits potent anti-cancer activity, it is cardiotoxic and also induces hypothyroidism. The current research aims to optimize the Sunitinib for cardio-toxicity and thyro-toxicity by scaffold hopping approach using the admetSAR server. The server has optimized the physico-chemical properties of Sunitinib, which were contributing to the cardiotoxicity and thyro-toxicity. The library of the optimized compounds was further screened by the molecular docking studies and results were validated by the MD simulation and DFT analysis for VEGFR-2 inhibition. Compounds 163 and 432 exhibited the highest affinity to VEGFR-2 receptor with minimal cardiotoxicity and thyro-toxicity. These two compounds could be the starting point for the further discovery of angiogenic inhibitors.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-022-00125-1.

3D-QSAR and scaffold hopping based designing of benzo[d]ox-azol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one derivatives as selective aldehyde dehydrogenase 1A1 inhibitors: Synthesis and biological evaluation

Aldehyde dehydrogenase 1 (ALDH1A1), an oxidoreductase class of enzymes, is overexpressed in various types of cancer cell lines and is the major cause of resistance to the Food and Drug Administration (FDA)-approved drug, cyclophosphamide (CP). In cancer conditions, CP undergoes a sequence of biotransformations to form an active metabolite, aldophosphamide, which further biotransforms to its putative cytotoxic metabolite, phosphoramide mustard. However, in resistant cancer conditions, aldophosphamide is converted into its inactive metabolite, carboxyphosphamide, via oxidation with ALDH1A1. Herein, to address the issue of ALDH1A1 mediated CP resistance, we report a series of benzo[d]oxazol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one derivatives as selective ALDH1A1 inhibitors. These inhibitors were designed using a validated 3D-quantitative structure activity relationship (3D-QSAR) model coupled with scaffold hopping. The 3D-QSAR model was developed using reported indole-2,3-diones based ALDH1A1 inhibitors, which provided field points in terms of electrostatic, van der Waals and hydrophobic potentials required for selectively inhibiting ALDH1A1. The most selective indole-2,3-diones-based compound, that is, cmp 3, was further considered for scaffold hopping. Two top-ranked bioisosteres, that is, benzo[d]oxazol-2(3H)-one and 2-oxazolo[4,5-b]pyridin-2(3H)-one, were selected for designing new inhibitors by considering the field pattern of 3D-QSAR. All designed molecules were mapped perfectly on the 3D-QSAR model and found to be predictive with good inhibitory potency (pIC₅₀ range: 7.5-6.8). Molecular docking was carried out for each designed molecule to identify key interactions that are required for ALDH1A1 inhibition and to authenticate the 3D-QSAR result. The top five inhibitor-ALDH1A1 complexes were also submitted for molecular dynamics simulations to access their stability. In vitro enzyme assays of 21 compounds suggested that these compounds are selective toward ALDH1A1 over the other two isoforms, that is, ALDH2 and ALDH3A1. All the compounds were found to be at least three and two times more selective toward ALDH1A1 over ALDH2 and ALDH3A1, respectively. All the compounds showed an IC₅₀ value in the range of 0.02-0.80 μM, which indicates the potential for these to be developed as adjuvant therapy for CP resistance.

Atomic level and structural understanding of natural ligands inhibiting Helicobacter pylori peptide deformylase through ligand and receptor based screening, SIFT, molecular dynamics and DFT - a structural computational approach

Helicobacter pylori is a Gram-negative microaerophilic gastric pathogen, responsible for the cause of peptic ulcer around half of the global population. Although several antibiotics and combination therapies have been employed for H. pylori-related gastric ulcer and cancer regiments, identifying potent inhibitors for specific targets of this bacterium will help assessing better treatment periodicity and methods to eradicate H. pylori. Herein, 1,000,000 natural compounds were virtually screened against Helicobacter pylori Peptide deformylase (HpPDF). Pharmacophore hypotheses were created using ligand and receptor-based pharmacophore modeling of GLIDE. Stringent HTVS and IFD docking protocol of GLIDE predicted leads with stable intermolecular bonds and scores. Molecular dynamics simulation of HpPDF was carried out for 100 ns using GROMACS. Hits ZINC00225109 and ZINC44896875 came up with a glide score of -9.967 kcal/mol and -12.114 kcal/mol whereas; reference compound actinonin produced a glide score of -9.730 kcal/mol. Binding energy values of these hits revealed the involvement of significant Van der Waals and Coulomb forces and the deduction of lipophilic forces that portray the deep hydrophobic residues in the S1pocket of H. pylori. The DFT analysis established the electron density-based features of the molecules and observed that the results correlate with intermolecular docking interactions. Analysis of the MD trajectories revealed the crucial residues involved in HpPDF - ligand binding and the conformational changes in the receptor. We have identified and deciphered the crucial features necessary for the potent ligand binding at catalytic site of HpPDF. The resulting ZINC natural compound hits from the study could be further employed for potent drug development.Communicated by Ramaswamy H. Sarma.

Substituted pteridinones, pyrimidines, pyrrolopyrimidines, and purines as p90 ribosomal S6 protein kinase-2 (RSK2) inhibitors: Pharmacophore modeling data

The RSK2 kinase is a downstream effector of the Ras/Raf/MEK/ERK pathway that is aberrantly active in a range of cancer types and has been recognized an anticancer target. The inhibition of RSK2 kinase activity would disrupt multiple pro-cancer processes; however, there are few RSK2 inhibitors. The data have been obtained for a series of pteridinone-, pyrimidine-, purine-, and pyrrolopyrimidine-based compounds, developed to establish a structure-activity relationship for RSK inhibition. The compounds were docked into the ATP-binding site of the N-terminal domain of the RSK2 kinase using Glide. The binding conformations of these molecules was then used to generate a set of pharmacophore models to determine the structural requirements for RSK2 inhibition. Through the combination of these models, common features (pharmacophores) can be identified that can inform the development of further small molecule RSK inhibitors. The synthesis and evaluation of the pteridinone- and pyrimidine-based compounds was reported in the related articles: Substituted pteridinones as p90 ribosomal S6 protein kinase (RSK) inhibitors: A structure-activity study (Casalvieri et al., 2020) and Molecular docking of substituted pteridinones and pyrimidines to the ATP-binding site of the N-terminal domain of RSK2 and associated MM/GBSA and molecular field datasets (Casalvieri et al., 2020). [1], [2]. The synthesis and evaluation of the purine- and pyrrolopyrimidine-based compounds was reported in the related research article: N-substituted pyrrolopyrimidines and purines as p90 ribosomal S6 protein kinase-2 (RSK2) inhibitors (Casalvieri et al., 2021) [3].

Novel thiomorpholine tethered isatin hydrazones as potential inhibitors of resistant Mycobacterium tuberculosis

Novel chemotherapeutic agents against multidrug resistant-tuberculosis (MDR-TB) are urgently needed at this juncture to save the life of TB-infected patients. In this work, we have synthesized and characterized novel isatin hydrazones 4(a-o) and their thiomorpholine tethered analogues 5(a-o). All the synthesized compounds were initially screened for their anti-mycobacterial activity against the H₃₇Rv strain of Mycobacterium tuberculosis (MTB) under level-I testing. Remarkably, five compounds 4f, 4h, 4n, 5f and 5m (IC₅₀ = 1.9 µM to 9.8 µM) were found to be most active, with 4f (IC₅₀ = 1.9 µM) indicating highest inhibition of H₃₇Rv. These compounds were further evaluated at level-II testing against the five drug-resistant strains such as isoniazid-resistant strains (INH-R1 and INH-R2), rifampicin-resistant strains (RIF-R1 and RIF-R2) and fluoroquinolone-resistant strain (FQ-R1) of MTB. Interestingly, 4f and 5f emerged as the most potent compounds with IC₅₀ of 3.6 µM and 1.9 µM against RIF-R1 MTB strain, followed by INH-R1 MTB strain with IC₅₀ of 3.5 µM and 3.4 µM, respectively. Against FQ-R1 MTB strain, the lead compounds 4f and 5f displayed excellent inhibition at IC₅₀ 5.9 µM and 4.9 µM, respectively indicating broad-spectrum of activity. Further, molecular docking, ADME pharmacokinetic and molecular dynamics simulations of the compounds were performed against the DNA gyrase B and obtained encouraging results.

N-Substituted pyrrolopyrimidines and purines as p90 ribosomal S6 protein kinase-2 (RSK2) inhibitors

The RSK2 kinase is the downstream effector of the Ras/Raf/MEK/ERK pathway, that is often aberrantly activated in acute myeloid leukemia (AML). Recently, we reported a structure-activity study for BI-D1870, the pan-RSK inhibitor, and identified pteridinones that inhibited cellular RSK2 activity that did not result in concomitant cytotoxicity. In the current study, we developed a series of pyrrolopyrimidines and purines to replace the pteridinone ring of BI-D1870, with a range of N-substituents that extend to the substrate binding site to probe complementary interactions, while retaining the 2,6-difluorophenol-4-amino group to maintain interactions with the hinge domain and the DFG motif. Several compounds inhibited cellular RSK2 activity, and we identified compounds that uncoupled cellular RSK2 inhibition from potent cytotoxicity in the MOLM-13 AML cell line. These N-substituted probes have revealed an opportunity to further examine substituents that extend from the ATP- to the substrate-binding site may confer improved RSK potency and selectivity.

Computational investigation of phytochemicals from Withania somnifera (Indian ginseng/ashwagandha) as plausible inhibitors of GluN2B-containing NMDA receptors

N-Methyl-d-aspartate receptor (NMDAR)-mediated excitotoxicity has been implicated in multi-neurodegenerative diseases. Owing to dearth of efficacy and adverse effects of NMDA receptor antagonists, search for herbal remedies acting like salutary agents is a dynamic expanse of investigation to contest neurodegenerative disease. Withania somnifera (W. somnifera) has been used since antiquity as a nerve tonic and nootropic agents in Ayurveda, an old Indian system of medicine. In the present study, we have explored phytochemicals from Ayurvedic herb W. somnifera as an inhibitor of NMDA receptor-mediated excitotoxicity through allosteric reticence of the GluN1-GluN2B encompassing NMDARs by dint of molecular docking and dynamics studies. Thus, steering and constraining GluN1-GluN2B may be effective in the management of neurodegenerative diseases including Alzheimer's disease. Out of the curtained phytochemicals, chlorogenic acid revealed significant docking scores of -8.856 and -8.645 kcal/mol and free binding energies of -49.84 and -50.67 kcal/mol in Chain AB and Chain CD of NMDARs, respectively. Chlorogenic acid in AB chain forms four hydrogen bonding with Glu110, Arg115, Leu135 and Asp136 amino acid residues and five hydrogen bond with Glu106, Ala107, Ile133, Ile335and Arg155 amino acid residues of CD chain. To further validate the interaction of top scored molecule chlorogenic acid, molecular dynamics study of 100 ns was carried out. It indicated that the protein-ligand complex was stable throughout the simulation period, and minimal backbone fluctuations have ensued in the system. In silico pharmacokinetic predictions of the screened phytochemicals were within the defined range described for human use.Communicated by Ramaswamy H. Sarma.

Optimizing Bedaquiline for cardiotoxicity by structure based virtual screening, DFT analysis and molecular dynamic simulation studies to identify selective MDR-TB inhibitors

Since the last 4 decades, Bedaquiline has been the first drug discovered as a new kind of anti-tubercular agent and received FDA approval in December 2012 to treat pulmonary multi-drug resistance tuberculosis (MDR-TB). It demonstrates excellent efficacy against MDR-TB by effectively inhibiting mycobacterial ATP synthase. In addition to these apparent assets of Bedaquiline, potential disadvantages of Bedaquiline include inhibition of the hERG (human Ether-à-go-related gene; KCNH2), potassium channel (concurrent risk of cardiac toxicity), and risk of phospholipidosis due to its more lipophilic nature. To assist the effective treatment of MDR-TB, highly active Bedaquiline analogs that display a better safety profile are urgently needed. A structure-based virtual screening approach was used to address the toxicity problems associated with Bedaquiline. Among the virtually screened compound, CID 15947587 had significant docking affinity (- 5.636 kcal/mol) and highest binding free energy (ΔG bind - 85.2703 kcal/mol) towards the Mycobacterial ATP synthase enzyme with insignificant cardiotoxicity and lipophilicity. During MD simulation studies (50 ns), the molecule optimizes its conformation to fit better the active receptor site justifying the binding affinity. The obtained results showed that CID15947587 could be a useful template for further optimizing the MDR-TB inhibitor.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00086-x.