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.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Photoassembly of the manganese cluster and oxygen evolution from monomeric and dimeric CP47 reaction center photosystem II complexes

Isolated subcomplexes of photosystem II from spinach (CP47RC), composed of D1, D2, cytochrome b559, CP47, and a number of hydrophobic small subunits but devoid of CP43 and the extrinsic proteins of the oxygen-evolving complex, were shown to reconstitute the Mn4Ca1Clx cluster of the water-splitting system and to evolve oxygen. The photoactivation process in CP47RC dimers proceeds by the same two-step mechanism as observed in PSII membranes and exhibits the same stoichiometry for Mn2+, but with a 10-fold lower affinity for Ca2+ and an increased susceptibility to photodamage. After the lower Ca2+ affinity and the 10-fold smaller absorption cross-section for photons in CP47 dimers is taken into account, the intrinsic rate constant for the rate-limiting calcium-dependent dark step is indistinguishable for the two systems. The monomeric form of CP47RC also showed capacity to photoactivate and catalyze water oxidation, but with lower activity than the dimeric form and increased susceptibility to photodamage. After optimization of the various parameters affecting the photoactivation process in dimeric CP47RC subcores, 18% of the complexes were functionally reconstituted and the quantum efficiency for oxygen production by reactivated centers approached 96% of that observed for reconstituted photosystem II-enriched membranes.

Ocular cicatricial pemphigoid antigen: Partial sequence and biochemical characterization

Ocular cicatricial pemphigoid (OCP) is an autoimmune disease that affects mainly conjunctiva and other squamous epithelia. OCP is histologically characterized by a separation of the epithelium from underlying tissues within the basement membrane zone. Immunopathological studies demonstrate the deposition of anti-basement membrane zone autoantibodies in vivo. Purified IgG from sera of patients with active OCP identified a cDNA clone from a human keratinocyte cDNA library that had complete homology with the cytoplasmic domain of β4-integrin. The sera recognized a 205-kDa protein in human epidermal, human conjunctiva, and tumor cell lysates that was identified as β4-integrin by its reaction with polyclonal and monoclonal antibodies to human β4-integrin. Sera from patients with bullous pemphigoid, pemphigus vulgaris, and cicatricial pemphigoid-like diseases did not recognize the 205-kDa protein, indicating the specificity of the binding. These data strongly implicate a role for human β4-integrin in the pathogenesis of OCP. It should be emphasized that multiple antigens in the basement membrane zone of squamous epithelia may serve as targets for a wide spectrum of autoantibodies observed in vesiculobullous diseases. Molecular definition of these autoantigens will facilitate the classification and characterization of subsets of cicatricial pemphigoid and help distinguishing them from bullous pemphigoid. This study highlights the function and importance of β4-integrin in maintaining the attachment of epithelial cells to the basement membrane.

Comparative genomics of the restriction-modification systems in Helicobacter pylori

Helicobacter pylori is a Gram-negative bacterial pathogen with a small genome of 1.64–1.67 Mb. More than 20 putative DNA restriction-modification (R-M) systems, comprising more than 4% of the total genome, have been identified in the two completely sequenced H. pylori strains, 26695 and J99, based on sequence similarities. In this study, we have investigated the biochemical activities of 14 Type II R-M systems in H. pylori 26695. Less than 30% of the Type II R-M systems in 26695 are fully functional, similar to the results obtained from strain J99. Although nearly 90% of the R-M genes are shared by the two H. pylori strains, different sets of these R-M genes are functionally active in each strain. Interestingly, all strain-specific R-M genes are active, whereas most shared genes are inactive. This agrees with the notion that strain-specific genes have been acquired more recently through horizontal transfer from other bacteria and selected for function. Thus, they are less likely to be impaired by random mutations. Our results also show that H. pylori has extremely diversified R-M systems in different strains, and that the diversity may be maintained by constantly acquiring new R-M systems and by inactivating and deleting the old ones.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.

Intrinsic noise in gene regulatory networks

Cells are intrinsically noisy biochemical reactors: low reactant numbers can lead to significant statistical fluctuations in molecule numbers and reaction rates. Here we use an analytic model to investigate the emergent noise properties of genetic systems. We find for a single gene that noise is essentially determined at the translational level, and that the mean and variance of protein concentration can be independently controlled. The noise strength immediately following single gene induction is almost twice the final steady-state value. We find that fluctuations in the concentrations of a regulatory protein can propagate through a genetic cascade; translational noise control could explain the inefficient translation rates observed for genes encoding such regulatory proteins. For an autoregulatory protein, we demonstrate that negative feedback efficiently decreases system noise. The model can be used to predict the noise characteristics of networks of arbitrary connectivity. The general procedure is further illustrated for an autocatalytic protein and a bistable genetic switch. The analysis of intrinsic noise reveals biological roles of gene network structures and can lead to a deeper understanding of their evolutionary origin.

Mutant LexA proteins with specific defects in autodigestion.

In self-processing biochemical reactions, a protein or RNA molecule specifically modifies its own structure. Many such reactions are regulated in response to the needs of the cell by an interaction with another effector molecule. In the system we study here, specific cleavage of the Escherichia coli LexA repressor, LexA cleaves itself in vitro at a slow rate, but in vivo cleavage requires interaction with an activated form of RecA protein. RecA acts indirectly as a coprotease to stimulate LexA autodigestion. We describe here a new class of lexA mutants, lexA (Adg-; for autodigestion-defective) mutants, termed Adg- for brevity. Adg- mutants specifically interfered with the ability of LexA to autodigest but left intact its ability to undergo RecA-mediated cleavage. The data are consistent with a conformational model in which RecA favors a reactive conformation capable of undergoing cleavage. To our knowledge, this is the first example of a mutation in a regulated self-processing reaction that impairs the rate of self-processing without markedly affecting the stimulated reaction. Had wild-type lexA carried such a substitution, discovery of its self-processing would have been difficult; we suggest that, in other systems, a slow rate of self-processing has prevented recognition that a reaction is of this nature.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

The Chlorella hexose/H+ symporter is a useful selectable marker and biochemical reagent when expressed in Volvox.

The multicellular obligately photoautotrophic alga Volvox is composed of only two types of cells, somatic and reproductive. Therefore, Volvox provides the simplest model system for the study of multicellularity. Metabolic labeling experiments using radioactive precursors are crucial for the detection of stage- and cell-type-specific proteins, glycoproteins, lipids, and carbohydrates. However, wild-type Volvox lacks import systems for sugars or amino acids. To circumvent this problem, the hexose/H+ symporter (HUP1) gene from the unicellular alga Chlorella was placed under the control of the constitutive Volvox beta-tubulin promoter. The corresponding transgenic Volvox strain synthesized the sugar transporter in a functional state and was able to efficiently incorporate 14C from labeled glucose or glucosamine. Sensitivity toward the toxic glucose/mannose analogue 2-deoxy-glucose increased by orders of magnitude in transformants. Thus we report the successful transformation of Volvox with a gene of heterologous origin. The chimeric gene may be selected for in either a positive or a negative manner, because transformants exhibit both prolonged survival in the dark in the presence of glucose and greatly increased sensitivity to the toxic sugar 2-deoxyglucose. The former trait may make the gene useful as a dominant selectable marker for use in transformation studies, whereas the latter trait may make it useful in development of a gene-targeting system.