Molecular Modeling Study for the Design of Novel Peroxisome Proliferator-Activated Receptor Gamma Agonists using 3D-QSAR and Molecular Docking

Identification of the structural features of quinazoline derivatives as EGFR inhibitors using 3D-QSAR modeling, molecular docking, molecular dynamics simulations and free energy calculations

Epidermal growth factor receptor (EGFR) is a promising target for the treatment of different types of malignant tumors. Therefore, a combined molecular modeling study was performed on a series of quinazoline derivatives as EGFR inhibitors. The optimum ligand-based CoMFA and CoMSIA models showed reliable and satisfactory predictability (with R²_cv=0.681, R²_ncv=0.844, R²_pred=0.8702 and R²_cv=0.643, R²_ncv=0.874, R²_pred=0.6423). The derived contour maps provide structural features to improve inhibitory activity. Furthermore, the contour maps, molecular docking, and molecular dynamics (MD) simulations have good consistency, illustrating that the derived models are reliable. In addition, MD simulations and binding free energy calculations were also carried out to understand the conformational fluctuations at the binding pocket of the receptor. The results indicate that hydrogen bond, hydrophobic and electrostatic interactions play significant roles on activity and selectivity. Furthermore, amino acids Val31, Lys50, Thr95, Leu149 and Asp160 are considered as essential residues to participate in the ligand-receptor interactions. Overall, this work would offer reliable theoretical basis for future structural modification, design and synthesis of novel EGFR inhibitors with good potency.Communicated by Ramaswamy H. Sarma.

Identification of novel inhibitors of high affinity iron permease (FTR1) through implementing pharmacokinetics index to fight against black fungus: An in silico approach

Mucormycosis is a life-threatening fungal infection, particularly in immunocompromised patients. Mucormycosis has been reported to show resistance to available antifungal drugs and was recently found in COVID-19 as a co-morbidity that demands new classes of drugs. In an attempt to find novel inhibitors against the high-affinity iron permease (FTR1), a novel target having fundamental importance on the pathogenesis of mucormycosis, 11,000 natural compounds were investigated in this study. Virtual screening and molecular docking identified two potent natural compounds [6',7,7,10',10',13'-hexamethylspiro[1,8-dihydropyrano[2,3-g]indole-3,11'-3,13-diazatetracyclo[5.5.2.01,9.03,7]tetradecane]-2,9,14'-trione and 5,7-dihydroxy-3-(2,2,8,8-tetramethylpyrano[2,3-f]chromen-6-yl)chromen-4-one] that effectively bind to the active cavity of FTR1 with a binding affinity of -9.9 kcal/mol. Multiple non-covalent interactions between the compounds and the active residues of this cavity were noticed, which is required for FTR1 inhibition. These compounds were found to have inhibitory nature and meet essential requirements to be drug-like compounds with a considerable absorption, distribution, metabolism, and excretion (ADME) profile with no toxicity probabilities. Molecular dynamics simulation confirms the structural compactness and less conformational variation of the drug-protein complexes maintaining structural stability and rigidity. MM-PBSA and post-simulation analysis predict binding stability of these compounds in the active cavity. This study hypothesizing that these compounds could be a potential inhibitor of FTR1 and will broaden the clinical prospects of mucormycosis.

Phytochemical Analysis of the Methanolic Extract and Essential Oil from Leaves of Industrial Hemp Futura 75 Cultivar: Isolation of a New Cannabinoid Derivative and Biological Profile Using Computational Approaches

Cannabis sativa L. is a plant belonging to the Cannabaceae family, cultivated for its psychoactive cannabinoid (Δ⁹-THC) concentration or for its fiber and nutrient content in industrial use. Industrial hemp shows a low Δ⁹-THC level and is a valuable source of phytochemicals, mainly represented by cannabinoids, flavones, terpenes, and alkaloids, with health-promoting effects. In the present study, we investigated the phytochemical composition of leaves of the industrial hemp cultivar Futura 75, a monoecious cultivar commercially used for food preparations or cosmetic purposes. Leaves are generally discarded, and represent waste products. We analyzed the methanol extract of Futura 75 leaves by HPLC and NMR spectroscopy and the essential oil by GC-MS. In addition, in order to compare the chemical constituents, we prepared the water infusion. One new cannabinoid derivative (1) and seven known components, namely, cannabidiol (2), cannabidiolic acid (3), β-cannabispirol (4), β-cannabispirol (5), canniprene (6), cannabiripsol (7), and cannflavin B (8) were identified. The content of CBD was highest in all preparations. In addition, we present the outcomes of a computational study focused on elucidating the role of 2α-hydroxy-Δ^3,7-cannabitriol (1), CBD (2), and CBDA (3) in inflammation and thrombogenesis.

Application of Various Molecular Modelling Methods in the Study of Estrogens and Xenoestrogens

In this review, applications of various molecular modelling methods in the study of estrogens and xenoestrogens are summarized. Selected biomolecules that are the most commonly chosen as molecular modelling objects in this field are presented. In most of the reviewed works, ligand docking using solely force field methods was performed, employing various molecular targets involved in metabolism and action of estrogens. Other molecular modelling methods such as molecular dynamics and combined quantum mechanics with molecular mechanics have also been successfully used to predict the properties of estrogens and xenoestrogens. Among published works, a great number also focused on the application of different types of quantitative structure-activity relationship (QSAR) analyses to examine estrogen's structures and activities. Although the interactions between estrogens and xenoestrogens with various proteins are the most commonly studied, other aspects such as penetration of estrogens through lipid bilayers or their ability to adsorb on different materials are also explored using theoretical calculations. Apart from molecular mechanics and statistical methods, quantum mechanics calculations are also employed in the studies of estrogens and xenoestrogens. Their applications include computation of spectroscopic properties, both vibrational and Nuclear Magnetic Resonance (NMR), and also in quantum molecular dynamics simulations and crystal structure prediction. The main aim of this review is to present the great potential and versatility of various molecular modelling methods in the studies on estrogens and xenoestrogens.

Structure-based drug design of peroxisome proliferator-activated receptor gamma inhibitors: ferulic acid and derivatives

Peroxisome proliferator-activated receptor gamma (PPARγ), ligand-activated transcription factor, is a key modulator of genes considered in diabetes development as well as treatment. Adipogenesis differentiation through PPARγ, CCAAT-enhancer protein alpha (C/EBPα) is identified as a critical mechanism in fat accumulation and weight gain. Polyphenols studied against adipocyte differentiation is taken up for consistent support and drug discovery. Structure-based drug design found useful to distinguish the underlying mechanism of receptor-ligand interaction and function. In this work, phenolic acids, ferulic acid and its derivatives are used as ligands. Molecular parameters have been set to filter and sort the 34 derivatives from ZINC and PubChem databases. Besides, for affinity and activity identification, troglitazone and resveratrol co-crystallized ligands have been studied. Absorption, distribution, metabolism, elimination and toxicity, density functional theory, highest occupied molecular orbital-lowest unoccupied molecular orbital values and docking scores define the drug candidate as a potential inhibitor. Residues Ser 342 and Arg 280 bind with the ligands by forming hydrogen bonds and hydrophobic contacts. Based on the docking score, pharmacophore properties and functional energy values of the top six compounds are chosen for molecular dynamics and simulation. Consistency and stability maintained throughout the simulation up to 50 ns were observed. Free binding energy values and standard deviation of receptor and ligand calculated using molecular mechanics-generalized Born and surface area solvation method (MM_GBSA) is found significant. Therefore, ferulic acid derivatives and phenolic acids could be a potential inhibitor for adipocyte differentiation and lipid accumulation.Communicated by Ramaswamy H. Sarma.

The recent use of Swietenia mahagoni (L.) Jacq. as antidiabetes type 2 phytomedicine: A systematic review

Background

Natural resources provide more efficient and safer alternatives in managing diabetes compare to the synthetic oral anti diabetes (OAD). The plants not only have hypoglycemic effect, but also prevent its complications; in which no synthetic drugs provide of both properties. Among antidiabetes plants, mahogany seed (Swietenia macrophylla) has been used as traditional medicine in Indonesia and India, though most popular utilized as timber wood.

Methods

The present study was performed of chemotaxonomic approach to review its phytochemical and anti-diabetic properties of Swietenia mahagoni (L.) Jacq seed/bark/leaves. The qualitative systematic review (SR) was carried out by analysing indexed journals and peer reviewed of Swietenia and Swietenia spp from Scopus, PubMed, Medline, Google Scholar, and Research Gate. Data selection criteria are accordance to botany, phytochemistry, in vitro, in vivo, and clinical test of related subject. The keywords used for the search in the databases were Swietenia, Swietenia mahagony, diabetes, and diabetes plants.

Results

Swietenia mahagoni (L.) Jacq. extracts have shown in vitro, in vivo and limited clinically test of its anti-diabetic properties. Ethanolic/methanolic/aqueous/petroleum/n-hexane extracts of mahagonis's seed/bark or leaves have anti-diabetic activities comparable to the synthetic drug and observed no to relatively mild toxic effect. The hypoglycemic mechanism suggested via reducing blood glucose level, restoring liver and β-cells islet function (might) blocking epinephrine function, inhibiting of α-amylase and β-glucosidase, antioxidant and antihiperlipidemia. Phytochemical compounds of S. mahagoni consist of the phenolics (flavonoids (swietemacrophyllanin, catechins and epichatechins) and tannins), triterpenoids and tetranortriterpenoids (limonoids: mahonin, secomahoganin, swietmanins, swiemahogins, swietenine and swietenolide), saponins and alkaloids which are known as anti-diabetic bioactive principles.

Conclusion

S. mahagoni was potentially used and developed as an antidiabetes source. To use it as an antidiabetic further, more extensive clinical trials and biomarkers of active compounds determination are needed.

Virtual identification of novel peroxisome proliferator-activated receptor (PPAR) α/δ dual antagonist by 3D-QSAR, molecule docking, and molecule dynamics simulation

The therapeutic potential of PPARs antagonists extends beyond diabetes. PPARs antagonists represent a new drug class that holds promise as a broadly applicable therapeutic approach for cancer treatment. Thus, there is a strong need to develop a rational design strategy for creating PPARs antagonists. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR) models of PPARα receptor (CoMFA-1, q ² = 0.636, r ² = 0.953; CoMSIA-1, q ² = 0.779, r ² = 0.999) and PPARδ receptor (CoMFA-2, q ² = 0.624, r ² = 0.906; CoMSIA-2, q ² = 0.627, r ² = 0.959) were successfully constructed using 35 triazolone ring derivatives. Contour map analysis revealed that the electrostatic and hydrophobic fields played vital roles in the bioactivity of dual antagonists. Molecular docking studies suggested that the hydrogen bonding, electrostatic and hydrophobic interactions all influenced the binding of receptor-ligand complex. Based on the information obtained above, we designed a series of compounds. The docking results were mutually validated with 3D-QSAR results. Three-dimensional-QSAR and absorption, distribution, metabolism, excretion and toxicity (ADMET) predictions indicated that 19 newly designed compounds possessed excellent biological activity and physicochemical properties. In summary, this research could provide theoretical guidance for the structural optimization of novel PPARα and δ dual antagonists. Communicated by Ramaswamy H. Sarma.

Discovery of a Natural Syk Inhibitor from Chinese Medicine through a Docking-Based Virtual Screening and Biological Assay Study

Spleen tyrosine kinase (Syk) is a critical target protein for treating immunoreceptor signalling-mediated allergies. In this study, a virtual screening of an in-house Chinese medicine database followed by biological assays was carried out to identify novel Syk inhibitors. A molecular docking method was employed to screen for compounds with potential Syk inhibitory activity. Then, an in vitro kinase inhibition assay was performed to verify the Syk inhibitory activity of the virtual screening hits. Subsequently, a β-hexosaminidase release assay was conducted to evaluate the anti-mast cell degranulation activity of the active compounds. Finally, tanshinone I was confirmed as a Syk inhibitor (IC₅₀ = 1.64 μM) and exhibited anti-mast cell degranulation activity in vitro (IC₅₀ = 2.76 μM). Docking studies showed that Pro455, Gln462, Leu377, and Lys458 were key amino acid residues for Syk inhibitory activity. This study demonstrated that tanshinone I is a Syk inhibitor with mast cell degranulation inhibitory activity. Tanshinone I may be a potential lead compound for developing effective and safe Syk-inhibiting drugs.

Identification of a Novel PPAR-γ Agonist through a Scaffold Tuning Approach

Peroxisome proliferator-activated receptors (PPARs) are important targets in metabolic diseases including obesity, metabolic syndrome, diabetes, and non-alcoholic fatty liver disease. Recently, they have been highlighted as attractive targets for the treatment of cardiovascular diseases and chronic myeloid leukemia. The PPAR agonist structure is consists of a polar head, a hydrophobic tail, and a linker. Each part interacts with PPARs through hydrogen bonds or hydrophobic interactions to stabilize target protein conformation, thus increasing its activity. Acidic head is essential for PPAR agonist activity. The aromatic linker plays an important role in making hydrophobic interactions with PPAR as well as adjusting the head-to-tail distance and conformation of the whole molecule. By tuning the scaffold of compound, the whole molecule could fit into the ligand-binding domain to achieve proper binding mode. We modified indol-3-ylacetic acid scaffold to (indol-1-ylmethyl)benzoic acid, whereas 2,4-dichloroanilide was fixed as the hydrophobic tail. We designed, synthesized, and assayed the in vitro activity of novel indole compounds with (indol-1-ylmethyl)benzoic acid scaffold. Compound 12 was a more potent PPAR-γ agonist than pioglitazone and our previous hit compound. Molecular docking studies may suggest the binding between compound 12 and PPAR-γ, rationalizing its high activity.