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.

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.

Unfolding a molecular trefoil derived from a zwitterionic metallopeptide to form self-assembled nanostructures.

While used extensively by nature to control the geometry of protein structures, and dynamics of proteins, such as self-organization, hydration forces and ionic interactions received less attention for controlling the behaviour of small molecules. Here we describe the synthesis and characterization of a novel zwitterionic metallopeptide consisting of a cationic core and three distal anionic groups linked by self-assembling peptide motifs. 2D NMR spectra, total correlated spectroscopy and nuclear Overhauser effect spectroscopy, show that the molecule exhibits a three-fold rotational symmetry and adopts a folded conformation in dimethyl sulfoxide due to Coulombic forces. When hydrated in water, the molecule unfolds to act as a self-assembling building block of supramolecular nanostructures. By combining ionic interactions with the unique geometry from metal complex and hydrophobic interactions from simple peptides, we demonstrate a new and effective way to design molecules for smart materials through mimicking a sophisticated biofunctional system using a conformational switch.

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.

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.

Spring-loading the active site of cytochrome P450cam.

A hydrogen bond network has been identified that adjusts protein-substrate contacts in cytochrome P450(cam) (CYP101A1). Replacing the native substrate camphor with adamantanone or norcamphor causes perturbations in NMR-detected NH correlations assigned to the network, which includes portions of a ß sheet and an adjacent helix that is remote from the active site. A mutation in this helix reduces enzyme efficiency and perturbs the extent of substrate-induced spin state changes at the haem iron that accompany substrate binding. In turn, the magnitude of the spin state changes induced by alternate substrate binding parallel the NMR-detected perturbations observed near the haem in the enzyme active site.

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.

Redox-dependent dynamics in cytochrome P450cam.

Local protein backbone dynamics of the camphor hydroxylase cytochrome P450(cam) (CYP101) depend upon the oxidation and ligation state of the heme iron. (1)H-(15)N correlation nuclear magnetic resonance experiments were used to compare backbone dynamics of oxidized and reduced forms of this 414-residue metalloenzyme via hydrogen-deuterium exchange kinetics (H-D exchange) and (15)N relaxation measurements, and these results are compared with previously published results obtained by H-D exchange mass spectrometry. In general, the reduced enzyme exhibits lower-amplitude motions of secondary structural features than the oxidized enzyme on all of the time scales accessible to these experiments, and these differences are more pronounced in regions of the enzyme involved in substrate access to the active site (B' helix and beta3 and beta5 sheets) and binding of putidaredoxin (C and L helices), the iron-sulfur protein that acts as the effector and reductant of CYP101 in vivo. These results are interpreted in terms of local structural effects of changes in the heme oxidation state, and the relevance of the observed effects to the enzyme mechanism is discussed.

Structural and dynamic implications of an effector-induced backbone amide cis-trans isomerization in cytochrome P450cam.

Experimental evidence has been provided for a functionally relevant cis-trans isomerization of the Ile88-Pro89 peptide bond in cytochrome P450(cam) (CYP101). The isomerization is proposed to be a key element of the structural reorganization leading to the catalytically competent form of CYP101 upon binding of the effector protein putidaredoxin (Pdx). A detailed comparison of the results of molecular dynamics simulations on the cis and trans conformations of substrate- and carbonmonoxy-bound ferrous CYP101 with sequence-specific Pdx-induced structural perturbations identified by nuclear magnetic resonance is presented, providing insight into the structural and dynamic consequences of the isomerization. The mechanical coupling between the Pdx binding site on the proximal face of CYP101 and the site of isomerization is described.

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.

Specific effects of potassium ion binding on wild-type and L358P cytochrome P450cam.

The camphor monoxygenase cytochrome P450cam (CYP101) requires potassium ion (K+) to drive formation of the characteristic high-spin state of the heme Fe+3 upon substrate binding. Amide 1H, 15N correlations in perdeuterated [U-15N] CYP101 were monitored as a function of K+ concentration by 2D-TROSY-HSQC in both camphor-bound oxidized (CYP-S) and camphor- and CO-bound reduced CYP101 (CYP-S-CO). In both forms, K+-induced spectral perturbations are detected in the vicinity of the K+ binding site proposed from crystallographic structures, but are larger and more widespread structurally in CYP-S than in CYP-S-CO. In CYP-S-CO, K+-induced perturbations occur primarily near the proposed K+ binding site in the B-B' loop and B' helix, which are also perturbed by binding of effector, putidaredoxin (Pdx). The spectral effects of K+ binding in CYP-S-CO oppose those observed upon Pdxr titration. However, Pdxr titration of CYP-S-CO in the absence of K+ results in multiple conformations. The spin-state equilibrium in the L358P mutant of CYP101 is more sensitive to K+ concentration than WT CYP101, consistent with a hypothesis that L358P preferentially populates conformations enforced by Pdx binding in WT CYP101. Thallium(I), a K+ mimic, minimizes the effects of Pdx titration on the NMR spectrum of CYP-S-CO, but is competent to replace K+ in driving the formation of high-spin CYP-S. These observations suggest that the role of K+ is to stabilize conformers of CYP-S that drive the spin-state change prior to the first electron transfer, and that K+ stabilizes the CYP-S-CO conformer that interacts with Pdx. However, upon binding of Pdx, further conformational changes occur that disfavor K+ binding.

One protein, two enzymes revisited: a structural entropy switch interconverts the two isoforms of acireductone dioxygenase.

Acireductone dioxygenase (ARD) catalyzes different reactions between O2 and 1,2-dihydroxy-3-oxo-5-(methylthio)pent-1-ene (acireductone) depending upon the metal bound in the active site. Ni2+ -ARD cleaves acireductone to formate, CO and methylthiopropionate. If Fe2+ is bound (ARD'), the same substrates yield methylthioketobutyrate and formate. The two forms differ in structure, and are chromatographically separable. Paramagnetism of Fe2+ renders the active site of ARD' inaccessible to standard NMR methods. The structure of ARD' has been determined using Fe2+ binding parameters determined by X-ray absorption spectroscopy and NMR restraints from H98S ARD, a metal-free diamagnetic protein that is isostructural with ARD'. ARD' retains the beta-sandwich fold of ARD, but a structural entropy switch increases order at one end of a two-helix system that bisects the beta-sandwich and decreases order at the other upon interconversion of ARD and ARD', causing loss of the C-terminal helix in ARD' and rearrangements of residues involved in substrate orientation in the active site.

Comparison of the complexes formed by cytochrome P450cam with cytochrome b5 and putidaredoxin, two effectors of camphor hydroxylase activity.

Structural perturbations in cytochrome P450cam (CYP101) induced by the soluble fragment of cytochrome b5, a nonphysiological effector of CYP101, were investigated by NMR spectroscopy and compared with the perturbations induced by the physiological reductant and effector putidaredoxin (Pdx). Chemical shifts of perdeuterated [U-15N]CYP101 backbone amide (NH) resonances were monitored as a function of cytochrome b5 concentration by 1H-15N TROSY-HSQC experiments. The association of cytochrome b5 with the reduced CYP101-camphor-carbon monoxide complex (CYP-S-CO) perturbs many of the same resonances that Pdx does, including regions of the CYP101 molecule implicated in substrate access and orientation. The perturbations are smaller in magnitude than those observed with Pdx(r) due to a lower binding affinity (a Kd of 13 +/- 3 mM, for the reduced cytochrome b5-CYP-S-CO complex compared to a Kd of 26 +/- 12 microM for the Pdx-CYP-S-CO complex). The results are in accord with our previous suggestion that the observed perturbations are related to effector activity and support the proposal that the primary role of the effector is to populate the active conformation of CYP101 to prevent uncoupling [Pochapsky, S. S., et al. (2003) Biochemistry 42, 5649-5656]. A titratable perturbation is observed at the 1H resonance of the 8-CH3 group of CYP101-bound camphor upon addition of cytochrome b5, a phenomenon also associated with the formation of the CYP101 x Pdx complex, albeit with larger perturbations [Wei, J. Y., et al. (2005) J. Am. Chem. Soc. 127, 6974-6976]. The effector activity of the particular rat cytochrome b5 construct used for NMR studies was confirmed by monitoring the enzymatic turnover that yielded 5-exo-hydroxycamphor using gas chromatography and mass spectrometry. Finally, the common features of the perturbations observed in the NMR spectra of the two complexes are discussed, and their relevance to effector activity is considered.

Detection of a high-barrier conformational change in the active site of cytochrome P450cam upon binding of putidaredoxin.

The orientation of the substrate camphor in the active site of reduced CO-bound cytochrome P450cam (CYP101) as a function of reduced putidaredoxin (Pdxr) addition has been examined by NMR using perdeuterated CYP101 and perdeuterated Pdx as well as isotopically labeled d-camphor. This permits the 1H resonances of CYP101-bound camphor to be observed without interference from the signals of CYP101 or Pdx and confirms assignments of the methyl signals of camphor in the bound form. The Cys4Fe2S2 ferredoxin Pdx is the physiological redox partner and effector of CYP101. The addition of Pdx to the reduced CYP101-camphor-CO complex results in a conformational selection that is slow on the chemical shift time scale with spectral effects observed primarily at the 8-CH3 group of the camphor. The camphor signals are ring current shifted by the heme, and for the 9- and 10-CH3 resonances, these shifts are reasonably well predicted by ring current calculations from the crystal structure of CO-bound CYP101. However, in the absence of Pdx, the 8-CH3 resonance of CYP101-bound camphor is observed at considerably higher field than predicted. Dynamic simulations using ring current shift restraints generated a structure with low chemical shift violations in which the hydrogen bond between the camphor carbonyl oxygen and the OH of Tyr96 is lost, and an expansion of the active site takes place that permits reorientation of the camphor within the active site.

A model for effector activity in a highly specific biological electron transfer complex: the cytochrome P450(cam)-putidaredoxin couple.

The camphor hydroxylase cytochrome P450(cam) (CYP101) catalyzes the 5-exo hydroxylation of camphor in the first step of camphor catabolism by Pseudomonas putida. CYP101 forms a specific electron transfer complex with its physiological reductant, the Cys(4)Fe(2)S(2) ferredoxin putidaredoxin (Pdx). Pdx, along with other proteins and small molecules, has also been shown to be an effector for turnover by CYP101. Multidimensional nuclear magnetic resonance (NMR) techniques have been used to make extensive sequential (1)H, (15)N, and (13)C resonance assignments in CYP101 that permit a more complete characterization of the complex formed by CYP101 and Pdx. NMR-detected perturbations in CYP101 upon Pdx binding encompass regions of the CYP101 remote from the putative Pdx binding site, including in particular a region of the CYP101 molecule that has been implicated in substrate access to the active site via dynamical processes. A model for effector activity is proposed in which the primary role of the effector is to prevent uncoupling (formation of reduced oxo species without formation of hydroxycamphor) by enforcing conformations of CYP101 that prevent loss of substrate and/or intermediates prior to turnover. A secondary role could also be to enforce conformations that permit efficient proton transfer into the active site for coupled proton/electron transfer.

Modeling and experiment yields the structure of acireductone dioxygenase from Klebsiella pneumoniae.

Here we report the structure of acireductone dioxygenase (ARD), the first determined for a new family of metalloenzymes. ARD represents a branch point in the methionine salvage pathway leading from methylthioadenosine to methionine and has been shown to catalyze different reactions depending on the type of metal ion bound in the active site. The solution structure of nickel-containing ARD (Ni-ARD) was determined using NMR methods. X-ray absorption spectroscopy, assignment of hyperfine shifted NMR resonances and conserved domain homology were used to model the metal-binding site because of the paramagnetism of the bound Ni2+. Although there is no structure in the Protein Data Bank within 3 A r.m.s deviation of that of Ni-ARD, the enzyme active site is located in a conserved double-stranded b-helix domain. Furthermore, the proposed Ni-ARD active site shows significant post-facto structural homology to the active sites of several metalloenzymes in the cupin superfamily.

Rapid recycle (13)C',(15)N and (13)C,(13)C' heteronuclear and homonuclear multiple quantum coherence detection for resonance assignments in paramagnetic proteins: example of Ni(2+)-containing acireductone dioxygenase.

NMR resonance assignments in the vicinity of paramagnetic metals in proteins are often difficult or impossible to make using conventional 1H detected 2-D and 3-D methods due to paramagnetic line broadening. The applicability of 13Calpha{13C'} and 13C'{15N} multiple quantum coherence methods for residue-specific assignments of resonances near paramagnetic centers is described, using the Ni2+-containing enzyme acireductone dioxygenase as an example.

Comparison of functional domains in vertebrate-type ferredoxins.

The vertebrate-type Cys(4)Fe(2)S(2) ferredoxins are a class of small acidic proteins that typically act as electron shuttles between NAD(P)H-dependent reductases and monoxygenases, particularly cytochromes P450. Nuclear magnetic resonance assignments and detailed analysis of nuclear Overhauser effects permit the direct comparison of the functional C-terminal domains of three vertebrate-type ferredoxins, the mammalian adrenodoxin (Adx) and the bacterial ferredoxins putidaredoxin (Pdx) and terpredoxin (Tdx). In particular, homologous hydrogen-bonding networks involving a conserved basic residue (His 49 in Pdx, His 56 in Adx, Arg 49 in Tdx) are detailed. This hydrogen bond network appears to play a role in the mechanical transmission of redox-dependent conformational and dynamic changes from the iron-sulfur binding loop to the C-terminal domain. Hydrogen/deuterium exchange measurements have been made in Adx as a function of oxidation state for comparison with previous studies of Pdx and Tdx. The results of these measurements highlight the importance of the conserved basic residue in the linkage between oxidation state and protein dynamics. Finally, a series of mutations have been made in the C-terminal domain of Pdx, including one, Y51F, that disrupts the proposed hydrogen-bonding network without perturbing steric and hydrophobic interactions in the functional domain. Although the mutant is considerably destabilized with respect to wild-type Pdx, relatively unperturbed chemical shifts for residues near the site of the mutation and NOEs between water and Phe 51 suggest that the network is reconstituted with a solvent water in place of the tyrosine hydroxyl group in this mutant.