Modelling the DFT structural and reactivity study of feverfew and evaluation of its potential antiviral activity against COVID-19 using molecular docking and MD simulations

A Comparison between Bulgarian Tanacetum parthenium Essential Oil from Two Different Locations

Tanacetum parthenium L. (Asteraceae) is a perennial herbaceous plant with a long-standing historical use in traditional medicine. Recently Tanacetum parthenium L. essential oil has been associated with a promising potential for future applications in the pharmaceutical industry, in the cosmetics industry, and in agriculture. Investigations on the essential oil (EO) have indicated antimicrobial, antioxidant, and repellent activity. The present study aimed to evaluate the chemical composition of Bulgarian T. parthenium essential oil from two different regions, to compare the results to those reported previously in the literature, and to point out some of its future applications. The essential oils of the air-dried flowering aerial parts were obtained by hydrodistillation using a Clevenger-type apparatus. The chemical composition was evaluated using gas chromatography with mass spectrometry (GC-MS). It was established that the oxygenated monoterpenes were the predominant terpene class, followed by the monoterpene hydrocarbons. Significant qualitative and quantitative differences between both samples were revealed. Camphor (50.90%), camphene (16.12%), and bornyl acetate (6.05%) were the major constituents in the feverfew EO from the western Rhodope Mountains, while in the EO from the central Balkan mountains camphor (45.54%), trans-chrysanthenyl acetate (13.87%), and camphene (13.03%) were the most abundant components.

Computational-based investigation of antioxidative potential polyphenolic compounds of Salvia officinalis L.: combined DFT and molecular docking approaches

CONTEXT

The antioxidant properties of the three polyphenolic compounds (carnosol, cirsiliol, and luteolin) of Salvia officinalis L. were investigated employing the density functional theory (DFT) calculations at the B3LYP of basis set at 6-311 +  + G (d, p) in order to evaluate their antioxidant activity. The enthalpies of reactions associated with the SET-PT, SPLET, and HAT mechanisms were analyzed in gas and in different solvents using the CPCM (conductor-like polarizable continuum) model. For all possible hydrogen donor sites, the corresponding parameters (BDE, AIP, PDE, PA, ETE, HOMOs, and LUMOs) and reactivity indices (IPE, EA, Χ, η, S, and ω) were also evaluated. The calculated results showed that derivatives 12-OH, 11-OH, 4'-OH, and 3'-OH had the lowest antioxidant activity. The results showed as well that carnosol, cirsiliol, and luteolin have higher reactivity compared to ascorbic acid and could be considered better antioxidants. According to research, the catechol group is crucial in influencing the studied compounds antioxidant activity. The theoretically predicted order of antioxidant efficiencies in this work agrees well with the QSAR (quantitative structure-activity relationship) data. The findings show that in the vacuum as well as benzene media. HAT would be the most effective mechanism; in contrast, the thermodynamic equilibrium approach in polar media is the SPLET mechanism. Likewise, the outcomes of the docking modeling confirm that the selected molecules have high inhibitory activity to glutathione-S-transferases (GSTs) receptors. Moreover, they have very important pharmacokinetic, chemical, and biological profiles. Finally, all the results show that the three natural molecules have good pharmacokinetic profiles, particularly the bioavailability and permeability toward biological membranes.

METHODS

The software packages used in this investigation are Gaussian 16, Discovery studio Visualizer, and AutoDock vina. The three compounds (carnosol, cirsiliol, and luteolin) of Salvia officinalis L. were optimized with DFT/B3LYP of basis set at 6-311 +  + G (d, p). The optimized structures were established via vibrational analysis (i.e., no imaginary frequencies in the frequency set). All enthalpies were zero-point (ZPE) corrected. Vibrational frequency calculations were performed at 298.15 K and 1 atmosphere pressure to determine the thermodynamic characteristics of the investigated reactions. The descriptors were associated with the antioxidant mechanisms for investigated molecules in vacuum and in various solvents. The molecular docking was used by AutoDock vina to estimate and evaluate the title compounds compatibility as potential antioxidant drugs utilizing appropriate receptor proteins. The solvation effect in the medium of benzene (ɛ = 2.27) and water (ɛ = 78.39) was taken into account. Furthermore, a methanol solvent (ɛ = 32.61) was also taken into consideration to compare with the empirical data.

Comprehensive quantum chemical study of the associative complex of para-aminobenzoic acid and 7-diethylamino 4-methyl coumarin by adsorption and aromatic bridges

CONTEXT

The present study accounts for the quantum chemical and nonlinear optical properties of the combination of para-aminobenzoic acid and 7-diethylamino 4-methyl coumarin. Three different complexes were designed, surface interaction (adsorption) and two by connecting both molecules with π-bridge benzene and biphenyl. The amino and carboxyl groups were observed to behave as strong donor and acceptor sites in all the complexes. The band gap of the adsorbed complex was found more suitable. The absorption wavelength and intensity both were seen to increase with the increase in the number of benzene rings in the π-bridge. The values of first- and second-order hyperpolarizability suggest the improved nonlinear optical responses of the introduced complexes. Additionally, the negative value of second-order hyperpolarizability suggests the possibility of the occurrence of reverse saturable absorption in these combinations. The reported work gives theoretical insights into the nonlinear optical properties of the combination of para-aminobenzoic acid and 7-diethylamino 4-methyl coumarin.

METHODS

The molecular modeling and the quantum chemical studies were performed with Gaussian software packages. The standard B3LYP-6-311++G(d,p) basis functionals were used for energy minimization and other spectral calculations. All the surface analyses reported here are obtained by employing Multiwfn software. The chemical reactivity was established by the global reactivity descriptors. The intramolecular interactions and charge localization were stated using inter-fragment charge transfer analysis.

Waste to Medicine: Evidence from Computational Studies on the Modulatory Role of Corn Silk on the Therapeutic Targets Implicated in Type 2 Diabetes Mellitus

Type 2 diabetes mellitus (T2DM) is characterized by insulin resistance and/or defective insulin production in the human body. Although the antidiabetic action of corn silk (CS) is well-established, the understanding of the mechanism of action (MoA) behind this potential is lacking. Hence, this study aimed to elucidate the MoA in different samples (raw and three extracts: aqueous, hydro-ethanolic, and ethanolic) as a therapeutic agent for the management of T2DM using metabolomic profiling and computational techniques. Ultra-performance liquid chromatography-mass spectrometry (UP-LCMS), in silico techniques, and density functional theory were used for compound identification and to predict the MoA. A total of 110 out of the 128 identified secondary metabolites passed the Lipinski's rule of five. The Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analysis revealed the cAMP pathway as the hub signaling pathway, in which ADORA1, HCAR2, and GABBR1 were identified as the key target genes implicated in the pathway. Since gallicynoic acid (-48.74 kcal/mol), dodecanedioc acid (-34.53 kcal/mol), and tetradecanedioc acid (-36.80 kcal/mol) interacted well with ADORA1, HCAR2, and GABBR1, respectively, and are thermodynamically stable in their formed compatible complexes, according to the post-molecular dynamics simulation results, they are suggested as potential drug candidates for T2DM therapy via the maintenance of normal glucose homeostasis and pancreatic β-cell function.

Revolutionizing anti-cancer drug discovery against breast cancer and lung cancer by modification of natural genistein: an advanced computational and drug design approach

Breast and lung cancer are two of the most lethal forms of cancer, responsible for a disproportionately high number of deaths worldwide. Both doctors and cancer patients express alarm about the rising incidence of the disease globally. Although targeted treatment has achieved enormous advancements, it is not without its drawbacks. Numerous medicines and chemotherapeutic drugs have been authorized by the FDA; nevertheless, they can be quite costly and often fall short of completely curing the condition. Therefore, this investigation has been conducted to identify a potential medication against breast and lung cancer through structural modification of genistein. Genistein is the active compound in Glycyrrhiza glabra (licorice), and it exhibits solid anticancer efficiency against various cancers, including breast cancer, lung cancer, and brain cancer. Hence, the design of its analogs with the interchange of five functional groups-COOH, NH2 and OCH3, Benzene, and NH-CH2-CH2-OH-have been employed to enhance affinities compared to primary genistein. Additionally, advanced computational studies such as PASS prediction, molecular docking, ADMET, and molecular dynamics simulation were conducted. Firstly, the PASS prediction spectrum was analyzed, revealing that the designed genistein analogs exhibit improved antineoplastic activity. In the prediction data, breast and lung cancer were selected as primary targets. Subsequently, other computational investigations were gradually conducted. The mentioned compounds have shown acceptable results for in silico ADME, AMES toxicity, and hepatotoxicity estimations, which are fundamental for their oral medication. It is noteworthy that the initial binding affinity was only -8.7 kcal/mol against the breast cancer targeted protein (PDB ID: 3HB5). However, after the modification of the functional group, when calculating the binding affinities, it becomes apparent that the binding affinities increase gradually, reaching a maximum of -11.0 and -10.0 kcal/mol. Similarly, the initial binding affinity was only -8.0 kcal/mol against lung cancer (PDB ID: 2P85), but after the addition of binding affinity, it reached -9.5 kcal/mol. Finally, a molecular dynamics simulation was conducted to study the molecular models over 100 ns and examine the stability of the docked complexes. The results indicate that the selected complexes remain highly stable throughout the 100-ns molecular dynamics simulation runs, displaying strong correlations with the binding of targeted ligands within the active site of the selected protein. It is important to further investigate and proceed to clinical or wet lab experiments to determine the practical value of the proposed compounds.

Cheminformatics bioprospection of selected medicinal plants metabolites against trypsin cleaved VP4 (spike protein) of rotavirus A

Rotaviruses have continued to be the primary cause of acute dehydrating diarrhoea in children under five years of age despite the global introduction of four World Health Organization (WHO) prequalified oral vaccines in over 106 countries. Currently, no medication is approved by the Food and Drug Administration (FDA) specifically for treating rotavirus A-induced diarrhoea. Consequently, it is important to focus on developing prophylactic and curative therapeutics to combat rotaviral infections. For the first time, this study computationally screened and identified metabolites from Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum as potential novel inhibitors with broad-spectrum activity against VP5* and VP8* (spike protein) of rotavirus A (RVA). The initial top 20 metabolites identified through molecular docking were further filtered using drug-likeness and pharmacokinetics parameters. The molecular properties of the resulting top-ranked compounds were predicted by conducting density functional theory (DFT) calculations, while molecular dynamics (MD) simulation revealed their thermodynamic compatibility with a significant affinity towards VP8* than VP5*. Except for ellagic acid (-11.78 kcal/mol), the lead compounds had higher binding free energy than the reference standard (VP5* (-11.81 kcal/mol), VP8* (-14.12 kcal/mol)) with 2SG (-20.98 kcal/mol) and apigenin-4'-glucoside (-23.56 kcal/mol) having the highest affinity towards VP5* and VP8*, respectively. Of all the top-ranked compounds, better broad-spectrum affinities for both VP5* and VP8* than tizoxanide were observed in 2SG (VP5* (-20.98 kcal/mol), VP8* (-20.13 kcal/mol)) and sericetin (VP5* (-20.46 kcal/mol), VP8* (-18.31 kcal/mol)). While the identified leads could be regarded as potential modulators of the investigated therapeutic targets for effective management of rotaviral infection, additional in vitro and in vivo evaluation is strongly recommended, and efforts are on-going in this regard.Communicated by Ramaswamy H. Sarma.

The therapeutic effect and mechanism of parthenolide in skeletal disease, cancers, and cytokine storm

Parthenolide (PTL or PAR) was first isolated from Magnolia grandiflora and identified as a small molecule cancer inhibitor. PTL has the chemical structure of C15H20O3 with characteristics of sesquiterpene lactones and exhibits the biological property of inhibiting DNA biosynthesis of cancer cells. In this review, we summarise the recent research progress of medicinal PTL, including the therapeutic effects on skeletal diseases, cancers, and inflammation-induced cytokine storm. Mechanistic investigations reveal that PTL predominantly inhibits NF-κB activation and other signalling pathways, such as reactive oxygen species. As an inhibitor of NF-κB, PTL appears to inhibit several cytokines, including RANKL, TNF-α, IL-1β, together with LPS induced activation of NF-κB and NF-κB -mediated specific gene expression such as IL-1β, TNF-α, COX-2, iNOS, IL-8, MCP-1, RANTES, ICAM-1, VCAM-1. It is also proposed that PTL could inhibit cytokine storms or hypercytokinemia triggered by COVID-19 via blocking the activation of NF-κB signalling. Understanding the pharmacologic properties of PTL will assist us in developing its therapeutic application for medical conditions, including arthritis, osteolysis, periodontal disease, cancers, and COVID-19-related disease.

Cyclohexane-1,3-dione Derivatives as Future Therapeutic Agents for NSCLC: QSAR Modeling, In Silico ADME-Tox Properties, and Structure-Based Drug Designing Approach

The abnormal expression of the c-Met tyrosine kinase has been linked to the proliferation of several human cancer cell lines, including non-small-cell lung cancer (NSCLC). In this context, the identification of new c-Met inhibitors based on heterocyclic small molecules could pave the way for the development of a new cancer therapeutic pathway. Using multiple linear regression (MLR)-quantitative structure-activity relationship (QSAR) and artificial neural network (ANN)-QSAR modeling techniques, we look at the quantitative relationship between the biological inhibitory activity of 40 small molecules derived from cyclohexane-1,3-dione and their topological, physicochemical, and electronic properties against NSCLC cells. In this regard, screening methods based on QSAR modeling with density-functional theory (DFT) computations, in silico pharmacokinetic/pharmacodynamic (ADME-Tox) modeling, and molecular docking with molecular electrostatic potential (MEP) and molecular mechanics-generalized Born surface area (MM-GBSA) computations were used. Using physicochemical (stretch-bend, hydrogen bond acceptor, Connolly molecular area, polar surface area, total connectivity) and electronic (total energy, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy levels) molecular descriptors, compound 6d is identified as the optimal scaffold for drug design based on in silico screening tests. The computer-aided modeling developed in this study allowed us to design, optimize, and screen a new class of 36 small molecules based on cyclohexane-1,3-dione as potential c-Met inhibitors against NSCLC cell growth. The in silico rational drug design approach used in this study led to the identification of nine lead compounds for NSCLC therapy via c-Met protein targeting. Finally, the findings are validated using a 100 ns series of molecular dynamics simulations in an aqueous environment on c-Met free and complexed with samples of the proposed lead compounds and Foretinib drug.

Study of nonlinear optical responses of phytochemicals of Clitoria ternatea by quantum mechanical approach and investigation of their anti-Alzheimer activity with in silico approach

Clitoria ternatea is a flowering plant with promising medicinal plants with a wide variety of active phytochemicals. The present study aimed at the computational investigation of the nonlinear optical (NLO) responses of the active phytochemicals of the Clitoria ternatea. The computational investigation of the NLO features was done by using the density functional theory (DFT) by B3LYP/6-311G + + (d, p) basis set. The structural parameters, mulliken charge distribution, and molecular electrostatic potential (MEP) surface clearly show the intramolecular charge transfer within Clitorin. The NLO properties were identified by computing the polarizability parameters. As the plant has high medicinal characteristics, the inhibiting properties of its phytochemicals were also investigated to combat Alzheimer disease (AD). The systematic in silico study identifies Clitorin as the most active and inhibiting phytochemicals of the plant. The results obtained from molecular dynamics (MD) simulation tell the stability of the complex and make it a fair selection as a drug-like molecule against AD. The cardio-toxicity analysis done for the Clitorin molecule verifies that it is harmless for the heart.

Supplementary information

The online version contains supplementary material available at 10.1007/s11224-022-01981-5.

Evaluation of apigenin-based biflavonoid derivatives as potential therapeutic agents against viral protease (3CLpro) of SARS-CoV-2 via molecular docking, molecular dynamics and quantum mechanics studies

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic COVID-19 disease that affects human respiratory function. Despite the scientific progression made in the development of the vaccine, there is an urgent need for the discovery of antiviral drugs for better performance at different stages of SARS-CoV-2 reproduction. The main protease (Mpro or 3CLpro) plays a pivotal role in the life cycle of the virus, making it an attractive target for the development of antiviral agents effective against the new strains of coronaviruses (CoVs). In this study, a series of apigenin-based natural biflavonoid derivatives as potential inhibitors of coronaviruses 3CLpro was investigated by in silico approaches. For this purpose, the molecular docking was performed to analyze the interaction of the natural biflavonoids with SARS-Cov-2 main protease and for further investigation, docking to the 3CLpro of SARS-CoV and MERS-CoV. Based on docking scores and comparison with the reference inhibitors (ritonavir and lopinavir), more than half of the biflavonoids had strong interactions with the residues of the binding pocket of the coronaviruses 3CLpro and exhibited better binding affinities toward the main protease than ritonavir and lopinavir. The top biflavonoids were further explored through molecular dynamics simulation, binding free energy calculation and residual energy contributions estimated by the MM-PBSA. Also, drug likeness property investigation by Swiss ADME tools and density functional theory (DFT) calculations were performed. The results confirmed that the 3CLpro-amentoflavone, 3CLpro-bilobetin, 3CLpro-ginkgetin, and 3CLpro-sotetsuflavone complexes possess a large amount of dynamic properties such as high stability, significant binding energy and fewer conformation fluctuations. Also, the pharmacokinetics and drug-likeness studies and HOMO-LUMO and DFT descriptor values indicated a promising result of the selected natural biflavonoids. Overall findings indicate that the apigenin-based biflavonoids may inhibit COVID-19 by significant interactions in the binding pocket and those results can pave the way in drug discovery although the effectiveness of these bioactive compounds should be further validated by in-vitro and in-vivo investigations. Communicated by Ramaswamy H. Sarma.