Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

Computational exploration of microsomal cytochrome P450 3A1 enzyme modulation by phytochemicals of Cichorium intybus L.: Insights into drug metabolism

Drug metabolism plays a crucial role in drug fate, including therapeutic inactivation or activation, as well as the formation of toxic compounds. This underscores the importance of understanding drug metabolism in drug discovery and development. Considering the substantial costs associated with traditional drug development methods, computational approaches have emerged as valuable tools for predicting the metabolic fate of drug candidates. With this in mind, the present study aimed to investigate the potential mechanisms underlying the modulation of microsomal cytochrome P450 3A1 (CYP3A1) enzyme activity by various phytochemicals found in Cichorium intybus L., commonly known as chicory. To achieve this goal, several in silico methods, including molecular docking and molecular dynamics (MD) simulation, were employed to explore computationally the microsomal CYP3A1 enzyme. Schrodinger software was utilized for the molecular docking study, which involved the interaction analysis between CYP3A1 and 28 phytoconstituents of Cichorium intybus. Virtual screening of 28 compounds from chicory led to the identification of the top five ranked compounds. These compounds were evaluated for drug-likeness properties, pharmacokinetic profiles, and predicted binding affinities to CYP3A1. Caffeoylshikimic acid and cichoric acid emerged as promising candidates due to their favorable characteristics, including good oral bioavailability and high binding affinities to CYP3A1. Molecular dynamics simulations were conducted to assess the stability of caffeoylshikimic acid within the CYP3A1 binding pocket. The results demonstrated that caffeoylshikimic acid maintained stable interactions with the enzyme throughout the simulation, suggesting its potential as an effective modulator of CYP3A1 activity. The findings of this study have the potential to provide valuable insights into the complex molecular mechanisms by which Cichorium intybus L. acts on hepatocytes and modulates CYP3A1 enzyme expression or activity. By elucidating the impact of these phytochemicals on drug metabolism, this research contributes to our understanding of how chicory may interact with drugs and influence their efficacy and safety profiles.

A cytochrome P450 insecticide detoxification mechanism is not conserved across the Megachilidae family of bees

Recent work has demonstrated that many bee species have specific cytochrome P450 enzymes (P450s) that can efficiently detoxify certain insecticides. The presence of these P450s, belonging or closely related to the CYP9Q subfamily (CYP9Q-related), is generally well conserved across the diversity of bees. However, the alfalfa leafcutter bee, Megachile rotundata, lacks CYP9Q-related P450s and is 170-2500 times more sensitive to certain insecticides than bee pollinators with these P450s. The extent to which these findings apply to other Megachilidae bee species remains uncertain. To address this knowledge gap, we sequenced the transcriptomes of four Megachile species and leveraged the data obtained, in combination with publicly available genomic data, to investigate the evolution and function of P450s in the Megachilidae. Our analyses reveal that several Megachilidae species, belonging to the Lithurgini, Megachilini and Anthidini tribes, including all species of the Megachile genus investigated, lack CYP9Q-related genes. In place of these genes Megachile species have evolved phylogenetically distinct CYP9 genes, the CYP9DM lineage. Functional expression of these P450s from M. rotundata reveal they lack the capacity to metabolize the neonicotinoid insecticides thiacloprid and imidacloprid. In contrast, species from the Osmiini and Dioxyini tribes of Megachilidae have CYP9Q-related P450s belonging to the CYP9BU subfamily that are able to detoxify thiacloprid. These findings provide new insight into the evolution of P450s that act as key determinants of insecticide sensitivity in bees and have important applied implications for pesticide risk assessment.

Design, synthesis, computational study and cytotoxic evaluation of some new quinazoline derivatives containing pyrimidine moiety

Quinazoline derivatives, as an important category of heterocyclic compounds, have received much attention for the design and development of new drugs due to their various pharmacological properties. Besides, there is a great deal of evidence showing pyrimidine analogs as anticancer agents. Thus, in the present study, for the design of new target compounds with cytotoxic activity, we focused on various quinazolinone and pyrimidine hybrids. A new series of quinazoline-pyrimidine hybrid derivatives (6a-6n) have been designed and synthesized as novel antiproliferative agents. All the synthesized compounds characterized based on their IR, NMR and Mass spectroscopic data. Antiproliferative activities of the new compounds were evaluated against three human cancer cell lines (MCF-7, A549, SW-480). The compounds were found to have appropriate potential with IC50 values ranging from 2.3 ± 5.91 to 176.5 ± 0.7 μM against the tested cell lines. Compound 6n exerted the highest antiproliferative activity with IC50 values of 5.9 ± 1.69 μM, 2.3 ± 5.91 μM and 5.65 ± 2.33 μM against A549, SW-480 and MCF-7 respectively. The results indicated that 6n could induce apoptosis in A549 cell line in a dose dependent manner and arrest in the S phase of cell cycle. Docking studies were also done to investigate the detailed binding pattern of the synthesized compounds against EGFR. Furthermore, molecular dynamic simulation and binding free energy calculation have been done to rescore initial docking pose of the synthesized compounds using ensemble-based MMGB/PBSA free energy method. According to the results, free energy calculation confirmed biological activity of compounds and also, Arg 817 and Lys 721 residues had the pivotal role in the high potency of 6n. Finally, the drug likeness and in silico ADME study were also predicted.

Toxic effects of Tripterygium glycoside tablets on the reproductive system of male rats by metabolomics, cytotoxicity, and molecular docking

BACKGROUND

Tripterygium glycoside tablets (TGT) is the most common preparation from Tripterygium wilfordii Hook F, which is widely used in clinical for treating rheumatoid arthritis (RA) and other autoimmune diseases. However, its serious reproductive toxicity limits its application.

PURPOSE

This study aimed to elucidate the toxic effects of TGT on the reproductive system of male RA rats and its potential toxic components and mechanism.

METHODS

Collagen-induced arthritis (CIA) rat model was established, and TGT suspension was given at low, medium, and high doses. Gonadal index, pathological changes, and the number of spermatogenic cells were used to evaluate the toxic effects of TGT on the reproductive system. Non-targeted metabolomics of testicular tissue was conducted by UHPLC-QTOF/MS. Combined with network toxicology, the key targets of TGT-induced reproductive toxicity were screened and RT-qPCR was used to validation. In vitro toxicity of 19 components of TGT was evaluated using TM3 and TM4 cell lines. Molecular docking was used to predict the interaction between toxic components and key targets.

RESULTS

TGT reduced testicular and epididymis weight. Pathology analysis showed a lot of deformed and atrophic spermatogenic tubules. The number of spermatogenic cells decreased significantly (P<0.0001). A total of 58 different metabolites including platelet-activating factor (PAF), lysophosphatidylcholine (Lyso PC), phosphatidylinositol (PI), glutathione (GSH), and adenosine monophosphate (AMP) were identified by testicular metabolomics. Glycerophospholipid metabolism, ether lipid metabolism, and glutathione metabolism were key pathways responsible for the reproductive toxicity of TGT. Ten key reproductive toxicity targets were screened by network toxicology. The cytotoxicity test showed that triptolide, triptonide, celastrol, and demethylzeylasteral could significantly reduce the viability of TM3 and TM4 cells. Alkaloids had no apparent toxic effects. Molecular docking showed that the four toxic components had a good affinity with 10 key targets. All binding energies were less than -7 kcal/mol. The RT-qPCR results showed the Cyp19a1 level was significantly up-regulated. Pik3ca and Pik3cg levels were significantly down-regulated.

CONCLUSION

Through testicular metabolomics, we found that TGT may cause reproductive toxicity through CYP19A1, PIK3CA, and PIK3CG three target, which was preliminarily revealed. This study laid the foundation for elucidating the toxicity mechanism of TGT and evaluating its safety and quality.

Study of tyramine-binding mechanism and insecticidal activity of oil extracted from Eucalyptus against Sitophilus oryzae

The rice weevil, Sitophilus oryzae (L.), is a major pest of stored grains throughout the world, which causes quantitative and qualitative losses of food commodities. Eucalyptus essential oils (EOs) possess insecticidal and repellent properties, which make them a potential option for insect control in stored grains with environmentally friendly properties. In the current study, the binding mechanism of tyramine (TA) as a control compound has been investigated by funnel metadynamics (FM) simulation toward the homology model of tyramine1 receptor (TyrR) to explore its binding mode and key residues involved in the binding mechanism. EO compounds have been extracted from the leaf and flower part of Eucalyptus camaldulensis and characterized by GC/MS, and their effectiveness has been evaluated by molecular docking and conventional molecular dynamic (CMD) simulation toward the TyrR model. The FM results suggested that Asp114 followed by Asp80, Asn91, and Asn427 are crucial residues in the binding and the functioning of TA toward TyrR in Sitophilus Oryzae. The GC/MS analysis confirmed a total of 54 and 31 constituents in leaf and flower, respectively, where most of the components (29) are common in both groups. This analysis also revealed the significant concentration of Eucalyptus and α-pinene in leaves and flower EOs. The docking followed by CMD was performed to find the most effective compound in Eucalyptus EOs. In this regard, butanoic acid, 3-methyl-, 3-methyl butyl ester (B12) and 2-Octen-1-ol, 3,7-dimethyl- (B23) from leaf and trans- β-Ocimene (G04) from flower showed the maximum dock score and binding free energy, making them the leading candidates to replace tyramine in TyrR. The MM-PB/GBSA and MD analysis proved that the B12 structure is the most effective compound in inhibition of TyrR.

Computer-Aided (In Silico) Modeling of Cytochrome P450-Mediated Food–Drug Interactions (FDI)

Modifications of the activity of Cytochrome 450 (CYP) enzymes by compounds in food might impair medical treatments. These CYP-mediated food–drug interactions (FDI) play a major role in drug clearance in the intestine and liver. Inter-individual variation in both CYP expression and structure is an important determinant of FDI. Traditional targeted approaches have highlighted a limited number of dietary inhibitors and single-nucleotide variations (SNVs), each determining personal CYP activity and inhibition. These approaches are costly in time, money and labor. Here, we review computational tools and databases that are already available and are relevant to predicting CYP-mediated FDIs. Computer-aided approaches such as protein–ligand interaction modeling and the virtual screening of big data narrow down hundreds of thousands of items in databanks to a few putative targets, to which the research resources could be further directed. Structure-based methods are used to explore the structural nature of the interaction between compounds and CYP enzymes. However, while collections of chemical, biochemical and genetic data are available today and call for the implementation of big-data approaches, ligand-based machine-learning approaches for virtual screening are still scarcely used for FDI studies. This review of CYP-mediated FDIs promises to attract scientists and the general public.

A Comprehensive Analysis of Human CYP3A4 Crystal Structures as a Potential Tool for Molecular Docking-Based Site of Metabolism and Enzyme Inhibition Studies

The notable ability of human liver cytochrome P450 3A4 (CYP3A4) to metabolize diverse xenobiotics encourages researchers to explore in-depth the mechanism of enzyme action. Numerous CYP3A4 protein crystal structures have been deposited in protein data bank (PDB) and are majorly used in molecular docking analysis. The quality of the molecular docking results depends on the three-dimensional CYP3A4 protein crystal structures from the PDB. Present review endeavors to provide a brief outline of some technical parameters of CYP3A4 PDB entries as valuable information for molecular docking research. PDB entries between 22 April 2004 and 2 June 2021 were compiled and the active sites were thoroughly observed. The present review identified 76 deposited PDB entries and described basic information that includes CYP3A4 from human genetic, Escherichia coli (E. coli) use for protein expression, crystal structure obtained from X-ray diffraction method, taxonomy ID 9606, Uniprot ID P08684, ligand–protein structure description, co-crystal ligand, protein site deposit and resolution ranges between 1.7[Formula: see text]Å and 2.95[Formula: see text]Å. The observation of protein–ligand interactions showed the various residues on the active site depending on the ligand. The residues Ala305, Ser119, Ala370, Phe304, Phe108, Phe213 and Phe215 have been found to frequently interact with ligands from CYP3A4 PDB. Literature surveys of 17 co-crystal ligands reveal multiple mechanisms that include competitive inhibition, noncompetitive inhibition, mixed-mode inhibition, mechanism-based inhibition, substrate with metabolite, inducer, or combination modes of action. This overview may help researchers choose a trustworthy CYP3A4 protein structure from the PDB database to apply the protein in molecular docking analysis for drug discovery.

Determination of Novel Anti-Cancer Agents by Targeting OGG1 Enzyme Using Integrated Bioinformatics Methods

The 8-oxoguanine DNA glycosylase (OGG1) enzyme is a key DNA glycosylase mediating the excision of 7,8-dihydro-8-oxoguanine (8-oxoG) from DNA molecule to the start base excision repair pathway. The OGG1 glycosylase function depletion has been seen to obstruct pathological conditions such as inflammation, A3 T-cell lymphoblastic acute leukemia growth, and neurodegenerative diseases, thus warranting OGG1 as an attractive anti-cancer enzyme. Herein, we employed several drug libraries intending to screen non-toxic inhibitory molecules against the active pocket of the enzyme that achieved stable binding mode in dynamics. Two anti-cancer compounds ([O-]C1=C(CC2=CC=CC=C2)SC(=[N+]1CC(=O)NC3=NC=C(CC4=CC=CC=C4)S3)S and CCCN(CCC)[S]-(=O)(=O)C1=CC=C(C=C1)C(=O)NNC2=NC3=CC=C(Br)C=C3C(=N2)C4=CC=CC=C4) from Selleckchem.com were identified to occupy the active pocket of OGG1 and bind with greater affinity than Control TH5487. The binding affinity of Top-1 is -11.6 kcal/mol while that of Top-2 is -10.7 kcal/mol in contrast to TH5487 Control (-9 kcal/mol). During molecular dynamic simulations versus time, the said compounds are tightly held by the enzyme with no minor structural deviations reported except flexible loops in particular those present at the N and C-terminal. Both the compounds produced extensive hydrophobic interactions with the enzyme along with stable hydrogen bonding. The docking and molecular dynamics simulations predictions were further validated by molecular mechanics with generalized Born and surface area solvation (MM/GBSA) and Poisson Boltzmann surface area (MM/PBSA), and WaterSwap binding energies that validated strong binding of the compounds to the enzyme. The MM/GBSA binding free energy for Top-1 complex is -28.10 kcal/mol, Top-2 complex is -50.14 kcal/mol) and Control is -46.91 kcal/mol while MM/PBSA value for Top-1, Top-2 and Control is -23.38 kcal/mol, -35.29 kcal/mol and -38.20 kcal/mol, respectively. Computational pharmacokinetics support good druglike candidacy of the compounds with acceptable profile of pharmacokinetics and very little toxicity. All these findings support the notion that the compounds can be used in experiments to test their anti-cancer activities.

Single Nucleotide Polymorphism Induces Divergent Dynamic Patterns in CYP3A5: A Microsecond Scale Biomolecular Simulation of Variants Identified in Sub-Saharan African Populations

Pharmacogenomics aims to reveal variants associated with drug response phenotypes. Genes whose roles involve the absorption, distribution, metabolism, and excretion of drugs, are highly polymorphic between populations. High coverage whole genome sequencing showed that a large proportion of the variants for these genes are rare in African populations. This study investigated the impact of such variants on protein structure to assess their functional importance. We used genetic data of CYP3A5 from 458 individuals from sub-Saharan Africa to conduct a structural bioinformatics analysis. Five missense variants were modeled and microsecond scale molecular dynamics simulations were conducted for each, as well as for the CYP3A5 wildtype and the Y53C variant, which has a known deleterious impact on enzyme activity. The binding of ritonavir and artemether to CYP3A5 variant structures was also evaluated. Our results showed different conformational characteristics between all the variants. No significant structural changes were noticed. However, the genetic variability seemed to act on the plasticity of the protein. The impact on drug binding might be drug dependant. We concluded that rare variants hold relevance in determining the pharmacogenomics properties of populations. This could have a significant impact on precision medicine applications in sub-Saharan Africa.

Computational Determination of Potential Multiprotein Targeting Natural Compounds for Rational Drug Design Against SARS-COV-2

SARS-CoV-2 caused the current COVID-19 pandemic and there is an urgent need to explore effective therapeutics that can inhibit enzymes that are imperative in virus reproduction. To this end, we computationally investigated the MPD3 phytochemical database along with the pool of reported natural antiviral compounds with potential to be used as anti-SARS-CoV-2. The docking results demonstrated glycyrrhizin followed by azadirachtanin, mycophenolic acid, kushenol-w and 6-azauridine, as potential candidates. Glycyrrhizin depicted very stable binding mode to the active pocket of the Mpro (binding energy, -8.7 kcal/mol), PLpro (binding energy, -7.9 kcal/mol), and Nucleocapsid (binding energy, -7.9 kcal/mol) enzymes. This compound showed binding with several key residues that are critical to natural substrate binding and functionality to all the receptors. To test docking prediction, the compound with each receptor was subjected to molecular dynamics simulation to characterize the molecule stability and decipher its possible mechanism of binding. Each complex concludes that the receptor dynamics are stable (Mpro (mean RMSD, 0.93 Å), PLpro (mean RMSD, 0.96 Å), and Nucleocapsid (mean RMSD, 3.48 Å)). Moreover, binding free energy analyses such as MMGB/PBSA and WaterSwap were run over selected trajectory snapshots to affirm intermolecular affinity in the complexes. Glycyrrhizin was rescored to form strong affinity complexes with the virus enzymes: Mpro (MMGBSA, -24.42 kcal/mol and MMPBSA, -10.80 kcal/mol), PLpro (MMGBSA, -48.69 kcal/mol and MMPBSA, -38.17 kcal/mol) and Nucleocapsid (MMGBSA, -30.05 kcal/mol and MMPBSA, -25.95 kcal/mol), were dominated mainly by vigorous van der Waals energy. Further affirmation was achieved by WaterSwap absolute binding free energy that concluded all the complexes in good equilibrium and stability (Mpro (mean, -22.44 kcal/mol), PLpro (mean, -25.46 kcal/mol), and Nucleocapsid (mean, -23.30 kcal/mol)). These promising findings substantially advance our understanding of how natural compounds could be shaped to counter SARS-CoV-2 infection.

Potential effective inhibitory compounds against Prostate Specific Membrane Antigen (PSMA): A molecular docking and molecular dynamics study

One of the most prevalent cancers in men is prostate cancer and could be managed with immunotoxins or antibody treatment. Because of the substantial rise of the Prostate-Specific Antigen and the Prostate-Specific Membrane Antigen (PSMA), cancer vaccination should be rendered with these antigens. Through pharmacodynamic experiments in a library of natural compounds from ZINC database, the current research sought to identify compounds that could suppress PSMA protein. To test the most productive compounds for further research, the Library has been scanned with Pharmacophore and ADMET analysis followed by molecular docking methods in the first phase. After selecting 15 ligands with the best pose related to docking results, to evaluate the stability of the ligand-protein bounds of the compounds, a molecular dynamics simulation considering the effect of the presence of zinc ions on the protein structure was performed. The measurement of ligand binding modes and free energy has shown that four compounds, including Z10, Z06, Z01, and Z03, have formed critical interactions with the active site's residues. Besides, multiple approaches were employed to determine their inhibition rating and describe the variables that facilitate the attachment of ligands to the protein active site. The results are obtained from the MMPBSA/GBSA analysis of four selected small molecules (Z10, Z06, Z01, and Z03), which are very close to the IC50 value of reference ligand (DCIBzl); they are -13.85 kcal/mol, -12.58 kcal/mol, -10.71 kcal/mol and -9.39 kcal/mol respectively. Finally, we evaluate the results obtained from selected ligands using hydrogen bond and decomposition analyzes. We have examined the effective interactions between ligands and S1/S1'pockets in protein. Our computational results illustrate the design of more efficient inhibitors of PSMA.

Conformational transition of Acinetobacter baumannii KdsC enzyme and the role of magnesium in binding: An insight from comparative molecular dynamics simulation and its implications in novel antibiotics design

The 3-deoxy-d-manno-octulosonate 8-phosphate phosphatase (KdsC) catalyzes the hydrolysis of 3-deoxy-d-manno-octulosonate 8-phosphate (KDO 8-P) to 3-deoxy-d-manno-octulosonate (KDO) and inorganic phosphate in KDO biosynthesis pathway of Gram-negative bacteria lipopolysaccharide (LPS) hydrophobic lipid-A core. The essentiality of KDO for bacterial cell viability presents the possibility of its targeting to develop broad-spectrum antibacterial agents. In this study, a receptor based virtually screening method was put forward to identify novel lead inhibitory molecules for KdsC enzyme. Dynamics evaluation in solution revealed three complexes: Asinex-1197, Asinex-1705, and Asinex-1710 from Asinex antibacterial library as highly stable, involving conformational transition from open to close upon lead molecules binding and eloquent role of active pocket magnesium towards inhibitors binding and movements. Interconversion of local secondary structure elements in sequence region of Asp192-Asp208 covering motif β-turn, β-hairpin, and β-sheets is seen recurrently that could be in all likelihood of the pressure excreted on this region during closing conformation event or magnesium driven inhibitor adjustments. The binding free energy estimation predicted gas phase energy for all the three complexes dominating with major contribution from van der Waals energy (in case of Asinex-1705 and Asinex-1710) and balanced contributions of both electrostatic and van der Waals (in case of Asinex-1197). Key residues-scanning shortlisted Leu45, Asp185, Gy188, Arg231, and Lys255 as vital in the interaction network of magnesium and inhibitors at the binding site. Their crucial roles in net binding energy were reaffirmed via in silico site directed alanine scanning method. The filtered hits might be useful to further scaffolds addition and structural optimization to yield high affinity binders of KdsC enzyme, whose inhibition, in turn, will disrupt the outer membrane synthesis.