Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of ACEIs Against SARS-CoV-2 Targeting the hACE2 Receptor

Computational analysis of zoanthamine alkaloids from Zoanthus sp. as potential DKK1 and GSK-3β inhibitors for osteoporosis therapy via Wnt signaling

Marine invertebrates are a rich source of structurally diverse secondary metabolites with broad biological activities, making them valuable for drug discovery. The genus Zoanthus is particularly noteworthy, producing numerous bioactive alkaloids, including the zoanthamines, which show promise in treating osteoporosis. Osteoporosis, a debilitating bone disease characterized by reduced bone mineral density and increased fracture risk, is linked to Wnt signaling pathway dysregulation. This highly conserved pathway maintains tissue homeostasis and is crucial for neurogenesis, synapse formation, and bone development. Dickkopf-1 (DKK1) and glycogen synthase kinase-3β (GSK-3β), key Wnt pathway regulators, are established therapeutic targets for osteoporosis. This study employed an integrated computational approach-combining molecular docking, extensive molecular dynamics (MD) simulations, and density functional theory (DFT) calculations-to assess the inhibitory potential of 69 zoanthamine-type alkaloids against DKK1 and GSK-3β. MD simulations, analyzing root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration, and free energy landscape, provided insights into protein-ligand complex stability and key interactions. Binding free energies were calculated using the MM-PBSA method combined with interaction entropy. DFT calculations further elucidated the electronic structure and reactivity of the most promising inhibitors (3α-hydroxyzoanthenamine, epioxyzoanthamine, 7α-hydroxykuroshine E, and norzoanthamine), which exhibited favorable binding interactions with key residues in target proteins. This integrative approach demonstrates the power of computational methods in drug discovery, highlighting the potential of zoanthamine alkaloids as lead compounds for innovative osteoporosis therapies.

Integrative approach to decipher pharmacological mechanism of Cinnamomum zeylanicum essential oil in prostate cancer

Prostate cancer has garnered much importance in recent years due to its rising incidence and mortality among men worldwide. The ineffectiveness of existing therapies and adverse events associated with conventional treatment have led patients to turn towards traditional medicine for the management of prostate cancer. Cinnamomum zeylanicum bark essential oil (CZEO) possesses promising anticancer properties, yet the exact mechanism of action of CZEO for the management of prostate cancer remains unclear. Therefore, the current study tried to elucidate the bioactive components and key potential targets through which CZEO may exert its anticancer effect for treating prostate cancer. Fifty-nine constituents were identified by GC-MS, of which 52 were drug-like constituents. A total of 2847 targets related to CZEO and 2283 targets related to prostate cancer were obtained from public databases and the GEO dataset. Twenty-three overlapping targets exist between CZEO and disease targets. Compound-disease-target network analysis revealed camphor, eugenol, methyl eugenol, trans farnesyl acetate and nerol as the core bioactive ingredients of CZEO. The topological screening of the PPI network revealed BCL2, TNF, NFKBIA, CREBBP and IL7R as potential hub targets. These hub targets were validated based on mRNA expression level, pathological stages, overall survival, immune infiltrate and genetic alteration analysis in prostate adenocarcinoma and normal patients. KEGG enrichment analysis proposed that CZEO exhibits its anticancer effect mainly by modulating the PI3-AKT and MAPK signalling pathway. Moreover, molecular docking and dynamics simulation studies revealed a good binding affinity of these core compounds with TNF, NFKBIA and BCL2. CZEO exhibited a remarkable anti-proliferative effect against PC-3 cells with an IC50 value of 13.56 µg/mL. CZEO promoted apoptosis and cell cycle arrest in the G2/M phase in PC-3 cells. CZEO-induced apoptosis was due to loss of mitochondrial membrane potential, increase in reactive oxygen species levels and activation of caspases (caspase 3, caspase 8 and caspase 9). RT-qPCR analysis revealed that CZEO modulated the mRNA expression level of hub genes (BCL2, TNF, NFKBIA, CREBBP, and IL7R, caspase 3, caspase 8 and caspase 9). The present study provides a mechanistic approach of Cinnamomum zeylanicum essential oil against prostate cancer.

Phytochemical, In Silico, In Vitro, and In Vivo Research on Piptadeniastrum africanum (Fabaceae) Unveiling Anti-Stereotypic, Anxiolytic, and Analgesic Effects in a Sodium Valproate-Induced Autistic Disorders Model

OBJECTIVE

Individuals with autistic spectrum disorders (ASD) primarily exhibit deficits in communication and social interaction, along with repetitive behaviors and restricted interests. This disorder is often associated with anxiety, nociceptive disorders, and pain. While medical treatment generally focuses on treating the symptoms rather than addressing the underlying causes, traditional medicine is sometimes used as an alternative. Piptadeniastrum africanum is used in Cameroonian medicinal folks to treat cognitive disorders. However, its effects and mechanisms of action regarding the inhibition of ASD-like symptoms remain unclear. The primary goal of the present study was to evaluate the anxiolytic and analgesic effects of the water extract of P. africanum on autistic triad induced in rats by sodium valproate.

MATERIAL AND METHODS

The study investigated the secondary metabolites in P. africanum extract using UHPLC-MS. DPPH, ABTS, and FRAP tests were performed to assess the extract's ability to neutralize free radicals. Molecular docking was utilized to evaluate the extract's binding to various receptors. For the experimental study, 33 pregnant female rats were divided into two groups after pregnancy was confirmed. One group was given distilled water orally at 10 mL/kg, while the other group received sodium valproate at 800 mg/kg on gestation days 11, 12, and 13. When the male offspring reached 3 weeks old, they were evaluated for anxiety, social interaction, and pain sensitivity, with those displaying any disorders selected for further study. The remaining rats were split into six groups of five and treated with either a vehicle, bumetanide, or P. africanum extract at 190 and 760 mg/kg. Behavioral assessments focusing on sociability, anxiety, and pain sensitivity were conducted on days 28 and 37 after weaning. In the end, biochemical markers related to GABA metabolism, serotonin levels, and oxidative status were analyzed in the cerebellum, prefrontal cortex, hippocampus, and amygdala alongside histopathological analyses in the brain.

RESULTS

UHPLC-MS allows us to identify several compounds. They bind to H3R (7F61) and HDAC2 through conventional hydrogen bonding. Findings showed that prenatal administration of sodium valproate induced in male offspring a deficit in social interaction (p < 0.001), anxiety disorders (p < 0.001), hypersensitivity to pain (p < 0.001), increased GABA and serotonin concentration (p < 0.001), disturbed oxidative status (p < 0.001), and neuronal loss (p < 0.001) as well as neuronal disorganization in the hippocampus, cerebellum and amygdala in young rats compared to neurotypical animals. P. africanum extract at doses used, like bumetanide, corrected these disorders and protected against neuronal loss. These results suggest that the extract has anxiolytic and anti-nociceptive effects. It has been found that the positive effects can be achieved by restoring GABAergic and serotonergic neurotransmission, coupled with antioxidant and neuromodulatory activity.

CONCLUSION

The current findings support that P. africanum induces anxiolytic and analgesic effects in a sodium valproate-induced autistic disorders model.

Molecular dynamics of SARS-CoV-2 omicron variants from Philippine isolates against hesperidin as spike protein inhibitor

SARS-CoV-2 remains a global threat with new sublineages posing challenges, particularly in the Philippines. Hesperidin (HD) is being studied as a potential prophylactic for COVID-19. However, the virus's rapid evolution could alter how HD binds to it, affecting its effectiveness. Here, we study the mutation-induced variabilities of HD dynamics and their effects on molecular energetics in SARS-CoV-2 spike receptor complex systems. We considered eight different point mutations present in the Omicron variant. Root-mean-square deviation and binding energy analysis showed that S477N and Omicron did not eject HD throughout the simulation. Hydrogen bond distribution analysis highlighted the involvement of hydrogen bonding in mutant-HD stabilization, especially for S477N and Omicron. Root-mean-square fluctuation analysis revealed evidence of Y505H destabilization on complex systems, while distal-end loop mutations increased loop flexibility for all models bearing the three mutations. Per-residue energy decomposition demonstrated that Q493R substitution increased HD interaction. Free energy landscape and essential dynamics through principal component analysis provided insights into the conformational subspace distribution of mutant model molecular dynamics trajectories. In conclusion, significant mutations contributed to the HD interaction in different ways. S477N has shown significant binding contributions through favorable ligand interaction, while other mutations contribute via conformational modifications, increased affinity due to sidechain mutations, and increased loop flexibility.

Impact of African-Specific ACE2 Polymorphisms on Omicron BA.4/5 RBD Binding and Allosteric Communication Within the ACE2-RBD Protein Complex

Severe acute respiratory symptom coronavirus 2 (SARS-CoV-2) infection occurs via the attachment of the spike (S) protein's receptor binding domain (RBD) to human ACE2 (hACE2). Natural polymorphisms in hACE2, particularly at the interface, may alter RBD-hACE2 interactions, potentially affecting viral infectivity across populations. This study identified the effects of six naturally occurring hACE2 polymorphisms with high allele frequency in the African population (S19P, K26R, M82I, K341R, N546D and D597Q) on the interaction with the S protein RBD of the BA.4/5 Omicron sub-lineage through post-molecular dynamics (MD), inter-protein interaction and dynamic residue network (DRN) analyses. Inter-protein interaction analysis suggested that the K26R variation, with the highest interactions, aligns with reports of enhanced RBD binding and increased SARS-CoV-2 susceptibility. Conversely, S19P, showing the fewest interactions and largest inter-protein distances, agrees with studies indicating it hinders RBD binding. The hACE2 M82I substitution destabilized RBD-hACE2 interactions, reducing contact frequency from 92 (WT) to 27. The K341R hACE2 variant, located distally, had allosteric effects that increased RBD-hACE2 contacts compared to WThACE2. This polymorphism has been linked to enhanced affinity for Alpha, Beta and Delta lineages. DRN analyses revealed that hACE2 polymorphisms may alter the interaction networks, especially in key residues involved in enzyme activity and RBD binding. Notably, S19P may weaken hACE2-RBD interactions, while M82I showed reduced centrality of zinc and chloride-coordinating residues, hinting at impaired communication pathways. Overall, our findings show that hACE2 polymorphisms affect S BA.4/5 RBD stability and modulate spike RBD-hACE2 interactions, potentially influencing SARS-CoV-2 infectivity-key insights for vaccine and therapeutic development.

High-throughput sequencing reveals the existence of a novel fucose binding lectin gene from the whole gut and fat body tissues of the grub of darkling beetle, Zophobas morio

Fucose binding lectins are a taxonomically, evolutionarily, and ecologically relevant class of lectins that have been identified among a wide range of taxa. These lectins possess a distinctive ability to identify and bind fucose-containing glycans, thereby lending them the ability to function as immune cell mitogens, diagnostic markers and anti-cancer agents. Our preliminary analysis revealed the existence of a single D-fucose binding lectin from the grub serum of the darkling beetle, Z. morio. Here, whole transcriptome analysis using entire gut and fat body tissues of grub of Z. morio with Illumina NovaSeq6000 sequencing platform revealed the existence of a novel fucose binding lectin (ZmFBL) displaying significant similarity to Diaphorina citri fucolectin-3 and possessing fucolectin-tachylectin4-pentraxin domain. Molecular docking between the structure of the predicted FBL transcript from Z. morio and D-fucose demonstrated the lowest binding energy of -5.4 kcal/mol. The MD simulation and MM/PBSA analysis furnished insights into the binding stability of D-fucose and the ZmFBL. A similar transcript was also identified from the fat body transcriptome. A 2.49-fold increase in the ZmFBL expression was observed in the fat body than the whole gut as evidenced from the relative quantitation using RT-PCR.

Self-assembly of differently charged trimesic based lithocholic amphiphiles and their assessment on antimicrobial and biostimulant properties

Biosurfactant based biostimulants plays a vital role in agriculture filed by enhancing the soil quality, promote plant growth, and eliminate plant pathogens, and increasing nutrient uptake. This manuscript describes the synthesis of trimesic based lithocholic ester functionalized amphiphiles (TMLCEA) with oppositely charged head groups using thiol-yne click chemistry, which is an effective and simple approach. The trimesic based lithocholic ester functionalized zwitterionic penicillamine (TMLCEPA), cationic cysteamine·HCl (TMLCECy), and anionic thiomalic acid (TMLCETM) exhibited hierarchically self-assembled microstructures from below to above the CMC. In below the CMC, TMLCEPA, TMLCECy, and TMLCETM showed a bundle of petals, flower-like morphology, and grass seed-like patterns respectively. The morphology of self-assembly was studied by FE-SEM, DLS, OPM, contact angle, and zeta potential measurements. Among these amphiphiles, TMLCECy exhibited potential antimicrobial activity at above the CMC. The biostimulant effect of different concentration of TMLCEA treated with maize and green gram seeds were evaluated under in vitro condition, wherein TMLCECy showed improved seed germination and seedling parameters at 750 µL/mL as compared to TMLCEPA, TMLCETM and untreated amphiphiles as control. Molecular docking and molecular dynamic simulations show that TMLCEPA and TMLCETM showed higher binding affinity for dengue methyltransferase protein. The result of the present study opens up new avenues for bile acid-based amphiphiles as bio-based and cost-effective biostimulants for sustainable agriculture.

Bidirectional approach of Punica granatum natural compounds: reduction in lung cancer and SARS-CoV-2 propagation

The spreading of COVID-19 has posed a risk to global health, especially for lung cancer patients. An investigation is needed to overcome the challenges of COVID-19 pathophysiology and lung cancer disease. This study was designed to evaluate the phytoconstituents in Punica granatum peel (PGP), its anti-lung cancer activity, and in silico evaluation for antiviral potential. GC-MS technique was used to detect the phytoconstituents. Cytotoxicity was analyzed using MTT dye, followed by apoptosis, ROS generation, and cell cycle phase detection in human lung cancer cells (A549). The glide module of Maestro software was used to investigate the molecular-docking interaction of the constituents against main protease (Mpro) and papain-like protease (PLpro) of SARS-CoV-2. GROMACS 2023.2 was utilized to evaluate the complex stability. A total of nineteen phytocomponents were detected in the PGP extract through GC-MS analysis. PGP has shown a potential to reduce lung cancer cell proliferation while evading normal cell death. PGP induced apoptosis by arresting cells in the G0/G1 phase and generating ROS. A total of six and eight phytocomponents had a high affinity for PLpro and Mpro proteins, respectively. The top docked complex, ethyl 5-oxo-2-pyrrolidinecarboxylate, with PLpro and Mpro proteins, showed likely stable interaction throughout 100 ns simulation. This finding raises the possibility of top-eight hits (docking score ≥ -1.0 kcal/mol) preventing SARS-CoV-2 severity. The phytoconstituents exhibited orally active drugs with no more than one violation and drug-likeness activity. The PGP phytoconstituents are suggested to be dual agents for lung cancer and SARS-CoV-2 pathogenesis.

Exploring the impact of deleterious missense nonsynonymous single nucleotide polymorphisms in the DRD4 gene using computational approaches

Dopamine receptor D4 (DRD4) plays a vital role in regulating various physiological functions, including attention, impulse control, and sleep, as well as being associated with various neurological diseases, including attention deficit hyperactivity disorder, novelty seeking, and so on. However, a comprehensive analysis of harmful nonsynonymous single nucleotide polymorphisms (nsSNPs) of the DRD4 gene and their effects remains unexplored. The aim of this study is to uncover novel damaging missense nsSNPs and their structural and functional effects on the DRD4 receptor. From the dbSNP database, we found 677 nsSNPs, and then we analyzed their functional consequences, disease associations, and effects on protein stability with fifteen in silico tools. Five variants, including L65ICL1P (rs1459150721), V1163.33D (rs761875546), I1293.46S (rs751467198), I1564.46T (rs757732258), and F2015.47S (rs199609858), were identified as the most deleterious mutations that were also present in the conserved region and showed lower interactions with neighboring residues. To comprehensively understand their impact, we docked agonist dopamine and antagonist nemonapride at the binding site of the receptor, followed by 200 ns molecular dynamics simulations. We identified the V116D and I129S mutations as the most damaging, followed by F201S in the dopamine-bound states. Both the V116D and I129S variants demonstrated significantly high RMSD, Rg, and SASA, and low thermodynamic stability. The F201S-dopamine complex exhibited lower compactness and higher motions, along with a significant loss of hydrogen bonds and active site interactions. By contrast, while interacting with nemonapride, the impact of the I156T and L65P mutations was highly deleterious; both showed lower stability, higher flexibility, and higher motions. Additionally, nemonapride significantly lost interactions with the active site, notably in the I156T variant. We also found the V116D-nemonapride complex as structurally damaging; however, the interaction patterns of nemonapride were less altered in the MMPBSA analysis. Overall, this study revealed five novel deleterious variants along with a comprehensive understanding of their effect in the presence of an agonist and antagonist, which could be helpful for understanding disease susceptibility, precision medicine, and developing potential drugs.

Evaluation of hesperidin as a potential larvicide against Culex pipiens with computational prediction of its mode of action via molecular docking

Hesperidin, a natural flavanone glycoside predominantly found in citrus fruits, has gained attention for its wide-ranging biological activities, including potential insecticidal properties. Culex pipiens, commonly known as the northern house mosquito, is a major vector of several human pathogens, such as the West Nile virus and filariasis, making it a key target in the fight against vector-borne diseases. In this study, we evaluated the larvicidal activity of Hesperidin against Culex pipiens larvae, assessing its potential as an alternative to chemical insecticides. Hesperidin demonstrated potent larvicidal effects, with a lethal concentration 50 (LC50) of 570.3 ± 0.04 µg/mL, outperforming the conventional insecticide Chlorpyrifos 588.3 ± 0.28 µg/mL in efficacy. Molecular docking simulations revealed a strong binding affinity between Hesperidin and crucial neuroreceptors in Culex pipiens, particularly acetylcholinesterase (AChE), a key enzyme involved in nerve signal transmission. The interaction between Hesperidin's hydroxyl groups and the AChE enzyme's active site suggests that AChE inhibition is the primary mechanism driving Hesperidin's insecticidal action. These findings position Hesperidin as a promising, environmentally friendly alternative to synthetic insecticides. However, further research is needed to assess its toxicity to non-target organisms and optimize its formulation for broader application in mosquito control.

Exploring the antifungal potential of Cannabis sativa-derived stilbenoids and cannabinoids against novel targets through in silico protein interaction profiling

Cannabinoid and stilbenoid compounds derived from Cannabis sativa were screened against eight specific fungal protein targets to identify potential antifungal agents. The proteins investigated included Glycosylphosphatidylinositol (GPI), Enolase, Mannitol-2-dehydrogenase, GMP synthase, Dihydroorotate dehydrogenase (DHODH), Heat shock protein 90 homolog (Hsp90), Chitin Synthase 2 (CaChs2), and Mannitol-1-phosphate 5-dehydrogenase (M1P5DH), all of which play crucial roles in fungal survival and pathogenicity. This research evaluates the binding affinities and interaction profiles of selected cannabinoids and stilbenoids with these eight proteins using molecular docking and molecular dynamics simulations. The ligands with the highest binding affinities were identified, and their pharmacokinetic profiles were analyzed using ADMET analysis. The results indicate that GMP synthase exhibited the highest binding affinity with Cannabistilbene I (-9.1 kcal/mol), suggesting hydrophobic solid interactions and multiple hydrogen bonds. Similarly, Chitin Synthase 2 demonstrated significant binding with Cannabistilbene I (-9.1 kcal/mol). In contrast, ligands such as Cannabinolic acid and 8-hydroxycannabinolic acid exhibited moderate binding affinities, underscoring the variability in interaction strengths among different proteins. Despite promising in silico results, experimental validation is necessary to confirm therapeutic potential. This research lays a crucial foundation for future studies, emphasizing the importance of evaluating binding affinities, pharmacokinetic properties, and multi-target interactions to identify promising antifungal agents.

Comprehensive In Silico Analysis of Uncaria Tomentosa Extract: Chemical Profiling, Antioxidant Assessment, and CLASP Protein Interaction for Drug Design in Neurodegenerative Diseases

BACKGROUND

Uncaria tomentosa is a traditional medicinal herb renowned for its anti-inflammatory, antioxidant, and immune-enhancing properties. In the realm of neurodegenerative diseases (NDDS), CLASP proteins, responsible for regulating microtubule dynamics in neurons, have emerged as critical players. Dysregulation of CLASP proteins is associated with NDDS, such as Alzheimer's, Parkinson's, and Huntington's diseases. Consequently, comprehending the role of CLASP proteins in NDDS holds promise for the development of innovative therapeutic interventions.

OBJECTIVES

The objectives of the research were to identify phytoconstituents in the hydroalcoholic extract of Uncaria tomentosa (HEUT), to evaluate its antioxidant potential through in vitro free radical scavenging assays and to explore its potential interaction with CLASP using in silico molecular docking studies.

METHODS

HPLC and LC-MS techniques were used to identify and quantify phytochemicals in HEUT. The antioxidant potential was assessed through DPPH, ferric reducing antioxidant power (FRAP), nitric oxide (NO) and superoxide (SO) free radical scavenging methods. Interactions between conventional quinovic acid, chlorogenic acid, epicatechin, corynoxeine, rhynchophylline and syringic acid and CLASP were studied through in silico molecular docking using Auto Dock 4.2.

RESULTS

The HEUT extract demonstrated the highest concentration of quinovic acid derivatives. HEUT exhibited strong free radical-scavenging activity with IC50 values of 0.113 μg/ml (DPPH) and 9.51 μM (FRAP). It also suppressed NO production by 47.1 ± 0.37% at 40 μg/ml and inhibited 77.3 ± 0.69% of SO generation. Additionally, molecular docking revealed the potential interaction of quinovic acid with CLASP for NDDS.

CONCLUSION

The strong antioxidant potential of HEUT and the interaction of quinovic acid with CLASP protein suggest a promising role in treating NDDS linked to CLASP protein dysregulation.

A comparative analysis of Marburg virus-infected bat and human models from public high-throughput sequencing data

Marburg virus (MARV) disease (MVD) is an uncommon yet serious viral hemorrhagic fever that impacts humans and non-human primates. In humans, infection by the MARV is marked by rapid onset, high transmissibility, and elevated mortality rates, presenting considerable obstacles to the development of vaccines and treatments. Bats, particularly Rousettus aegyptiacus, are suspected to be natural hosts of MARV. Previous research reported asymptomatic MARV infection in bats, in stark contrast to the severe responses observed in humans and other primates. Recent MARV outbreaks highlight significant public health concerns, underscoring the need for gene expression studies during MARV progression. To investigate this, we employed two models from the Gene Expression Omnibus, including kidney cells from Rousettus aegyptiacus and primary proximal tubular cells from Homo sapiens. These models were chosen to identify changes in gene expression profiles and to examine co-regulated genes and pathways involved in MARV disease progression. Our analysis of differentially expressed genes (DEGs) revealed that these genes are mainly associated with pathways related to the complement system, innate immune response via interferons (IFNs), Wnt/β-catenin signaling, and Hedgehog signaling, which played crucial roles in MARV infection across both models. Furthermore, we also identified several potential compounds that may be useful against MARV infection. These findings offer valuable insights into the mechanisms underlying MARV's pathophysiology and suggest potential strategies for preventing transmission, managing post-infection effects, and developing future vaccines.

Neuroprotective Efficacy of Phytoconstituents of Methanolic Shoots Extract of Calligonum polygonoides L. in Hypercholesterolemia-associated Neurodegenerations

BACKGROUND

Small molecule phytocompounds can potentially ameliorate degenerative changes in cerebral tissues. Thus, the current study aimed to evaluate the neuroprotective efficacy of phytocompounds of methanolic shoots extract of Calligonum polygonoides L. (MSECP) in hypercholesterolemia-associated neurodegenerations.

METHODS

Phytochemical screening of the extract was made by LCMS/MS and validated by a repository of the chemical library. The hypercholesterolemia was induced through the intraperitoneal administration of poloxamer-407 with a high-fat diet. The in silico assessments were accomplished by following the molecular docking, ADME and molecular dynamics. MMPBSA and PCA (Principal Component Analysis) analyzed the molecular dynamics simulations. Consequently, in-vivo studies were examined by lipid metabolism, free radical scavenging capabilities and histopathology of brain tissues (cortex and hippocampus).

RESULTS

22 leading phytocompounds were exhibited in the test extract, as revealed by LCMS/ MS scrutiny. Molecular docking evaluated significant interactions of apigenin triacetate with target proteins (HMGCR (HMG-CoA reductase), (AChE-Acetylcholinesterase) and (BuChE- Butyrylcholinesterase). Molecular dynamics examined the interactions through assessments of the radius of gyration, RSMD, RSMF and SASA at 100 ns, which were further analyzed by MMPBSA (Molecular Mechanics Poisson-Boltzmann) and PCA (Principal Component Analysis). Accordingly, the treatment of test extract caused significant alterations in lipid profile, dyslipidemia indices, antioxidant levels and histopathology of brain tissues.

CONCLUSION

It can be concluded that apigenin triacetate is a potent phytoconstituent of MSEPC and can interact with HMGCR, AChE, and BuChE, which resulted in improved hypercholesterolemia along with neuroprotective ameliorations in the cortex and hippocampus.

Computational-driven discovery of AI-2 quorum sensing inhibitor targeting the 5'- methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) to combat drug-resistant Helicobacter pylori

MTAN is an attainable therapeutic target for H. pylori because it may minimize virulence production, limit resistance, and impair quorum sensing without affecting gut flora. Here, 457 compounds with anti-H. pylori activity were methodically analyzed, revealing a diverse array of chemical classes and unique compounds. Molecular docking studies identified three potential compounds with high binding affinities, Dehydrocostus lactone, keramamine B, and ZINC00013531409, each having binding affinity of -7.9, -9.2, and -8.3 kcal/mol, respectively. Molecular dynamics simulations of the ZINC00013531409-MTAN interactions in comparison with Apo-MTAN demonstrated stability and interactions of 300 ns, with key residues involved in protein-ligand binding illuminated. Analysis of hydrogen bonds (Ile52, Met174, and Arg194) and secondary structure variations further elucidated the binding interactions and conformational changes within the complex. Binding free energy calculations shed light on the energetics and interactions governing the complex formation of the ZINC00013531409-MTAN complex. PCA elucidated the dominant modes of motion, along with FEL revealed the energetically favorable states and then DCCM shed light on the correlated motions between residues. Overall, this study offers a detailed computational evaluation of ZINC00013531409 with anti-H. pylori activity, highlighting toxicity profile, conformational stability, and binding interactions, providing a foundation for further drug development efforts toward bacterial resistance.

Design of new Mcl-1 inhibitors for cancer using fragments hybridization, molecular docking, and molecular dynamics studies

Apoptosis is a critical process that regulates cell survival and death and plays an essential role in cancer development. The Bcl-2 protein family, including myeloid leukemia 1 (Mcl-1), is a key regulator of the intrinsic apoptosis pathway, and its overexpression in many human cancers has prompted efforts to develop Mcl-1 inhibitors as potential anticancer agents. In this study, we aimed to design new Mcl-1 inhibitors using various computational techniques. First, we used the Mcl-1 receptor-ligand complex to build an e-pharmacophore hypothesis and screened a library of 567,000 fragments from the Enamine database. We obtained 410 fragments and used them to design 92,384 novel compounds, which we then docked into the Mcl-1 binding cavity using HTVS, SP, and XP docking modes of Glide. To assess their suitability as drug candidates, we conducted MM-GBSA calculations and ADME prediction, leading to the identification of 10 compounds with excellent binding affinity and favorable pharmacokinetic properties. To further investigate the interaction strength, we performed molecular dynamics simulations on the top three Mcl-1 receptor-ligand complexes to study their interaction stability. Overall, our findings suggest that these compounds have promising potential as anticancer agents, pending further experimental validation such as Mcl-1 apoptosis Assay. By combining experimental methods with various in silico approaches, these techniques prove to be invaluable for identifying novel drug candidates with distinct therapeutic applications using fragment-based drug design. This methodology has the potential to expedite the drug discovery process while also reducing its costs.Communicated by Ramaswamy H. Sarma.

Elucidating the potential of bioactive of Trichoderma sp.. in combating pathogenesis by Fusarium sp.. by targeting pectin lyases: a bioinformatics approach

Pectin lyase is an industrially important enzyme, predominately used in fruit juice clarification and retting of fibers. It also promotes pathogenesis via the degradation of the pectin. The phytopathogen, Fusarium infects various crops and causes several diseases. Trichoderma sp. is a promising biocontrol agent that is vital in maintaining plant health and disease prevention. In the current study, a computational approach utilizing structure prediction, molecular docking, molecular dynamics, and MM-PBSA analysis was used to analyze the potential role of bioactive compounds secreted by Trichoderma sp. in inhibiting the pectin lyase enzyme from Fusarium proliferatum, F. fujikuroi, F. graminearum, F. oxysporum and F. verticillioides. Molecular docking with secondary metabolites revealed that Viridiofungin A secreted by Trichoderma harzianum and Virone secreted by T. virens are bioactive compounds with immense potential to inhibit PNLs of Fusarium species. Further, the rigidity of the structure and stability of the docked complex were confirmed via Molecular dynamic simulations assessed through multiple parameters from the simulation trajectory data. Dual culture assay of T. harzianum and T. virens with F. proliferatum, F. fujikuroi, F. graminearum, F. oxysporum, and F. verticillioides showed variable mycelial inhibition. The research provides insight into the potential of the bioactive compounds secreted by Trichoderma species as an effective agent for the inhibition of pectin lyases produced by phytopathogens, especially Fusarium species. The proposed research can be used to develop bioformulations that function as biopesticides, offering a sustainable replacement for chemical products.

In silico vaccine design for Yersinia enterocolitica: A comprehensive approach to enhanced immunogenicity, efficacy and protection

Yersinia enterocolitica, a foodborne pathogen, has emerged as a significant public health concern due to its increased prevalence and multidrug resistance. This study employed reverse vaccinology to identify novel vaccine candidates against Y. enterocolitica through comprehensive in silico analyses. The core genome's conserved protein translocase subunit SecY was selected as the target, and potential B-cell, MHC class I, and MHC class II epitopes were mapped. 3B-cell epitopes, 3 MHCI and 11 MHCII epitopes were acquired. A multi-epitope vaccine construct was designed by incorporating the identified epitopes, TLR4 Agonist was used as adjuvants to enhance the immunogenic response. EAAAK, CPGPG and AYY linkers were used to form a vaccine construct, followed by extensive computational evaluations. The vaccine exhibited desirable physicochemical properties, stable secondary and tertiary structures as evaluated by PDBSum and trRosetta. Moreover, favorable interactions with the human Toll-like receptor 4 (TLR4) was observed by ClusPro. Population coverage analysis estimated the vaccine's applicability across 99.74 % in diverse populations. In addition, molecular dynamics simulations and normal mode analysis confirmed the vaccine's structural stability and dynamics in a simulated biological environment. Furthermore, codon optimization and in silico cloning facilitated the evaluation of the vaccine's expression potential in E. coli and pET-28a was used a recombinant plasmid. This study provides a promising foundation for the development of an efficacious vaccine against Y. enterocolitica infections.

Elucidating the bioremediation potential of laccase and peroxidase enzymes from Bacillus ligniniphilus L1 in antibiotic degradation: A computationally guided study

This study showcased the antibiotic degradation abilities of laccase and catalase-peroxidase from Bacillus ligniniphilus L1, an extremophile, against 18 common antibiotics using computationally guided approach. Molecular docking and simulation identified six enzyme-antibiotic complexes for laccase and four for catalase-peroxidase, demonstrating significant binding affinity and stability. Enzyme activity assays corroborated computational results, indicating both enzymes could degrade all tested antibiotics with varying efficiencies. L1 laccase outperformed commercial laccase against five antibiotics, notably vancomycin, levofloxacin, tobramycin, linezolid, and rifamycin, with enhanced degradation potential upon ABTS addition. Catalase-peroxidase from L1 exhibited superior degradation efficiency over commercial peroxidase against vancomycin, linezolid, tobramycin, and clindamycin. Overall, this study underscores the computational approach's utility in understanding enzyme-mediated antibiotic degradation, offering insights into environmental contaminant remediation.

Muti-target rationale design of novel substituted N-phenyl-2-((6-phenylpyridazin-3-yl)thio)acetamide candidates as telomerase/JAK1/STAT3/TLR4 inhibitors: In vitro and in vivo investigations

In this work, additional effort was applied to design new BIBR1532-based analogues with potential inhibitory activity against telomerase and acting as multitarget antitumor candidates to overcome the resistance problem. Therefore, novel substituted N-phenyl-2-((6-phenylpyridazin-3-yl)thio)acetamide candidates (4a-n) were synthesized. Applying the lead optimization strategy of the previously designed compound 8e; compound 4l showed an improved telomerase inhibition of 64.95 % and a superior growth inhibition of 79 % suggesting its potential use as a successful "multitarget-directed drug" for cancer therapy. Accordingly, compound 4l was further selected to evaluate its additional JAK1/STAT3/TLR4 inhibitory potentials. Compound 4l represented a very promising JAK1 inhibitory potential with a 0.46-fold change, compared to that of pacritinib reference standard (0.33-fold change). Besides, it showed a superior STAT3-inhibitory potential with a 0.22-fold change compared to sorafenib (0.33-fold change). Additionally, compound 4l downregulated TLR4 protein expression by 0.81-fold change compared to that of resatorvid (0.29-fold change). Also, molecular docking was performed to investigate the binding mode and affinity of the superior candidate 4l towards the four target receptors (telomerase, JAK1, STAT3, and TLR4). Furthermore, the therapeutic potential of compound 4l as an antitumor agent was additionally explored through in vivo studies involving female mice implanted with Solid Ehrlich Carcinoma (SEC). Remarkably, compound 4l led to prominent reductions in tumor size and mass. Concurrent enhancements in biochemical, hematologic, histopathologic, and immunohistochemical parameters further confirmed the suppression of angiogenesis and inflammation, elucidating additional mechanisms by which compound 4l exerts its anticancer effects.

SLNP-based CDK4- targeted nanotherapy against glioblastoma

Introduction

Glioblastoma is a grade IV solid brain tumor and has a 15-month survival rate even after treatment. Glioblastoma development is heavily influenced by retinoblastoma protein (pRB) pathway changes. The blood-brain barrier, drug resistance, and severe toxicity of Temozolamide are key obstacles in treating glioblastoma. Innovative treatments targeting the pRB pathway with efficient delivery vehicles are required to treat glioblastoma.

Methods

For this purpose, a library of 691 plant extracts previously tested in vitro for anti-cancerous, anti inflammatory, and anti-proliferative characteristics was created after thorough literature investigations. Compounds were docked against pRB pathway protein ligands using molecular operating environment and chimera. Their nuclear structure and drug-like properties were predicted through Lipinski rule and density functional theory analysis. Physio-chemical characterizations of naked and drug-encapsulated SLNPs assessed size, stability, entrapment efficiency, and drug release rate. Anti-cancer potential of drug and drug- loaded SLNPs was evaluated using U87, U251, and HEK cell lines. Formulations were tested for cancer cell metastatic potential using cell migration assays.

Results

Silymarin (Sil) was identified as the most potent compound against CDK4, which was then encapsulated in stearic acid solid lipid nanoparticles (SLNP-Sil). Sil showed decreased cell viability 72 h after treatment against both U87 and U251 cell lines but had negligible cytotoxic effect on HEK-293. IC50 value of Sil was 155.14 µM for U87 and 195.93 µM for U251. Sil and SLNP-Sil effectively inhibited U87 and U251 cell migration 24 h after treatment.

Discussion

Our results indicated that Sil and SLNP-Sil are promising therapeutic approaches against glioblastoma and merit in vivo experimental verification using orthotropic xenograft mouse models against glioblastoma.

Molecular hybridization modification improves the stability and immunomodulatory activity of TP5 peptide

Thymopentin (TP5) plays an important role in host immunomodulation, yet its bioavailability is significantly limited by its short half-life. YW12D is a peptide with strong stability but relatively weak immunoactivity. Tuning the physicochemical properties of such molecules may yield synthetic molecules displaying optimal stability, safety and enhanced immunological activity. Here, natural peptides were modified to improve their activity by hybridization strategies. A hybrid peptide YW12D-TP5 (YTP) that combines TP5 and YW12D is designed. The half-life of YTP in plasma is significantly longer than that of YW12D and TP5. YTP also displays an improved ability to protect the host from CTX-induced weight loss and thymus and spleen indices decrease than YW12D and TP5. In addition, YTP promotes dendritic cell maturation and increases the expression of cytokines IL-1β, IL-6, TNF-α and immunoglobulins IgA, IgG, and IgM. A combination of antibody-specific blocking assay, SPR, molecular dynamics simulations and western blotting suggest that the immunomodulatory effect of YTP is associated with its activation of the TLR2-NF-кB signaling axis. In sum, we demonstrate that peptide hybridization is an effective strategy for redirecting biological activity to generate novel bioactive molecules with desired properties.

Discovery of novel MLK4 inhibitors against colorectal cancer through computational approaches

Colorectal cancer (CRC) is a significant health issue globally, affecting approximately 10 % of the world's population. The prevalence of CRC highlights the need for effective treatments and prevention strategies. The current therapeutic option, such as chemotherapy, has significant side effects. Thus, this study investigated the anticancer properties of Sanguinarine derivatives, an alkaloid found in traditional herbs via chemoinformatic approaches. Six Sanguinarine derivatives were discovered through virtual screening and molecular docking to determine their binding affinities against the mixed lineage kinase (MLK4) protein which is responsible for CRC. All the compounds were found to be more effective than standard drug used for colorectal cancer treatment, with Sanguinarine derivative 11 showing the highest affinity. The stability of the drug was confirmed through molecular dynamics simulations at 500 ns. This suggests that compound 11 has a higher chance of replacing 5-Fluorouracil, which is currently a widely used chemotherapy drug. Before molecular dynamics simulations, the pharmacokinetic and chemical properties of Sanguinarine derivatives were determined using pkCSM server and DFT method, respectively. The results support that compound 11 is a good drug candidate, as evidenced by Lipinski's Rule of Five. Therefore, compound 11 is recommended for further analysis via in vivo and in vitro studies to confirm its efficacy and safety.

Peptide TaY Attenuates Inflammatory Responses by Interacting with Myeloid Differentiation 2 and Inhibiting NF-κB Signaling Pathway

A balanced inflammatory response is crucial for the organism to defend against external infections, however, an exaggerated response may lead to detrimental effects, including tissue damage and even the onset of disease. Therefore, anti-inflammatory drugs are essential for the rational control of inflammation. In this study, we found that a previously screened peptide TaY (KEKKEVVEYGPSSYGYG) was able to inhibit the LPS-induced RAW264.7 inflammatory response by decreasing a series of proinflammatory cytokines, such as TNF-α, IL-6, and nitric oxide (NO). To elucidate the underlying mechanism, we conducted further investigations. Western blot analysis showed that TaY reduced the phosphorylation of key proteins (IKK-α/β, IκB-α,NF-κB (P65)) in the TLR4-NF-κB signaling pathway and inhibited the inflammatory response. Furthermore, molecular docking and molecular dynamic simulations suggested that TaY binds to the hydrophobic pocket of MD2 through hydrogen bonding and hydrophobic interactions, potentially competing with LPS for MD2 binding. Collectively, TaY is a promising candidate for the development of novel therapeutic strategies against inflammatory disorders.

Structure-Based Discovery of Phytocompounds from Azadirachta indica as Potential Inhibitors of Thioredoxin Glutathione Reductase in Schistosoma mansoni

Schistosomiasis, a parasitic disease caused by Schistosoma species such as S. haematobium, S. mansoni, and S. japonicum, poses a significant global health burden. The thioredoxin glutathione reductase (TGR) enzyme, crucial for maintaining the parasite's redox balance and preventing oxidative stress, has been identified as a promising target for anti-schistosomal drug development. This study aims to identify potential TGR inhibitors from Azadirachta indica phytochemicals using molecular modeling approaches. We screened 60 compounds derived from A. indica bark and leaves through molecular docking to assess their binding affinity, followed by the evaluation of binding-free energies for the most promising candidates. Drug-likeness and pharmacokinetic properties were assessed, and molecular dynamics simulations were conducted to explore the conformational stability of the protein-ligand complexes. Our findings revealed that several A. indica compounds exhibited significantly lower docking scores (up to -9.669 kcal/mol) compared to the standard drug praziquantel (-4.349 kcal/mol). Notably, Isorhamnetin, Isomargolonone, Nimbaflavone, Quercetin, and Nimbionol demonstrated strong interactions with TGR, although Isorhamnetin showed potential mutagenicity. Further binding free energy calculations and molecular dynamics simulations confirmed the stability of Isomargolonone, Nimbionol, and Quercetin as potential TGR inhibitors. In conclusion, these findings suggest that Isomargolonone, Nimbionol, and Quercetin warrant further experimental validation as promising candidates for anti-schistosomal therapy.

Insights into the mechanism of crotamine and potential targets involved in obesity-related metabolic pathways

Crotamine (Ctm) is a peptide isolated from Crotalus durissus terrificus venom. This molecule has been demonstrated to diminish body weight gain and enhance browning in adipose tissue, glucose tolerance, and insulin sensitivity; hence, it has been postulated as an anti-obesogenic peptide. However, the mechanism to elicit the anti-obesogenic effects has yet to be elucidated. Thus, we investigated the possible interaction of Ctm with receptors involved in obesity-related metabolic pathways through protein-protein docking and molecular dynamics refinement. To test the anti-obesogenic mechanism of Ctm, we selected and retrieved 18 targets involved in obesity-related drug discovery from Protein Data Bank. Then, we performed protein-protein dockings. The best three Ctm-target models were selected and refined by molecular dynamics simulations. Molecular docking demonstrated that Ctm was able to interact with 13 of the 18 targets tested. Having a better docking score with glucagon-like peptide-1 receptor (GLP-1R) (-1430.2 kcal/mol), DPP-IV (dipeptidyl peptidase-IV) (-1781.7 kcal/mol) and α-glucosidase (-1232.3 kcal/mol). These three models were refined by molecular dynamics. Ctm demonstrated a higher affinity for GLP-1R (ΔG: -41.886 ± 2.289 kcal/mol). However, Ctm interaction was more stable with DPP-IV (RMSD: 0.360 ± 0.015 nm, Radius of gyration: 2.781 ± 0.009 nm). Moreover, the number of interactions and the molecular mechanics energies of Ctm residues suggest that the interaction of Ctm with these receptors is mainly mediated by basic-hydrophobic dyads Y1-K2, W31-R32, and W33-R34. Together, all these results allow elucidating a possible molecular mechanism behind the previously described anti-obesogenic effects.

Alisol C 23-acetate might be a lead compound of potential lipase inhibitor from Alismatis Rhizoma: Screening, identification and molecular dynamics simulation

Alismatis Rhizoma (AR), a traditional Chinese medicine for treating obesity in traditional Chinese medicine clinic, is recognized as a promising source of lead compounds of lipase inhibitors. Ultrafiltration centrifugal combined with liquid chromatography-mass spectrometry (UF-LC-MS) was used for screening potential lipase inhibitors from AR, and the result indicated the binding capacity between compound 7 and lipase (92.3 ± 1.28 %) was significantly higher than other triterpenoids, and was identified as alisol C 23-acetate. It exhibited a mixed-type inhibitory behavior with an IC50 value of 84.88 ± 1.03 μM. Subsequently, the binding pockets of alisol C 23-acetate to lipase were predicted, and their binding mechanism was explored with molecular simulation. Pocket 1 (active center) and pocket 4 might be the orthosteric and allosteric binding sites of alisol C 23-acetate to lipase, respectively. The interaction between alisol C 23-acetate and lipase was identified to involve key amino acid residues such as GLY-77, PHE-78, TYR-115, LEU-154, PRO-181, PHE-216, LEU-264, ASP-278, GLN-306, ARG-313, and VAL-426. Meanwhile, alisol C 23-acetate remained stable during the intestinal digestive but degraded in the gastric digestion. Overall, alisol C 23-acetate is expected to be the lead compound of lipase inhibitors for treating obesity.

Computational targeting of iron uptake proteins in Covid-19 induced mucormycosis to identify inhibitors via molecular dynamics, molecular mechanics and density function theory studies

Mucormycosis is a concerning invasive fungal infection with difficult diagnosis, high mortality rates, and limited treatment options. Iron availability is crucial for fungal growth that causes this disease. This study aimed to computationally target iron uptake proteins in Rhizopus arrhizus, Lichtheimia corymbifera, and Mucor circinelloides to identify inhibitors, thereby halting fungal growth and intervening in mucormycosis pathogenesis. Seven important iron uptake proteins were identified, modeled, and validated using Ramachandran plots. An in-house antifungal library of ~ 15,401 compounds was screened in molecular docking studies with these proteins. The best small molecule-protein complexes were simulated at 100 ns using Maestro, Schrodinger. Toxicity predictions suggested all six molecules, identified as the best binding compounds to seven proteins, belonged to lower toxicity levels per GHS classification. A molecular mechanics GBSA study for all seven complexes indicated low standard deviations after calculating free binding energies every 10 ns of the 100 ns trajectory. Density functional theory via quantum mechanics approaches highlighted the HOMO, LUMO, and other properties of the six best-bound molecules, revealing their binding capabilities and behaviour. This study sheds light on the molecular mechanisms and protein-ligand interactions, providing a multi-dimensional view towards the use of FDBD01920, FDBD01923, and FDBD01848 as stable antifungal ligands.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00264-7.

Computational Screening of FDA‐Approved Hepatitis C Drugs for Inhibition of VEGFR2 in Liver Cancer

Liver cancer (LC) is one of the most common tumours and the leading cause of cancer‐related death globally. Amidst the problems associated with existing treatments, such as hepatotoxicity, recurrence, drug resistance, and other adverse effects, researchers are under pressure to find alternatives. Towards a comprehensive rationalisation of the search for new anti‐LC drugs among approved ones, we employed an in‐silico approach to accelerate the selection of the most efficacious LC drugs. The FDA‐approved hepatitis C virus (HCV) drugs were docked with the LC protein using the AutoDock Vina software. Compared to the control compound, two FDA‐approved HCV drugs (DB09102 and DB09027) were selected based on their binding energies and interactions with the target protein, which showed comparable binding energies. Furthermore, these compounds were then subjected to molecular dynamic simulation, principle component analysis, and MMGBSA using the AMBER20 software, and the results showed stable complexes compared to the control complex. All things considered, this study will help the scientific community and society find a novel drug to treat LC.

Computational analysis of spike protein of SARS-CoV-2 (Omicron variant) for development of peptide-based therapeutics and diagnostics

In the last few years, the worldwide population has suffered from the SARS-CoV-2 pandemic. The WHO dashboard indicated that around 504,079,039 people were infected and 6,204,155 died from COVID-19 caused by different variants of SARS-CoV-2. Recently, a new variant of SARS-CoV-2 (B.1.1.529) was reported by South Africa known as Omicron. The high transmissibility rate and resistance towards available anti-SARS-CoV-2 drugs/vaccines/monoclonal antibodies, make Omicron a variant of concern. Because of various mutations in spike protein, available diagnostic and therapeutic treatments are not reliable. Therefore, the present study explored the development of some therapeutic peptides that can inhibit the SARS-CoV-2 virus interaction with host ACE2 receptors and can also be used for diagnostic purposes. The screened linear B cell epitopes derived from receptor-binding domain of spike protein of Omicron variant were evaluated as peptide inhibitor/vaccine candidates through different bioinformatics tools including molecular docking and simulation to analyze the interaction between Omicron peptide and human ACE2 receptor. Overall, in-silico studies revealed that Omicron peptides OP1-P12, OP14, OP20, OP23, OP24, OP25, OP26, OP27, OP28, OP29, and OP30 have the potential to inhibit Omicron interaction with ACE2 receptor. Moreover, Omicron peptides OP20, OP22, OP23, OP24, OP25, OP26, OP27, and OP30 have shown potential antigenic and immunogenic properties that can be used in design and development vaccines against Omicron. Although the in-silico validation was performed by comparative analysis with the control peptide inhibitor, further validation through wet lab experimentation is required before its use as therapeutic peptides.Communicated by Ramaswamy H. Sarma.

Employing comparative QSAR techniques for the recognition of dibenzofuran and dibenzothiophene derivatives toward MMP-12 inhibition

Among various matrix metalloproteinases (MMPs), MMP-12 is one of the potential targets for cancer and other diseases. However, none of the MMP-12 inhibitors has passed the clinical trials to date. Therefore, designing potential MMP-12 inhibitors as new drug molecules can provide effective therapeutic strategies for several diseases. In this study, a series of dibenzofuran and dibenzothiophene derivatives were subjected to different 2D and 3D-QSAR techniques to point out the crucial structural contributions highly influential toward the MMP-12 inhibitory activity. These techniques identified some structural attributes of these compounds that are responsible for influencing their MMP-12 inhibition. The carboxylic group may enhance proper binding with catalytic Zn2+ ion at the MMP-12 active site. Again, the i-propyl sulfonamido carboxylic acid function contributed positively toward MMP-12 inhibition. Moreover, the dibenzofuran moiety conferred stable binding at the S1' pocket for higher MMP-12 inhibition. The steric and hydrophobic groups were found favourable near the furan ring substituted at the dibenzofuran moiety. Besides these ligand-based approaches, molecular docking and molecular dynamic (MD) simulation studies not only elucidated the importance of several aspects of these MMP-12 inhibitors while disclosing the significance of the finding of these QSAR studies and their influences toward MMP-12 inhibition. The MD simulation study also revealed stable and compact binding between such compounds at the MMP-12 active site. Therefore, the findings of these validated ligand-based and structure-based molecular modeling studies can aid the development of selective and potent lead molecules that can be used for the treatment of MMP-12-associated diseases.Communicated by Ramaswamy H. Sarma.

A review of biophysical strategies to investigate protein-ligand binding: What have we employed?

The protein-ligand binding frequently occurs in living organisms and plays a crucial role in the execution of the functions of proteins and drugs. It is also an indispensable part of drug discovery and screening. While the methods for investigating protein-ligand binding are diverse, each has its own objectives, strengths, and limitations, which all influence the choice of method. Many studies concentrate on one or a few specific methods, suggesting that comprehensive summaries are lacking. Therefore in this review, these methods are comprehensively summarized and are discussed in detail: prediction and simulation methods, thermal and thermodynamic methods, spectroscopic methods, methods of determining three-dimensional structures of the complex, mass spectrometry-based methods and others. It is also important to integrate these methods based on the specific objectives of the research. With the aim of advancing pharmaceutical research, this review seeks to deepen the understanding of the protein-ligand binding process.

GK, KT. N, Selvaraj C, K. L. Explication of Pharmacological Proficiency of Phytoconstituents from Adansonia digitata Bark: An In Vitro and In Silico Approaches

Compared to other drug discovery sources, traditional medicine has significantly contributed to developing innovative therapeutic molecules for preventive and curative medicine. The Baobab tree, also known as Adansonia digitata L., is significant in Africa due to its multitude of benefits and various parts that serve different purposes, providing economic support to rural communities. The analysis of a plant sample using Fourier transform infrared (FT-IR) spectroscopy detected multiple functional groups, such as carboxyl and aromatic groups. Additionally, gas chromatography-mass spectroscopy (GC-MS) was utilized to identify various compounds present in the sample, including tetrachloroethylene and octyl ester. The results of different assays, such as α-diphenyl-β-picrylhydrazyl (DPPH), superoxide, nitric oxide scavenging assays, and total antioxidant by thiobarbituric acid method (TBA) and ferric thiocyanate (FTC) method, demonstrated a substantial scavenging of free radicals and an effective antioxidant efficacy. The bark's antimicrobial activity was tested through agar diffusion, resulting in a range of zone of inhibition from 10.1 ± 0.36 mm to 20.85 ± 0.76 mm. The minimum inhibitory concentration (MIC) value was observed to be approximately 0.625 µg/mL. The biofilm inhibition percentage ranged from 9.89% to 57.92%, with the highest percentage being 57.92%. The GC-MS and FT-IR studies revealed phytocompounds, which were then analyzed for their potential therapeutic properties. Computational studies were conducted on the phytocompounds against Pseudomonas aeruginosa and C2 kinase (antioxidant). The study concluded that the Adansonia digitata bark extract and its phytocompound have potential therapeutic efficacy against the target proteins. The best docking scores were about -7.053 kcal/mol and -7.573 kcal/mol for Pseudomonas aeruginosa and C2 kinase (antioxidant), respectively. The interaction patterns with the crucial amino acid residues elucidate the inhibitory efficacy of the phytocompounds.

Novel Insights into the Antimicrobial and Antibiofilm Activity of Pyrroloquinoline Quinone (PQQ); In Vitro, In Silico, and Shotgun Proteomic Studies

Microbial infections pose a significant global health threat, affecting millions of individuals and leading to substantial mortality rates. The increasing resistance of microorganisms to conventional treatments requires the development of novel antimicrobial agents. Pyrroloquinoline quinone (PQQ), a natural medicinal drug involved in various cellular processes, holds promise as a potential antimicrobial agent. In the present study, our aim was, for the first time, to explore the antimicrobial activity of PQQ against 29 pathogenic microbes, including 13 fungal strains, 8 Gram-positive bacteria, and 8 Gram-negative bacteria. Our findings revealed potent antifungal properties of PQQ, particularly against Syncephalastrum racemosum, Talaromyces marneffei, Candida lipolytica, and Trichophyton rubrum. The MIC values varied between fungal strains, and T. marneffei exhibited a lower MIC, indicating a greater susceptibility to PQQ. In addition, PQQ exhibited notable antibacterial activity against Gram-positive and -negative bacteria, with a prominent inhibition observed against Staphylococcus epidermidis, Proteus vulgaris, and MRSA strains. Remarkably, PQQ demonstrated considerable biofilm inhibition against the MRSA, S. epidermidis, and P. vulgaris strains. Transmission electron microscopy (TEM) studies revealed that PQQ caused structural damage and disrupted cell metabolism in bacterial cells, leading to aberrant morphology, compromised cell membrane integrity, and leakage of cytoplasmic contents. These findings were further affirmed by shotgun proteomic analysis, which revealed that PQQ targets several important cellular processes in bacteria, including membrane proteins, ATP metabolic processes, DNA repair processes, metal-binding proteins, and stress response. Finally, detailed molecular modeling investigations indicated that PQQ exhibits a substantial binding affinity score for key microbial targets, including the mannoprotein Mp1P, the transcriptional regulator TcaR, and the endonuclease PvuRTs1I. Taken together, our study underscores the effectiveness of PQQ as a broad-spectrum antimicrobial agent capable of combating pathogenic fungi and bacteria, while also inhibiting biofilm formation and targeting several critical biological processes, making it a promising therapeutic option for biofilm-related infections.

Recent advances in COVID-19-induced liver injury: causes, diagnosis, and management

Since the start of the pandemic, considerable advancements have been made in our understanding of the effects of SARS-CoV-2 infection and the associated COVID-19 on the hepatic system. There is a broad range of clinical symptoms for COVID-19. It affects multiple systems and has a dominant lung illness depending on complications. The progression of COVID-19 in people with pre-existing chronic liver disease (CLD) has also been studied in large multinational groups. Notably, SARS-CoV-2 infection is associated with a higher risk of hepatic decompensation and death in patients with cirrhosis. In this review, the source, composition, mechanisms, transmission characteristics, clinical characteristics, therapy, and prevention of SARS-CoV-2 were clarified and discussed, as well as the evolution and variations of the virus. This review briefly discusses the causes and effects of SARS-CoV-2 infection in patients with CLD. As part of COVID-19, In addition, we assess the potential of liver biochemistry as a diagnostic tool examine the data on direct viral infection of liver cells, and investigate potential pathways driving SARS-CoV-2-related liver damage. Finally, we explore how the pandemic has had a significant impact on patient behaviors and hepatology services, which may increase the prevalence and severity of liver disease in the future. The topics encompassed in this review encompass the intricate relationships between SARS-CoV-2, liver health, and broader health management strategies, providing valuable insights for both current clinical practice and future research directions.

Diterpenoid and C20 diterpenoid alkaloid as a potent inhibitor of SARS-CoV-2 main protease (Mpro): from Piper barberi Gamble, an endemic and endangered species of Southern Western Ghats

The present study investigated the phytochemicals and in silico anti-nCoV properties of Piper barberi, an endangered and endemic species of Southern Western Ghats. Using conventional soxhlet extraction method, the leaf and stem were extracted separately with methanol (PBLM and PBSM). The bioactive compounds from the extracts were identified using HR-LCMS/MS-qTOF analysis. These compounds were subjected to various in silico analyses to identify potential drug candidates against nCoV. The HR LCMS/MS analysis of PBLM and PBSM revealed the presence of phenols, flavonoids, alkaloids, and terpenoids in it and this is the first report of the phytoconstituents present in the species P. barberi. All the identified bioactive compounds were subjected to predict ADMET. Out of 49 identified compounds, only 31 passed drug-likeness properties and toxicity tests. Molecular interaction studies were conducted using the AutoDockTools 4.2.6., which showed that only 13 compounds exhibited acceptable binding affinity with the nCoV target Mpro. Structural stability and binding free energy analyses of the five compounds with the higher binding affinity indicated that the bioactive compounds Hetisine and Ajaconine are stable with both hydrogen bonds and hydrophobic interactions. Hetisine shows stable binding among these two compounds with two hydrogen bond interactions with the crucial catalytic dyad residue (His41). Thus, this study concludes that these compounds might potentially be used as an alternative drug candidate for managing nCoV. However, further experimental validation, including in vitro and in vivo assays, is required to substantiate the results.Communicated by Ramaswamy H. Sarma.

A computational approach to developing a multi-epitope vaccine for combating Pseudomonas aeruginosa-induced pneumonia and sepsis

Pseudomonas aeruginosa is a complex nosocomial infectious agent responsible for numerous illnesses, with its growing resistance variations complicating treatment development. Studies have emphasized the importance of virulence factors OprE and OprF in pathogenesis, highlighting their potential as vaccine candidates. In this study, B-cell, MHC-I, and MHC-II epitopes were identified, and molecular linkers were active to join these epitopes with an appropriate adjuvant to construct a vaccine. Computational tools were employed to forecast the tertiary framework, characteristics, and also to confirm the vaccine's composition. The potency was weighed through population coverage analysis and immune simulation. This project aims to create a multi-epitope vaccine to reduce P. aeruginosa-related illness and mortality using immunoinformatics resources. The ultimate complex has been determined to be stable, soluble, antigenic, and non-allergenic upon inspection of its physicochemical and immunological properties. Additionally, the protein exhibited acidic and hydrophilic characteristics. The Ramachandran plot, ProSA-web, ERRAT, and Verify3D were employed to ensure the final model's authenticity once the protein's three-dimensional structure had been established and refined. The vaccine model showed a significant binding score and stability when interacting with MHC receptors. Population coverage analysis indicated a global coverage rate of 83.40%, with the USA having the highest coverage rate, exceeding 90%. Moreover, the vaccine sequence underwent codon optimization before being cloned into the Escherichia coli plasmid vector pET-28a (+) at the EcoRI and EcoRV restriction sites. Our research has developed a vaccine against P. aeruginosa that has strong binding affinity and worldwide coverage, offering an acceptable way to mitigate nosocomial infections.

Interaction of systemic drugs causing ocular toxicity with organic cation transporter: an artificial intelligence prediction

Chronic disease patients (cancer, arthritis, cardiovascular diseases) undergo long-term systemic drug treatment. Membrane transporters in ocular barriers could falsely recognize these drugs and allow their trafficking into the eye from systemic circulation. Hence, despite their pharmacological activity, these drugs accumulate and cause toxicity at the non-target site, such as the eye. Since around 40% of clinically used drugs are organic cation in nature, it is essential to understand the role of organic cation transporter (OCT1) in ocular barriers to facilitate the entry of systemic drugs into the eye. We applied machine learning techniques and computer simulation models (molecular dynamics and metadynamics) in the current study to predict the potential OCT1 substrates. Artificial intelligence models were developed using a training dataset of a known substrates and non-substrates of OCT1 and predicted the potential OCT1 substrates from various systemic drugs causing ocular toxicity. Computer simulation studies was performed by developing the OCT1 homology model. Molecular dynamic simulations equilibrated the docked protein-ligand complex. And metadynamics revealed the movement of substrates across the transporter with minimum free energy near the binding pocket. The machine learning model showed an accuracy of about 80% and predicted the potential substrates for OCT1 among systemic drugs causing ocular toxicity - not known earlier, such as cyclophosphamide, bupivacaine, bortezomib, sulphanilamide, tosufloxacin, topiramate, and many more. However, further invitro and invivo studies are required to confirm these predictions.Communicated by Ramaswamy H. Sarma.

Exploring the effectiveness of flavone derivatives for treating liver diseases: Utilizing DFT, molecular docking, and molecular dynamics techniques

In exploring nature's potential in addressing liver-related conditions, this study investigates the therapeutic capabilities of flavonoids. Utilizing in silico methodologies, we focus on flavone and its analogs (1-14) to assess their therapeutic potential in treating liver diseases. Molecular change calculations using density functional theory (DFT) were conducted on these compounds, accompanied by an evaluation of each analog's physiochemical and biochemical properties. The study further assesses these flavonoids' binding effectiveness and locations through molecular docking studies against six target proteins associated with human cancer. Tropoflavin and taxifolin served as reference drugs. The structurally modified flavone analogs (1-14) displayed a broad range of binding affinities, ranging from -7.0 to -9.4 kcal mol⁻¹, surpassing the reference drugs. Notably, flavonoid (7) exhibited significantly higher binding affinities with proteins Nrf2 (PDB:1 × 2 J) and DCK (PDB:1 × 2 J) (-9.4 and -8.1 kcal mol⁻¹) compared to tropoflavin (-9.3 and -8.0 kcal mol⁻¹) and taxifolin (-9.4 and -7.1 kcal mol⁻¹), respectively. Molecular dynamics (MD) simulations revealed that the docked complexes had a root mean square deviation (RMSD) value ranging from 0.05 to 0.2 nm and a root mean square fluctuation (RMSF) value between 0.35 and 1.3 nm during perturbation. The study concludes that 5,7-dihydroxyflavone (7) shows substantial promise as a potential therapeutic agent for liver-related conditions. However, further validation through in vitro and in vivo studies is necessary. Key insights from this study include:•Screening of flavanols and their derivatives to determine pharmacological and bioactive properties using ADMET, molinspiration, and pass prediction analysis.•Docking of shortlisted flavone derivatives with proteins having essential functions.•Analysis of the best protein-flavonoid docked complexes using molecular dynamics simulation to determine the flavonoid's efficiency and stability within a system.

Identification of FDA-approved drugs with triple targeting mode of action for the treatment of monkeypox: a high throughput virtual screening study

According to the Center for Disease Control and Prevention, as of August 23, 94 countries had confirmed 42,954 Monkeypox Virus cases. As specific monkeypox drugs are not yet developed, the treatment depends on repurposed FDA-approved drugs. According to a recent study, the Monkeypox outbreak is caused by a strain with a unique mutation, raising the likelihood that the virus will develop resistance to current drugs by acquiring mutations in the targets of currently used drugs. The probability of multiple mutations in two or more drug targets at a time is always low than mutation in a single drug target. Therefore, we identified 15 triple-targeting FDA-approved drugs that can inhibit three viral targets, including topoisomerase1, p37, and thymidylate kinase, using high throughput virtual screening approach. Further, the molecular dynamics simulation analysis of the top hits such as Naldemedine and Saquinavir with their respective targets reveals the formation of stable conformational changes of the ligand-protein complexes inside the dynamic biological environment. We suggest further research on these triple-targeting molecules to develop an effective therapy for the currently spreading Monkeypox.

Exploring different classification-dependent QSAR modelling strategies for HDAC3 inhibitors in search of meaningful structural contributors

Histone deacetylase 3 (HDAC3), a Zn2+-dependent class I HDACs, contributes to numerous disorders such as neurodegenerative disorders, diabetes, cardiovascular disease, kidney disease and several types of cancers. Therefore, the development of novel and selective HDAC3 inhibitors might be promising to combat such diseases. Here, different classification-based molecular modelling studies such as Bayesian classification, recursive partitioning (RP), SARpy and linear discriminant analysis (LDA) were conducted on a set of HDAC3 inhibitors to pinpoint essential structural requirements contributing to HDAC3 inhibition followed by molecular docking study and molecular dynamics (MD) simulation analyses. The current study revealed the importance of hydroxamate function for Zn2+ chelation as well as hydrogen bonding interaction with Tyr298 residue. The importance of hydroxamate function for higher HDAC3 inhibition was noticed in the case of Bayesian classification, recursive partitioning and SARpy models. Also, the importance of substituted thiazole ring was revealed, whereas the presence of linear alkyl groups with carboxylic acid function, any type of ester function, benzodiazepine moiety and methoxy group in the molecular structure can be detrimental to HDAC3 inhibition. Therefore, this study can aid in the design and discovery of effective novel HDAC3 inhibitors in the future.

Powerful Approach for New Drugs as Antibacterial Agents via Molecular Docking and In Vitro Studies of Some New Cyclic Imides and Quinazoline-2,5-diones

We generated novel elven 1,2,3,6-tetrahydrophthalimides and tetrahydroquinazoline derivatives from 1,2,3,6-tetrahydrophthalic anhydride (1) in response to our interest in using the anhydrides to produce heterocyclic nitrogen compounds. The elemental and spectral analyses of the produced compounds validated the recommended configurations and MOE 2014.09 (Molecular Operating Environment) computations were used to perform their in silico analysis. The synthesized compounds have been analyzed and put through various experiments, including in vitro and in silico methods to assess their biological activity against Escherichia coli Penicillin-Binding Protein 3 (PBP3) and Staphylococcus aureus Penicillin-Binding Protein 2 (PBP2), among these compounds showing promising data as antibacterial drugs.

Repurposing Drugs to Modulate Sortilin: Structure-Guided Strategies Against Atherogenesis, Coronary Artery Disease, and Neurological Disorders

Sortilin (SORT1) is a multifunctional protein intricately involved in atherogenesis, coronary artery disease (CAD), and various neurological disorders. It has materialized as a potential pharmacological target for therapeutic development due to its diverse biological roles in pathological processes. Despite its central role under these conditions, effective therapeutic strategies targeting SORT1 remain challenging. In this study, we introduce a drug repurposing strategy guided by structural insights to identify potent SORT1 inhibitors with broad therapeutic potential. Our approach combines molecular docking, virtual screening, and molecular dynamics (MD) simulations, enabling the systematic evaluation of 3648 FDA-approved drugs for their potential to modulate SORT1. The investigation reveals a subset of repurposed drugs exhibiting highly favorable binding profiles and stable interactions within the binding site of SORT1. Notably, two hits, ergotamine and digitoxin, were carefully chosen based on their drug profiles and subjected to analyze their interactions with SORT1 and stability assessment via all-atom MD simulations spanning 300 ns (ns). The structural analyses uncover the complex binding interactions between these identified compounds and SORT1, offering essential mechanistic insights. Additionally, we explore the clinical implications of repurposing these compounds as potential therapeutic agents, emphasizing their significance in addressing atherogenesis, CAD, and neurological disorders. Overall, this study highlights the efficacy of structure-guided drug repurposing and provides a solid foundation for future research endeavors aimed at the development of effective therapies targeting SORT1 under diverse pathological conditions.

Identification of new potent NLRP3 inhibitors by multi-level in-silico approaches

Nod-like receptor protein 3 (NLRP-3), is an intracellular sensor that is involved in inflammasome activation, and the aberrant expression of NLRP3 is responsible for diabetes mellitus, its complications, and many other inflammatory diseases. NLRP3 is considered a promising drug target for novel drug design. Here, a pharmacophore model was generated from the most potent inhibitor, and its validation was performed by the Gunner-Henry scoring method. The validated pharmacophore was used to screen selected compounds databases. As a result, 646 compounds were mapped on the pharmacophore model. After applying Lipinski's rule of five, 391 hits were obtained. All the hits were docked into the binding pocket of target protein. Based on docking scores and interactions with binding site residues, six compounds were selected potential hits. To check the stability of these compounds, 100 ns molecular dynamic (MD) simulations were performed. The RMSD, RMSF, DCCM and hydrogen bond analysis showed that all the six compounds formed stable complex with NLRP3. The binding free energy with the MM-PBSA approach suggested that electrostatic force, and van der Waals interactions, played a significant role in the binding pattern of these compounds. Thus, the outcomes of the current study could provide insights into the identification of new potential NLRP3 inflammasome inhibitors against diabetes and its related disorders.

In silico exploration of phenolics as modulators of penicillin binding protein (PBP) 2× of Streptococcus pneumoniae

Infections caused by multidrug-resistant Streptococcus pneumoniae remain the leading cause of pneumonia-related deaths in children < 5 years globally, and mutations in penicillin-binding protein (PBP) 2 × have been identified as the major cause of resistance in the organism to beta-lactams. Thus, the development of new modulators with enhanced binding of PBP2x is highly encouraged. In this study, phenolics, due to their reported antibacterial activities, were screened against the active site of PBP2x using structure-based pharmacophore and molecular docking techniques, and the ability of the top-hit phenolics to inhibit the active and allosteric sites of PBP2x was refined through 120 ns molecular dynamic simulation. Except for gallocatechin gallate and lysidicichin, respectively, at the active and allosteric sites of PBP2x, the top-hit phenolics had higher negative binding free energy (ΔGbind) than amoxicillin [active site (- 19.23 kcal/mol), allosteric site (- 33.75 kcal/mol)]. Although silicristin had the best broad-spectrum effects at the active (- 38.41 kcal/mol) and allosteric (- 50.54 kcal/mol) sites of PBP2x, the high thermodynamic entropy (4.90 Å) of the resulting complex might suggest the need for its possible structural refinement for enhanced potency. Interestingly, silicristin had a predicted synthetic feasibility score of < 5 and quantum calculations using the DFT B3LYP/6-31G+ (dp) revealed that silicristin is less stable and more reactive than amoxicillin. These findings point to the possible benefits of the top-hit phenolics, and most especially silicristin, in the direct and synergistic treatment of infections caused by S. pneumoniae. Accordingly, silicristin is currently the subject of further confirmatory in vitro research.

Computational screening of natural MtbDXR inhibitors for novel anti-tuberculosis compound discovery

DXR (1-deoxy-d-xylulose-5-phosphate reductoisomerase) is an essential enzyme in the Methylerythritol 4-phosphate (MEP) pathway, which is used by M. tuberculosis and a few other pathogens. This essential enzyme in the isoprenoid synthesis pathway has been previously reported as an important target for antibiotic drug design. However, till now, there is no record of any drug-like safe molecule to inhibit MtbDXR. Numerous plant species have been traditionally used for tuberculosis therapies. In this study, we selected six plant species with anti-tubercular properties. The chemoinformatic screening was performed on 352 phytochemicals from those plants against the MtbDXR protein. After molecular docking analysis, we filtered the top five compounds, CID: 5280443 (Apigenin), CID: 3220 (Emodin), CID: 5280863 (Kaempferol), CID: 5280445 (Luteolin), and CID: 6101979 (beta-Hydroxychalcone), based on binding affinity. Molecular dynamics simulations disclosed the stability of the compounds at the active site of the proteins. Finally, in silico ADME and toxicity evaluations confirmed the compounds to be effective and safe for oral administration. Thus, our findings identified three drug-like safe molecules- Apigenin, Kaempferol, and beta-Hydroxychalcone, that showed good stability in the protein's active site. The results of this computational approach may act as an initial instruction for future in vitro and in vivo testing to identify natural drug-like compounds to treat tuberculosis.Communicated by Ramaswamy H. Sarma.

Study of the binding interaction of salmon sperm DNA with nintedanib, a tyrosine kinase inhibitor using multi-spectroscopic, thermodynamic, and in silico approaches

The study of the intermolecular binding interaction of small molecules with DNA can guide the rational drug design with greater efficacy and improved or more selective activity. In the current study, nintedanib's binding interaction with salmon sperm DNA (ssDNA) was thoroughly investigated using UV-vis spectrophotometry, spectrofluorimetry, ionic strength measurements, viscosity measurements, thermodynamics, molecular docking, and molecular dynamic simulation techniques under physiologically simulated conditions (pH 7.4). The obtained experimental results showed that nintedanib and ssDNA had an apparent binding interaction. Nintedanib's binding constant (Kb) with ssDNA, as determined using the Benesi-Hildebrand plot, was 7.9 × 104 M-1 at 298 K, indicating a moderate binding affinity. The primary binding contact forces were hydrophobic and hydrogen bonding interactions, as verified by the enthalpy and entropy changes (ΔH0 and ΔS0), which were - 16.25 kJ.mol-1 and 39.30 J mol-1 K-1, respectively. According to the results of UV-vis spectrophotometry, viscosity assays, and competitive binding interactions with ethidium bromide or rhodamine B, the binding mode of nintedanib to ssDNA was minor groove. Molecular docking and molecular dynamic simulation studies showed that nintedanib fitted into the B-DNA minor groove's AT-rich region with high stability. This study can contribute to further understanding of nintedanib's molecular mechanisms and pharmacological effects.

Repurposing immune boosting and anti-viral efficacy of Parkia bioactive entities as multi-target directed therapeutic approach for SARS-CoV-2: exploration of lead drugs by drug likeness, molecular docking and molecular dynamics simulation methods

The COVID-19 pandemic has caused adverse health (severe respiratory, enteric and systemic infections) and environmental impacts that have threatened public health and the economy worldwide. Drug repurposing and small molecule multi-target directed herbal medicine therapeutic approaches are the most appropriate exploration strategies for SARS-CoV-2 drug discovery. This study identified potential multi-target-directed Parkia bioactive entities against SARS-CoV-2 receptors (S-protein, ACE2, TMPRSS2, RBD/ACE2, RdRp, MPro, and PLPro) using ADMET, drug-likeness, molecular docking (AutoDock, FireDock and HDOCK), molecular dynamics simulation and MM-PBSA tools. One thousand Parkia bioactive entities were screened out by virtual screening and forty-five bioactive phytomolecules were selected based on favorable binding affinity and acceptable pharmacokinetic and pharmacodynamics properties. The binding affinity values of Parkia phyto-ligands (AutoDock: -6.00--10.40 kcal/mol; FireDock: -31.00--62.02 kcal/mol; and HDOCK: -150.0--294.93 kcal/mol) were observed to be higher than the reference antiviral drugs (AutoDock: -5.90--9.10 kcal/mol; FireDock: -35.64--59.35 kcal/mol; and HDOCK: -132.82--211.87 kcal/mol), suggesting a potent modulatory action of Parkia bioactive entities against the SARS-CoV-2. Didymin, rutin, epigallocatechin gallate, epicatechin-3-0-gallate, hyperin, ursolic acid, lupeol, stigmasta-5,24(28)-diene-3-ol, ellagic acid, apigenin, stigmasterol, and campesterol strongly bound with the multiple targets of the SARS-CoV-2 receptors, inhibiting viral entry, attachment, binding, replication, transcription, maturation, packaging and spread. Furthermore, ACE2, TMPRSS2, and MPro receptors possess significant molecular dynamic properties, including stability, compactness, flexibility and total binding energy. Residues GLU-589, and LEU-95 of ACE2, GLN-350, HIS-186, and ASP-257 of TMPRSS2, and GLU-14, MET-49, and GLN-189 of MPro receptors contributed to the formation of hydrogen bonds and binding interactions, playing vital roles in inhibiting the activity of the receptors. Promising results were achieved by developing multi-targeted antiviral Parkia bioactive entities as lead and prospective candidates under a small molecule strategy against SARS-CoV-2 pathogenesis. The antiviral activity of Parkia bioactive entities needs to be further validated by pre-clinical and clinical trials.

Molecular Dynamics as a Tool for Virtual Ligand Screening

Rational drug design is essential for new drugs to emerge, especially when the structure of a target protein or nucleic acid is known. To that purpose, high-throughput virtual ligand screening campaigns aim at discovering computationally new binding molecules or fragments to modulate particular biomolecular interactions or biological activities, related to a disease process. The structure-based virtual ligand screening process primarily relies on docking methods which allow predicting the binding of a molecule to a biological target structure with a correct conformation and the best possible affinity. The docking method itself is not sufficient as it suffers from several and crucial limitations (lack of full protein flexibility information, no solvation and ion effects, poor scoring functions, and unreliable molecular affinity estimation).At the interface of computer techniques and drug discovery, molecular dynamics (MD) allows introducing protein flexibility before or after a docking protocol, refining the structure of protein-drug complexes in the presence of water, ions, and even in membrane-like environments, describing more precisely the temporal evolution of the biological complex and ranking these complexes with more accurate binding energy calculations. In this chapter, we describe the up-to-date MD, which plays the role of supporting tools in the virtual ligand screening (VS) process.Without a doubt, using docking in combination with MD is an attractive approach in structure-based drug discovery protocols nowadays. It has proved its efficiency through many examples in the literature and is a powerful method to significantly reduce the amount of required wet experimentations (Tarcsay et al, J Chem Inf Model 53:2990-2999, 2013; Barakat et al, PLoS One 7:e51329, 2012; De Vivo et al, J Med Chem 59:4035-4061, 2016; Durrant, McCammon, BMC Biol 9:71-79, 2011; Galeazzi, Curr Comput Aided Drug Des 5:225-240, 2009; Hospital et al, Adv Appl Bioinforma Chem 8:37-47, 2015; Jiang et al, Molecules 20:12769-12786, 2015; Kundu et al, J Mol Graph Model 61:160-174, 2015; Mirza et al, J Mol Graph Model 66:99-107, 2016; Moroy et al, Future Med Chem 7:2317-2331, 2015; Naresh et al, J Mol Graph Model 61:272-280, 2015; Nichols et al, J Chem Inf Model 51:1439-1446, 2011; Nichols et al, Methods Mol Biol 819:93-103, 2012; Okimoto et al, PLoS Comput Biol 5:e1000528, 2009; Rodriguez-Bussey et al, Biopolymers 105:35-42, 2016; Sliwoski et al, Pharmacol Rev 66:334-395, 2014).