Working at the membrane interface: Ligand-induced changes in dynamic conformation and oligomeric structure in human aromatase

CYP116B5-SOX: An artificial peroxygenase for drug metabolites production and bioremediation

CYP116B5 is a class VII P450 in which the heme domain is linked to a FMN and 2Fe2S-binding reductase. Our laboratory has proved that the CYP116B5 heme domain (CYP116B5-hd) is capable of catalyzing the oxidation of substrates using H2O2. Recently, the Molecular Lego approach was applied to join the heme domain of CYP116B5 to sarcosine oxidase (SOX), which provides H2O2 in-situ by the sarcosine oxidation. In this work, the chimeric self-sufficient fusion enzyme CYP116B5-SOX was heterologously expressed, purified, and characterized for its functionality by absorbance and fluorescence spectroscopy. Differential scanning calorimetry (DSC) experiments revealed a TM of 48.4 ± 0.04 and 58.3 ± 0.02°C and a enthalpy value of 175,500 ± 1850 and 120,500 ± 1350 cal mol-1 for the CYP116B5 and SOX domains respectively. The fusion enzyme showed an outstanding chemical stability in presence of up to 200 mM sarcosine or 5 mM H2O2 (4.4 ± 0.8 and 11.0 ± 2.6% heme leakage respectively). Thanks to the in-situ H2O2 generation, an improved kcat/KM for the p-nitrophenol conversion was observed (kcat of 20.1 ± 0.6 min-1 and KM of 0.23 ± 0.03 mM), corresponding to 4 times the kcat/KM of the CYP116B5-hd. The aim of this work is the development of an engineered biocatalyst to be exploited in bioremediation. In order to tackle this challenge, an E. coli strain expressing CYP116B5-SOX was employed to exploit this biocatalyst for the oxidation of the wastewater contaminating-drug tamoxifen. Data show a 12-fold increase in tamoxifen N-oxide production-herein detected for the first time as CYP116B5 metabolite-compared to the direct H2O2 supply, equal to the 25% of the total drug conversion.

Albumin/Mitotane Interaction Affects Drug Activity in Adrenocortical Carcinoma Cells: Smoke and Mirrors on Mitotane Effect with Possible Implications for Patients' Management

BACKGROUND

Mitotane is the only drug approved for the treatment of adrenocortical carcinoma (ACC). Although it has been used for many years, its mechanism of action remains elusive. H295R cells are, in ACC, an essential tool to evaluate drug mechanisms, although they often lead to conflicting results.

METHODS

Using different in vitro biomolecular technologies and biochemical/biophysical experiments, we evaluated how the presence of "confounding factors" in culture media and patient sera could reduce the pharmacological effect of mitotane and its metabolites.

RESULTS

We discovered that albumin, the most abundant protein in the blood, was able to bind mitotane. This interaction altered the effect of the drug by blocking its biological activity. This blocking effect was independent of the albumin source or methodology used and altered the assessment of drug sensitivity of the cell lines.

CONCLUSIONS

In conclusion, we have for the first time demonstrated that albumin does not only act as an inert drug carrier when mitotane or its metabolites are present. Indeed, our experiments clearly indicated that both albumin and human serum were able to suppress the pharmacological effect of mitotane in vitro. These experiments could represent a first step towards the individualization of mitotane treatment in this rare tumor.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Molecular Lego of Human Cytochrome P450: The Key Role of Heme Domain Flexibility for the Activity of the Chimeric Proteins

The cytochrome P450 superfamily are heme-thiolate enzymes able to carry out monooxygenase reactions. Several studies have demonstrated the feasibility of using a soluble bacterial reductase from Bacillus megaterium, BMR, as an artificial electron transfer partner fused to the human P450 domain in a single polypeptide chain in an approach known as ‘molecular Lego’. The 3A4-BMR chimera has been deeply characterized biochemically for its activity, coupling efficiency, and flexibility by many different biophysical techniques leading to the conclusion that an extension of five glycines in the loop that connects the two domains improves all the catalytic parameters due to improved flexibility of the system. In this work, we extend the characterization of 3A4-BMR chimeras using differential scanning calorimetry to evaluate stabilizing role of BMR. We apply the ‘molecular Lego’ approach also to CYP19A1 (aromatase) and the data show that the activity of the chimeras is very low (<0.003 min−1) for all the constructs tested with a different linker loop length: ARO-BMR, ARO-BMR-3GLY, and ARO-BMR-5GLY. Nevertheless, the fusion to BMR shows a remarkable effect on thermal stability studied by differential scanning calorimetry as indicated by the increase in Tonset by 10 °C and the presence of a cooperative unfolding process driven by the BMR protein domain. Previously characterized 3A4-BMR constructs show the same behavior of ARO-BMR constructs in terms of thermal stabilization but a higher activity as a function of the loop length. A comparison of the ARO-BMR system to 3A4-BMR indicates that the design of each P450-BMR chimera should be carefully evaluated not only in terms of electron transfer, but also for the biophysical constraints that cannot always be overcome by chimerization.

Polymorphism on human aromatase affects protein dynamics and substrate binding: spectroscopic evidence

Human aromatase is a member of the cytochrome P450 superfamily, involved in steroid hormones biosynthesis. In particular, it converts androgen into estrogens being therefore responsible for the correct sex steroids balance. Due to its capacity in producing estrogens it has also been considered as a promising target for breast cancer therapy. Two single-nucleotide polymorphisms (R264C and R264H) have been shown to alter aromatase activity and they have been associated to an increased or decreased risk for estrogen-dependent pathologies. Here, the effect of these mutations on the protein dynamics is investigated by UV/FTIR and time resolved fluorescence spectroscopy. H/D exchange rates were measured by FTIR for the three proteins in the ligand-free, substrate- and inhibitor-bound forms and the data indicate that the wild-type enzyme undergoes a conformational change leading to a more compact tertiary structure upon substrate or inhibitor binding. Indeed, the H/D exchange rates are decreased when a ligand is present. In the variants, the exchange rates in the ligand-free and -bound forms are similar, indicating that a structural change is lacking, despite the single amino acid substitution is located in the peripheral shell of the protein molecule. Moreover, the fluorescence lifetimes data show that the quenching effect on tryptophan-224 observed upon ligand binding in the wild-type, is absent in both variants. Since this residue is located in the catalytic pocket, these findings suggest that substrate entrance and/or retention in the active site is partially compromised in both mutants. A contact network analysis demonstrates that the protein structure is organized in two main clusters, whose connectivity is altered by ligand binding, especially in correspondence of helix-G, where the amino acid substitutions occur. Our findings demonstrate that SNPs resulting in mutations on aromatase surface modify the protein flexibility that is required for substrate binding and catalysis. The cluster analysis provides a rationale for such effect, suggesting helix G as a possible target for aromatase inhibition.

Molecular determinant of substrate binding and specificity of cytochrome P450 2J2

Cytochrome P450 2J2 (CYP2J2) is responsible for the epoxidation of endogenous arachidonic acid, and is involved in the metabolism of exogenous drugs. To date, no crystal structure of CYP2J2 is available, and the proposed structural basis for the substrate recognition and specificity in CYP2J2 varies with the structural models developed using different computational protocols. In this study, we developed a new structural model of CYP2J2, and explored its sensitivity to substrate binding by molecular dynamics simulations of the interactions with chemically similar fluorescent probes. Our results showed that the induced-fit binding of these probes led to the preferred active poses ready for the catalysis by CYP2J2. Divergent conformational dynamics of CYP2J2 due to the binding of each probe were observed. However, a stable hydrophobic clamp composed of residues I127, F310, A311, V380, and I487 was identified to restrict any substrate access to the active site of CYP2J2. Molecular docking of a series of compounds including amiodarone, astemizole, danazol, ebastine, ketoconazole, terfenadine, terfenadone, and arachidonic acid to CYP2J2 confirmed the role of those residues in determining substrate binding and specificity of CYP2J2. In addition to the flexibility of CYP2J2, the present work also identified other factors such as electrostatic potential in the vicinity of the active site, and substrate strain energy and property that have implications for the interpretation of CYP2J2 metabolism.

Ligand stabilization and effect on unfolding by polymorphism in human flavin-containing monooxygenase 3

Pharmacogenomics is a powerful tool to prevent adverse reactions caused by different response of individuals to drug administration. Single nucleotide polymorphisms (SNPs) represent up to 90% of genetic variations among individuals. Drug metabolizing enzymes are highly polymorphic therefore the kinetic parameters of their catalytic reactions can be significantly influenced. This work reports on the unfolding process of a phase I drug metabolizing enzyme, human flavin-containing monooxygenase 3 (hFMO3) and its single nucleotide polymorphic variants (SNPs) V257M, E158K and E308G. Differential scanning calorimetry (DSC) indicates that the thermal denaturation of the enzyme is irreversible. The melting temperature (T_m) for the (Wild Type) WT and its polymorphic variants is found to be in a range from 46 °C to 50 °C. Also the activation energies of unfolding (E_a) show no significant differences among all proteins investigated (290-328 KJ/mol), except for the E308G variant that showed a significantly higher E_a of 412 KJ/mol. The presence of the bound NADP⁺ cofactor is found to stabilize all the variants by shifting the main T_m by 4-5 °C for all the proteins, exception made for E308G where no changes are observed. Isothermal titration calorimetry (ITC) was used to characterize the interaction of the protein with NADP⁺ in terms of dissociation constant (K_d), enthalpy (ΔH) and entropy (ΔS). K_d values of 1.6 and 0.7 μM, ΔH of -13.9 Kcal/mol and -16.8 Kcal/mol, ΔS of -20.5 cal/mol/deg, and -28.5 cal/mol/deg were found for V257M and E158K respectively. E308G was found to be unable to bind the NADP⁺ cofactor, a result that is in line with the T_m results. Circular dichroism also confirmed an overall lower stability of E308G, while NADP⁺ was found to give a strong positive shift of the T_m stabilizing the structure of E158K (46.2 to 50.6 °C). Previous data highlighted significant differences in terms of activity among the SNPs of hFMO3. In this work a minor impact of the SNPs was found on the stability of the enzyme in the ligand free form, except for E308G, whereas the binding of NADP⁺ reveals major differences among WT and polymorphic variants that are all measurable in terms of heat capacity, enthalpy and secondary structure content. These data provide the first direct evidence of ligand stabilization effects on hFMO3 that can explain the differences observed in catalytic efficiencies and serve as the starting point for the development of inhibitors of this enzyme.

Effector role of cytochrome P450 reductase for androstenedione binding to human aromatase

Cytochromes P450 constitute a large superfamily of monooxygenases involved in many metabolic pathways. Most of them are not self-sufficient and need a reductase protein to provide the electrons necessary for catalysis. It was shown that the redox partner plays a role in the modulation of the structure and function of some bacterial P450 enzymes. Here, the effect of NADPH-cytochrome reductase (CPR) on human aromatase (Aro) is studied for what concerns its role in substrate binding. Pre-steady-state kinetic experiments indicate that both the substrate binding rates and the percentage of spin shift detected for aromatase are increased when CPR is present. Moreover, aromatase binds the substrate through a conformational selection mechanism, suggesting a possible effector role of CPR. The thermodynamic parameters for the formation of the CPR-Aro complex were studied by isothermal titration calorimetry. The dissociation constant of the complex formation is 4.5 folds lower for substrate-free compared to the substrate-bound enzyme. The enthalpy change observed when the CPR-Aro complex forms in the absence of the substrate are higher than in its presence, indicating that more interactions are formed/broken in the former case. Taken together, our data confirm that CPR has a role in promoting aromatase conformation optimal for substrate binding.

Estrogen receptor-positive (ER+) breast cancer treatment: Are multi-target compounds the next promising approach?

Endocrine therapy is currently the main therapeutic approach for estrogen receptor-positive (ER⁺) breast cancer, the most frequent subtype of breast cancer in women worldwide. For this subtype of tumors, the current clinical treatment includes aromatase inhibitors (AIs) and anti-estrogenic compounds, such as Tamoxifen and Fulvestrant, being AIs the first-line treatment option for post-menopausal women. Moreover, the recent guidelines also suggest the use of these compounds by pre-menopausal women after suppressing ovaries function. However, besides its therapeutic efficacy, the prolonged use of this type of therapies may lead to the development of several adverse effects, as well as, endocrine resistance, limiting the effectiveness of such treatments. In order to surpass this issues and clinical concerns, during the last years, several studies have been suggesting alternative therapeutic approaches, considering the function of aromatase, ERα and ERβ. Here, we review the structural and functional features of these three targets and their importance in ER⁺ breast cancer treatment, as well as, the current treatment strategies used in clinic, emphasizing the importance of the development of multi-target compounds able to simultaneously modulate these key targets, as a novel and promising therapeutic strategy for this type of cancer.

Differential effects of variations in human P450 oxidoreductase on the aromatase activity of CYP19A1 polymorphisms R264C and R264H

Aromatase (CYP19A1) converts androgens into estrogens and is required for female sexual development and growth and development in both sexes. CYP19A1 is a member of cytochrome P450 family of heme-thiolate monooxygenases located in the endoplasmic reticulum and depends on reducing equivalents from the reduced nicotinamide adenine dinucleotide phosphate via the cytochrome P450 oxidoreductase coded by POR. Both the CYP19A1 and POR genes are highly polymorphic, and mutations in both these genes are linked to disorders of steroid biosynthesis. We have previously shown that R264C and R264H mutations in CYP19A1, as well as mutations in POR, reduce CYP19A1 activity. The R264C is a common polymorphic variant of CYP19A1, with high frequency in Asian and African populations. Polymorphic alleles of POR are found in all populations studied so far and, therefore, may influence activities of CYP19A1 allelic variants. So far, the effects of variations in POR on enzymatic activities of allelic variants of CYP19A1 or any other steroid metabolizing cytochrome P450 proteins have not been studied. Here we are reporting the effects of three POR variants on the aromatase activities of two CYP19A1 variants, R264C, and R264H. We used bacterially expressed and purified preparations of WT and variant forms of CYP19A1 and POR and constructed liposomes with embedded CYP19A1 and POR proteins and assayed the CYP19A1 activities using radiolabeled androstenedione as a substrate. With the WT-POR as a redox partner, the R264C-CYP19A1 showed only 15% of aromatase activity, but the R264H had 87% of aromatase activity compared to WT-CYP19A1. With P284L-POR as a redox partner, R264C-CYP19A1 lost all activity but retained 6.7% of activity when P284T-POR was used as a redox partner. The R264H-CYP19A1 showed low activities with both the POR-P284 L as well as the POR-P284 T. When the POR-Y607C was used as a redox partner, the R264C-CYP19A1 retained approximately 5% of CYP19A1 activity. Remarkably, The R264H-CYP19A1 had more than three-fold higher activity compared to WT-CYP19A1 when the POR-Y607C was used as the redox partner, pointing toward a beneficial effect. The slight increase in activity of R264C-CYP19A1 with the P284T-POR and the three-fold increase in activity of the R264H-CYP19A1 with the Y607C-POR point toward a conformational effect and role of protein-protein interaction governed by the R264C and R264H substitutions in the CYP19A1 as well as P284 L, P284 T and Y607C variants of POR. These studies demonstrate that the allelic variants of P450 when present with a variant form of POR may show different activities, and combined effects of variations in the P450 enzymes as well as in the POR should be considered when genetic data are available. Recent trends in the whole-exome and whole-genome sequencing as diagnostic tools will permit combined evaluation of variations in multiple genes that are interdependent and may guide treatment options by adjusting therapeutic interventions based on laboratory analysis.

Impact of peripheral mutations on the access channels of human cytochrome P450 1A2

Communicated by Ramaswamy H. Sarma.

Crystal structure of bacterial CYP116B5 heme domain: New insights on class VII P450s structural flexibility and peroxygenase activity

Class VII cytochromes P450 are self-sufficient enzymes carrying a phthalate family oxygenase-like reductase domain and a P450 domain fused in a single polypeptide chain. The biocatalytic applications of CYP116B members are limited by the need of the NADPH cofactor and the lack of crystal structures as a starting point for protein engineering. Nevertheless, we demonstrated that the heme domain of CYP116B5 can use hydrogen peroxide as electron donor bypassing the need of NADPH. Here, we report the crystal structure of CYP116B5 heme domain in complex with histidine at 2.6 Å of resolution. The structure reveals the typical P450 fold and a closed conformation with an active site cavity of 284 ų in volume, accommodating a histidine molecule forming a hydrogen bond with the water molecule present as 6th heme iron ligand. MD simulations in the absence of any ligand revealed the opening of a tunnel connecting the active site to the protein surface through the movement of F-, G- and H-helices. A structural alignment with bacterial cytochromes P450 allowed the identification of amino acids in the proximal heme site potentially involved in peroxygenase activity. The availability of the crystal structure provides the bases for the structure-guided design of new biocatalysts.

Human cytochrome P450 enzymes 5-51 as targets of drugs and natural and environmental compounds: mechanisms, induction, and inhibition - toxic effects and benefits

Cytochrome P450 (P450, CYP) enzymes have long been of interest due to their roles in the metabolism of drugs, pesticides, pro-carcinogens, and other xenobiotic chemicals. They have also been of interest due to their very critical roles in the biosynthesis and metabolism of steroids, vitamins, and certain eicosanoids. This review covers the 22 (of the total of 57) human P450s in Families 5-51 and their substrate selectivity. Furthermore, included is information and references regarding inducibility, inhibition, and (in some cases) stimulation by chemicals. We update and discuss important aspects of each of these 22 P450s and questions that remain open.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.