Cameroonian medicinal plants as potential candidates of SARS-CoV-2 inhibitors

In silico discovery of potent and selective Janus kinase 3 (JAK3) inhibitors through 3D-QSAR, covalent docking, ADMET analysis, molecular dynamics simulations, and binding free energy of pyrazolopyrimidine derivatives

Rheumatoid arthritis is a prevalent and debilitating chronic disease worldwide. Targeting Janus kinase 3 (JAK3) has emerged as a crucial molecular strategy to treat this condition. In this study, we employed a comprehensive theoretical approach that included 3D-QSAR, covalent docking, ADMET, and molecular dynamics to propose and optimize new anti-JAK3 compounds. We investigated a series of 28 1H-pyrazolo[3.4-d]pyrimidin-4-amino inhibitors and developed a highly accurate 3D-QSAR model using comparative molecular similarity index analysis (COMSIA). The model predicted with Q2 = 0.59, R2 = 0.96, and R2(Pred) = 0.89, was validated using Y-randomization and external validation methods. Our covalent docking studies identified T3 and T5 as highly potent inhibitors of JAK3 compared to the reference ligand 17. Additionally, we evaluated the ADMET properties and drug similarity of our newly developed compounds and reference ligand, providing critical insights for further optimization of anti-JAK3 medications. Furthermore, MM-GBSA analysis showed promising results for the designed compounds. Finally, we validated our docking results using molecular dynamics simulations, which confirmed the stability of hydrogen bonding contacts with key residues required to block JAK3 activity. Our findings offer new chemical scaffolds and insights that could lead to the development of novel and effective JAK3 therapeutic targets for treating rheumatoid arthritis.Communicated by Ramaswamy H. Sarma.

Identification of bioactive compounds of Zanthoxylum armatum as potential inhibitor of pyruvate kinase M2 (PKM2): Computational and virtual screening approaches

PKM2 (Pyruvate kinase M2) is the isoform of pyruvate kinase which is known to catalyse the last step of glycolysis that is responsible for energy production. This specific isoform is known to be highly expressed in certain cancerous conditions. Considering the role of this protein in various cancer conditions, we used PKM2 as a target protein to identify the potential compounds against this target. In this study, we have examined 96 compounds of Zanthoxylum armatum using an array of computational and in silico tools. The compounds were assessed for toxicity then their anticancer potential was predicted. The virtual screening was done with molecular docking followed by a detailed examination using molecular dynamics simulation. The majority of the compounds showed a higher probability of being antineoplastic. Based on toxicity, predicted anticancer potential, binding affinity, and binding site, three compounds (nevadensin, asarinin, and kaempferol) were selected as hit compounds. The binding energy of these compounds with PKM2 ranged from -7.7 to -8.3 kcal/mol and all hit compounds interact at the active site of the protein. The selected hit compounds formed a stable complex with PKM2 when simulated under physiological conditions. The dynamic analysis showed that these compounds remained attached to the active site till the completion of molecular simulation. MM-PBSA analysis showed that nevadensin exhibited a higher affinity towards PKM2 compared to asarinin and kaempferol. These compounds need to be assessed properties in vivo and in vitro to validate their efficacy.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule's binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R2 = 0.8521, Q2LOO = 0.8438, and Q2LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: -8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies.Communicated by Ramaswamy H. Sarma.

Insight into the interaction of isochroman with bovine serum albumin: extensive experimental and computational investigations

The way therapeutic compounds interact with serum protein provides valuable information on their pharmacokinetics, toxicity, effectiveness, and even their structural-related information. Isochroman (IC) is a phytochemical compound obtained from the leaves of Olea europea plant. The derivatives of IC have various pharmacological properties including antidepressants, antihistamines, antiinflammation, anticonvulsants, appetite depressants, etc. The binding of small molecules to bovine serum albumin (BSA) is useful to ensure their efficacy. Thus, in this study, we have found out the binding mode of IC with BSA using several spectroscopic and in silico studies. UV and fluorescence spectroscopy suggested the complex formation between IC and BSA with a binding constant of 103 M-1. IC resulted in fluorescence quenching in BSA through static mechanism. The microenvironmental and conformational changes in BSA were confirmed using synchronous and three-dimensional studies. Site marker experiment revealed the IC binding in site-III of BSA. The influence of vitamins, metals and β-cyclodextrin (β-CD) on binding constant of IC-BSA complex was also examined. Circular dichroism spectra showed that α-helical of BSA decreased upon interaction with IC. Computational and experimental results were complimentary with one another and assisted in determining the binding sites, nature of bonds and amino acids included in the IC-BSA complex formation.Communicated by Ramaswamy H. Sarma.

Imaging metabolic mechanisms and the binding behavior of nutrients/transporters of edible Matricaria flowers VOCs

To promote the consumption of flowers and to utilize the nutritional value of proteins, the efficacy of the beneficial components of flowers has been intensively studied. Anthemis nobilis was used as the study object, and all its volatile components (VOCs) were fingerprinted using headspace solid-phase micro-extraction gas-mass spectrometry (HS-SPME/GC-MS). GC-MS fingerprints of five parts of Anthemis nobilis were established using three proteins, bovine lactoferrin (BLF), bovine lactoglobulin (β-Lg), and human serum albumin (HSA), as nutrient transporters. The interactions between the volatile components from different parts of the mother chrysanthemum plant and the nutrient/transport proteins were investigated. The results of fingerprinting showed that the flavor components were dominated by alkenes. In addition, this study revealed that among the three nutrient transporters, the strongest binding to the adsorbed volatile components was HSA, followed by BLF, and β-Lg was second. In addition, a characteristic molecule, camphene, was screened. Integrated molecular simulation using fluorescence spectroscopy was used to validate the results of the interaction of the nutrient/transport proteins systems with characteristic molecule. The properties of the characteristic molecules such as absorption, distribution, metabolism, excretion and toxicity in vivo were analyzed using ADMET to provide a theoretical basis for the preparation of flower-flavored dairy products.

Phytocompounds screening of Nigella sativa in terms of human cancer by targeting sphingosine kinase-1 and pyruvate kinase-M2: a study based on in silico analysis

Cancer is a multifactorial disease that can cause morbidity and mortality in humans. An altered gene expression in cancer leads to a change in the overall activity of the human cell. Overexpression of cancer protein may give a piece of wide information about the specific type of tumor. Sphingosine kinase-1 (SK-1) is a metabolic enzyme that is mainly overexpressed in several types of cancer and other inflammatory diseases. Similarly, pyruvate kinase-M2 (PK-M2) is an important oncogenic ATP-producing glycolytic enzyme that is upregulated in most cancer cells. The phytocompound of medicinal plants such as Nigella sativa contains a variety of micronutrients that inhibit the proliferation and activity of tumor cells. In this study, the role of phytocompounds in combating cancer was studied against the model kinase proteins, that is, PK-M2 and SK-1. In silico tool like the PASS-Way2Drug server was used to predict the anticancer properties of phytocompounds. Moreover, the CLC-Pred web server provided the cytotoxicity prediction of chemical compounds against several human cancer cell lines. The pharmacokinetics and toxicity profiles were predicted by the SwissADME and pkCSM software. The binding energies were obtained by molecular docking to confirm the intermolecular interaction of selected phytocompounds with proteins. Consequently, molecular dynamics (MD) simulation confirmed the stability, conformational changes, and dynamic behavior of the kinase proteins complexed with the lead phytocompounds, that is, epicatechin, apigenin, and kaempferol.Communicated by Ramaswamy H. Sarma.

Application of in-silico drug discovery techniques to discover a novel hit for target-specific inhibition of SARS-CoV-2 Mpro's revealed allosteric binding with MAO-B receptor: A theoretical study to find a cure for post-covid neurological disorder

Several studies have revealed that SARS-CoV-2 damages brain function and produces significant neurological disability. The SARS-CoV-2 coronavirus, which causes COVID-19, may infect the heart, kidneys, and brain. Recent research suggests that monoamine oxidase B (MAO-B) may be involved in metabolomics variations in delirium-prone individuals and severe SARS-CoV-2 infection. In light of this situation, we have employed a variety of computational to develop suitable QSAR model using PyDescriptor and genetic algorithm-multilinear regression (GA-MLR) models (R2 = 0.800-793, Q2LOO = 0.734-0.727, and so on) on the data set of 106 molecules whose anti-SARS-CoV-2 activity was empirically determined. QSAR models generated follow OECD standards and are predictive. QSAR model descriptors were also observed in x-ray-resolved structures. After developing a QSAR model, we did a QSAR-based virtual screening on an in-house database of 200 compounds and found a potential hit molecule. The new hit's docking score (-8.208 kcal/mol) and PIC50 (7.85 M) demonstrated a significant affinity for SARS-CoV-2's main protease. Based on post-covid neurodegenerative episodes in Alzheimer's and Parkinson's-like disorders and MAO-B's role in neurodegeneration, the initially disclosed hit for the SARS-CoV-2 main protease was repurposed against the MAO-B receptor using receptor-based molecular docking, which yielded a docking score of -12.0 kcal/mol. This shows that the compound that inhibits SARS-CoV-2's primary protease may bind allosterically to the MAO-B receptor. We then did molecular dynamic simulations and MMGBSA tests to confirm molecular docking analyses and quantify binding free energy. The drug-receptor complex was stable during the 150-ns MD simulation. The first computational effort to show in-silico inhibition of SARS-CoV-2 Mpro and allosteric interaction of novel inhibitors with MAO-B in post-covid neurodegenerative symptoms and other disorders. The current study seeks a novel compound that inhibits SAR's COV-2 Mpro and perhaps binds MAO-B allosterically. Thus, this study will enable scientists design a new SARS-CoV-2 Mpro that inhibits the MAO-B receptor to treat post-covid neurological illness.

Inhibitory effect of baicalein against glycation in HSA: an in vitro approach

Hyperglycaemia accelerates the aging process significantly. Diabetes problems can be mitigated by inhibiting glycation. To learn more about glycation and antiglycation mediated by methyl glyoxal and baicalein, we studied human serum albumin as a model protein. A Methylglyoxal (MGO) incubation period of seven days at 37 degrees Celsius induced glycation of Human Serum Albumin.s Hyperchromicity, decreased tryptophan and intrinsic fluorescence, increased AGE-specific fluorescence, and reduced mobility were all seen in glycated human serum albumin (MGO-HSA) in sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE). Fourier transform infrared spectroscopy (FT-IR) and then far ultraviolet dichroism were used to detect secondary and tertiary structural perturbations (CD). The Congo red assay (CR), scanning electron microscopy (SEM), and transmission electron microscopy (TEM) all verified the presence of amyloid-like clumps. Structure (carbonyl groups on ketoamine moieties) (CO), physiological problems including diabetes mellitus, and cardiovascular disease, etc. are linked to the structural and functional changes in glycated HSA, as proven by these studies.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening, ADMET analysis, and molecular dynamics simulation of Moroccan natural compounds as candidates for the SARS-CoV-2 inhibitors

The lack of treatments and vaccines effective against SARS-CoV-2 has forced us to explore natural compounds that could potentially inhibit this virus. Additionally, Morocco is renowned for its rich plant diversity and traditional medicinal uses, which inspires us to leverage our cultural heritage and the abundance of natural resources in our country for therapeutic purposes. In this study, an extensive investigation was conducted to gather a collection of phytoconstituents extracted from Moroccan plants, aiming to evaluate their ability to inhibit the proliferation of the SARS-CoV-2 virus. Molecular docking of the studied compounds was performed at the active sites of the main protease (6lu7) and spike (6m0j) proteins to assess their binding affinity to these target proteins. Compounds exhibiting high affinity to the proteins underwent further evaluation based on Lipinski's rule and ADME-Tox analysis to gain insights into their oral bioavailability and safety. The results revealed that the two compounds demonstrated strong binding affinity to the target proteins, making them potential candidates for oral antiviral drugs against SARS-CoV-2. The molecular dynamics results from this computational analysis supported the overall stability of the resulting complex.

Preparation of nanoformulation of 5-fluorouracil to improve anticancer efficacy: integrated spectroscopic, docking, and MD simulation approaches

Nanoformulations (NFs) can be used as a novel drug delivery system to treat all cancer types. One of the major drawbacks of conventional anticancer drugs is that they have poor specificity and higher toxicity towards normal cells. 5-fluorouracil (5-FU) is a well-studied anticancer drug that has a significant role in various cancers, specifically colorectal cancer therapy. This study was performed to determine the functional groups, particle size, surface charge, heterogeneity, and stability of the NF. The NFs of 5-FU were prepared through the ultrasonication technique by increasing the surfactant (Tween-80) concentrations. Among all three NFs, nanoformulated 5-FU (n5-FU) showed the most effective particle size (10.72 nm) with a zeta potential of (-4.57 mV). The cytotoxicity and apoptosis profiles confirmed that n5-FU enhanced the anticancer effect of the pure drug in HCT-116 cells, as evident from MTT assay, fluorescence microscopy, and FACS analysis. In HCT-116 cells, the IC50 values of pure and n5-FU were obtained as 41.3 μM and 18.8 μM, respectively, indicating that n5-FU was more effective against the cancer cell line. The cellular uptake study was performed to check the intake of NF in cancer cells. However, the microtubule-affinity regulating kinase-4 (MARK-4), a cancer-target protein, was purified to study the inhibition and interaction studies. The inhibition assay confirmed the inhibitory potential of 5-FU against MARK-4 protein. the multi-spectroscopic, molecular docking and MD simulation studies were performed to analyse the conformational changes, binding studies, intermolecular interactions, and stability of MARK-4 protein upon binding 5-FU. This demonstrates that NF can enhance the effectiveness of anticancer drugs.Communicated by Ramaswamy H. Sarma.

In vitro and in silico binding studies of phytochemical isochroman with calf thymus DNA using multi-spectroscopic and computational modelling techniques

A wide range of therapeutic molecules uses deoxyribonucleic acid (DNA) as an intracellular target. The interaction of small molecules to DNA is a key feature in pharmacology and plays a vital role in the development of novel and more efficient drugs with increased selective activity and enhanced therapeutic effectiveness. Isochroman (IC) is a constituent of Olea europea plant, which has been shown to exhibit several beneficial pharmacological activities. At present, its interaction studies using calf thymus DNA (ct-DNA) have not been explained. A set of multi-spectroscopic techniques has been performed to determine the interaction mechanism of isochroman with ct-DNA. Absorption spectra and quenching in fluorescence studies show that isochroman and ct-DNA form a complex. The static mode of quenching was determined by the Stern-Volmer plot. The value of binding constant, Kb = 4.0 × 103 M-1 revealed moderate type of binding. Effects of single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) and ionic strength were studied to examine the isochroman binding to ct-DNA. Potassium iodide (KI) quenching effects and competitive binding studies clearly showed that isochroman binds in the minor groove of ct-DNA. Circular dichroic and DNA melting experiments also confirmed these results. The experimental outputs were further corroborated via in silico computational modelling studies. Lipinski's rule of 5 and SwissADME showed drug-likeness and oral bioavailability scores. Protox ІІ online software predicts oral and organ toxicity.Communicated by Ramaswamy H. Sarma.

Identification of a novel Sorcin isoform with a different C-terminal but intact dimerization property

Sorcin (Sri), a member of penta EF-hand protein family plays a diverse role in maintaining calcium homeostasis, cell cycle and vesicular trafficking. Sri is highly conserved amongst mammals and consists of N-terminal glycine rich domain and C-terminal calcium binding domain that mediates its dimerization and interacts with different compounds. In the present study, with the help of combination of computational and molecular biology techniques, we have identified a novel isoform (Sri-N) in mouse which differs only in the C-terminal domain with that of Sri reported earlier. The novel isoform contains a new last exon that is different from the one present in the reported transcript (Sri). The presence of the novel isoform was further validated in different tissues by RT-PCR and DNA sequencing. The transcript was conceptually translated and subjected to in-silico analysis using different bioinformatics tools. The novel transcript variant encodes for a longer protein isoform without any change in the sub-cellular localization as predicted by PSORT-II online tool. Molecular modelling was performed to compare the structural changes in Sri-N and Sri isoforms. The structural characterization of the novel isoform using MD simulation depicted its overall stability under the physiological conditions. The molecular docking of proteins with various chemotherapeutic drugs revealed that their binding affinity is more for Sri-N as compared to that for the previously reported transcript Sri.

Therapeutic Promises of Plant Metabolites against Monkeypox Virus: An In Silico Study

The monkeypox virus was still spreading in May 2022, with the first case identified in a person with travel ties to Nigeria. Using molecular docking-based techniques, we evaluated the efficiency of different bioactive chemicals obtained from plants against the monkeypox virus. A total of 56 plant compounds were evaluated for antimonekypox capabilities, with the top four candidates having a higher binding affinity than the control. We targeted the monkeypox profilin-like protein, which plays a key role in viral replication and assembly. Among the metabolites, curcumin showed the strongest binding affinity with a value of -37.43 kcal/mol, followed by gedunin (-34.89 kcal/mol), piperine (-34.58 kcal/mol), and coumadin (-34.14 kcal/mol). Based on ADME and toxicity assessments, the top four substances had no negative impacts. Furthermore, four compounds demonstrated resistance to deformability, which was corroborated by normal mode analysis. According to the bioactivity prediction study, the top compound target class was an enzyme, membrane receptor, and oxidoreductase. Furthermore, the study discovered that wortmannin, a gedunin analogue, can behave as an orthopoxvirus. The study found that these bioactive natural drug candidates could potentially work as monkeypox virus inhibitors. We recommended further experimental validation to confirm the promising findings of the study.

Studying the Mechanism of Interaction of Doxofylline with Human Lysozyme: A Biophysical and In Silico Approach

In this study, multiple spectroscopic and computational methods were utilized to investigate the binding mechanism of doxofylline with lysozyme. The in vitro methods were used to obtain the binding kinetics and thermodynamics. UV-vis spectroscopy indicated the formation of complex between doxofylline and lysozyme. The Gibb's free energy and binding constant from UV-vis data was obtained as -7.20 kcal M^-1 and 1.929 × 10⁵ M^-1, respectively. Doxofylline successfully quenched the fluorescence of lysozyme, confirming the formation of complex. The k_q and K_sv values for the quenching of lysozyme's fluorescence by doxofylline were 5.74 × 10¹¹ M^-1 s^-1 and 3.32 × 10³ M^-1, respectively. These values signified a moderate binding affinity between doxofylline and lysozyme. In synchronous spectroscopy, red shifts were observed for indicating the changes in microenvironment of lysozyme following the binding of doxofylline. The secondary structural analysis was determined using circular dichroism (CD) which revealed an increase in % α-helical as a result of doxofylline interaction. The binding affinity and flexibility of lysozyme upon complexation have been revealed via molecular docking and molecular dynamic (MD) simulations, respectively. According to the many parameters of the MD simulation, the lysozyme-doxofylline complex was stable under physiological conditions. All during the simulation time, hydrogen bonds were continuously present. The MM-PBSA binding energy for lysozyme and doxofylline binding was found to be -30.55 kcal mol^-1.

Calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) inhibitors: a novel approach in small molecule discovery

The calcium/calmodulin dependent protein kinase kinase 2 (CAMKK2) plays a key role in regulation of intracellular calcium levels and signaling pathways. It is involved in activation of downstream signaling pathways that regulate various cellular processes. Dysregulation of CAMKK2 activity has been linked to various diseases including cancer, suggesting that CAMKK2 inhibitors might be beneficial in oncological, metabolic and inflammatory indications. The most pressing issues in small molecule discovery are synthesis feasibility, novel chemical structure and desired biological characteristics. To circumvent this constraint, we employed 'DrugspaceX' for rapid lead identification, followed by repositioning seven FDA-approved drugs for CAMKK2 inhibition. Further, first-level transformation (Set1 analogues) was performed in 'DrugspaceX', followed by virtual screening. The t-SNE visualization revealed that the transformations surrounding Rucaparib, Treprostinil and Canagliflozin are more promising for developing CAMKK2 inhibitors. Second, using the top-ranked Set1 analogues, Set2 analogues were generated, and virtual screening revealed the top-ranked five analogues. Among the top five Set2 analogues, DE273038_5 had the lowest docking score of -11.034 kcal/mol and SA score of 2.59, retaining the essential interactions with Hotspot residues LYS194 and VAL270 across 250 ns simulation period. When compared to the other four compounds, the ligand effectiveness score was 0.409, and the number of rotatable penalties was only three. Further, DE273038_5 after two rounds of transformations was discovered to be novel and had not been previously described in other databases. These data suggest that the new candidate DE273038_5 is likely to have inhibitory activity at the CAMKK2 active site, implying potential therapeutic use.Communicated by Ramaswamy H. Sarma.

Analysis and prediction pathways of natural products and their cytotoxicity against HeLa cell line protein using docking, molecular dynamics and ADMET

Natural product such as flavonoids and their derivatives have a discernible capability to inhibit tumor formation and the growth of cancer cell, which have a vital link between diet and chronic disease prevention. Several plants and spices that contain flavonoid derivatives have been used in traditional medicine as disease preventative and therapeutic agents. Therefore, flavonoids could be used as chemotherapeutic drugs, indicating their potential clinical utility in cancer treatment. The purpose of this research was to discover and produce innovative pharmaceuticals from natural sources by introducing structural changes into flavonoids' backbones and changing their structures to improve biological activity and anticancer effects. In the current study, it was expected that the percent unbound values for the 15 compounds in human plasma would be low, ranging between 0.188 and 0.391. However, all compounds have a safe range and are not toxic to the brain. Compounds 2, 10, and 13 were shown to be permeable to the CNS (log PS > -3), but all other compounds had difficulty penetrating the CNS. Furthermore, all compounds had a low total clearance, ranging from 0.038 to 1.216 ml/min/kg, indicating that these compounds have a long half-life. None of the compounds caused skin sensitization (SS), and only compounds 1, 11, and 12 are expected to be AMES-positive, suggesting that the other compounds are not mutagenic. The result of the study showed based on the Drug-likeness and ADMET studies, only 3 compounds, including 3, 4, and 15, have a good pharmacokinetics propriety, the lowest toxicity, and good binding affinity towards Caspase 3 V266APDB (ID: 5I9B) as potential inhibitor candidates for the HeLa cell line, they have a low total clearance property and no AMES mutagenicity or hERG inhibition properties. These compounds (3,4,15) were examined to act as new cytotoxic drug candidates and would have an interest as starting point for designing compounds against the HeLa cell line.Communicated by Ramaswamy H. Sarma.

An in-silico pharmacophore-based molecular docking study to evaluate the inhibitory potentials of novel fungal triterpenoid Astrakurkurone analogues against a hypothetical mutated main protease of SARS-CoV-2 virus

BACKGROUND

The main protease is an important structural protein of SARS-CoV-2, essential for its survivability inside a human host. Considering current vaccines' limitations and the absence of approved therapeutic targets, Mpro may be regarded as the potential candidate drug target. Novel fungal phytocompound Astrakurkurone may be studied as the potential Mpro inhibitor, considering its medicinal properties reported elsewhere.

METHODS

In silico molecular docking was performed with Astrakurkurone and its twenty pharmacophore-based analogues against the native Mpro protein. A hypothetical Mpro was also constructed with seven mutations and targeted by Astrakurkurone and its analogues. Furthermore, multiple parameters such as statistical analysis (Principal Component Analysis), pharmacophore alignment, and drug likeness evaluation were performed to understand the mechanism of protein-ligand molecular interaction. Finally, molecular dynamic simulation was done for the top-ranking ligands to validate the result.

RESULT

We identified twenty Astrakurkurone analogues through pharmacophore screening methodology. Among these twenty compounds, two analogues namely, ZINC89341287 and ZINC12128321 showed the highest inhibitory potentials against native and our hypothetical mutant Mpro, respectively (-7.7 and -7.3 kcal mol-1) when compared with the control drug Telaprevir (-5.9 and -6.0 kcal mol-1). Finally, we observed that functional groups of ligands namely two aromatic and one acceptor groups were responsible for the residual interaction with the target proteins. The molecular dynamic simulation further revealed that these compounds could make a stable complex with their respective protein targets in the near-native physiological condition.

CONCLUSION

To conclude, Astrakurkurone analogues ZINC89341287 and ZINC12128321 can be potential therapeutic agents against the highly infectious SARS-CoV-2 virus.

Combined molecular docking and dynamics simulations studies of natural compounds as potent inhibitors against SARS-CoV-2 main protease

Main protease (Mpro) of SARS-CoV-2 is a key CoV enzyme that plays a pivotal role in mediating viral replication and transcription, making it an attractive drug target for SARS-CoV-2 the new strain of coronavirus. In this study, we evaluated biologically active compounds present in medicinal plants as potential SARS-CoV-2 Mpro inhibitors, using a molecular docking study with Autodock Vina software. Top seven compounds Afzelin, Phloroglucinol, Myricetin-3-O- rutinosid Tricin 7-neohesperidoside, Silybin, Kaempferol and Silychristin among 50 molecules of natural Origin (Algerian Medicinal plants) were selected which had better and significantly low binding energy as compared to the reference molecule with binding affinities of -9.3, -9.3, -9, -8.9, -8.5, 8.3 and -8.3 kcal mol^-1 respectively. Then, we analyzed the ADME properties of the best 7 ligands using the Web server SwissADME. Two of small molecules have been shown to be the ideal candidates for further drug development. Finally, the stability of the both compounds complexed with Mpro was validated through molecular dynamics (MD) simulation, they displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations, moreover, Silybin could form more stable complex with Mpro than Silychristin.Communicated by Ramaswamy H. Sarma.

Combined computational approaches for developing new anti-Alzheimer drug candidates: 3D-QSAR, molecular docking and molecular dynamics studies of liquiritigenin derivatives

Butyrylcholinesterase is an acetylcholine-degrading enzyme involved in the memorization process, which is becoming an interesting target for the symptomatic treatment of Alzheimer's disease. In the present investigation, the structure-activity relationship of a set of Liquiritigenin derivatives recently revealed to be Butyrylcholinesterase inhibitors was studied basing on comparative field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMISA). As a result, performant models with high predictive capability have been developed (CoMFA model: R² = 0.91, Q² = 0.62, R² _pred = 0.85; CoMISA model: R² = 0.92, Q² = 0.59, R² _pred = 0.83) and implemented to design new Liquiritigenin derivatives with improved activity. Besides, the affinity of the designed derivatives towards the active site of Butyrylcholinesterase, was confirmed by molecular docking and molecular dynamics studies. Moreover, they exhibited good pharmacokinetics properties. Accordingly, the outcomes of the present investigations can provide important direction for the development of new anti-Alzheimer's drug candidates.

Unsymmetrical aromatic disulfides as SARS-CoV-2 Mpro inhibitors: Molecular docking, molecular dynamics, and ADME scoring investigations

COVID-19 pandemic caused by very severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) agent is an ongoing major global health concern. The disease has caused more than 452 million affected cases and more than 6 million death worldwide. Hence, there is an urgency to search for possible medications and drug treatments. There are no approved drugs available to treat COVID-19 yet, although several vaccine candidates are already available and some of them are listed for emergency use by the world health organization (WHO). Identifying a potential drug candidate may make a significant contribution to control the expansion of COVID-19. The in vitro biological activity of asymmetric disulfides against coronavirus through the inhibition of SARS-CoV-2 main protease (Mpro) protein was reported. Due to the lack of convincing evidence those asymmetric disulfides have favorable pharmacological properties for the clinical treatment of Coronavirus, in silico evaluation should be performed to assess the potential of these compounds to inhibit the SARS-CoV-2 Mpro.

In this context, we report herein the molecular docking for a series of 40 unsymmetrical aromatic disulfides as SARS-CoV-2 Mpro inhibitor. The optimal binding features of disulfides within the binding pocket of SARS-CoV-2 endoribonuclease protein (Protein Data Bank [PDB]: 6LU7) was described. Studied compounds were ranked for potential effectiveness, and those have shown high molecular docking scores were proposed as novel drug candidates against SARS-CoV-2. Moreover, the outcomes of drug similarity and ADME (Absorption, Distribution, Metabolism, and Excretion) analyses have may have the effectiveness of acting as medicines, and would be of interest as promising starting point for designing compounds against SARS-CoV-2. Finally, the stability of these three compounds in the complex with Mpro was validated through molecular dynamics (MD) simulation, in which they displayed stable trajectory and molecular properties with a consistent interaction profile.

In-Silico Analysis of Phytocompounds of Olea europaea as Potential Anti-Cancer Agents to Target PKM2 Protein

Globally, cancer is the second leading cause of mortality and morbidity. The growth and development of cancer are extremely complex. It is caused by a variety of pathways and involves various types of enzymes. Pyruvate kinase M2 (PKM2) is an isoform of pyruvate kinase, that catalyses the last steps of glycolysis to produce energy. PKM2 is relatively more expressed in tumour cells where it tends to exist in a dimer form. Various medicinal plants are available that contain a variety of micronutrients to combat against different cancers. The phytocompounds of the olive tree (Olea europaea) leaves play an important role in inhibiting the proliferation of several cancers. In this study, the phytocompounds of olive leaf extract (OLE) were studied using various in silico tools, such as pkCSM software to predict ADMET properties and PASS Online software to predict anticancer activity. However, the molecular docking study provided the binding energies and inhibition constant and confirmed the interaction between PKM2 and the ligands. The dynamic behaviour, conformational changes, and stability between PKM2 and the top three hit compounds (Verbascoside (Ver), Rutin (Rut), and Luteolin_7_O_glucoside (Lut)) are studied by MD simulations.

Rational identification of small molecules derived from 9,10-dihydrophenanthrene as potential inhibitors of 3CLpro enzyme for COVID-19 therapy: a computer-aided drug design approach

Small molecules such as 9,10-dihydrophenanthrene derivatives have remarkable activity toward inhibition of SARS-CoV-2 3CL^pro and COVID-19 proliferation, which show a strong correlation between their structures and bioactivity. Therefore, these small compounds could be suitable for clinical pharmaceutical use against COVID-19. The objective of this study was to remodel the structures of 9,10-dihydrophenanthrene derivatives to achieve a powerful biological activity against 3CL^pro and favorable pharmacokinetic properties for drug design and discovery. Therefore, by the use of bioinformatics techniques, we developed robust 3D-QSAR models that are capable of describing the structure–activity relationship for 46 molecules based on 9,10-dihydrophenanthrene derivatives using CoMFA/SE (R² = 0.97, Q² = 0.81, R²_pred = 0.95, ^cR²_p = 0.71) and CoMSIA/SEHDA (R² = 0.94, Q² = 0.76, R²_pred = 0.91, ^cR²_p = 0.65) techniques. Accordingly, 96 lead compounds were generated based on a template molecule that showed the highest observed activity in vitro (T40, pIC₅₀ = 5.81) and predicted their activities and bioavailability in silico. The rational screening outputs of 3D-QSAR, Molecular docking, ADMET, and MM-GBSA led to the identification of 9 novel modeled molecules as potent noncovalent drugs against SARS-CoV-2-3CL^pro. Finally, by molecular dynamics simulations, the stability and structural dynamics of 3CL^pro free and complex (PDB code: 6LU7) were discussed in the presence of samples of 9,10-dihydrophenanthrene derivative in an aqueous environment. Overall, the retrosynthesis of the proposed drug compounds in this study and the evaluation of their bioactivity in vitro and in vivo may be interesting for designing and discovering a new drug effective against COVID-19.

In silico investigation of phytoconstituents from Cameroonian medicinal plants towards COVID-19 treatment

In silico studies performed on the metabolites of four Cameroonian medicinal plants with a view to propose potential molecules to fight against COVID-19 were carried out. At first, molecular docking was performed for a set of 84 selected phytochemicals with SARS-CoV-2 main protease (PDB ID: 6lu7) protein. It was further followed by assessing the pharmacokinetics and pharmacological abilities of 15 compounds, which showed low binding energy values. As the screening criteria for their ADMET properties were performed, only two compounds have shown suitable pharmacological properties for human administration which were shortlisted. Furthermore, the stability of binding of these compounds was assessed by performing molecular dynamics (MD) simulations. Based on further analysis through molecular dynamics simulations and reactivity studies, it was concluded that only the Pycnanthuquinone C (17) and the Pycnanthuquinone A (18) extracted from the Pycnanthus angolensis could be considered as candidate inhibitors for targeted protein. Indeed, we expect that these compounds could show excellent in vitro and in vivo activity against SARS-CoV-2.

3D-QSAR, ADME-Tox, and molecular docking of semisynthetic triterpene derivatives as antibacterial and insecticide agents

In the present work, 27 triterpene derivatives have been subjected to 3D-QSAR, ADME-Tox, and molecular docking for their insecticidal activity. The selected derivatives are previously semi-synthesized based on compounds obtained from Euphorbia resinifera and Euphorbia officinarum latex. The in silico studies were used to predict and to evaluate the antibacterial and insecticidal properties of the 3D structure of triterpene derivatives. The 3D-QSAR models are developed using CoMFA and CoMSIA techniques, and they have showed excellent statistical results (R² = 0.99; Q² = 0.672; R²_pred = 0.91 for CoMFA and R² = 0.97; Q² = 0.61; R²_pred = 0.94 for CoMSIA). The results indicate that the built models are able to describe the relationship between the structure of triterpene derivatives and the pLD₅₀ bioactivity. Based on contour maps obtained from CoMFA and CoMSIA models, 38 new molecules are designed and their pLD₅₀ activities are predicted. The drug-like and ADME-Tox properties of the molecule designed are examined and led to the selection of four molecules (55, 56, 59, 64) as promising antibacterial and insecticidal agents. Compounds 55, 56, 59, and 64 are able to inhibit the MurE (PDB code: 1E8C) and EcR (PDB code: 1R20) proteins involved in the process of antibacterial and insecticidal activities. This hypothesis is confirmed by the implementation of a molecular docking test. This test predicted the most important referential interactions that occur between the structure of triterpene derivatives and the targeted receptors. Among the four docked molecules, three molecules (55, 56, and 59) showed greater stability than the reference molecule 16 inside the MurE and EcR receptors pocket. Therefore, the structure of the three new triterpene derivatives can be adopted as reference for the synthesis of antibacterial drugs and also in the development of insecticides.

In silico screening and in vitro validation of phytocompounds as multidrug efflux pump inhibitor against E. coli

Multiple drug resistance (MDR) in bacteria has increased globally in recent times. This has reduced the efficacy of antibiotics and increasing the rate of therapeutic failure. Targeting efflux pump by natural and synthetic compounds is one of the strategies to develop an ideal broad-spectrum resistance-modifying agent. Very few inhibitors of AcrB from natural sources have been reported till date. In the current study, 19 phytocompounds were screened for efflux pump inhibitory activity against AcrB protein of E. coli TG1 using molecular docking studies. The molecular dynamics simulation provided stability the protein (AcrB) and its complex with chlorogenic acid under physiological conditions. Moreover, the detailed molecular insights of the binding were also explored. The Lipinski rule of 5 and the drug-likeness prediction was determined using Swiss ADME server, while toxicity prediction was done using admetSAR and PROTOX-II webservers. Chlorogenic acid showed the highest binding affinity (-9.1 kcal mol^-1) with AcrB protein among all screened phytocompounds. Consequently, all the phytocompounds that accede to Lipinski's rule, demonstrated a high LD50 value indicating that they are non-toxic except the phytocompound reserpine. Chlorogenic acid and capsaicin are filtered out based on the synergy with tetracycline having FIC index of 0.25 and 0.28. The percentage increase of EtBr fluorescence by chlorogenic acid was 36.6% followed by piperine (24.2%). Chlorogenic acid may be a promising efflux pump inhibitor that might be employed in combination therapy with tetracycline against E. coli, based on the above relationship between in silico screening and in vitro positive efflux inhibitory activity.Communicated by Ramaswamy H. Sarma.

Antioxidant, antibacterial, and antimutagenic activity of Piper nigrum seeds extracts

Piper nigrum is a widely used plant in traditional remedies and known for its numerous biological properties. However, fraction-based antioxidant activity and their antimutagenic potential are not yet fully investigated. Different extracts of the seeds P. nigrum were obtained by sequential extraction in different solvents. All extracts were evaluated for antibacterial and antioxidant activities using different methods. The most active fraction was analyzed for antimutagenic activity using the Ames Salmonella test. The antibacterial activity against methicillin-resistant Staphylococcus aureus (MRSA) was found to be more prominent compared to ESβL producing Klebsiella pneumoniae isolates. The MIC values were found to be lower against MRSA than K. pneumoniae. The extract showing highest antioxidant activity (methanol extract) was further tested for antimutagenic activity both against direct and indirect-acting mutagens. A varying level of antimutagenic activity was shown by methanol extract at highest tested concentration (200 µg/plate). Alkaloids, phenols, and flavonoids were detected as major class of compounds in methanol extract. Gas chromatography-mass spectrometry (GC–MS) analysis showed the presence of various phytocompounds. Based on molecular docking of two major active phytocompounds (piperine and copaene), they were found to interact at the minor groove of DNA. Molecular dynamics (MD) simulation revealed that both the ligands were quite stable with DNA under physiological conditions. The ability of phytocompounds to interact with DNA might be reducing the interaction of mutagens and could be one of the possible mechanism of anti-mutagenic activity of P. nigrum extract. This study highlights the antioxidant and antimutagenic potential of Piper nigrum. The role of phytocompounds present in the bioactive extract is needed to be explored further for herbal drug research.

QSAR, molecular docking and ADMET properties in silico studies of novel 4,5,6,7-tetrahydrobenzo[D]-thiazol-2-Yl derivatives derived from dimedone as potent anti-tumor agents through inhibition of C-Met receptor tyrosine kinase

A quantitative structure-activity relationship (QSAR) study is performed on 48 novel 4,5,6,7-tetrahydrobenzo[D]-thiazol-2 derivatives as anticancer agents capable of inhibiting c-Met receptor tyrosine kinase. The present study is conducted using multiple linear regression, multiple nonlinear regression and artificial neural networks. Three QSAR models are developed after partitioning the database into two sets (training and test) via the k-means method. The obtained values of the correlation coefficients by the three developed QSAR models are 0.90, 0.91 and 0.92, respectively. The resulting models are validated by using the external validation, leave-one-out cross-validation, Y-randomization test, and applicability domain methods. Moreover, we evaluated the drug-likeness properties of seven selected molecules based on their observed high activity to inhibit the c-Met receptor. The results of the evaluation showed that three of the seven compounds present drug-like characteristics. In order to identify the important active sites for the inhibition of the c-Met receptor responsible for the development of cancer cell lines, the crystallized form of the Crizotinib-c-Met complex (PDB code: 2WGJ) is used. These sites are used as references in the molecular docking test of the three selected molecules to identify the most suitable molecule for use as a new c-Met inhibitor. A comparative study is conducted based on the evaluation of the predicted properties of ADMET in silico between the candidate molecule and the Crizotinib inhibitor. The comparison results show that the selected molecule can be used as new anticancer drug candidates.