Characterization of the CYP3A4 active site by homology modeling

Magnetic chitosan-copper nanocomposite: A plant assembled catalyst for the synthesis of amino- and N-sulfonyl tetrazoles in eco-friendly media

A greener, cost efficient and simple method is described to prepare copper nanoparticles (NPs) immobilized on the magnetic chitosan (one of the more versatile polysaccharides) using Euphorbia falcata leaf extract as reducing/stabilizing agent. The prepared catalyst (Cu NPs@Fe₃O₄-chitosan) was authenticated by field emission scanning electron microscope (FESEM), X-ray diffraction (XRD), transmission electron microscopy (TEM), Scanning transmission electron microscopy (STEM), energy-dispersive X-ray spectroscopy (EDS), thermogravimetry/derivative thermogravimetry (TG/DTG), Vibrating sample magnetometer (VSM), and elemental mapping. TEM analysis indicates that Cu NPs with average sizes in 5-10 nm range is formed on magnetic chitosan with the spherical morphology. The Cu NPs@Fe₃O₄-chitosan was employed as a new catalyst for the synthesis of different tetrazoles by the reaction of various secondary or tertiary cyanamides with sodium azide in water under reflux conditions. Easy separation by external magnetic field, mild reaction conditions, low cost and the reusability are some of the beneficial features of this catalyst.

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.

Characterization of the active site properties of CYP4F12

Cytochrome P450 4F12 is a drug-metabolizing enzyme that is primarily expressed in the liver, kidney, colon, small intestine, and heart. The properties of CYP4F12 that may impart an increased catalytic selectivity (decreased promiscuity) were explored through in vitro metabolite elucidation, kinetic isotope effect experiments, and computational modeling of the CYP4F12 active site. By using astemizole as a probe substrate for CYP4F12 and CYP3A4, it was observed that although CYP4F12 favored astemizole O-demethylation as the primary route of metabolism, CYP3A4 was capable of metabolizing astemizole at multiple sites on the molecule. Deuteration of astemizole at the site of O-demethylation resulted in an isotope effect of 7.1 as well as an 8.3-fold decrease in the rate of clearance for astemizole by CYP4F12. Conversely, although an isotope effect of 3.8 was observed for the formation of the O-desmethyl metabolite when deuterated astemizole was metabolized by CYP3A4, there was no decrease in the clearance of astemizole. Development of a homology model of CYP4F12 based on the crystal structure of cytochrome P450 BM3 predicted an active site volume for CYP4F12 that was approximately 76% of the active site volume of CYP3A4. As predicted, multiple favorable binding orientations were available for astemizole docked into the active site of CYP3A4, but only a single binding orientation with the site of O-demethylation oriented toward the heme was identified for CYP4F12. Overall, it appears that although CYP4F12 may be capable of binding similar ligands to other cytochrome P450 enzymes such as CYP3A4, the ability to achieve catalytically favorable orientations may be inherently more difficult because of the increased steric constraints of the CYP4F12 active site.

The effect of chlorpyrifos-oxon and other xenobiotics on the human cytochrome P450-dependent metabolism of naphthalene and deet

Chlorpyrifos-oxon (CPO), a metabolite of chlorpyrifos, is a potent inhibitor of acetylcholinesterase and, although the neurotoxicological impact of this organophosphorus compound has been broadly studied both in vitro and in vivo, there are few studies of metabolic interactions of CPO with other xenobiotics. CPO significantly activated the production of 1-naphthol (5-fold), 2-naphthol (10-fold), trans-1,2-dihydro-1,2-naphthalenediol (1.5-fold), and 1,4-naphthoquinone from naphthalene by human liver microsomes (HLM). It was further demonstrated that the production of naphthalene metabolites by CYP2C8, 2C9*(1), 2C19, 2D6*(1), 3A4, 3A5, and 3A7 was activated by CPO, while the production of naphthalene metabolites by CYP1A1, 1A2, 1B1, and 2B6 was inhibited by CPO. CPO inhibited CYP1A2 production of naphthalene metabolites, while activating their production by CYP3A4. Similarly, CPO inhibited the production of N,N-diethyl-m-hydroxymethylbenzamide (BALC) from DEET by human liver microsomes, but activated the production of N-ethyl-m-toluamide (ET) from this substrate. CYP2B6, the most efficient isoform for BALC production, was inhibited by CPO, while CYP3A4, the most efficient isoform for ET production, was activated by CPO. CPO inhibited CYP2B6 production of both BALC and ET from DEET, but activated CYP3A4 production of ET, while inhibiting CYP3A4 BALC production. CPO appears to facilitate the binding of naphthalene to CYP3A4. This metabolic activation is independent of cytochrome b5, suggesting that activation of CYP3A4 by CPO is associated with a conformational change of the isoform rather than facilitating electron transfer.

Pharmacokinetics in drug discovery

The aim of this current review is to summarize the present status of pharmacokinetics in Drug Discovery. The review is structured into four sections. The first section is a general overview of what we understand by pharmacokinetics and the different LADMET aspects: Liberation, Absorption, Distribution, Metabolism, Excretion, and Toxicity. The second section highlights the different computational or in silico approaches to estimate/predict one or several aspects of the pharmacokinetic profile of a discovery lead compound. The third section discusses the most commonly used in vitro methodologies. The fourth and last section examines the various approaches employed towards the pharmacokinetic assessment of discovery molecules; including all the LADME processes, discussing the different mathematical methodologies available to establish the PK profile of a test compound; what the main differences are and what should be the criteria for using one or another mathematical approach. The major conclusion of this review is that the use of the appropriate preclinical assays has a key role in the long-term viability of a pharmaceutical company since applying the right tools early in discovery will play a key role in determining the company's ability to discover novel safe and effective therapeutics to patients as quickly as possible.

Construction of 3D models of the CYP11B family as a tool to predict ligand binding characteristics

Aldosterone is synthesised by aldosterone synthase (CYP11B2). CYP11B2 has a highly homologous isoform, steroid 11beta-hydroxylase (CYP11B1), which is responsible for the biosynthesis of aldosterone precursors and glucocorticoids. To investigate aldosterone biosynthesis and facilitate the search for selective CYP11B2 inhibitors, we constructed three-dimensional models for CYP11B1 and CYP11B2 for both human and rat. The models were constructed based on the crystal structure of Pseudomonas Putida CYP101 and Oryctolagus Cuniculus CYP2C5. Small steric active site differences between the isoforms were found to be the most important determinants for the regioselective steroid synthesis. A possible explanation for these steric differences for the selective synthesis of aldosterone by CYP11B2 is presented. The activities of the known CYP11B inhibitors metyrapone, R-etomidate, R-fadrazole and S-fadrazole were determined using assays of V79MZ cells that express human CYP11B1 and CYP11B2, respectively. By investigating the inhibitors in the human CYP11B models using molecular docking and molecular dynamics simulations we were able to predict a similar trend in potency for the inhibitors as found in the in vitro assays. Importantly, based on the docking and dynamics simulations it is possible to understand the enantioselectivity of the human enzymes for the inhibitor fadrazole, the R-enantiomer being selective for CYP11B2 and the S-enantiomer being selective for CYP11B1.

Engineering of Cytochrome P450 3A4 for Enhanced Peroxide-Mediated Substrate Oxidation Using Directed Evolution and Site-Directed Mutagenesis

CYP3A4 has been subjected to random and site-directed mutagenesis to enhance peroxide-supported metabolism of several substrates. Initially, a high-throughput screening method using whole cell suspensions was developed for H2O2-supported oxidation of 7-benzyloxyquinoline. Random mutagenesis by error-prone polymerase chain reaction and activity screening yielded several CYP3A4 mutants with enhanced activity. L216W and F228I showed a 3-fold decrease in Km, HOOH and a 2.5-fold increase in kcat/Km, HOOH compared with CYP3A4. Subsequently, T309V and T309A were created based on the observation that T309V in CYP2D6 has enhanced cumene hydroperoxide (CuOOH)-supported activity. T309V and T309A showed a > 6- and 5-fold higher kcat/Km, CuOOH than CYP3A4, respectively. Interestingly, L216W and F228I also exhibited, respectively, a > 4- and a > 3-fold higher kcat/Km, CuOOH than CYP3A4. Therefore, several multiple mutants were constructed from rationally designed and randomly isolated mutants; among them, F228I/T309A showed an 11-fold higher kcat/Km, CuOOH than CYP3A4. Addition of cytochrome b5, which is known to stimulate peroxide-supported activity, enhanced the kcat/Km, CuOOH of CYP3A4 by 4- to 7-fold. When the mutants were tested with other substrates, T309V and T433S showed enhanced kcat/Km, CuOOH with 7-benzyloxy-4-(trifluoromethyl)coumarin and testosterone, respectively, compared with CYP3A4. In addition, in the presence of cytochrome b5, T433S has the potential to produce milligram quantities of 6beta-hydroxytestosterone through peroxide-supported oxidation. In conclusion, a combination of random and site-directed mutagenesis approaches yielded CYP3A4 enzymes with enhanced peroxide-supported metabolism of several substrates.

Designing better drugs: predicting cytochrome P450 metabolism

Many 3D ligand-based and structure-based computational approaches have been used to predict, and thus help explain, the metabolism catalyzed by the enzymes of the cytochrome P450 superfamily (P450s). P450s are responsible for >90% of the metabolism of all drugs, so the computational prediction of metabolism can help to design out drug-drug interactions in the early phases of the drug discovery process. Computational methodologies have focused on a few P450s that are directly involved in drug metabolism. The recently derived crystal structures for human P450s enable better 3D modelling of these important metabolizing enzymes. Models derived for P450s have evolved from simple comparisons of known substrates to more-elaborate experiments that require considerable computer power involving 3D overlaps and docking experiments. These models help to explain and, more importantly, predict the involvement of P450s in the metabolism of specific compounds and guide the drug-design process.

Structure of microsomal cytochrome P450 2B4 complexed with the antifungal drug bifonazole: insight into P450 conformational plasticity and membrane interaction

To better understand ligand-induced structural transitions in cytochrome P450 2B4, protein-ligand interactions were investigated using a bulky inhibitor. Bifonazole, a broad spectrum antifungal agent, inhibits monooxygenase activity and induces a type II binding spectrum in 2B4dH(H226Y), a modified enzyme previously crystallized in the presence of 4-(4-chlorophenyl)imidazole (CPI). Isothermal titration calorimetry and tryptophan fluorescence quenching indicate no significant burial of protein apolar surface nor altered accessibility of Trp-121 upon bifonazole binding, in contrast to recent results with CPI. A 2.3 A crystal structure of 2B4-bifonazole reveals a novel open conformation with ligand bound in the active site, which is significantly different from either the U-shaped cleft of ligand-free 2B4 or the small active site pocket of 2B4-CPI. The O-shaped active site cleft of 2B4-bifonazole is widely open in the middle but narrow at the top. A bifonazole molecule occupies the bottom of the active site cleft, where helix I is bent approximately 15 degrees to accommodate the bulky ligand. The structure also defines unanticipated interactions between helix C residues and bifonazole, suggesting an important role of helix C in azole recognition by mammalian P450s. Comparison of the ligand-free 2B4 structure, the 2B4-CPI structure, and the 2B4-bifonazole structure identifies structurally plastic regions that undergo correlated conformational changes in response to ligand binding. The most plastic regions are putative membrane-binding motifs involved in substrate access or substrate binding. The results allow us to model the membrane-associated state of P450 and provide insight into how lipophilic substrates access the buried active site.

Structural diversities of active site in clinical azole-bound forms between sterol 14alpha-demethylases (CYP51s) from human and Mycobacterium tuberculosis

To gain insights into the molecular basis of the design for the selective azole anti-fungals, we compared the binding properties of azole-based inhibitors for cytochrome P450 sterol 14alpha-demethylase (CYP51) from human (HuCYP51) and Mycobacterium tuberculosis (MtCYP51). Spectroscopic titration of azoles to the CYP51s revealed that HuCYP51 has higher affinity for ketoconazole (KET), an azole derivative that has long lipophilic groups, than MtCYP51, but the affinity for fluconazole (FLU), which is a member of the anti-fungal armamentarium, was lower in HuCYP51. The affinity for 4-phenylimidazole (4-PhIm) to MtCYP51 was quite low compared with that to HuCYP51. In the resonance Raman spectra for HuCYP51, the FLU binding induced only minor spectral changes, whereas the prominent high frequency shift of the bending mode of the heme vinyl group was detected in the KET- or 4-PhIm-bound forms. On the other hand, the bending mode of the heme propionate group for the FLU-bound form of MtCYP51 was shifted to high frequency as found for the KET-bound form, but that for 4-PhIm was shifted to low frequency. The EPR spectra for 4-PhIm-bound MtCYP51 and FLU-bound HuCYP51 gave multiple g values, showing heterogeneous binding of the azoles, whereas the single gx and gz values were observed for other azole-bound forms. Together with the alignment of the amino acid sequence, these spectroscopic differences suggest that the region between the B' and C helices, particularly the hydrophobicity of the C helix, in CYP51s plays primary roles in determining strength of interactions with azoles; this differentiates the binding specificity of azoles to CYP51s.

Characterization of the CYP2C8 Active Site by Homology Modeling

To compare the features of the active sites of CYP2C8, CYP2C9, and CYP2C19, homology modeling was performed based on the crystallographic coordinates of mammalian CYP2C5. It was found that CYP2C8 has a much larger pocket than the other forms due to the existence of an additional pocket. The approach to the additional pocket is comprised of Ile102, Ser114, Leu208, Val366, and Ile476, and the side chains of Ser114, Val366, and Ile476, which are smaller than the corresponding residues in the other CYPs, enable access to the pocket. The general features of the active site in the CYP2C8 model are similar to those of the previously constructed CYP3A4 model, which may account for the 2 CYPs sharing some of their substrates. The CYP2C8 model was validated by examining the bound orientation of paclitaxel and showing that it is consistent with the formation of the 6-beta hydroxylated derivative during metabolism. Docked paclitaxel was found to form a hydrogen bond with the side chain of Asn 99, which is a characteristic residue of CYP2C8 and is located in the additional pocket. Descriptors for CYP2C8 and CYP2C9 substrates were also examined with the molecular operating environment (MOE). The descriptor by which CYP2C8 and CYP2C9 substrates were classified most distinctly was found to be molar refractivity, which might be related to the longer shape and more polar nature of the active site of CYP2C8 in the CYP2C subfamily.