Small angle scattering reveals the orientation of cytochrome P450 19A1 in lipoprotein nanodiscs.

Human aromatase (CYP19A1), the cytochrome P450 enzyme responsible for conversion of androgens to estrogens, was incorporated into lipoprotein nanodiscs (NDs) and interrogated by small angle X-ray and neutron scattering (SAXS/SANS). CYP19A1 was associated with the surface and centered at the edge of the long axis of the ND membrane. In the absence of the N-terminal anchor, the amphipathic A'- and G'-helices were predominately buried in the lipid head groups, with the possibly that their hydrophobic side chains protrude into the hydrophobic, aliphatic tails. The prediction is like that for CYP3A4 based on SAXS employing a similar modeling approach. The orientation of CYP19A1 in a ND is consistent with our previous predictions based on molecular dynamics simulations and lends additional credibility to the notion that CYP19A1 captures substrates from the membrane.

HCV Antiviral Drugs Have the Potential to Adversely Perturb the Fetal-Maternal Communication Axis through Inhibition of CYP3A7 DHEA-S Oxidation.

The hepatitis C virus (HCV) poses a great risk to pregnant people and their developing fetus, yet no HCV antiviral treatment guidelines have been established. While there has been a substantial increase in the development of HCV antivirals, the effect they have on the developing fetus remains poorly defined. Many of these drugs are metabolized through the cytochrome P450 CYP3A pathway, which is mediated by cytochrome P450 3A7 (CYP3A7) in the fetus and developing infant. In this study, we sought to investigate the effect HCV antivirals have on CYP3A7 metabolism, as this CYP enzyme plays a vital role in proper fetal and neonatal development. Of the 13 HCV antivirals we investigated, 8 (∼62%) inhibited CYP3A7 metabolic activity by 50% or more at a concentration of 20 µM. Furthermore, paritaprevir, asunaprevir, simeprevir, danoprevir, and glecaprevir all had observed half-maximal inhibitory concentrations between the range of 10 and 20 µM, which is physiologically relevant in comparison with the Km of dehydroepiandrosterone-sulfate (DHEA-S) oxidation (reported to be between 5 and 20 µM). We also discovered that paritaprevir is a time-dependent inhibitor of CYP3A7, which shifts the IC50 ∼twofold from 11 µM to 5 µM. Upon further characterization, paritaprevir inactivates DHEA-S metabolism by CYP3A7, with KI and Kinact values of 4.66 µM and 0.00954 minute-1, respectively. Depending on treatment plan and off-label drug use, HCV treatment could adversely affect the fetal-maternal communication axis by blocking fetal CYP3A7 metabolism of important endogenous hormones. SIGNIFICANCE STATEMENT: The prevalence of HCV in pregnant people is estimated at between 1% and 8% of the global population, yet little to no information exists about the risk antiviral treatment poses to the developing fetus. There is a potential risk of drugs adversely affecting mother-fetal communication by inhibiting fetal hepatic CYP3A7, an integral enzyme for estriol production. We discovered that five HCV antivirals inhibited DHEA-S metabolism by CYP3A7, and paritaprevir inactivated the enzyme. Our studies demonstrate the potential threat these drugs pose to proper fetal development.

Induced Fit Describes Ligand Binding to Membrane-Associated Cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant P450 involved in human xenobiotic metabolism. Competition for CYP3A4 therefore underlies several adverse drug-drug interactions. Despite its clinical significance, the mechanisms CYP3A4 uses to bind diverse ligands remain poorly understood. Highly monodisperse CYP3A4 embedded in anionic lipoprotein nanodiscs containing an equal mixture of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) were used to determine which of the limiting kinetic schemes that include protein conformational change, conformational selection (CS) or induced fit (IF), best described the binding of four known irreversible inhibitors. Azamulin, retapamulin, pleuromutilin, and mibrefadil binding to CYP3A4 nanodiscs conformed to a single-site binding model. Exponential fits of stopped-flow UV-visible absorption spectroscopy data supported multiple-step binding mechanisms. Trends in the rates of relaxation to equilibrium with increasing ligand concentrations were ambiguous as to whether IF or CS was involved; however, global fitting and consideration of the rate constants favored an IF mechanism. In the case of mibrefadil, a transient complex was observed in the stopped-flow UV-visible experiment, definitively assigning the presence of IF in ligand binding. While these studies only consider a small region of CYP3A4's vast ligand space, they provide kinetic evidence that CYP3A4 can use an IF mechanism. SIGNIFICANCE STATEMENT: CYP3A4 is capable of oxidizing numerous xenobiotics, including many drugs. Such promiscuity could not be achieved without conformational changes to accommodate diverse substrates. It is unknown whether conformational heterogeneity is present before (conformational selection) or after (induced fit) ligand binding. Stopped-flow measurements of suicide inhibitors binding to nanodisc-embedded CYP3A4 combined with sophisticated numerical analyses support that induced fit better describes ligand binding to this important enzyme.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements.

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-ß1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Structures and Dynamics of Anionic Lipoprotein Nanodiscs.

Nanolipoprotein particles known as nanodiscs (NDs) have emerged as versatile and powerful tools for the stabilization of membrane proteins permitting a plethora of structural and biophysical studies. Part of their allure is their flexibility to accommodate many types of lipids and precise control of the composition. However, little is known about how variations in lipid composition impact their structures and dynamics. Herein, we investigate how the introduction of the anionic lipid POPG into POPC NDs impacts these features. Small-angle X-ray and neutron scattering (SAXS and SANS) of variable-composition NDs are complemented with molecular dynamics simulations to interrogate how increasing the concern of POPG impacts the ND shape, structure of the lipid core, and the dynamics of the popular membrane scaffold protein, MSP1D1(-). A convenient benefit of including POPG is that it eliminates D2O-induced aggregation observed in pure POPC NDs, permitting studies by SANS at multiple contrasts. SAXS and SANS data could be globally fit to a stacked elliptical cylinder model as well as an extension of the model that accounts for membrane curvature. Fitting to both models supports that the introduction of POPG results in strongly elliptical NDs; however, MD simulations predict the curvature of the membrane, thereby supporting the use of the latter model. Trends in the model-independent parameters suggest that increases in POPG reduce the conformational heterogeneity of the MSP1D1(-), which is in agreement with MD simulations that show that the incorporation of sufficient POPG suppresses disengagement of the N-terminal helix from the lipid core. These studies highlight novel structural changes in NDs in response to an anionic lipid and will inform the interpretation of future structural studies of membrane proteins embedded in NDs of mixed lipid composition.

PSI relieves the pressure of membrane fusion.

Some plant proteases contain a latent sequence known as the plant-specific insert (PSI) that, upon release from the full protease sequence, initiates membrane fusion to defend from pathogens. However, the mechanism by which it exerts its effects has been unclear. Zhao et al. report an elegant integration of biophysical experiments and molecular dynamics simulations to reveal events leading up to PSI-mediated membrane fusion. Their results demonstrate a pH-dependent monomer-to-dimer transition, clear evidence of membrane association, and probable structures of prefusion intermediates. These data expand our understanding of the elusive PSIs and may provide new directions for antimicrobial development.

Dynamics and Mechanism of Binding of Androstenedione to Membrane-Associated Aromatase.

Aromatase (CYP19A1) catalyzes the synthesis of estrogens from androgens and is an invaluable target of pharmacotherapy for estrogen-dependent cancers. CYP19A1 is also one of the most primordial human CYPs and, to the extent that its fundamental dynamics are conserved, is highly relevant to understanding those of the more recently evolved and promiscuous enzymes. A complementary approach employing molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry (HDX-MS) was employed to interrogate the changes in CYP19A1 dynamics coupled to binding androstenedione (ASD). Gaussian-accelerated molecular dynamics and HDX-MS agree that ASD globally suppresses CYP19A1 dynamics. Bimodal HDX patterns of the B'-C loop potentially arising from at least two conformations are present in free 19A1 only, supporting the possibility that conformational selection is operative. Random-acceleration molecular dynamics and adaptive biasing force simulations illuminate ASD's binding pathway, predicting ASD capture in the lipid headgroups and a pathway to the active site shielded from solvent. Intriguingly, the predicted access channel in 19A1 aligns well with the steroid binding sites of other human sterol-oxidizing CYPs.

Conformational selection is present in ligand binding to cytochrome P450 19A1 lipoprotein nanodiscs.

Cytochromes P450 (CYPs) display remarkable plasticity in their ability to bind substrates and catalyze a broad array of chemical reactions. Herein we evaluate binding of androstenedione, testosterone, and 7-hydroxyflavone to CYP19A1, also known as aromatase, in phospholipid nanodiscs by stopped-flow UV-vis spectroscopy. Exponential fitting of the kinetic traces supports the possibility of a multi-step binding mechanism. Subsequent global fitting of the data to the solutions of the coupled differential equations describing the fundamental mechanisms of induced fit and conformational selection, consistently support presence of the latter. To our knowledge, this is the first discrimination of conformational selection from induced fit for a mono-disperse CYP in a native-like membrane environment. In addition, 7-hydroxyflavone binds to CYP19A1 nanodiscs with comparable affinity to the substrates and induces an unusual spectral response likely attributable to hydrogen bonding to, rather than displacement of the heme-coordinated water molecule.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs.

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the ß1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

Noncovalent interactions dominate dynamic heme distortion in cytochrome P450 4B1.

Cytochrome P450 4B1 (4B1) functions in both xenobiotic and endobiotic metabolism. An ester linkage between Glu-310 in 4B1 and the 5-methyl group of heme facilitates preferential hydroxylation of terminal (ω) methyl groups of hydrocarbons (HCs) and fatty acids compared with ω-1 sites bearing weaker C-H bonds. This preference is retained albeit diminished 4-fold for the E310A mutant, but the reason for this is unclear. Here, a crystal structure of the E310A-octane complex disclosed that noncovalent interactions maintain heme deformation in the absence of the ester linkage. Consistent with the lower symmetry of the heme, resonance Raman (RR) spectroscopy revealed large enhancements of RR peaks for high-spin HC complexes of 4B1 and the E310A mutant relative to P450 3A4. Whereas these enhancements were diminished in RR spectra of a low-spin 4B1-N-hydroxy-N'-(4-butyl-2-methylphenyl)formamidine complex, a crystal structure indicated that this inhibitor does not alter heme ruffling. RR spectra of Fe2+-CO HC complexes revealed larger effects of HC length in E310A than in 4B1, suggesting that reduced rigidity probably underlies increased E310A-catalyzed (ω-1)-hydroxylation. Diminished effects of the HC on the position of the Fe-CO stretching mode in 4B1 suggested that the ester linkage limits substrate access to the CO. Heme ruffling probably facilitates autocatalytic ester formation by reducing inhibitory coordination of Glu-310 with the heme iron. This also positions the 5-methyl for a reaction with the proposed glutamyl radical intermediate and potentially enhances oxo-ferryl intermediate reactivity for generation of the glutamyl radical to initiate ester bond formation and ω-hydroxylation.

Biophysical characterization of Aptenodytes forsteri cytochrome P450 aromatase.

Cytochrome P450 19 (CYP19, aromatase) catalyzes the conversion of androgens to estrogens in a sequence of three reactions that each depend on NADPH and O2. Aromatase is a phylogenetically-ancient enzyme and its breadth of expression in other species has highlighted distinct physiological functions. In songbirds, estrogen production is required for programming the neural circuits controlling song and in the determination of sex in fish and reptiles. This work describes the expression, purification, and biophysical characterization of Aptenodytes forsteri (Emperor penguin, af) aromatase. Using human cytochrome P450 reductase as a redox partner, afCYP19 displayed similar substrate turnover and LC/MS/MS confirmed that afCYP19 catalyzes the transformations through the intermediates 19-hydroxy- and 19-oxo-androstenedione. Androstenedione and anastrozole had the highest affinity for the enzyme and were followed closely by 19-hydroxyandrostenedione and testosterone. The affinity of 19-oxo-androstenedione for afCYP19 was ten-fold lower. The time-dependent changes in the Soret bands observed in stopped-flow mixing experiments of the steroidal ligands and the inhibitor anastrozole with afCYP19 were best described by a two-step binding mechanism. In summary, these studies describe the first biophysical characterization of an avian aromatase that displays strikingly similar enzyme kinetics and ligand binding properties to the human enzyme and could serve as a convenient model system for studies of the enigmatic transformation of androgens to estrogens.

Membrane-embedded substrate recognition by cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic-metabolizing enzyme in the liver and intestine and is involved in the disposition of more than 50% of drugs. Because of its ability to bind multiple substrates, its reaction kinetics are complex, and its association with the microsomal membrane confounds our understanding of how this enzyme recognizes and recruits diverse substrates. Testosterone (TST) hydroxylation is the prototypical CYP3A4 reaction, displaying positive homotropic cooperativity with three binding sites. Here, exploiting the capability of accelerated molecular dynamics (aMD) to sample events in the millisecond regime, I performed >25-µs aMD simulations in the presence of three TST molecules. These simulations identified high-occupancy surface-binding sites as well as a pathway for TST ingress into the CYP3A4 active site originating in the membrane. Adaptive biasing force analysis of the latter pathway revealed a metastable intermediate that could constitute a third binding site at high TST concentrations. Prompted by the observation that interactions between TST and the G'-helix mobilize the ligand into the active site, a free-energy analysis of TST distribution in the membrane was conducted and revealed that the depth of the G'-helix is optimal for extracting TST. In summary, these simulations confirm separate, but adjacent substrate-binding sites within the enzyme and the existence of an auxiliary TST-binding site. The broader impact of these simulations is that they support a mechanism in which cytochromes P450 directly recruit membrane-solubilized substrates.

Engineered P450 biocatalysts show improved activity and regio-promiscuity in aromatic nitration.

Nitroaromatics are among the most important and commonly used chemicals but their production often suffers from multiple unsolved challenges. We have previously described the development of biocatalytic nitration processes driven by an engineered P450 TxtE fusion construct. Herein we report the creation of improved nitration biocatalysts through constructing and characterizing fusion proteins of TxtE with the reductase domain of CYP102A1 (P450BM3, BM3R). The majority of constructs contained variable linker length while one was rationally designed for optimizing protein-protein interactions. Detailed biochemical characterization identified multiple active chimeras that showed improved nitration activity, increased coupling efficiency and higher total turnover numbers compared with TxtE. Substrate promiscuity of the most active chimera was further assessed with a substrate library. Finally, a biocatalytic nitration process was developed to nitrate 4-Me-DL-Trp. The production of both 4-Me-5-NO2-L-Trp and 4-Me-7-NO2-L-Trp uncovered remarkable regio-promiscuity of nitration biocatalysts.

N-Heterocyclic Carbene Capture by Cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant P450 enzyme involved in human drug metabolism, and its inhibition may result in adverse interactions or, conversely, favorably reduce the systemic elimination rates of poorly bioavailable drugs. Herein we describe a spectroscopic investigation of the interaction of CYP3A4 with N-methylritonavir, an analog of ritonavir, widely used as a pharmacoenhancer. In contrast to ritonavir, the binding affinity of N-methylritonavir for CYP3A4 is pH-dependent. At pH <7.4, the spectra are definitively type I, whereas at pH ≥7.4 the spectra have split Soret bands, including a red-shifted component characteristic of a P450-carbene complex. Variable-pH UV-visible spectroscopy binding studies with molecular fragments narrows the source of this pH dependence to its N-methylthiazolium fragment. The C2 proton of this group is acidic, and variable-pH resonance Raman spectroscopy tentatively assigns it a pKa of 7.4. Hence, this fragment of N-methylritonavir is expected to be readily deprotonated under physiologic conditions to yield a thiazol-2-ylidene, which is an N-heterocyclic carbene that has high-affinity for and is presumed to be subsequently captured by the heme iron. This mechanism is supported by time-dependent density functional theory with an active site model that accurately reproduces distinguishing features of the experimental UV-visible spectra of N-methylritonavir bound to CYP3A4. Finally, density functional theory calculations support that this novel interaction is as strong as the tightest-binding azaheterocycles found in P450 inhibitors and could offer new avenues for inhibitor development.

Photoswitchable anticancer activity via trans-cis isomerization of a combretastatin A-4 analog.

Combretastatin A-4 (CA4) is highly potent anticancer drug that acts as an inhibitor of tubulin polymerization. The core of the CA4 structure contains a cis-stilbene, and it is known that the trans isomer is significantly less potent. We prepared an azobenzene analog of CA4 (Azo-CA4) that shows 13-35 fold enhancement in potency upon illumination. EC50 values in the light were in the mid nM range. Due to its ability to thermally revert to less toxic trans form, Azo-CA4 also has the ability to automatically turn its activity off with time. Azo-CA4 is less potent than CA-4 because it degrades in the presence of glutathione as evidenced by UV-Vis spectroscopy and ESI-MS. Nevertheless, Azo-CA4 represents a promising strategy for switchable potency for treatment of cancer.

Ab initio dynamics of the cytochrome P450 hydroxylation reaction.

The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.

Evidence for an elevated aspartate pK(a) in the active site of human aromatase.

Aromatase (CYP19A1), the enzyme that converts androgens to estrogens, is of significant mechanistic and therapeutic interest. Crystal structures and computational studies of this enzyme shed light on the critical role of Asp(309) in substrate binding and catalysis. These studies predicted an elevated pK(a) for Asp(309) and proposed that protonation of this residue was required for function. In this study, UV-visible absorption, circular dichroism, resonance Raman spectroscopy, and enzyme kinetics were used to study the impact of pH on aromatase structure and androstenedione binding. Spectroscopic studies demonstrate that androstenedione binding is pH-dependent, whereas, in contrast, the D309N mutant retains its ability to bind to androstenedione across the entire pH range studied. Neither pH nor mutation perturbed the secondary structure or heme environment. The origin of the observed pH dependence was further narrowed to the protonation equilibria of Asp(309) with a parallel set of spectroscopic studies using exemestane and anastrozole. Because exemestane interacts with Asp(309) based on its co-crystal structure with the enzyme, its binding is pH-dependent. Aromatase binding to anastrozole is pH-independent, consistent with the hypothesis that this ligand exploits a distinct set of interactions in the active site. In summary, we assign the apparent pK(a) of 8.2 observed for androstenedione binding to the side chain of Asp(309). To our knowledge, this work represents the first experimental assignment of a pK(a) value to a residue in a cytochrome P450. This value is in agreement with theoretical calculations (7.7-8.1) despite the reliance of the computational methods on the conformational snapshots provided by crystal structures.

Spin equilibrium and O2-binding kinetics of Mycobacterium tuberculosis CYP51 with mutations in the histidine-threonine dyad.

The acidic residues of the "acid-alcohol pair" in CYP51 enzymes are uniformly replaced with histidine. Herein, we adopt the Mycobacterium tuberculosis (mt) enzyme as a model system to investigate these residues' roles in finely tuning the heme conformation, iron spin state, and formation and decay of the oxyferrous enzyme. Properties of the mtCYP51 and the T260A, T260V, and H259A mutants were interrogated using UV-Vis and resonance Raman spectroscopies. Evidence supports that these mutations induce comprehensive changes in the heme environment. The heme iron spin states are differentially sensitive to the binding of the substrate, dihydrolanosterol (DHL). DHL and clotrimazole perturb the local environments of the heme vinyl and propionate substituents. Molecular dynamics (MD) simulations of the DHL-enzyme complexes support that the observed perturbations are attributable to changes in the DHL binding mode. Furthermore, the rates of the oxyferrous formation were measured using stopped-flow methods. These studies demonstrate that both HT mutations and DHL modulate the rates of oxyferrous formation. Paradoxically, the binding rate to the H259A mutant-DHL complex was approximately four-fold that of mtCYP51, a phenomenon that is predicted to result from the creation of an additional diffusion channel from loss of the H259-E173 ion pair in the mutant. Oxyferrous enzyme auto-oxidation rates were relatively constant, with the exception of the T260V-DHL complex. MD simulations lead us to speculate that this behavior may be attributed to the distortion of the heme macrocycle by the substrate.