Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Exploring natural products library to identify potential inhibitors targeting isoniazid-resistant tuberculosis

Mycobacterium tuberculosis (MTB) causing tuberculosis (TB) infection is a leading source of illness and death in developing nations, and the emergence of drug-resistant TB remains a significant global threat and a challenge in treating the disease. Mutations in the inhA and katG genes are connected to the principal molecular mechanism of isoniazid (INH) resistance, and continuous treatment of INH for more than a decade led to the evolution of INH resistant-TB (inhR-TB). Structure-based drug discovery approaches on traditional natural compounds are the contemporary source to identify significant lead molecules. This work focuses on discovering effective small compounds from natural compound libraries and applying pharmacophore-based virtual screening to filter out the molecules. The best-identified hit complexes were used for molecular dynamics simulations (MDS) to observe their stability and compactness. A three-dimensional e-pharmacophore hypothesis and screening generated 62 hits based on phase fitness scores from the pharmacophore-based virtual screening. Molecular docking experiments in Maestro's GLIDE module indicated that ZINC000002383126 and ASN22022 may be potential inhibitors of inhA and katG (native, inhA mutants S94A, Y158A, Y158F and Y158S and D137S, Y229F, S315T, W321F, and R418L mutants of katG). In addition, MDS analysis indicated that the native and mutant docked complexes of inhA and katG had good stability and remained compact in the binding pocket of the targets. In vitro studies can further validate the compounds that can act as INH competitive inhibitors.Communicated by Ramaswamy H. Sarma.

Role of Mitochondrial Cytochrome P450 2E1 in Healthy and Diseased Liver

Cytochrome P450 2E1 (CYP2E1) is pivotal in hepatotoxicity induced by alcohol abuse and different xenobiotics. In this setting, CYP2E1 generates reactive metabolites inducing oxidative stress, mitochondrial dysfunction and cell death. In addition, this enzyme appears to play a role in the progression of obesity-related fatty liver to nonalcoholic steatohepatitis. Indeed, increased CYP2E1 activity in nonalcoholic fatty liver disease (NAFLD) is deemed to induce reactive oxygen species overproduction, which in turn triggers oxidative stress, necroinflammation and fibrosis. In 1997, Avadhani’s group reported for the first time the presence of CYP2E1 in rat liver mitochondria, and subsequent investigations by other groups confirmed that mitochondrial CYP2E1 (mtCYP2E1) could be found in different experimental models. In this review, we first recall the main features of CYP2E1 including its role in the biotransformation of endogenous and exogenous molecules, the regulation of its expression and activity and its involvement in different liver diseases. Then, we present the current knowledge on the physiological role of mtCYP2E1, its contribution to xenobiotic biotransformation as well as the mechanism and regulation of CYP2E1 targeting to mitochondria. Finally, we discuss experimental investigations suggesting that mtCYP2E1 could have a role in alcohol-associated liver disease, xenobiotic-induced hepatotoxicity and NAFLD.

Prediction of potential inhibitors of the dimeric SARS-CoV2 main proteinase through the MM/GBSA approach

Since the emergence of SARS-CoV2, to date, no effective antiviral drug has been approved to treat the disease, and no vaccine against SARS-CoV2 is available. Under this scenario, the combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has attracted attention since they have been previously employed against the SARS-CoV main proteinase (Mpro) and exhibited some signs of effectiveness. Recently, the 3D structure of SARS-CoV2 Mpro was constructed based on the monomeric SARS-CoV Mpro and employed to identify potential approved small inhibitors against SARS-CoV2 Mpro, allowing the selection of 15 drugs among 1903 approved drugs to be employed. In this study, we performed docking of these 15 approved drugs against the recently solved X-ray crystallography structure of SARS-CoV2 Mpro in the monomeric and dimeric states; the latter is the functional state that was determined in a biological context, and these were submitted to molecular dynamics (MD) simulations coupled with the molecular mechanics generalized Born surface area (MM/GBSA) approach to obtain insight into the inhibitory activity of these compounds. Similar studies were performed with lopinavir and ritonavir coupled to monomeric and dimeric SARS-CoV Mpro and SARS-CoV2 Mpro to compare the inhibitory differences. Our study provides the structural and energetic basis of the inhibitory properties of lopinavir and ritonavir on SARS-CoV Mpro and SARS-CoV2 Mpro, allowing us to identify two FDA-approved drugs that can be used against SARS-CoV2 Mpro. This study also demonstrated that drug discovery requires the dimeric state to obtain good results.

Identification of saquinavir as a potent inhibitor of dimeric SARS-CoV2 main protease through MM/GBSA

Among targets selected for studies aimed at identifying potential inhibitors against COVID-19, SARS-CoV2 main proteinase (M^pro) is highlighted. M^pro is indispensable for virus replication and is a promising target of potential inhibitors of COVID-19. Recently, monomeric SARS-CoV2 M^pro, drug repurposing, and docking methods have facilitated the identification of several potential inhibitors. Results were refined through the assessment of dimeric SARS-CoV2 M^pro, which represents the functional state of enzyme. Docking and molecular dynamics (MD) simulations combined with molecular mechanics/generalized Born surface area (MM/GBSA) studies indicated that dimeric M^pro most significantly impacts binding affinity tendency compared with the monomeric state, which suggests that dimeric state is most useful when performing studies aimed at identifying drugs targeting M^pro. In this study, we extend previous research by performing docking and MD simulation studies coupled with an MM/GBSA approach to assess binding of dimeric SARS-CoV2 M^pro to 12 FDA-approved drugs (darunavir, indinavir, saquinavir, tipranavir, diosmin, hesperidin, rutin, raltegravir, velpatasvir, ledipasvir, rosuvastatin, and bortezomib), which were identified as the best candidates for the treatment of COVID-19 in some previous dockings studies involving monomeric SARS-CoV2 M^pro. This analysis identified saquinavir as a potent inhibitor of dimeric SARS-CoV2 M^pro; therefore, the compound may have clinical utility against COVID-19.

Exploring the biotransformation of N-(2-hydroxyphenyl)-2-propylpentanamide (an aryl valproic acid derivative) by CYP2C11, using in silico predictions and in vitro studies

OBJECTIVES

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA), a derivative of valproic acid (VPA), has been proposed as a potential anticancer agent due to its improved antiproliferative effects in some cancer cell lines. Although there is evidence that VPA is metabolized by cytochrome P450 2C11 rat isoform, HO-AAVPA CYP-mediated metabolism has not yet been fully explored. Therefore, in this work, the biotransformation of HO-AAVPA by CYP2C11 was investigated.

METHODS

Kinetic parameters and spectral interaction between HO-AAVPA and CYP were evaluated using rat liver microsomes. The participation of CYP2C11 in metabolism of HO-AAVPA was confirmed by cimetidine (CIM) inhibition assay. Docking and molecular dynamics simulations coupled to MMGBSA methods were used in theoretical study.

KEY FINDINGS

HO-AAVPA is metabolized by CYP enzymes (K_M  = 38.94 µm), yielding a hydroxylated metabolite according to its HPLC retention time (5.4 min) and MS analysis (252.2 m/z). In addition, CIM inhibition in rat liver microsomes (K_i  = 59.23 µm) confirmed that CYP2C11 is mainly involved in HO-AAVPA metabolism. Furthermore, HO-AAVPA interacts with CYP2C11 as a type I ligand. HO-AAVPA is stabilized at the CYP2C11 ligand recognition site through a map of interactions similar to other typical CYP2C11 substrates.

CONCLUSION

Therefore, rat liver CYP2C11 isoform is able to metabolize HO-AAVPA.

Exploring the inhibitory activity of valproic acid against the HDAC family using an MMGBSA approach

Valproic acid (VPA) is a compound currently used in clinical practice for the treatment of epilepsy as well as bipolar and mood disorders. VPA targets histone deacetylases (HDACs), which participate in the removal of acetyl groups from lysine in several proteins, regulating a wide variety of functions within the organism. An imbalance or malfunction of these enzymes is associated with the development and progression of several diseases, such as cancer and neurodegenerative diseases. HDACs are divided into four classes, but VPA only targets Class I (HDAC1-3 and 8) and Class IIa (HDAC4-5, 7 and 9) HDACs; however, structural and energetic information regarding the manner by which VPA inhibits these HDACs is lacking. Here, the structural and energetic features that determine this recognition were studied using molecular docking and molecular dynamics (MD) simulation. It was found that VPA reaches the catalytic site in HDAC1-3 and 7, whereas in HDAC6, VPA only reaches the catalytic tunnel. In HDAC4, VPA was bound adjacent to L1 and L2, a zone that participates in corepressor binding, and in HDAC8, VPA was bound to the hydrophobic active site channel (HASC), in line with previous reports.

Metatox - Web application for generation of metabolic pathways and toxicity estimation

Xenobiotics biotransformation in humans is a process of the chemical modifications, which may lead to the formation of toxic metabolites. The prediction of such metabolites is very important for drug development and ecotoxicology studies. We created the web-application MetaTox ( http://way2drug.com/mg ) for the generation of xenobiotics metabolic pathways in the human organism. For each generated metabolite, the estimations of the acute toxicity (based on GUSAR software prediction), organ-specific carcinogenicity and adverse effects (based on PASS software prediction) are performed. Generation of metabolites by MetaTox is based on the fragments datasets, which describe transformations of substrates structures to a metabolites structure. We added three new classes of biotransformation reactions: Dehydrogenation, Glutathionation, and Hydrolysis, and now metabolite generation for 15 most frequent classes of xenobiotic's biotransformation reactions are available. MetaTox calculates the probability of formation of generated metabolite - it is the integrated assessment of the biotransformation reactions probabilities and their sites using the algorithm of PASS ( http://way2drug.com/passonline ). The prediction accuracy estimated by the leave-one-out cross-validation (LOO-CV) procedure calculated separately for the probabilities of biotransformation reactions and their sites is about 0.9 on the average for all reactions.

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

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

Impact of tetramerization on the ligand recognition of N1 influenza neuraminidase via MMGBSA approach

Influenza virus neuraminidase (NA) is a homotetrameric surface protein that, in contrast to other non-influenza NAs, requires a quaternary assembly to exhibit enzymatic activity, suggesting that the oligomeric state significantly impacts the active site of influenza NA. Nevertheless, most structure-based drug design studies have been reported by employing the monomeric state in the closed or open-loop due to the computational cost of employing the tetrameric NA. In this work, we present MD simulations coupled to the MMGBSA approach of avian N1 type NA in its monomeric and tetrameric closed and open-loop state both with and without the inhibitor oseltamivir and its natural substrate, sialic acid. Structural and energetic analyses revealed that the tetrameric state impacts flexibility as well as the map of interactions participating in stabilizing the protein-ligand complexes with respect to the monomeric state. It was observed that the tetrameric state exerts dissimilar effects in binding affinity, characteristic of positive and negative cooperativity for oseltamivir and sialic acid, respectively. Based on our results, to perform a confident structure-based drug design, as well as to evaluate the impact of key mutations through MD simulations, it is important to consider the tetrameric state closed-loop state.

Molecular recognition between potential natural inhibitors of the Keap1-Nrf2 complex

Disrupting the Keap1-Nrf2 pathway enhances Nrf2 activity, which has been identified as an important approach for the prevention of different chronic diseases in which oxidative stress and inflammation are present, such as cancer, diabetes, Alzheimer's and Parkinson's. Based on the high potential to modulate antioxidant, anti-inflammatory and anticancer properties that the discovery of Keap1-Nrf2 protein-protein interaction inhibitors would represent, the utilization of some natural compounds has emerged as a promising strategy to identify new drugs. To gain insight into the structural and energetic basis of the molecular recognition between some natural inhibitors that could work as inhibitors of the Keap1-Nrf2 complex, we evaluated the binding properties between four natural compounds present in the extract of Geranium schiedeanum (Gs): 3-O-a-L arabinofuranoside-7-O-a-l-rhamnopyranoside of kaempferol (KAM), gallic acid (GAL), ellagic acid (ELL) and geranium acetonitrile (ACE), which based on experimental findings have been proposed as possible Keap1-Nrf2 PPI inhibitors. Computational studies combining docking and MD simulations accompanied by the MMGBSA approach revealed that KAM and ACE directly interact with residues in the Kelch domain that participate in the molecular recognition of Nrf2, indicating that both natural compounds could act as activators of Nrf2, whereas GAL and ELL are possible free radical scavengers.

P450 catalysed dehydrogenation

Cytochrome P450s belong to a superfamily of enzymes that catalyse a wide variety of oxidative transformations. Hydroxylation is one the most thoroughly investigated of all identified P450-catalysed reactions whilst dehydrogenation has been relatively much less explored to date. P450-catalysed dehydrogenation is often found to occur with hydroxylation and thus, it was initially suspected to be a stepwise process consisting of hydroxylation and subsequent dehydration to yield the final olefin product. This theory has been proven to be invalid and the olefin was shown to be the direct product of a P450-catalysed reaction. This interesting reaction plays a vital role in the metabolism of xenobiotics and the biosynthesis of endogenous compounds, including a number of steroids. A number of well-known examples of P450 mediated dehydrogenation, including those in the metabolism of valproic acid, capsaicin and 3-methylindole and those in the biosynthesis of plant and fungal sterols are discussed in this review.

Molecular Dynamics Assisted Mechanistic Study of Isoniazid-Resistance against Mycobacterium tuberculosis InhA

Examination of InhA mutants I16T, I21V, I47T, S94A, and I95P showed that direct and water mediated H-bond interactions between NADH and binding site residues reduced drastically. It allowed conformational flexibility to NADH, particularly at the pyrophosphate region, leading to weakening of its binding at dinucleotide binding site. The highly scattered distribution of pyrophosphate dihedral angles and chi1 side chain dihedral angles of corresponding active site residues therein confirmed weak bonding between InhA and NADH. The average direct and water mediated bridged H-bond interactions between NADH and mutants were observed weaker as compared to wild type. Further, estimated NADH binding free energy in mutants supported the observed weakening of InhA-NADH interactions. Similarly, per residue contribution to NADH binding was also found little less as compared to corresponding residues in wild type. This investigation clearly depicted and supported the effect of mutations on NADH binding and can be accounted for isoniazid resistance as suggested by previous biochemical and mutagenic studies. Further, structural analysis of InhA provided the crucial points to enhance the NADH binding affinity towards InhA mutants in the presence of direct InhA inhibitors to combat isoniazid drug resistance. This combination could be a potential alternative for treatment of drug resistant tuberculosis.

Interaction mechanism exploration of R-bicalutamide/S-1 with WT/W741L AR using molecular dynamics simulations

R-Bicalutamide is a first generation antiandrogen used to treat prostate cancer, which inhibits androgen action by competitively binding to the androgen receptor (AR). However, R-bicalutamide was discovered to exhibit some agonistic properties in clinical application. According to reports, the W741L AR mutation may lead to resistance towards R-bicalutamide. But the mechanism of the R-bicalutamide switch from an antagonist to an agonist due to the mutation of AR W741L is still not so clear. Another molecule, S-1, owing to a very similar structure to R-bicalutamide, is always agonistic to both the wild type and W741L AR. The main difference between these two chemicals is that S-1 has an ether linkage while R-bicalutamide has a sulfonyl group. To study the drug-resistant mechanism caused by W741L mutation and the opposite effects arising from subtle structure differences, molecular dynamics (MD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) calculations were employed to explore the interaction mechanisms between R-bicalutamide/S-1 and WT/W741L AR. The calculated binding free energies are in accordance with the reported experimental values. The obtained results indicate that M895 and W741 are vital amino acids in the antagonism of R-bicalutamide. The bulkier substitution of sulfonyl and tryptophan push aside M895, together with helix 12 (H12), to expose the ligand-binding domain resulting in the antagonistic conformation of the AR. If W741 is mutated to L741, the B-ring of these two chemicals would shift toward L741. At the same time, M895 dragging helix H12, would also move closer to L741. So H12 tends to cover the AR ligand-binding domain to a certain degree, changing the androgen receptor from an antagonistic to an agonistic conformation, which may explain the agonism of R-bicalutamide to the mutant W741L AR.

A Combined Molecular Docking/Dynamics Approach to Probe the Binding Mode of Cancer Drugs with Cytochrome P450 3A4

Cytarabine, daunorubicin, doxorubicin and vincristine are clinically used for combinatorial therapies of cancers in different combinations. However, the knowledge about the interaction of these drugs with the metabolizing enzyme cytochrome P450 is limited. Therefore, we utilized computational methods to predict and assess the drug-binding modes. In this study, we performed docking, MD simulations and free energy landscape analysis to understand the drug-enzyme interactions, protein domain motions and the most populated free energy minimum conformations of the docked protein-drug complexes, respectively. The outcome of docking and MD simulations predicted the productive, as well as the non-productive binding modes of the selected drugs. Based on these interaction studies, we observed that S119, R212 and R372 are the major drug-binding residues in CYP3A4. The molecular mechanics Poisson–Boltzmann surface area analysis revealed the dominance of hydrophobic forces in the CYP3A4-drug association. Further analyses predicted the residues that may contain favorable drug-specific interactions. The probable binding modes of the cancer drugs from this study may extend the knowledge of the protein-drug interaction and pave the way to design analogs with reduced toxicity. In addition, they also provide valuable insights into the metabolism of the cancer drugs.