The role of water molecules in the association of cytochrome P450cam with putidaredoxin. An osmotic pressure study

Osmotic pressure effects identify dehydration upon cytochrome c-cytochrome c oxidase complex formation contributing to a specific electron pathway formation

In the electron transfer (ET) reaction from cytochrome c (Cyt c) to cytochrome c oxidase (CcO), we determined the number and sites of the hydration water released from the protein surface upon the formation of the ET complex by evaluating the osmotic pressure dependence of kinetics for the ET from Cyt c to CcO. We identified that ∼20 water molecules were dehydrated in complex formation under turnover conditions, and systematic Cyt c mutations in the interaction site for CcO revealed that nearly half of the released hydration water during the complexation were located around Ile81, one of the hydrophobic amino acid residues near the exposed heme periphery of Cyt c. Such a dehydration dominantly compensates for the entropy decrease due to the association of Cyt c with CcO, resulting in the entropy-driven ET reaction. The energetic analysis of the interprotein interactions in the ET complex predicted by the docking simulation suggested the formation of hydrophobic interaction sites surrounding the exposed heme periphery of Cyt c in the Cyt c-CcO interface (a 'molecular breakwater'). Such sites would contribute to the formation of the hydrophobic ET pathway from Cyt c to CcO by blocking water access from the bulk water phase.

Effects of Osmolytes on Ligand Binding to Dihydropteroate Synthase from Bacillus anthracis

Osmolyte interactions with ligands can affect their affinity for proteins and are dependent upon the cosolute and the functional groups of the ligand. Here, we explored ligand binding to Bacillus anthracis dihydropteroate synthase (BaDHPS) under osmotic stress conditions. Osmolyte effects were specific to the cosolute and ligand, suggesting interaction of the osmolytes with the free ligands in solution. The association rates of pterin pyrophosphate were mostly unaffected by the osmolytes, except for a 2-fold decrease in the presence of 1 M trehalose, while the dissociation rates decreased in most osmolyte solutions. The viscosity and dielectric constant of the solution did not correlate with the effects of the osmolytes. Experimental results were compared with predicted preferential interaction coefficients (Δμ₂₃/RT) between the osmolytes and ligands. The Δμ₂₃/RT were able to predict the experimental data for most of the osmolytes. Trehalose and proline effects did not correlate with the predicted values, indicating that these two osmolytes may affect binding in more complex ways than simple preferential interactions. Additionally, osmolytes weakly interacted with the sulfa drug sulfathiazole, which altered its affinity for BaDHPS, suggesting that these types of weak interactions can also impact drug binding. As osmolytes affect ligands binding to two different folate cycle enzymes (DHFRs and DHPS), we predicted how ligand binding to other folate cycle enzymes will be altered by the presence of osmolytes.

Photosynthetic fuel for heterologous enzymes: the role of electron carrier proteins

Plants, cyanobacteria, and algae generate a surplus of redox power through photosynthesis, which makes them attractive for biotechnological exploitations. While central metabolism consumes most of the energy, pathways introduced through metabolic engineering can also tap into this source of reducing power. Recent work on the metabolic engineering of photosynthetic organisms has shown that the electron carriers such as ferredoxin and flavodoxin can be used to couple heterologous enzymes to photosynthetic reducing power. Because these proteins have a plethora of interaction partners and rely on electrostatically steered complex formation, they form productive electron transfer complexes with non-native enzymes. A handful of examples demonstrate channeling of photosynthetic electrons to drive the activity of heterologous enzymes, and these focus mainly on hydrogenases and cytochrome P450s. However, competition from native pathways and inefficient electron transfer rates present major obstacles, which limit the productivity of heterologous reactions coupled to photosynthesis. We discuss specific approaches to address these bottlenecks and ensure high productivity of such enzymes in a photosynthetic context.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Thermodynamics and solvent linkage of macromolecule-ligand interactions

Binding involves two steps, desolvation and association. While water is ubiquitous and occurs at high concentration, it is typically ignored. In vitro experiments typically use infinite dilution conditions, while in vivo, the concentration of water is decreased due to the presence of high concentrations of molecules in the cellular milieu. This review discusses isothermal titration calorimetry approaches that address the role of water in binding. For example, use of D2O allows the contribution of solvent reorganization to the enthalpy component to be assessed. Further, the addition of osmolytes will decrease the water activity of a solution and allow effects on Ka to be determined. In most cases, binding becomes tighter in the presence of osmolytes as the desolvation penalty associated with binding is minimized. In other cases, the osmolytes prefer to interact with the ligand or protein, and if their removal is more difficult than shedding water, then binding can be weakened. These complicating layers can be discerned by different slopes in ln(Ka) vs osmolality plots and by differential scanning calorimetry in the presence of the osmolyte.

Investigation of osmolyte effects on FolM: comparison with other dihydrofolate reductases

A weak association between osmolytes and dihydrofolate (DHF) decreases the affinity of the substrate for the Escherichia coli chromosomal and R67 plasmid dihydrofolate reductase (DHFR) enzymes. To test whether the osmolyte-DHF association also interferes with binding of DHF to FolM, an E. coli enzyme that possesses weak DHFR activity, ligand binding was monitored in the presence of osmolytes. The affinity of FolM for DHF, measured by kcat/Km(DHF), was decreased by the addition of an osmolyte. Additionally, binding of the antifolate drug, methotrexate, to FolM was weakened by the addition of an osmolyte. The changes in ligand binding with water activity were unique for each osmolyte, indicating preferential interaction between the osmolyte and folate and its derivatives; however, additional evidence provided support for further interactions between FolM and osmolytes. Binding of the reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor to FolM was monitored by isothermal titration calorimetry as a control for protein-osmolyte association. In the presence of betaine (proposed to be the osmolyte most excluded from protein surfaces), the NADPH Kd decreased, consistent with dehydration effects. However, other osmolytes did not tighten binding to the cofactor. Rather, dimethyl sulfoxide (DMSO) had no effect on the NADPH Kd, while ethylene glycol and polyethylene glycol 400 weakened cofactor binding. Differential scanning calorimetry of FolM in the presence of osmolytes showed that both DMSO and ethylene glycol decreased the stability of FolM, while betaine increased the stability of the protein. These results suggest that some osmolytes can destabilize FolM by preferentially interacting with the protein. Further, these weak attractions can impede ligand binding. These various contributions have to be considered when interpreting osmotic pressure results.

Weak interactions between folate and osmolytes in solution

Previous osmotic stress studies on the role of solvent in two structurally unrelated dihydrofolate reductases (DHFRs) found weaker binding of dihydrofolate (DHF) to either enzyme in the presence of osmolytes. To explain these unusual results, weak interactions between DHF and osmolytes were proposed, with a competition between osmolyte and DHFR for DHF. High osmolyte concentrations will inhibit binding of the cognate pair. To evaluate this hypothesis, we devised a small molecule approach. Dimerization of folate, monitored by nuclear magnetic resonance, was weakened 2-3-fold upon addition of betaine or dimethyl sulfoxide (DMSO), supporting preferential interaction of either osmolyte with the monomer (as it possesses a larger surface area). Nuclear Overhauser effect (NOE) spectroscopy experiments found a positive NOE for the interaction of the C3'/C5' benzoyl ring protons with the C9 proton in buffer; however, a negative NOE was observed upon addition of betaine or DMSO. This change indicated a decreased tumbling rate, consistent with osmolyte interaction. Osmotic stress experiments also showed that betaine, DMSO, and sucrose preferentially interact with folate. Further, studies with the folate fragments, p-aminobenzoic acid and pterin 6-carboxylate, revealed interactions for both model compounds with betaine and sucrose. In contrast, DMSO was strongly excluded from the pterin ring but preferentially interacted with the p-aminobenzoyl moiety. These interactions are likely to be important in vivo because of the crowded conditions of the cell where weak contacts can more readily compete with specific binding interactions.

Thermodynamics and solvent effects on substrate and cofactor binding in Escherichia coli chromosomal dihydrofolate reductase

Chromosomal dihydrofolate reductase from Escherichia coli catalyzes the reduction of dihydrofolate to tetrahydrofolate using NADPH as a cofactor. The thermodynamics of ligand binding were examined using an isothermal titration calorimetry approach. Using buffers with different heats of ionization, zero to a small, fractional proton release was observed for dihydrofolate binding, while a proton was released upon NADP(+) binding. The role of water in binding was additionally monitored using a number of different osmolytes. Binding of NADP(+) is accompanied by the net release of ∼5-24 water molecules, with a dependence on the identity of the osmolyte. In contrast, binding of dihydrofolate is weakened in the presence of osmolytes, consistent with "water uptake". Different effects are observed depending on the identity of the osmolyte. The net uptake of water upon dihydrofolate binding was previously observed in the nonhomologous R67-encoded dihydrofolate reductase (dfrB or type II enzyme) [Chopra, S., et al. (2008) J. Biol. Chem. 283, 4690-4698]. As R67 dihydrofolate reductase possesses a nonhomologous sequence and forms a tetrameric structure with a single active site pore, the observation of weaker DHF binding in the presence of osmolytes in both enzymes implicates cosolvent effects on free dihydrofolate. Consistent with this analysis, stopped flow experiments find betaine mostly affects DHF binding via changes in k(on), while betaine mostly affects NADPH binding via changes in k(off). Finally, nonadditive enthalpy terms when binary and ternary cofactor binding events are compared suggest the presence of long-lived conformational transitions that are not included in a simple thermodynamic cycle.

Investigation of the low frequency dynamics of heme proteins: native and mutant cytochrome P450(cam) and redox partner complexes

Vibrational coherence spectroscopy (VCS) is used to investigate the low-frequency dynamics of camphor-free and camphor-bound cytochrome P450(cam) (CYP 101) and its L358P mutant. The low-frequency heme vibrations are found to be perturbed upon binding to the electron transfer partner putidaredoxin (Pdx). A strong correlation between the "detuned" vibrational coherence spectrum, which monitors frequencies between 100 and 400 cm(-1), and the lower frequency part of the Raman spectrum is also demonstrated. The very low frequency region ≤200 cm(-1), uniquely accessed by open-band VCS measurements, reveals a mode near 103 cm(-1) in P450(cam) when camphor is not present in the distal pocket. This reflects the presence of a specific heme distortion, such as saddling or ruffling, in the substrate-free state where water is coordinated to the low-spin iron atom. Such distortions are likely to retard the rate of electron transfer to the substrate-free protein. The presence of strong mode near ∼33 cm(-1) in the camphor-bound form suggests a significant heme-doming distortion, which is supported by analysis using normal coordinate structural decomposition. Pdx also displays a strong coherent vibration near 30 cm(-1) that in principle could be involved in vibrational resonance with its electron transfer target. A splitting of the 33 cm(-1) feature and intensification of a mode near 78 cm(-1) appear when the P450(cam)/Pdx complex is formed. These observations are consistent with vibrational mixing and heme geometric distortions upon Pdx binding that are coincident with the increased thiolate electron donation to the heme. The appearance of a mode near 65 cm(-1) in the coherence spectra of the L358P mutant is comparable to the mode at 78 cm(-1) seen in the P450(cam)/Pdx complex and is consistent with the view that the heme and its environment in the L358P mutant are similar to the Pdx-bound native protein. Resonance Raman spectra are presented for both P450(cam) and the L358P mutant and the changes are correlated with an increased amount of thiolate electron donation to the heme in the mutant sample.

P450cam visits an open conformation in the absence of substrate

P450cam from Pseudomonas putida is the best characterized member of the vast family of cytochrome P450s, and it has long been believed to have a more rigid and closed active site relative to other P450s. Here we report X-ray structures of P450cam crystallized in the absence of substrate and at high and low [K(+)]. The camphor-free structures are observed in a distinct open conformation characterized by a water-filled channel created by the retraction of the F and G helices, disorder of the B' helix, and loss of the K(+) binding site. Crystallization in the presence of K(+) alone does not alter the open conformation, while crystallization with camphor alone is sufficient for closure of the channel. Soaking crystals of the open conformation in excess camphor does not promote camphor binding or closure, suggesting resistance to conformational change by the crystal lattice. This open conformation is remarkably similar to that seen upon binding large tethered substrates, showing that it is not the result of a perturbation by the ligand. Redissolved crystals of the open conformation are observed as a mixture of P420 and P450 forms, which is converted to the P450 form upon addition of camphor and K(+). These data reveal that P450cam can dynamically visit an open conformation that allows access to the deeply buried active site without being induced by substrate or ligand.

Oxygen activation by cytochrome P450 monooxygenase

Unlike photosystem II (PSII) that catalyzes formation of the O-O bond, the cytochromes P450 (P450), members of a superfamily of hemoproteins, catalyze the scission of the O-O bond of dioxygen molecules and insert a single oxygen atom into unactivated hydrocarbons through a hydrogen abstraction-oxygen rebound mechanism. Hydroxylation of the unactivated hydrocarbons at physiological temperatures is vital for many cellar processes such as the biosynthesis of many endogenous compounds and the detoxification of xenobiotics in humans and plants. Even though it carries out the opposite of the water splitting reaction, P450 may share similarities to PSII in proton delivery networks, oxygen and water access channels, and consecutive electron transfer processes. In this article, we review recent advances in understanding the molecular mechanisms by which P450 activates dioxygen.

Solution NMR structure of putidaredoxin-cytochrome P450cam complex via a combined residual dipolar coupling-spin labeling approach suggests a role for Trp106 of putidaredoxin in complex formation

The 58-kDa complex formed between the [2Fe-2S] ferredoxin, putidaredoxin (Pdx), and cytochrome P450cam (CYP101) from the bacterium Pseudomonas putida has been investigated by high-resolution solution NMR spectroscopy. Pdx serves as both the physiological reductant and effector for CYP101 in the enzymatic reaction involving conversion of substrate camphor to 5-exo-hydroxycamphor. In order to obtain an experimental structure for the oxidized Pdx-CYP101 complex, a combined approach using orientational data on the two proteins derived from residual dipolar couplings and distance restraints from site-specific spin labeling of Pdx has been applied. Spectral changes for residues in and near the paramagnetic metal cluster region of Pdx in complex with CYP101 have also been mapped for the first time using (15)N and (13)C NMR spectroscopy, leading to direct identification of the residues strongly affected by CYP101 binding. The new NMR structure of the Pdx-CYP101 complex agrees well with results from previous mutagenesis and biophysical studies involving residues at the binding interface such as formation of a salt bridge between Asp38 of Pdx and Arg112 of CYP101, while at the same time identifying key features different from those of earlier modeling studies. Analysis of the binding interface of the complex reveals that the side chain of Trp106, the C-terminal residue of Pdx and critical for binding to CYP101, is located across from the heme-binding loop of CYP101 and forms non-polar contacts with several residues in the vicinity of the heme group on CYP101, pointing to a potentially important role in complex formation.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

FTIR studies of the redox partner interaction in cytochrome P450: The Pdx–P450cam couple

Recently we have developed a new approach to study protein-protein interactions using Fourier transform infrared spectroscopy in combination with titration experiments and principal component analysis (FTIR-TPCA). In the present paper we review the FTIR-TPCA results obtained for the interaction between cytochrome P450 and the redox partner protein in two P450 systems, the Pseudomonas putida P450cam (CYP101) with putidaredoxin (P450cam-Pdx), and the Bacillus megaterium P450BM-3 (CYP102) heme domain with the FMN domain (P450BMP-FMND). Both P450 systems reveal similarities in the structural changes that occur upon redox partner complex formation. These involve an increase in beta-sheets and alpha-helix content, a decrease in the population of random coil/3(10)-helix structure, a redistribution of turn structures within the interacting proteins and changes in the protonation states or hydrogen-bonding of amino acid carboxylic side chains. We discuss in detail the P450cam-Pdx interaction in comparison with literature data and conclusions drawn from experiments obtained by other spectroscopic techniques. The results are also interpreted in the context of a 3D structural model of the Pdx-P450cam complex.

Putidaredoxin-to-cytochrome P450cam electron transfer: differences between the two reductive steps required for catalysis

Cytochrome P450cam (P450cam) is the terminal monooxygenase in a three-component camphor-hydroxylating system from Pseudomonas putida. The reaction cycle requires two distinct electron transfer (ET) processes from the [2Fe-2S] containing putidaredoxin (Pdx) to P450cam. Even though the mechanism of interaction and ET between the two proteins has been under investigation for over 30 years, the second reductive step and the effector role of Pdx are not fully understood. We utilized mutagenesis, kinetic, and computer modeling approaches to better understand differences between the two Pdx-to-P450cam ET events. Our results indicate that interacting residues and the ET pathways in the complexes formed between reduced Pdx (Pdx(r)) and the ferric and ferrous dioxygen-bound forms of P450cam (oxy-P450cam) are different. Pdx Asp38 and Trp106 were found to be key players in both reductive steps. Compared to the wild-type Pdx, the D38A, W106A, and delta106 mutants exhibited considerably higher Kd values for ferric P450cam and retained ca. 20% of the first electron transferring ability. In contrast, the binding affinity of the mutants for oxy-P450cam was not substantially altered while the second ET rates were <1%. On the basis of the kinetic and modeling data we conclude that (i) P450cam-Pdx interaction is highly specific in part because it is guided/controlled by the redox state of both partners; (ii) there are alternative ET routes from Pdx(r) to ferric P450cam and a unique pathway to oxy-P450cam involving Asp38; (iii) Pdx Trp106 is a key structural element that couples the second ET event to product formation possibly via its "push" effect on the heme-binding loop.

NMR study on the structural changes of cytochrome P450cam upon the complex formation with putidaredoxin. Functional significance of the putidaredoxin-induced structural changes

We investigated putidaredoxin-induced structural changes in carbonmonoxy P450cam by using NMR spectroscopy. The resonance from the beta-proton of the axial cysteine was upfield shifted by 0.12 ppm upon the putidaredoxin binding, indicating that the axial cysteine approaches to the heme-iron by about 0.1 A. The approach of the axial cysteine to the heme-iron would enhance the electronic donation from the axial thiolate to the heme-iron, resulting in the enhanced heterolysis of the dioxygen bond. In addition to the structural perturbation on the axial ligand, the structural changes in the substrate and ligand binding site were observed. The resonances from the 5-exo- and 9-methyl-protons of d-camphor, which were newly identified in this study, were upfield shifted by 1.28 and 0.20 ppm, respectively, implying that d-camphor moves to the heme-iron by 0.15-0.7 A. Based on the radical rebound mechanism, the approach of d-camphor to the heme-iron could promote the oxygen transfer reaction. On the other hand, the downfield shift of the resonance from the gamma-methyl group of Thr-252 reflects the movement of the side chain away from the heme-iron by approximately 0.25 A. Because Thr-252 regulates the heterolysis of the dioxygen bond, the positional rearrangement of Thr-252 might assist the scission of the dioxygen bond. We, therefore, conclude that putidaredoxin induces the specific heme environmental changes of P450cam, which would facilitate the oxygen activation and the oxygen transfer reaction.

Crystal Structure of Putidaredoxin, the [2Fe–2S] Component of the P450cam Monooxygenase System from Pseudomonas putida

Stability of the [2Fe-2S]-containing putidaredoxin (Pdx), the electron donor to cytochrome P450cam in Pseudomonas putida, was improved by mutating non-ligating cysteine residues, Cys73 and Cys85, to serine singly and in combination. The increasing order of stability is Cys73Ser/Cys85Ser>Cys73Ser>Cys85Ser>WT Pdx. Crystal structures of Cys73Ser/Cys85Ser and Cys73Ser mutants of Pdx, solved by single-wavelength anomalous dispersion phasing using the [2Fe-2S] iron atoms to 1.47 A and 1.65 A resolution, respectively, are nearly identical and very similar to those of bovine adrenodoxin (Adx) and Escherichia coli ferredoxin. However, unlike the Adx structure, no motion between the core and interaction domains of Pdx is observed. This higher conformational stability of Pdx might be due to the presence of a more extensive hydrogen bonding network at the interface between the two structural domains around the conserved His49. In particular, formation of a hydrogen bond between the side-chain of Tyr51 and the carbonyl oxygen atom of Glu77 and the presence of two well-ordered water molecules linking the interaction domain and the C-terminal peptide to the core of the molecule are unique to Pdx. The folding topology of the NMR model is similar to that of the X-ray structure of Pdx. The overall rmsd of Calpha positions between the two models is 1.59 A. The largest positional differences are observed for residues 18-21 and 33-37 in the loop regions and the C terminus. The latter two peptides display conformational heterogeneity in the crystal structures. Owing to flexibility, the aromatic ring of the C-terminal Trp106 can closely approach the side-chains of Asp38 and Thr47 (3.2-3.9 A) or move away and leave the active site solvent-exposed. Therefore, Trp106, previously shown to be important in the Pdr-to-Pdx and Pdx-to-P450cam electron transfer reactions is in a position to regulate and/or mediate electron transfer to or from the [2Fe-2S] center of Pdx.

Reactivities of organic phase biosensors: 6. Square-wave and differential pulse studies of genetically engineered cytochrome P450(cam) (CYP101) bioelectrodes in selected solvents

Cytochrome P450(cam) (CYP101) bioelectrodes suitable for application in organic phases were prepared from genetically engineered CYP101 and vesicular dispersions of didodecyldimethylammonium bromide. The amperometric biosensor system was characterised under anaerobic conditions by cyclic and square-wave voltammetric methods. Cyclic- and square-wave-voltammetry studies showed that the biosensors exhibited direct reversible electron transfer between the haem iron atom and the glassy carbon electrode surface. The formal redox potential estimated for the electrode in acetonitrile was -380 mV/Ag-AgCl. The formal potential shifted anodically as the organic phase biosensor responded irreversibly to substrate (camphor) under anaerobic and aerobic conditions in acetonitrile. Differential pulse analysis of the reactivities of the CYP101 enzyme electrode confirmed the square-wave voltammetry result, which showed that the binding of substrate decreased the redox potential necessary for initiating the monooxygenation reaction of cytochrome P450(cam).

The effects of osmotic and hydrostatic pressures on macromolecular systems

Osmotic pressure and hydrostatic pressure can be used effectively to probe the behavior of biologically important macromolecules and their complexes. Using the two techniques requires a theoretical framework as well as knowledge of the more common pitfalls. Both are discussed in this review in the context of several examples.