Uracil-DNA glycosylase–DNA substrate and product structures: Conformational strain promotes catalytic efficiency by coupled stereoelectronic effects

Enzymatic transformations of macromolecular substrates such as DNA repair enzyme/DNA transformations are commonly interpreted primarily by active-site functional-group chemistry that ignores their extensive interfaces. Yet human uracil–DNA glycosylase (UDG), an archetypical enzyme that initiates DNA base-excision repair, efficiently excises the damaged base uracil resulting from cytosine deamination even when active-site functional groups are deleted by mutagenesis. The 1.8-Å resolution substrate analogue and 2.0-Å resolution cleaved product cocrystal structures of UDG bound to double-stranded DNA suggest enzyme–DNA substrate-binding energy from the macromolecular interface is funneled into catalytic power at the active site. The architecturally stabilized closing of UDG enforces distortions of the uracil and deoxyribose in the flipped-out nucleotide substrate that are relieved by glycosylic bond cleavage in the product complex. This experimentally defined substrate stereochemistry implies the enzyme alters the orientation of three orthogonal electron orbitals to favor electron transpositions for glycosylic bond cleavage. By revealing the coupling of this anomeric effect to a delocalization of the glycosylic bond electrons into the uracil aromatic system, this structurally implicated mechanism resolves apparent paradoxes concerning the transpositions of electrons among orthogonal orbitals and the retention of catalytic efficiency despite mutational removal of active-site functional groups. These UDG/DNA structures and their implied dissociative excision chemistry suggest biology favors a chemistry for base-excision repair initiation that optimizes pathway coordination by product binding to avoid the release of cytotoxic and mutagenic intermediates. Similar excision chemistry may apply to other biological reaction pathways requiring the coordination of complex multistep chemical transformations.

The chemistry of the S-nitrosoglutathione/glutathione system

S-Nitrosothiols have generated considerable interest due to their ability to act as nitric oxide (NO) donors and due to their possible involvement in bioregulatory systems—e.g., NO transfer reactions. Elucidation of the reaction pathways involved in the modification of the thiol group by S-nitrosothiols is important for understanding the role of S-nitroso compounds in vivo. The modification of glutathione (GSH) in the presence of S-nitrosoglutathione (GSNO) was examined as a model reaction. Incubation of GSNO (1 mM) with GSH at various concentrations (1–10 mM) in phosphate buffer (pH 7.4) yielded oxidized glutathione, nitrite, nitrous oxide, and ammonia as end products. The product yields were dependent on the concentrations of GSH and oxygen. Transient signals corresponding to GSH conjugates, which increased by one mass unit when the reaction was carried out with 15N-labeled GSNO, were identified by electrospray ionization mass spectrometry. When morpholine was present in the reaction system, N-nitrosomorpholine was formed. Increasing concentrations of either phosphate or GSH led to lower yields of N-nitrosomorpholine. The inhibitory effect of phosphate may be due to reaction with the nitrosating agent, nitrous anhydride (N2O3), formed by oxidation of NO. This supports the release of NO during the reaction of GSNO with GSH. The products noted above account quantitatively for virtually all of the GSNO nitrogen consumed during the reaction, and it is now possible to construct a complete set of pathways for the complex transformations arising from GSNO + GSH.

Relative role of heme nitrosylation and β-cysteine 93 nitrosation in the transport and metabolism of nitric oxide by hemoglobin in the human circulation

To quantify the reactions of nitric oxide (NO) with hemoglobin under physiological conditions and to test models of NO transport on hemoglobin, we have developed an assay to measure NO–hemoglobin reaction products in normal volunteers, under basal conditions and during NO inhalation. NO inhalation markedly raised total nitrosylated hemoglobin levels, with a significant arterial–venous gradient, supporting a role for hemoglobin in the transport and delivery of NO. The predominant species accounting for this arterial–venous gradient is nitrosyl(heme)hemoglobin. NO breathing increases S-nitrosation of hemoglobin β-chain cysteine 93, however only to a fraction of the level of nitrosyl(heme)hemoglobin and without a detectable arterial–venous gradient. A strong correlation between methemoglobin and plasma nitrate formation was observed, suggesting that NO metabolism is a primary physiological cause of hemoglobin oxidation. Our results demonstrate that NO–heme reaction pathways predominate in vivo, NO binding to heme groups is a rapidly reversible process, and S-nitrosohemoglobin formation is probably not a primary transport mechanism for NO but may facilitate NO release from heme.

Multiple pathways for ultrafast transduction of light energy in the photosynthetic reaction center of Rhodobacter sphaeroides

A pathway of electron transfer is described that operates in the wild-type reaction center (RC) of the photosynthetic bacterium Rhodobacter sphaeroides. The pathway does not involve the excited state of the special pair dimer of bacteriochlorophylls (P*), but instead is driven by the excited state of the monomeric bacteriochlorophyll (BA*) present in the active branch of pigments along which electron transfer occurs. Pump-probe experiments were performed at 77 K on membrane-bound RCs by using different excitation wavelengths, to investigate the formation of the charge separated state P+HA−. In experiments in which P or BA was selectively excited at 880 nm or 796 nm, respectively, the formation of P+HA− was associated with similar time constants of 1.5 ps and 1.7 ps. However, the spectral changes associated with the two time constants are very different. Global analysis of the transient spectra shows that a mixture of P+BA− and P* is formed in parallel from BA* on a subpicosecond time scale. In contrast, excitation of the inactive branch monomeric bacteriochlorophyll (BB) and the high exciton component of P (P+) resulted in electron transfer only after relaxation to P*. The multiple pathways for primary electron transfer in the bacterial RC are discussed with regard to the mechanism of charge separation in the RC of photosystem II from higher plants.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Apg7p/Cvt2p Is Required for the Cytoplasm-to-Vacuole Targeting, Macroautophagy, and Peroxisome Degradation Pathways

Proper functioning of organelles necessitates efficient protein targeting to the appropriate subcellular locations. For example, degradation in the fungal vacuole relies on an array of targeting mechanisms for both resident hydrolases and their substrates. The particular processes that are used vary depending on the available nutrients. Under starvation conditions, macroautophagy is the primary method by which bulk cytosol is sequestered into autophagic vesicles (autophagosomes) destined for this organelle. Molecular genetic, morphological, and biochemical evidence indicates that macroautophagy shares much of the same cellular machinery as a biosynthetic pathway for the delivery of the vacuolar hydrolase, aminopeptidase I, via the cytoplasm-to-vacuole targeting (Cvt) pathway. The machinery required in both pathways includes a novel protein modification system involving the conjugation of two autophagy proteins, Apg12p and Apg5p. The conjugation reaction was demonstrated to be dependent on Apg7p, which shares homology with the E1 family of ubiquitin-activating enzymes. In this study, we demonstrate that Apg7p functions at the sequestration step in the formation of Cvt vesicles and autophagosomes. The subcellular localization of Apg7p fused to green fluorescent protein (GFP) indicates that a subpopulation of Apg7pGFP becomes membrane associated in an Apg12p-dependent manner. Subcellular fractionation experiments also indicate that a portion of the Apg7p pool is pelletable under starvation conditions. Finally, we demonstrate that the Pichia pastoris homologue Gsa7p that is required for peroxisome degradation is functionally similar to Apg7p, indicating that this novel conjugation system may represent a general nonclassical targeting mechanism that is conserved across species.

Determination of the binding sites of the proton transfer inhibitors Cd2+ and Zn2+ in bacterial reaction centers

The reaction center (RC) from Rhodobacter sphaeroides couples light-driven electron transfer to protonation of a bound quinone acceptor molecule, QB, within the RC. The binding of Cd2+ or Zn2+ has been previously shown to inhibit the rate of reduction and protonation of QB. We report here on the metal binding site, determined by x-ray diffraction at 2.5-Å resolution, obtained from RC crystals that were soaked in the presence of the metal. The structures were refined to R factors of 23% and 24% for the Cd2+ and Zn2+ complexes, respectively. Both metals bind to the same location, coordinating to Asp-H124, His-H126, and His-H128. The rate of electron transfer from QA− to QB was measured in the Cd2+-soaked crystal and found to be the same as in solution in the presence of Cd2+. In addition to the changes in the kinetics, a structural effect of Cd2+ on Glu-H173 was observed. This residue was well resolved in the x-ray structure—i.e., ordered—with Cd2+ bound to the RC, in contrast to its disordered state in the absence of Cd2+, which suggests that the mobility of Glu-H173 plays an important role in the rate of reduction of QB. The position of the Cd2+ and Zn2+ localizes the proton entry into the RC near Asp-H124, His-H126, and His-H128. Based on the location of the metal, likely pathways of proton transfer from the aqueous surface to QB⨪ are proposed.

Identification of the proton pathway in bacterial reaction centers: Replacement of Asp-M17 and Asp-L210 with Asn reduces the proton transfer rate in the presence of Cd2+

The reaction center (RC) from Rhodobacter sphaeroides converts light into chemical energy through the reduction and protonation of a bound quinone molecule QB (the secondary quinone electron acceptor). We investigated the proton transfer pathway by measuring the proton-coupled electron transfer, kAB(2) [QA⨪QB⨪ + H+ → QA(QBH)−] in native and mutant RCs in the absence and presence of Cd2+. Previous work has shown that the binding of Cd2+ decreases kAB(2) in native RCs ≈100-fold. The preceding paper shows that bound Cd2+ binds to Asp-H124, His-H126, and His-H128. This region represents the entry point for protons. In this work we investigated the proton transfer pathway connecting the entry point with QB⨪ by searching for mutations that greatly affect kAB(2) (≳10-fold) in the presence of Cd2+, where kAB(2) is limited by the proton transfer rate (kH). Upon mutation of Asp-L210 or Asp-M17 to Asn, kH decreased from ≈60 s−1 to ≈7 s−1, which shows the important role that Asp-L210 and Asp-M17 play in the proton transfer chain. By comparing the rate of proton transfer in the mutants (kH ≈ 7 s−1) with that in native RCs in the absence of Cd2+ (kH ≥ 104 s−1), we conclude that alternate proton transfer pathways, which have been postulated, are at least 103-fold less effective.

Staphylococcal protein A: Unfolding pathways, unfolded states, and differences between the B and E domains

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein—contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the Cα root-mean-square deviation from the starting structure is less than 1.4 Å. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

Phytochrome B binds with greater apparent affinity than phytochrome A to the basic helix–loop–helix factor PIF3 in a reaction requiring the PAS domain of PIF3

The signaling pathways by which the phytochrome (phy) family of photoreceptors transmits sensory information to light-regulated genes remain to be fully defined. Evidence for a relatively direct pathway has been provided by the binding of one member of the family, phyB, to a promoter-element-bound, basic helix–loop–helix protein, PIF3, specifically upon light-induced conversion of the photoreceptor molecule to its biologically active conformer (Pfr). Here, we show that phyA also binds selectively and reversibly to PIF3 upon photoconversion to Pfr, but that the apparent affinity of PIF3 for phyA is 10-fold lower than for phyB. This result is consistent with previous in vivo data from PIF3-deficient Arabidopsis, indicating that PIF3 has a major role in phyB signaling, but a more minor role in phyA signaling. We also show that phyB binds stoichiometrically to PIF3 at an equimolar ratio, suggesting that the resultant complex is the unit active in transcriptional regulation at target promoters. Deletion mapping suggests that a 37-aa segment present at the N terminus of phyB, but absent from phyA, contributes strongly to the high binding affinity of phyB for PIF3. Conversely, deletion mapping and point mutation analysis of PIF3 for determinants involved in recognition of phyB indicates that the PAS domain of PIF3 is a major contributor to this interaction, but that a second determinant in the C-terminal domain is also necessary.

Directed evolution of a (βα)8-barrel enzyme to catalyze related reactions in two different metabolic pathways

Enzymes participating in different metabolic pathways often have similar catalytic mechanisms and structures, suggesting their evolution from a common ancestral precursor enzyme. We sought to create a precursor-like enzyme for N′-[(5′-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (ProFAR) isomerase (HisA; EC 5.3.1.16) and phosphoribosylanthranilate (PRA) isomerase (TrpF; EC 5.3.1.24), which catalyze similar reactions in the biosynthesis of the amino acids histidine and tryptophan and have a similar (βα)8-barrel structure. Using random mutagenesis and selection, we generated several HisA variants that catalyze the TrpF reaction both in vivo and in vitro, and one of these variants retained significant HisA activity. A more detailed analysis revealed that a single amino acid exchange could establish TrpF activity on the HisA scaffold. These findings suggest that HisA and TrpF may have evolved from an ancestral enzyme of broader substrate specificity and underscore that (βα)8-barrel enzymes are very suitable for the design of new catalytic activities.