Tracking protein-protein interactions by NMR: conformational selection in human steroidogenic cytochrome P450 CYP17A1 induced by cytochrome b5.

The human steroidogenic cytochrome P450 CYP17A1 catalyzes two types of reactions in the biosynthetic pathway leading from pregnenolone to testosterone and several other steroid hormones. The first is the hydroxylation of pregnenolone or progesterone to the corresponding 17α-hydroxy steroid, followed by a lyase reaction that converts these 17α-hydroxy intermediates to the androgens dehydroepiandrosterone and androstenedione, respectively. cytochrome b5 (cytb5) is known to act as both an effector and electron donor for the lyase oxidations, markedly stimulating the rate of the lyase reaction in its presence relative to the rate in its absence. Extensive sequential backbone 1H,15N and 13C nuclear magnetic resonance assignments have now been made for oxidized CYP17A1 bound to the prostate cancer drug and inhibitor abiraterone. This is the first eukaryotic P450 for which such assignments are now available. These assignments allow more complete interpretation of the structural perturbations observed upon cytb5 addition. Possible mechanism(s) for the effector activity of cytb5 are discussed in light of this new information.

Hydroxylation Regiochemistry Is Robust to Active Site Mutations in Cytochrome P450cam (CYP101A1).

Cytochrome P450cam (CYP101A1) catalyzes the hydroxylation of d-camphor by molecular oxygen. The enzyme-catalyzed hydroxylation exhibits a high degree of regioselectivity and stereoselectivity, with a single major product, d-5-exo-hydroxycamphor, suggesting that the substrate is oriented to facilitate this specificity. In previous work, we used an elastic network model and perturbation response scanning to show that normal deformation modes of the enzyme structure are highly responsive not only to the presence of a substrate but also to the substrate orientation. This work examines the effects of mutations near the active site on substrate localization and orientation. The investigated mutations were designed to promote a change in substrate orientation and/or location that might give rise to different hydroxylation products, while maintaining the same carbon and oxygen atom balances as in the wild type (WT) enzyme. Computational experiments and parallel in vitro site-directed mutations of CYP101A1 were used to examine reaction products and enzyme activity. 1H-15N TROSY-HSQC correlation maps were used to compare the computational results with detectable perturbations in the enzyme structure and dynamics. We found that all of the mutant enzymes retained the same regio- and stereospecificity of hydroxylation as the WT enzyme, with varying degrees of efficiency, which suggests that large portions of the enzyme have been subjected to evolutionary pressure to arrive at the appropriate sequence-structure combination for efficient 5-exo hydroxylation of camphor.

A Model for the Solution Structure of Human Fe(II)-Bound Acireductone Dioxygenase and Interactions with the Regulatory Domain of Matrix Metalloproteinase I (MMP-I).

The metalloenzyme acireductone dioxygenase (ARD) shows metal-dependent physical and enzymatic activities depending upon the metal bound in the active site. The Fe(II)-bound enzyme catalyzes the penultimate step of the methionine salvage pathway, converting 1,2-dihydroxy-5-(methylthio)pent-1-en-3-one (acireductone) into formate and the ketoacid precursor of methionine, 2-keto-4-thiomethyl-2-oxobutanoate, using O2 as the oxidant. If Ni(II) is bound, an off-pathway shunt occurs, producing 3-methylthiopropionate, formate, and carbon monoxide from the same acireductone substrate. The solution structure of the Fe(II)-bound human enzyme, HsARD, is described and compared with the structures of Ni-bound forms of the closely related mouse enzyme, MmARD. Potential rationales for the different reactivities of the two isoforms are discussed. The human enzyme has been found to regulate the activity of matrix metalloproteinase I (MMP-I), which is involved in tumor metastasis, by binding the cytoplasmic transmembrane tail peptide of MMP-I. Nuclear magnetic resonance titration of HsARD with the MMP-I tail peptide permits identification of the peptide binding site on HsARD, a cleft anterior to the metal binding site adjacent to a dynamic proline-rich loop.

What Your Crystal Structure Will Not Tell You about Enzyme Function.

Enzyme function requires that enzyme structures be dynamic. Substrate binding, product release, and transition state stabilization typically involve different enzyme conformers. Furthermore, in multistep enzyme-catalyzed reactions, more than one enzyme conformation may be important for stabilizing different transition states. While X-ray crystallography provides the most detailed structural information of any current methodology, X-ray crystal structures of enzymes capture only those conformations that fit into the crystal lattice, which may or may not be relevant to function. Solution nuclear magnetic resonance (NMR) methods can provide an alternative approach to characterizing enzymes under nonperturbing and controllable conditions, allowing one to identify and localize dynamic processes that are important to function. However, many enzymes are too large for standard approaches to making sequential resonance assignments, a critical first step in analyzing and interpreting the wealth of information inherent in NMR spectra. This Account describes our long-standing NMR-based research into structural and dynamic aspects of function in the cytochrome P450 monooxygenase superfamily. These heme-containing enzymes typically catalyze the oxidation of unactivated C-H and C═C bonds in a multitude of substrates, often with complete regio- and stereospecificity. Over 600 000 genes in GenBank have been assigned to P450s, yet all known P450 structures exhibit a highly conserved and unique fold. This combination of functional and structural conservation with a vast substrate clientele, each substrate having multiple possible sites for oxidation, makes the P450s a unique target for understanding the role of enzyme structure and dynamics in determining a particular substrate-product combination. P450s are large by solution NMR standards, requiring us to develop specialized approaches for making sequential resonance assignments and interpreting the spectral changes that occur as a function of changing conditions (e.g., oxidation and spin state changes, ligand, substrate or effector binding). Solution conformations are characterized by the fitting of residual dipolar couplings (RDCs) measured for sequence-specifically assigned amide N-H correlations to alignment tensors optimized in the course of restrained molecular dynamics (MD) simulations. The conformational ensembles obtained by such RDC-restrained simulations, which we call "soft annealing", are then tested by site-directed mutation and spectroscopic and activity assays for relevance. These efforts have gained us insights into cryptic conformational changes associated with substrate and redox partner binding that were not suspected from crystal structures, but were shown by subsequent work to be relevant to function. Furthermore, it appears that many of these changes can be generalized to P450s besides those that we have characterized, providing guidance for enzyme engineering efforts. While past research was primarily directed at the more tractable prokaryotic P450s, our current efforts are aimed at medically relevant human enzymes, including CYP17A1, CYP2D6, and CYP3A4.

The Metal Drives the Chemistry: Dual Functions of Acireductone Dioxygenase.

Acireductone dioxygenase (ARD) from the methionine salvage pathway (MSP) is a unique enzyme that exhibits dual chemistry determined solely by the identity of the divalent transition-metal ion (Fe2+ or Ni2+) in the active site. The Fe2+-containing isozyme catalyzes the on-pathway reaction using substrates 1,2-dihydroxy-3-keto-5-methylthiopent-1-ene (acireductone) and dioxygen to generate formate and the ketoacid precursor of methionine, 2-keto-4-methylthiobutyrate, whereas the Ni2+-containing isozyme catalyzes an off-pathway shunt with the same substrates, generating methylthiopropionate, carbon monoxide, and formate. The dual chemistry of ARD was originally discovered in the bacterium Klebsiella oxytoca, but it has recently been shown that mammalian ARD enzymes (mouse and human) are also capable of catalyzing metal-dependent dual chemistry in vitro. This is particularly interesting, since carbon monoxide, one of the products of off-pathway reaction, has been identified as an antiapoptotic molecule in mammals. In addition, several biochemical and genetic studies have indicated an inhibitory role of human ARD in cancer. This comprehensive review describes the biochemical and structural characterization of the ARD family, the proposed experimental and theoretical approaches to establishing mechanisms for the dual chemistry, insights into the mechanism based on comparison with structurally and functionally similar enzymes, and the applications of this research to the field of artificial metalloenzymes and synthetic biology.

NADH reduction of nitroaromatics as a probe for residual ferric form high-spin in a cytochrome P450.

The existence of a substrate-sensitive equilibrium between high spin (S=5/2) and low spin (S=1/2) ferric iron is a well-established phenomenon in the cytochrome P450 (CYP) superfamily, although its origins are still a subject of discussion. A series of mutations that strongly perturb the spin state equilibrium in the camphor hydroxylase CYP101A1 were recently described (Colthart et al., Sci. Rep. 6, 22035 (2016)). Wild type CYP101A1 as well as some CYP101A1 mutants are herein shown to be capable of catalyzing the reduction of nitroacetophenones by NADH to the corresponding anilino compounds (nitroreductase or NRase activity). The distinguishing characteristic between those mutants that catalyze the reduction and those that cannot appears to be the extent to which residual high spin form exists in the absence of the native substrate d-camphor, with those showing the largest spin state shifts upon camphor binding also exhibiting NRase activity. Optical and EPR spectroscopy was used to further examine these phenomena. These results suggest that reduction of nitroaromatics may provide a useful probe of residual high spin states in the CYP superfamily. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Solution Conformations and Dynamics of Substrate-Bound Cytochrome P450 MycG.

MycG is a P450 monooxygenase that catalyzes the sequential hydroxylation and epoxidation of mycinamicin IV (M-IV), the last two steps in the biosynthesis of mycinamicin II, a macrolide antibiotic isolated from Micromonospora griseorubida. The crystal structure of MycG with M-IV bound was previously determined but showed the bound substrate in an orientation that did not rationalize the observed regiochemistry of M-IV hydroxylation. Nuclear magnetic resonance paramagnetic relaxation enhancements provided evidence of an orientation of M-IV in the MycG active site more compatible with the observed chemistry, but substrate-induced changes in the enzyme structure were not characterized. We now describe the use of amide 1H-15N residual dipolar couplings as experimental restraints in solvated "soft annealing" molecular dynamics simulations to generate solution structural ensembles of M-IV-bound MycG. Chemical shift perturbations, hydrogen-deuterium exchange, and 15N relaxation behavior provide insight into the dynamic and electronic perturbations in the MycG structure in response to M-IV binding. The solution and crystallographic structures are compared, and the possibility that the crystallographic orientation of bound M-IV represents an inhibitory mode is discussed.

Metal-Dependent Function of a Mammalian Acireductone Dioxygenase.

The two acireductone dioxygenase (ARD) isozymes from the methionine salvage pathway of Klebsiella oxytoca are the only known pair of naturally occurring metalloenzymes with distinct chemical and physical properties determined solely by the identity of the divalent transition metal ion (Fe(2+) or Ni(2+)) in the active site. We now show that this dual chemistry can also occur in mammals. ARD from Mus musculus (MmARD) was studied to relate the metal ion identity and three-dimensional structure to enzyme function. The iron-containing isozyme catalyzes the cleavage of 1,2-dihydroxy-3-keto-5-(thiomethyl)pent-1-ene (acireductone) by O2 to formate and the ketoacid precursor of methionine, which is the penultimate step in methionine salvage. The nickel-bound form of ARD catalyzes an off-pathway reaction resulting in formate, carbon monoxide (CO), and 3-(thiomethyl) propionate. Recombinant MmARD was expressed and purified to obtain a homogeneous enzyme with a single transition metal ion bound. The Fe(2+)-bound protein, which shows about 10-fold higher activity than that of others, catalyzes on-pathway chemistry, whereas the Ni(2+), Co(2+), or Mn(2+) forms exhibit off-pathway chemistry, as has been seen with ARD from Klebsiella. Thermal stability of the isozymes is strongly affected by the metal ion identity, with Ni(2+)-bound MmARD being the most stable, followed by Co(2+) and Fe(2+), and Mn(2+)-bound ARD being the least stable. Ni(2+)- and Co(2+)-bound MmARD were crystallized, and the structures of the two proteins found to be similar. Enzyme-ligand complexes provide insight into substrate binding, metal coordination, and the catalytic mechanism.

A soluble α-synuclein construct forms a dynamic tetramer.

A heterologously expressed form of the human Parkinson disease-associated protein α-synuclein with a 10-residue N-terminal extension is shown to form a stable tetramer in the absence of lipid bilayers or micelles. Sequential NMR assignments, intramonomer nuclear Overhauser effects, and circular dichroism spectra are consistent with transient formation of α-helices in the first 100 N-terminal residues of the 140-residue α-synuclein sequence. Total phosphorus analysis indicates that phospholipids are not associated with the tetramer as isolated, and chemical cross-linking experiments confirm that the tetramer is the highest-order oligomer present at NMR sample concentrations. Image reconstruction from electron micrographs indicates that a symmetric oligomer is present, with three- or fourfold symmetry. Thermal unfolding experiments indicate that a hydrophobic core is present in the tetramer. A dynamic model for the tetramer structure is proposed, based on expected close association of the amphipathic central helices observed in the previously described micelle-associated "hairpin" structure of α-synuclein.

Substrate recognition by the multifunctional cytochrome P450 MycG in mycinamicin hydroxylation and epoxidation reactions.

The majority of characterized cytochrome P450 enzymes in actinomycete secondary metabolic pathways are strictly substrate-, regio-, and stereo-specific. Examples of multifunctional biosynthetic cytochromes P450 with broader substrate and regio-specificity are growing in number and are of particular interest for biosynthetic and chemoenzymatic applications. MycG is among the first P450 monooxygenases characterized that catalyzes both hydroxylation and epoxidation reactions in the final biosynthetic steps, leading to oxidative tailoring of the 16-membered ring macrolide antibiotic mycinamicin II in the actinomycete Micromonospora griseorubida. The ordering of steps to complete the biosynthetic process involves a complex substrate recognition pattern by the enzyme and interplay between three tailoring modifications as follows: glycosylation, methylation, and oxidation. To understand the catalytic properties of MycG, we structurally characterized the ligand-free enzyme and its complexes with three native metabolites. These include substrates mycinamicin IV and V and their biosynthetic precursor mycinamicin III, which carries the monomethoxy sugar javose instead of the dimethoxylated sugar mycinose. The two methoxy groups of mycinose serve as sensors that mediate initial recognition to discriminate between closely related substrates in the post-polyketide oxidative tailoring of mycinamicin metabolites. Because x-ray structures alone did not explain the mechanisms of macrolide hydroxylation and epoxidation, paramagnetic NMR relaxation measurements were conducted. Molecular modeling based on these data indicates that in solution substrate may penetrate the active site sufficiently to place the abstracted hydrogen atom of mycinamicin IV within 6 Å of the heme iron and ~4 Å of the oxygen of iron-ligated water.

Conformational heterogeneity suggests multiple substrate binding modes in CYP106A2.

CYP106A2 (cytochrome P450meg) is a bacterial enzyme originally isolated from B. megaterium, and has been shown to hydroxylate a wide variety of substrates, including steroids. The regio- and stereochemistry of CYP106A2 hydroxylation has been shown to be dependent on a variety of factors, and hydroxylation often occurs at more than one site and/or with lack of stereospecificity for some substrates. Comprehensive backbone 15N, 1H and 13C resonance assignments based on multidimensional nuclear magnetic resonance (NMR) experiments performed with uniform and selective isotopically labeled CYP106A2 samples are reported herein, and broadening and splitting of resonances assigned to regions of the enzyme shown to be affected by substrate binding in other P450 enzymes indicate that substrate binding does not reduce structural heterogeneity as has been observed previously in P450 enzymes CYP101A1 and MycG. Paramagnetic relaxation enhancement (PRE) due to proximity between substrate protons and the heme iron were measured for three different substrates, and the relatively uniform nature of the PREs support the proposal that multiple substrate binding modes are occupied at saturating substrate concentrations.

Selective steroidogenic cytochrome P450 haem iron ligation by steroid-derived isonitriles.

Alkyl isonitriles, R-NC, have previously been shown to ligate the heme (haem) iron of cytochromes P450 in both accessible oxidation states (ferrous, Fe2+, and ferric, Fe3+). Herein, the preparation of four steroid-derived isonitriles and their interactions with several P450s, including the steroidogenic CYP17A1 and CYP106A2, as well as the more promiscuous drug metabolizers CYP3A4 and CYP2D6, is described. It was found that successful ligation of the heme iron by the isonitrile functionality for a given P450 depends on both the position and stereochemistry of the isonitrile on the steroid skeleton. Spectral studies indicate that isonitrile ligation of the ferric heme is stable upon reduction to the ferrous form, with reoxidation resulting in the original complex. A crystallographic structure of CYP17A1 with an isonitrile derived from pregnanalone further confirmed the interaction and identified the absolute stereochemistry of the bound species.

Solution structural ensembles of substrate-free cytochrome P450(cam).

Removal of substrate (+)-camphor from the active site of cytochrome P450(cam) (CYP101A1) results in nuclear magnetic resonance-detected perturbations in multiple regions of the enzyme. The (1)H-(15)N correlation map of substrate-free diamagnetic Fe(II) CO-bound CYP101A permits these perturbations to be mapped onto the solution structure of the enzyme. Residual dipolar couplings (RDCs) were measured for (15)N-(1)H amide pairs in two independent alignment media for the substrate-free enzyme and used as restraints in solvated molecular dynamics (MD) simulations to generate an ensemble of best-fit structures of the substrate-free enzyme in solution. Nuclear magnetic resonance-detected chemical shift perturbations reflect changes in the electronic environment of the NH pairs, such as hydrogen bonding and ring current shifts, and are observed for residues in the active site as well as in hinge regions between secondary structural features. RDCs provide information about relative orientations of secondary structures, and RDC-restrained MD simulations indicate that portions of a ß-rich region adjacent to the active site shift so as to partially occupy the vacancy left by removal of the substrate. The accessible volume of the active site is reduced in the substrate-free enzyme relative to the substrate-bound structure calculated using the same methods. Both symmetric and asymmetric broadening of multiple resonances observed upon substrate removal as well as localized increased errors in RDC fits suggest that an ensemble of enzyme conformations are present in the substrate-free form.

Acireductone dioxygenase 1 (ARD1) is an effector of the heterotrimeric G protein beta subunit in Arabidopsis.

Heterotrimeric G protein complexes are conserved from plants to mammals, but the complexity of each system varies. Arabidopsis thaliana contains one Gα, one Gß (AGB1), and at least three Gγ subunits, allowing it to form three versions of the heterotrimer. This plant model is ideal for genetic studies because mammalian systems contain hundreds of unique heterotrimers. The activation of these complexes promotes interactions between both the Gα subunit and the Gßγ dimer with enzymes and scaffolds to propagate signaling to the cytoplasm. However, although effectors of Gα and Gß are known in mammals, no Gß effectors were previously known in plants. Toward identifying AGB1 effectors, we genetically screened for dominant mutations that suppress Gß-null mutant (agb1-2) phenotypes. We found that overexpression of acireductone dioxygenase 1 (ARD1) suppresses the 2-day-old etiolated phenotype of agb1-2. ARD1 is homologous to prokaryotic and eukaryotic ARD proteins; one function of ARDs is to operate in the methionine salvage pathway. We show here that ARD1 is an active metalloenzyme, and AGB1 and ARD1 both control embryonic hypocotyl length by modulating cell division; they also may contribute to the production of ethylene, a product of the methionine salvage pathway. ARD1 physically interacts with AGB1, and ARD enzymatic activity is stimulated by AGB1 in vitro. The binding interface on AGB1 was deduced using a comparative evolutionary approach and tested using recombinant AGB1 mutants. A possible mechanism for AGB1 activation of ARD1 activity was tested using directed mutations in a loop near the substrate-binding site.

Experimentally restrained molecular dynamics simulations for characterizing the open states of cytochrome P450cam.

Residual dipolar couplings (RDCs) were used as restraints in fully solvated molecular dynamics simulations of reduced substrate- and carbonmonoxy-bound cytochrome P450(cam) (CYP101A1), a 414-residue soluble monomeric heme-containing camphor monooxygenase from the soil bacterium Pseudomonas putida. The (1)D(NH) residual dipolar couplings used as restraints were measured in two independent alignment media. A soft annealing protocol was used to heat the starting structures while incorporating the RDC restraints. After production dynamics, structures with the lowest total violation energies for RDC restraints were extracted to identify ensembles of conformers accessible to the enzyme in solution. The simulations result in substrate orientations different from that seen in crystallographic structures and a more open and accessible enzyme active site and largely support previously reported differences between the open and closed states of CYP101A1.

A dynamic understanding of cytochrome P450 structure and function through solution NMR.

Many economically important biosyntheses incorporate regiospecific and stereospecific oxidations at unactivated carbons. Such oxidations are commonly catalyzed by cytochrome P450 monooxygenases, heme-containing enzymes that activate molecular oxygen while selectively binding and orienting the substrate for reaction. Despite the plethora of P450-catalyzed reactions, the P450 fold is highly conserved, and static structures are often insufficient for characterizing conformational states that contribute to specificity. High-resolution solution nuclear magnetic resonance (NMR) offers insights into dynamic processes and conformational changes that are required of a P450 in order to attain the combination of specificity and efficiency required for these reactions.

Communication between the zinc and nickel sites in dimeric HypA: metal recognition and pH sensing.

Helicobacter pylori , a pathogen that colonizes the human stomach, requires the nickel-containing metalloenzymes urease and NiFe-hydrogenase to survive this low pH environment. The maturation of both enzymes depends on the metallochaperone, HypA. HypA contains two metal sites, an intrinsic zinc site and a low-affinity nickel binding site. X-ray absorption spectroscopy (XAS) shows that the structure of the intrinsic zinc site of HypA is dynamic and able to sense both nickel loading and pH changes. At pH 6.3, an internal pH that occurs during acid shock, the zinc site undergoes unprecedented ligand substitutions to convert from a Zn(Cys)(4) site to a Zn(His)(2)(Cys)(2) site. NMR spectroscopy shows that binding of Ni(II) to HypA results in paramagnetic broadening of resonances near the N-terminus. NOEs between the beta-CH(2) protons of Zn cysteinyl ligands are consistent with a strand-swapped HypA dimer. Addition of nickel causes resonances from the zinc binding motif and other regions to double, indicating more than one conformation can exist in solution. Although the structure of the high-spin, 5-6 coordinate Ni(II) site is relatively unaffected by pH, the nickel binding stoichiometry is decreased from one per monomer to one per dimer at pH = 6.3. Mutation of any cysteine residue in the zinc binding motif results in a zinc site structure similar to that found for holo-WT-HypA at low pH and is unperturbed by the addition of nickel. Mutation of the histidines that flank the CXXC motifs results in a zinc site structure that is similar to holo-WT-HypA at neutral pH (Zn(Cys)(4)) and is no longer responsive to nickel binding or pH changes. Using an in vitro urease activity assay, it is shown that the recombinant protein is sufficient for recovery of urease activity in cell lysate from a HypA deletion mutant, and that mutations in the zinc-binding motif result in a decrease in recovered urease activity. The results are interpreted in terms of a model wherein HypA controls the flow of nickel traffic in the cell in response to nickel availability and pH.

A functional proline switch in cytochrome P450cam.

The two-protein complex between putidaredoxin (Pdx) and cytochrome P450(cam) (CYP101) is the catalytically competent species for camphor hydroxylation by CYP101. We detected a conformational change in CYP101 upon binding of Pdx that reorients bound camphor appropriately for hydroxylation. Experimental evidence shows that binding of Pdx converts a single X-proline amide bond in CYP101 from trans or distorted trans to cis. Mutation of proline 89 to isoleucine yields a mixture of both bound camphor orientations, that seen in Pdx-free and that seen in Pdx-bound CYP101. A mutation in CYP101 that destabilizes the cis conformer of the Ile 88-Pro 89 amide bond results in weaker binding of Pdx. This work provides direct experimental evidence for involvement of X-proline isomerization in enzyme function.