Semi-continuum electrostatic calculations of redox potentials in photosystem I

The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill.     

Calculations of free-energy contributions to protein-RNA complex stabilization

The problem of calculating binding affinities of protein-RNA complexes is addressed by analyzing a computational strategy of modeling electrostatic free energies based on a nonlinear Poisson-Boltzmann (NLPB) model and linear response approximation (LRA). The underlying idea is to treat binding as a two-step process. Solutions to the NLPB equation calculate free energies arising from electronic polarizability and the LRA is constructed from molecular dynamics simulations to model reorganization free energies due to conformational transitions. By implementing a consistency condition of requiring the NLPB model to reproduce the solute-solvent free-energy transitions determined by the LRA, a "macromolecule dielectric constant" (epsilon(m)) for treating reorganization is obtained. The applicability of this hybrid approach was evaluated by calculating the absolute free energy of binding and free-energy changes for amino acid substitutions in the complex between the U1A spliceosomal protein and its cognate RNA hairpin. Depending on the residue substitution, epsilon(m) varied from 3 to 18, and reflected dipolar reorientation not included in the polarization modeled by epsilon(m) = 2. Although the changes in binding affinities from substitutions modeled strictly at the implicit level by the NLPB equation with epsilon(m) = 4 reproduced the experimental values with good overall agreement, substitutions problematic to this simple treatment showed significant improvement when solved by the NLPB-LRA approach.     

Coupling of electron and proton transfer in oxidative water cleavage in photosynthesis

This minireview addresses questions on the mechanism of oxidative water cleavage with special emphasis on the coupling of electron (ET) and proton transfer (PT) of each individual redox step of the reaction sequence and on the mode of O-O bond formation. The following topics are discussed: (1) the multiphasic kinetics of Y(Z)(ox) formation by P680(+*) originate from three different types of rate limitations: (i) nonadiabatic electron transfer for the "fast" ns reaction, (ii) local "dielectric" relaxation for the "slow" ns reaction, and (iii) "large-scale" proton shift for the micros kinetics; (2) the ET/PT-coupling mode of the individual redox transitions within the water oxidizing complex (WOC) driven by Y(Z)(ox) is assumed to depend on the redox state S(i): the oxidation steps of S(0) and S(1) comprise separate ET and PT pathways while those of S(2) and S(3) take place via proton-coupled electron transfer (PCET) analogous to Jerry Babcock's hydrogen atom abstractor model [Biochim. Biophys. Acta, 1458 (2000) 199]; (3) S(3) is postulated to be a multistate redox level of the WOC with fast dynamic equilibria of both redox isomerism and proton tautomerism. The primary event in the essential O-O bond formation is the population of a state S(3)(P) characterized by an electronic configuration and nuclear geometry that corresponds with a complexed hydrogen peroxide; (4) the peroxidic type S(3)(P) is the entatic state for formation of complexed molecular oxygen through S(3) oxidation by Y(Z)(ox); and (5) the protein matrix itself is proposed to exert catalytic activity by functioning as "PCET director". The WOC is envisaged as a supermolecule that is especially tailored for oxidative water cleavage and acts as a molecular machine.     

Asp-99 donates a hydrogen bond not to Tyr-14 but to the steroid directly in the catalytic mechanism of Delta 5-3-ketosteroid isomerase from Pseudomonas putida biotype B

Delta 5-3-ketosteroid isomerase (KSI) catalyzes the allylic isomerization of Delta 5-3-ketosteroids at a rate approaching the diffusion limit by an intramolecular transfer of a proton. Despite the extensive studies on the catalytic mechanism, it still remains controversial whether the catalytic residue Asp-99 donates a hydrogen bond to the steroid or to Tyr-14. To clarify the role of Asp-99 in the catalysis, two single mutants of D99E and D99L and three double mutants of Y14F/D99E, Y14F/D99N, and Y14F/D99L have been prepared by site-directed mutagenesis. The D99E mutant whose side chain at position 99 is longer by an additional methylene group exhibits nearly the same kcat as the wild-type while the D99L mutant exhibits ca. 125-fold lower kcat than that of the wild-type. The mutations made at positions 14 and 99 exert synergistic or partially additive effect on kcat in the double mutants, which is inconsistent with the mechanism based on the hydrogen-bonded catalytic dyad, Asp-99 COOH...Tyr-14 OH...C3-O of the steroid. The crystal structure of D99E/D38N complexed with equilenin, an intermediate analogue, at 1.9 A resolution reveals that the distance between Tyr-14 O eta and Glu-99 O epsilon is ca. 4.2 A, which is beyond the range for a hydrogen bond, and that the distance between Glu-99 O epsilon and C3-O of the steroid is maintained to be ca. 2.4 A, short enough for a hydrogen bond to be formed. Taken together, these results strongly support the idea that Asp-99 contributes to the catalysis by donating a hydrogen bond directly to the intermediate.     

Catalysis of electron transfer by selenocysteine

Selenium is an essential element that is involved in biological redox processes. The electrode potentials of the selenocysteine half-reactions RSe(*) + e(-) --> RSe-, (RSeSeR)(*)(-) + e(-) --> 2 RSe(-), and RSeSeR + 2 e(-) --> 2 RSe(-) [E degrees' (pH 7)] are +0.43, +0.18, and -0.38 V, respectively, at pH 7. The spectra of RSe(*) and (RSeSeR)(*)(-) are characterized by absorption maxima at 460 nm (epsilon = 560 M(-)(1) cm(-)(1)) and 455 nm (epsilon = 7100 M(-)(1) cm(-)(1)), respectively. The bond dissociation energy of RSe-H has been calculated, and the value of 310 kJ/mol is in agreement with literature values. In comparison with the sulfur analogue cysteine, the more facile accessibility of the radical oxidation state is striking and may have biological implications, such as in mediation of one-electron- and two-electron-transfer processes, as illustrated by catalysis by selenocysteine of the electron transfer between dithiothreitol and benzyl viologen.     

Regulation of the properties of the heme-NO complexes in nitric-oxide synthase by hydrogen bonding to the proximal cysteine

Nitric-oxide synthase (NOS) catalyzes the formation of NO and citrulline from l-arginine and oxygen. However, the NO so formed has been found to auto-inhibit the enzymatic activity significantly. We hypothesized that the NO reactivity is in part controlled by hydrogen bonding between the conserved tryptophan residue (position 409 in the neuronal isoform of NOS (nNOS)) and the cysteine residue that forms the proximal bond to the heme. By using resonance Raman spectroscopy and NO as a probe of the heme environment, we show that in the W409F and W409Y mutants of the oxygenase domain of the neuronal enzyme (nNOSox), the Fe-NO bond in the Fe3+NO complex is weaker than in the wild type enzyme, consistent with the loss of a hydrogen bond on the sulfur atom of the proximal cysteine residue. The weaker Fe-NO bond in the W409F and W409Y mutants might result in a faster rate of NO dissociation from the ferric heme in the Trp-409 mutants as compared with the wild type enzyme, which could contribute to the lower accumulation of the inhibitory NO-bound complexes observed during catalysis with the Trp-409 mutants (Adak, S., Crooks, C., Wang, Q., Crane, B. R., Tainer, J. A., Getzoff, E. D., and Stuehr, D. J. (1999) J. Biol. Chem. 274, 26907-26911). The optical and resonance Raman spectra of the Fe2+NO complexes of the Trp-409 mutants differ from those of the wild type enzyme and indicate that a significant population of a five-coordinate Fe2+NO complex is present. These data show that the hydrogen bond provided by the Trp-409 residue is necessary to maintain the thiolate coordination when NO binds to the ferrous heme. Taken together our results indicate that the heme environment on the proximal side of nNOS is critical for the formation of a stable iron-cysteine bond and for the control of the electronic properties of heme-NO complexes.     

Critical role of the neuronal nitric-oxide synthase heme proximal side residue, Arg418, in catalysis and electron transfer.

Nitric oxide (NO) is synthesized from L-Arg in the P450-type heme active site of nitric-oxide synthase (NOS). The internal axial ligand of the heme, Cys415, may hydrogen-bond to the side chain of the conserved Arg418 residue in neuronal NOS (nNOS). To understand the role of Arg418, we generated the nNOS mutants, Arg418Ala and Arg418Leu. NO formation activities with the mutants using both L-Arg and NHA as substrates were less than 0.1 nmol/min/nmol heme, in contrast to rates of 34-35 nmol/min/nmol heme with the wild-type enzyme. The heme reduction rate of the mutants was very slow, less than 10(-2) min(-1), in contrast with that (more than 10 min(-1)) of the wild type. The backbone amide group of Arg418 interacts with the Cys415 thiolate through van der Waals contact, whereas the carbonyl oxygen of Cys415 and the guanidino N(epsilon) atom of Arg418 form a tight hydrogen bond. The results suggest that Arg418 is critical in preserving the heme proximal structure and thus, is indirectly involved in both catalysis and electron transfer from the reductase domain to the heme.     

Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

Harmonic and anharmonic dynamics of Fe-CO and Fe-O(2) in heme models

We present density-functional molecular dynamics simulations of FeP(Im)(AB) heme models (AB = CO, O(2), Im = imidazole) as a way of sketching the dynamic motion of the axial ligands at room temperature. The FeP(Im)(CO) model is characterized by an essentially upright FeCO unit, undergoing small deviations with respect to its linear equilibrium structure (bending and tilting up to 10 degrees and 7 degrees, often occur). The motion of the carbon monoxide ligand is found to be quite complex and fast, its projection on the porphyrin plane sampling all the porphyrin quadrants in a short time ( approximately 0.5 ps). Simultaneously, the imidazole ligand rotates slowly around the Fe-N(epsilon) bond. In contrast to carbon monoxide, the oxygen ligand in FeP(Im)(O(2)) prefers a conformation where the projection of the O-O axis on the porphyrin plane bisects one of the porphyrin quadrants. A transition to other quadrants takes place through an O-O/Fe-N(p) overlapping conformation, within 4-6 ps. Further details of these mechanisms and their implications are discussed.     

Chirality inversion in the bilirubin molecular exciton

The bichromophoric pigment bilirubin acts as a molecular exciton in its UV-visible and circular dichroism (CD) spectroscopy. In both polar and nonpolar solvents, an optically active analog, (beta R,beta 'R)-dimethylmesobilirubin-XIII alpha (1), exhibits intense bisignate CD Cotton effects in the region of its long wavelength UV-vis absorption near 400 nm: Delta epsilon(434)(max) + 337, Delta epsilon(389)(max) - 186 (CHCl(3)), and Delta epsilon(431)(max) + 285, Delta epsilon(386)(max) - 177 (CH(3)OH). However, introduction of an amine into a CHCl(3) solution of 1 causes the Cotton effect signs to become inverted, e.g., after addition of NH(3), Delta epsilon(433)(max) - 345, Delta epsilon(389)(max) + 243, and after addition of ethylene diamine, Delta epsilon(435)(max) - 420, Delta epsilon(390)(max) + 299. The sign inversions imply inversion of molecular chirality of the bilirubin and the phenomenon appears to be general for amines, including alpha,omega-diamines. 1,8-Diaminooctane was found to be more effective than longer or shorter chain analogs in producing CD sign inversion.     

Probing electrostatic interactions in cytochrome c using site-directed chemical modification.

This communication reports the generation of an electrostatic probe using chemical modification of methionine side chains. The alkylation of methionine by iodoacetamide was achieved in a set of Saccharomyces cerevisiae iso-1-cytochrome c mutants, introducing the nontitratable, nondelocalized positive charge of a carboxyamidomethylmethionine sulfonium (CAMMS) ion at five surface and one buried site in the protein. Changes in redox potential and its variation with temperature were used to calculate microscopic effective dielectric constants operating between the probe and the heme iron. Dielectric constants (epsilon) derived from deltadeltaG values were not useful due to entropic effects, but epsilon(deltadeltaH) gave results that supported the theory. The effect on biological activity of surface derivatization was interpreted in terms of protein-protein interactions. The introduction of an electrostatic probe in cytochrome c often resulted in marked effects on activity with one of two physiological partners: cytochrome c reductase, especially if introduced at position 65, and cytochrome c oxidase, if at position 28.     

Role of residues 143 and 278 of the human nuclear Vitamin D receptor in the full-length and Delta165-215 deletion mutant

Most of the actions of 1,25-dihydroxyvitamin D(3) [1,25(OH)(2)D(3)] are mediated by binding to the Vitamin D nuclear receptor (VDR). The crystal structure of a deletion mutant (Delta165-215) of the VDR ligand-binding domain (LBD) bound to 1,25(OH)(2)D(3) indicates that amino acid residues tyrosine-143 and serine-278 form hydrogen bonding interactions with the 3-hydroxyl group of 1,25(OH)(2)D(3). Studies of VDR and three mutants (Y143F, S278A, and Y143F/S278A) did not indicate any differences in the binding affinity between the variant receptors and the wild-type receptor. This might indicate that the 3-hydroxyl group binds differently to the full-length VDR than the to deletion mutant. To further investigate, four deletion VDR mutants were constructed: VDR(Delta165-215), VDR(Delta165-215) (Y143F), VDR(Delta165-215) (S278A), VDR(Delta165-215) (Y143F/S278A). There were no significant differences in binding affinity between the wild-type receptor and the deletion mutants except for VDR(Delta165-215) (Y143F/S278A). In gene activation assays, VDR constructs with the single mutation Y143F and the double mutation Y143F/S278A, but not the single mutation S278A required higher doses of 1,25(OH)(2)D(3) for half-maximal response. This suggests that there are some minor structural and functional differences between the wild-type VDR and the Delta165-215 deletion mutant and that Y143 residue is more important for receptor function than residue S278.     

Reaction of human myoglobin and nitric oxide. Heme iron or protein sulfhydryl (s) nitrosation dependence on the absence or presence of oxygen

The amino acid sequence of human myoglobin (Mb) is similar to other mammalian Mb except for a unique cysteine residue at position 110 (Cys(110)). Anaerobic treatment of ferrous forms of wild-type human Mb, the C110A variant of human Mb or horse heart Mb, with either authentic NO or chemically derived NO in vitro yields heme-NO complexes as detected by electron paramagnetic resonance spectroscopy (EPR). By contrast, no EPR-detectable heme-NO complex was observed from the aerobic reactions of NO and either the ferric or oxy-Mb forms of wild-type human or horse heart myoglobins. Mass analyses of wild-type human Mb treated aerobically with NO indicated a mass increase of approximately 30 atomic mass units (i.e., NO/Mb = 1 mol/mol). Mass analyses of the corresponding apoprotein after heme removal showed that NO was associated with the apoprotein fraction. New electronic maxima were detected at A(333 nm) (epsilon = 3665 +/- 90 mol(-)(1) cm(-)(1); mean +/- S.D.) and A(545 nm) (epsilon = 44 +/- 3 mol(-)(1) cm(-)(1)) in solutions of S-nitrosated wild-type human Mb (similar to S-nitrosoglutathione). Importantly, the sulfhydryl S-H stretch vibration for Cys(110) measured by Fourier transform infrared (nu approximately 2552 cm(-)(1)) was absent for both holo- and apo- forms of the wild-type human protein after aerobic treatment of the protein with NO. Together, these data indicate that the reaction of wild-type human Mb and NO yields either heme-NO or a novel S-nitrosated protein dependent on the oxidation state of the heme iron and the presence or absence of dioxygen.     

Cytochrome b5 reductase: role of the si-face residues, proline 92 and tyrosine 93, in structure and catalysis

The conserved sequence motif "RxY(T)(S)xx(S)(N)" coordinates flavin binding in NADH:cytochrome b(5) reductase (cb(5)r) and other members of the flavin transhydrogenase superfamily of oxidoreductases. To investigate the roles of Y93, the third and only aromatic residue of the "RxY(T)(S)xx(S)(N)" motif, that stacks against the si-face of the flavin isoalloxazine ring, and P92, the second residue in the motif that is also in close proximity to the FAD moiety, a series of rat cb(5)r variants were produced with substitutions at either P92 or Y93, respectively. The proline mutants P92A, G, and S together with the tyrosine mutants Y93A, D, F, H, S, and W were recombinantly expressed in E. coli and purified to homogeneity. Each mutant protein was found to bind FAD in a 1:1 cofactor:protein stoichiometry while UV CD spectra suggested similar secondary structure organization among all nine variants. The tyrosine variants Y93A, D, F, H, and S exhibited varying degrees of blue-shift in the flavin visible absorption maxima while visible CD spectra of the Y93A, D, H, S, and W mutants exhibited similar blue-shifted maxima together with changes in absorption intensity. Intrinsic flavin fluorescence was quenched in the wild type, P92S and A, and Y93H and W mutants while Y93A, D, F, and S mutants exhibited increased fluorescence when compared to free FAD. The tyrosine variants Y93A, D, F, and S also exhibited greater thermolability of FAD binding. The specificity constant (k(cat)/K(m)(NADH)) for NADH:FR activity decreased in the order wild type > P92S > P92A > P92G > Y93F > Y93S > Y93A > Y93D > Y93H > Y93W with the Y93W variant retaining only 0.5% of wild-type efficiency. Both K(s)(H4NAD) and K(s)(NAD+) values suggested that Y93A, F, and W mutants had compromised NADH and NAD(+) binding. Thermodynamic measurements of the midpoint potential (E degrees ', n = 2) of the FAD/FADH(2) redox couple revealed that the potentials of the Y93A and S variants were approximately 30 mV more positive than that of wild-type cb(5)r (E degrees ' = -268 mV) while that of Y93H was approximately 30 mV more negative. These results indicate that neither P92 nor Y93 are critical for flavin incorporation in cb(5)r and that an aromatic side chain is not essential at position 93, but they demonstrate that Y93 forms contacts with the FAD that effectively modulate the spectroscopic, catalytic, and thermodynamic properties of the bound cofactor.     

The effect of mutation at valine-45 on the stability and redox potentials of trypsin-cleaved cytochrome b5

In an attempt to elucidate the determinants of redox potential and protein stability in cytochrome b5, three mutants at a highly conserved residue Val45, which is a member of heme hydrophobic pocket residues have been characterized. The V45Y mutant was designed to introduce a bulkier residue and a hydroxyl group to the heme pocket. The mutants V45H and V45E were constructed to test the effect of positive and negative charge on the stability and redox potential of proteins. The influence of these mutants on the protein stability towards thermal, urea, acid, ethanol and on the redox potential were studied. It is concluded that the decrease of hydrophobic free energy and the larger volume of the tyrosine make the phenylhydroxyl group of tyrosine still sitting inside the hydrophobic pocket, while the side chain of the mutant V45E and V45H shift away from the heme pocket. The redox potentials of mutants V45Y, V45H, V45E and wild-type of cytochrome b5 are -35 mV, 8 mV, -26 mV and -3 mV, respectively. The bigger change of the V45Y on redox potential is due to the close contact between the hydroxyl group and the heme, while the changes of the V45E and V45H result from the alteration of charge density and distribution around the heme. Different relative stability of these mutants towards heat have been observed with the order: WT > V45Y-V45H > V45E being both in the oxidized and reduced state. The relative stability induced by addition of urea decreases in the order: WT > V45Y > V45H > V45E. These results suggest that the difference in the hydrophobic free energy is a major factor contributing to the stability of the Val45 mutants. Also the loose of the helix III in the mutant V45E makes it more unstable. These results indicate that residue Val45 plays an important role in the stability and redox potential of the protein.     

Neutrophil-bead collision assay: pharmacologically induced changes in membrane mechanics regulate the PSGL-1/P-selectin adhesion lifetime

Visualization of flowing neutrophils colliding with adherent 1-mum-diameter beads presenting P-selectin allowed the simultaneous measurement of collision efficiency (epsilon), membrane tethering fraction (f), membrane tether growth dynamics, and PSGL-1/P-selectin binding lifetime. For 1391 collisions analyzed over venous wall shear rates from 25 to 200 s(-1), epsilon decreased from 0.17 to 0.004, whereas f increased from 0.15 to 0.70, and the average projected membrane tether length, L(tether)(m), increased from 0.35 mum to approximately 2.0 mum over this shear range. At all shear rates tested, adhesive collisions lacking membrane tethers had average bond lifetimes less than those observed for collisions with tethers. For adhesive collisions that failed to form membrane tethers, the regressed Bell parameters (consistent with single bond Monte Carlo simulation) were zero-stress off-rate, k(off)(0) = 0.56 s(-1) and reactive compliance, r = 0.10 nm, similar to published atomic force microscopy (AFM) measurements. For all adhesion events (+/- tethers), the bond lifetime distributions were more similar to those obtained by rolling assay and best simulated by Monte Carlo with the above Bell parameters and an average of 1.48 bonds (n = 1 bond (67%), n = 2 (22%), and n = 3-5 (11%)). For collisions at 100 s(-1), pretreatment of neutrophils with actin depolymerizing agents, latrunculin or cytochalasin D, had no effect on epsilon, but increased L(tether)(m) by 1.74- or 2.65-fold and prolonged the average tether lifetime by 1.41- or 1.65-fold, respectively. Jasplakinolide, an actin polymerizing agent known to cause blebbing, yielded results similar to the depolymerizing agents. Conversely, cholesterol-depletion with methyl-beta-cyclodextrin or formaldehyde fixation had no effect on epsilon, but reduced L(tether)(m) by 66% or 97% and reduced the average tether lifetime by 30% or 42%, respectively. The neutrophil-bead collision assay combines advantages of atomic force microscopy (small contact zone), aggregometry (discrete interactions), micropipette manipulation (tether visualization), and rolling assays (physiologic flow loading). Membrane tether growth can be enhanced or reduced pharmacologically with consequent effects on PSGL-1/P-selectin lifetimes.     

Characterization of heme-regulated eIF2alpha kinase: roles of the N-terminal domain in the oligomeric state, heme binding, catalysis, and inhibition

Heme-regulated eIF2alpha kinase [heme-regulated inhibitor (HRI)] plays a critical role in the regulation of protein synthesis by heme iron. The kinase active site is located in the C-terminal domain, whereas the N-terminal domain is suggested to regulate catalysis in response to heme binding. Here, we found that the rate of dissociation for Fe(III)-protoporphyrin IX was much higher for full-length HRI (1.5 x 10(-)(3) s(-)(1)) than for myoglobin (8.4 x 10(-)(7) s(-)(1)) or the alpha-subunit of hemoglobin (7.1 x 10(-)(6) s(-)(1)), demonstrating the heme-sensing character of HRI. Because the role of the N-terminal domain in the structure and catalysis of HRI has not been clear, we generated N-terminal truncated mutants of HRI and examined their oligomeric state, heme binding, axial ligands, substrate interactions, and inhibition by heme derivatives. Multiangle light scattering indicated that the full-length enzyme is a hexamer, whereas truncated mutants (truncations of residues 1-127 and 1-145) are mainly trimers. In addition, we found that one molecule of heme is bound to the full-length and truncated mutant proteins. Optical absorption and electron spin resonance spectra suggested that Cys and water/OH(-) are the heme axial ligands in the N-terminal domain-truncated mutant complex. We also found that HRI has a moderate affinity for heme, allowing it to sense the heme concentration in the cell. Study of the kinetics showed that the HRI kinase reaction follows classical Michaelis-Menten kinetics with respect to ATP but sigmoidal kinetics and positive cooperativity between subunits with respect to the protein substrate (eIF2alpha). Removal of the N-terminal domain decreased this cooperativity between subunits and affected the other kinetic parameters including inhibition by Fe(III)-protoporphyrin IX, Fe(II)-protoporphyrin IX, and protoporphyrin IX. Finally, we found that HRI is inhibited by bilirubin at physiological/pathological levels (IC(50) = 20 microM). The roles of the N-terminal domain and the binding of heme in the structural and functional properties of HRI are discussed.     

Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar

The contributions of electrostatic interactions to the binding stability of barnase and barstar were studied by the Poisson-Boltzmann model with three different protocols: a), the dielectric boundary specified as the van der Waals (vdW) surface of the protein along with a protein dielectric constant (epsilon (p)) of 4; b), the dielectric boundary specified as the molecular (i.e., solvent-exclusion (SE)) surface along with epsilon (p) = 4; and c), "SE + epsilon (p) = 20." The "vdW + epsilon (p) = 4" and "SE + epsilon (p) = 20" protocols predicted an overall electrostatic stabilization whereas the "SE + epsilon (p) = 4" protocol predicted an overall electrostatic destabilization. The "vdW + epsilon (p) = 4" protocol was most consistent with experiment. It quantitatively reproduced the observed effects of 17 mutations neutralizing charged residues lining the binding interface and the measured coupling energies of six charge pairs across the interface and reasonably rationalized the experimental ionic strength and pH dependences of the binding constant. In contrast, the "SE + epsilon (p) = 4" protocol predicted significantly larger coupling energies of charge pairs whereas the "SE + epsilon (p) = 20" protocol did not predict any pH dependence. This study calls for further scrutiny of the different Poisson-Boltzmann protocols and demonstrates potential danger in drawing conclusions on electrostatic contributions based on a particular calculation protocol.     

Cytochrome c reconstituted from two peptide fragments displays native-like redox properties.

Recombination of two fragments of horse cytochrome c (the heme-containing N-fragment, residues 1-56, and the C-fragment, residues 57-104), which are substantially unstructured at neutral pH, gives rise to a 1:1 fragment complex with a compact conformation, in which the alpha helical structure and the native Met80-Fe(III) axial bond are recovered. With respect to the native protein, the ferric complex shows a less rigid atomic packing and a decreased stability [Delta(DeltaG(o))D = 14.7 kJ.mol(-1)], ascribed to perturbations involving the Trp59 microenvironment and, to a lower extent, the heme pocket region. The redox potential, E1/2 = 234 +/- 5 mV vs. normal hydrogen electrode at 25 degrees C, is close to that of the intact protein, consistent with recovery of the native Met80-heme Fe(III) axial bond. Furthermore, the fragment complex shows reactivity similar to intact cytochrome c, in the reaction with cytochrome c oxidase. We conclude that the absence in the complex of some native cross-links and interlocked packing important for protein rigidity and stability is not as relevant for maintaining the native redox properties of the protein, provided that some structural requirements (i.e. recovering of the native-like alpha helical structure) are fulfilled and coordination of Met80 to the heme-iron is restored.     

Water oxidation by photosystem II: H(2)O-D(2)O exchange and the influence of pH support formation of an intermediate by removal of a proton before dioxygen creation

Understanding the chemistry of photosynthetic water oxidation requires deeper insight into the interrelation between electron transfer (ET) and proton relocations. In photosystem II membrane particles, the redox transitions of the water-oxidizing Mn complex were initiated by nanosecond laser flashes and monitored by absorption spectroscopy at 360 nm (A(360)). In the oxygen evolution transition (S(3) + hν → S(0) + O(2)), an exponential decrease in A(360) (τ(O(2)) = 1.6 ms) can be assigned to Mn reduction and O(2) formation. The corresponding rate-determining step is the ET from the Mn complex to a tyrosine radical (Y(Z)(ox)). We find that this A(360) decrease is preceded by a lag phase with a duration of 170 ± 40 μs (τ(lag) at pH 6.2), indicating formation of an intermediate before ET and O-O bond formation and corroborating results obtained by time-resolved X-ray spectroscopy. Whereas τ(O(2)) exhibits a minor kinetic isotope effect and negligible pH dependence, formation of the intermediate is slowed significantly both in D(2)O (τ(lag) increase of ～140% in D(2)O) and at low pH (τ(lag) of 30 ± 20 μs at pH 7.0 vs τ(lag) of 470 ± 80 μs at pH 5.5). These findings support the fact that in the oxygen evolution transition an intermediate is created by deprotonation and removal of a proton from the Mn complex, after Y(Z)(ox) formation but before the onset of electron transfer and O-O bond formation.