Identification of secondary metabolites from Crescentia cujete as promising antibacterial therapeutics targeting type 2A topoisomerases through molecular dynamics simulation

Cheminformatics identification of phenolics as modulators of penicillin-binding protein-3 of Pseudomonas aeruginosa towards interventive antibacterial therapy

Antibacterial resistance to β-lactams in microorganisms has been attributed majorly to alterations in penicillin-binding proteins (PBPs) coupled with β-lactams' inactivation by β-lactamase. Consequently, the identification of a novel class of therapeutics with improved modulatory action on the PBPs is imperative and plant secondary metabolites, including phenolics, have found relevance in this regard. For the first time in this study, the over 10,000 phenolics currently known were computationally evaluated against PBP3 of Pseudomonas aeruginosa, a superbug implicated in several nosocomial infections. In doing this, a library of phenolics with an affinity for PBP3 of P. aeruginosa was screened using structure-activity relationship-based pharmacophore and molecular docking approaches. Subsequent thermodynamic screening of the top five phenolics with higher docking scores, more drug-likeness attributes, and feasible synthetic accessibility was achieved through a 120 ns molecular dynamic (MD) simulation. Four of the top five hits had higher binding free energy than cefotaxime (-18.72 kcal/mol), with catechin-3-rhamside having the highest affinity (-28.99 kcal/mol). All the hits were stable at the active site of the PBP3, with catechin-3-rhamside being the most stable (2.14 Å), and established important interactions with Ser294, implicated in the catalytic activity of PBP3. Also, PBP3 became more compact with less fluctuation of the active site amino acid residues following the binding of the hits. These observations are indicative of the potential of the test compounds as PBP3 inhibitors, with catechin-3-rhamside being the most prominent of the compounds that could be further improved for enhanced druggability against PBP3 in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

Cheminformatics identification of modulators of key carbohydrate-metabolizing enzymes from C. cujete for type-2 diabetes mellitus intervention

Purpose

The therapeutic use of oral hypoglycaemic agents in the management of type-2 diabetes mellitus (T2DM) is without adverse effects; thus, calls for alternative and novel candidates from natural products in medicinal plants.

Method

The study explored molecular docking and molecular dynamics (MD) simulation approaches to identify key antidiabetic metabolites from Crescentia cujete.

Results

Molecular docking results identified four and/or five best compounds against each target enzyme (alpha-glucosidase, dipeptidyl peptidase-IV, aldose reductase, and protein tyrosine phosphatase-1B (PTP-1B)) implicated in diabetes. The resulting complexes (except against PTP-1B) had higher docking scores above respective standards (acarbose, Diprotin A, ranirestat). The MD simulation results revealed compounds such as benzoic acid (-48.414 kcal/mol) and phytol (-45.112 kcal/mol) as well as chlorogenic acid (-42.978 kcal/mol) and naringenin (-31.292 kcal/mol) had higher binding affinities than the standards [acarbose (-28.248 kcal/mol), ranirestat (-21.042 kcal/mol)] against alpha-glucosidase and aldose reductase, respectively while Diprotin A (-45.112 kcal/mol) and ursolic acid (-18.740 kcal/mol) presented superior binding affinities than the compounds [luteolin (-41.957 kcal/mol and naringenin (-16.518 kcal/mol)] against DPP-IV and PTP-1B respectively.

Conclusion

While isoflavone (alpha-glucosidase), xylocaine (DPP-IV), luteolin (aldose reductase,) and chlorogenic acid (PTP-1B) were affirmed as the best inhibitors of respective enzyme targets, luteolin, and chlorogenic acid may be suggested and proposed as probable candidates against T2DM and related retinopathy complication based on their structural stability, compactness and affinity for three (DPP-IV, aldose reductase, and PTP-1B) of the four targets investigated. Further studies are warranted in vitro and in vivo on the antihyperglycaemic effects of these drug candidates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40200-023-01249-7.

Selected phytochemicals of Momordica charantia L. as potential anti-DENV-2 through the docking, DFT and molecular dynamic simulation

Dengue fever is now one of the major global health concerns particularly for tropical and sub-tropical countries. However, there has been no FDA approved medication to treat dengue fever. Researchers are looking into DENV NS5 RdRp protease as a potential therapeutic target for discovering effective anti-dengue agents. The aim of this study to discover dengue virus inhibitor from a set of five compounds from Momordica charantia L. using a series of in-silico approaches. The compounds were docked into the active area of the DENV-2 NS5 RdRp protease to obtain the hit compounds. The successful compounds underwent additional testing for a study on drug-likeness similarity. Our study obtained Momordicoside-I as a lead compound which was further exposed to the Cytochrome P450 (CYP450) toxicity analysis to determine the toxicity based on docking scores and drug-likeness studies. Moreover, DFT studies were carried out to calculate the thermodynamic, molecular orbital and electrostatic potential properties for the lead compound. Moreover, the lead compound was next subjected to molecular dynamic simulation for 200 ns in order to confirm the stability of the docked complex and the binding posture discovered during docking experiment. Overall, the lead compound has demonstrated good medication like qualities, non-toxicity, and significant binding affinity towards the DENV-2 RdRp enzyme.Communicated by Ramaswamy H. Sarma.

Application of molecular dynamics simulation for exploring the roles of plant biomolecules in promoting environmental health

Understanding the dynamic changes of plant biomolecules is vital for exploring their mechanisms in the environment. Molecular dynamics (MD) simulation has been widely used to study structural evolution and corresponding properties of plant biomolecules at the microscopic scale. Here, this review (i) outlines structural properties of plant biomolecules, and the crucial role of MD simulation in advancing studies of the biomolecules; (ii) describes the development of MD simulation in plant biomolecules, determinants of simulation, and analysis parameters; (iii) introduces the applications of MD simulation in plant biomolecules, including the response of the biomolecules to multiple stresses, their roles in corrosive environments, and their contributions in improving environmental health; (iv) reviews techniques integrated with MD simulation, such as molecular biology, quantum mechanics, molecular docking, and machine learning modeling, which bridge gaps in MD simulation. Finally, we make suggestions on determination of force field types, investigation of plant biomolecule mechanisms, and use of MD simulation in combination with other techniques. This review provides comprehensive summaries of the mechanisms of plant biomolecules in the environment revealed by MD simulation and validates it as an applicable tool for bridging gaps between macroscopic and microscopic behavior, providing insights into the wide application of MD simulation in plant biomolecules.

Molecular Modeling Identification of Key Secondary Metabolites from Xylopia aethiopica as Promising Therapeutics Targeting Essential Measles Viral Proteins

This study computationally screened three key compounds (vanillin (VAN), oxophoebine (OPB), and dihydrochalcone (DHC)) derived from Xylopia aethiopica (Guinea pepper), a medicinal plant with known antiviral activity, against key druggable measles virus (MV) proteins (fusion protein (FUP), haemagglutinin protein (HMG), and phosphoprotein (PSP)). Each molecular species was subjected to a 100 ns molecular dynamics (MD) simulation following docking, and a range of postdynamic parameters including free binding energy and pharmacokinetic properties were determined. The docking scores of the resulting OPB-FUP (-5.4 kcal/mol), OPB-HMG (-8.1 kcal/mol), and OPB-PSP (-8.0 kcal/mol) complexes were consistent with their respective binding energy values (-25.37, -28.74, and -40.68 kcal/mol), and higher than that of the reference standard, ribavirin (RBV) in each case. Furthermore, all the investigated compounds were thermodynamically compact and stable, especially HMG of MV, and this observation could be attributed to the resulting intermolecular interactions in each system. Overall, OPB may possess inhibitory properties against MV glycoproteins (FUP and HMG) and PSP that play important roles in the replication of MV and measles pathogenesis. While OPB could serve as a scaffold for the development of novel MV fusion and entry inhibitors, further in vitro and in vivo evaluation is highly recommended.

Network Pharmacology- and Molecular Dynamics Simulation-Based Bioprospection of Aspalathus linearis for Type-2 Diabetes Care

The medicinal herb Aspalathus linearis (rooibos) is globally recognized in type-2 diabetes mellitus (T2DM) treatment due to its known and distinctive compounds. This work utilized network pharmacology (NP) coupled with molecular dynamics simulation in gaining new insight into the anti-diabetic molecular mechanism of action of rooibos teas. It looked at the interactions between rooibos constituents with various relevant protein receptors and signaling routes associated with T2DM progression. The initial analysis revealed 197 intersecting gene targets and 13 bioactive rooibos constituents linked to T2DM. The interactions between proteins and compounds to the target matrix were generated with the Cystoscope platform and STRING database. These analyses revealed intersecting nodes active in T2DM and hypoxia-inducible factor 1 (HIF-1) as an integral receptors target. In addition, KEGG analysis identified 11 other pathways besides the hub HIF-1 signaling route which may also be targeted in T2DM progression. In final molecular docking and dynamics simulation analysis, a significant binding affinity was confirmed for key compound-protein matrices. As such, the identified rooibos moieties could serve as putative drug candidates for T2DM control and therapy. This study shows rooibos constituents' interaction with T2DM-linked signaling pathways and target receptors and proposes vitexin, esculin and isovitexin as well as apigenin and kaempferol as respective pharmacologically active rooibos compounds for the modulation of EGFR and IGF1R in the HIF-1 signaling pathway to maintain normal homeostasis and function of the pancreas and pancreatic β-cells in diabetics.

Cheminformatics Identification of Phenolics as Modulators of Penicillin-Binding Protein 2a of Staphylococcus aureus: A Structure–Activity-Relationship-Based Study

The acquisition of penicillin-binding protein (PBP) 2a in resistant strains of Staphylococcus aureus allows for the continuous production of cell walls even after the inactivation of intrinsic PBPs. Thus, the discovery of novel therapeutics with enhanced modulatory activity on PBP2a is crucial, and plant secondary metabolites, such as phenolics, have found relevance in this regard. In this study, using computational techniques, phenolics were screened against the active site of PBP2a, and the ability of the lead phenolics to modulate PBP2a’s active and allosteric sites was studied. The top-five phenolics (leads) identified through structure–activity-based screening, pharmacokinetics and synthetic feasibility evaluations were subjected to molecular dynamics simulations. Except for propan-2-one at the active site, the leads had a higher binding free energy at both the active and allosteric sites of PBP2a than amoxicillin. The leads, while promoting the thermodynamic stability of PBP2a, showed a more promising affinity at the allosteric site than the active site, with silicristin (−25.61 kcal/mol) and epicatechin gallate (−47.65 kcal/mol) having the best affinity at the active and allosteric sites, respectively. Interestingly, the modulation of Tyr446, the active site gatekeeper residue in PBP2a, was noted to correlate with the affinity of the leads at the allosteric site. Overall, these observations point to the leads’ ability to inhibit PBP2a, either directly or through allosteric modulation with conventional drugs. Further confirmatory in vitro studies on the leads are underway.

Identification of Flavonoid C-Glycosides as Promising Antidiabetics Targeting Protein Tyrosine Phosphatase 1B

Protein tyrosine phosphatase 1B (PTP1B), a negative regulator of the insulin signaling pathway, has gained attention as a validated druggable target in the management of type 2 diabetes mellitus (T2DM). The lack of clinically approved PTP1B inhibitors has continued to prompt research in plant-derived therapeutics possibly due to their relatively lesser toxicity profiles. Flavonoid C-glycosides are one of the plant-derived metabolites gaining increased relevance as antidiabetic agents, but their possible mechanism of action remains largely unknown. This study investigates the antidiabetic potential of flavonoid C-glycosides against PTP1B in silico and in vitro. Of the seven flavonoid C-glycosides docked against the enzyme, three compounds (apigenin, vitexin, and orientin) had the best affinity for the enzyme with a binding score of -7.3 kcal/mol each, relative to -7.4 kcal/mol for the reference standard, ursolic acid. A further probe (in terms of stability, flexibility, and compactness) of the complexes over a molecular dynamics time study of 100 ns for the three compounds suggested orientin as the most outstanding inhibitor of PTP1B owing to its overall -34.47 kcal/mol binding energy score compared to ursolic acid (-19.24 kcal/mol). This observation was in accordance with the in vitro evaluation result, where orientin had a half maximal inhibitory concentration (IC₅₀) of 0.18 mg/ml relative to 0.13 mg/ml for the reference standard. The kinetics of inhibition of PTP1B by orientin was mixed-type with V _max and K _m values of 0.004 μM/s and 0.515 μM. Put together, the results suggest orientin as a potential PTP1B inhibitor and could therefore be further explored in the management T2DM as a promising therapeutic agent.