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.

Separating Growth from Elastic Deformation during Cell Enlargement1

Plants change size by deforming reversibly (elastically) whenever turgor pressure changes, and by growing. The elastic deformation is independent of growth because it occurs in nongrowing cells. Its occurrence with growth has prevented growth from being observed alone. We investigated whether the two processes could be separated in internode cells of Chara corallina Klien ex Willd., em R.D.W. by injecting or removing cell solution with a pressure probe to change turgor while the cell length was continuously measured. Cell size changed immediately when turgor changed, and growth rates appeared to be altered. Low temperature eliminated growth but did not alter the elastic effects. This allowed elastic deformation measured at low temperature to be subtracted from elongation at warm temperature in the same cell. After the subtraction, growth alone could be observed for the first time. Alterations in turgor caused growth to change rapidly to a new, steady rate with no evidence of rapid adjustments in wall properties. This turgor response, together with the marked sensitivity of growth to temperature, suggested that the growth rate was not controlled by inert polymer extension but rather by biochemical reactions that include a turgor-sensitive step.

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.

Kinetic discrimination in T-cell activation.

We propose a quantitative model for T-cell activation in which the rate of dissociation of ligand from T-cell receptors determines the agonist and antagonist properties of the ligand. The ligands are molecular complexes between antigenic peptides and proteins of the major histocompatibility complex on the surfaces of antigen-presenting cells. Binding of ligand to receptor triggers a series of biochemical reactions in the T cell. If the ligand dissociates after these reactions are complete, the T cell receives a positive activation signal. However, dissociation of ligand after completion of the first reaction but prior to generation of the final products results in partial T-cell activation, which acts to suppress a positive response. Such a negative signal is brought about by T-cell ligands containing the variants of antigenic peptides referred to as T-cell receptor antagonists. Results of recent experiments with altered peptide ligands compare favorably with T-cell responses predicted by this model.

Biochemical evolution II: Origin of life in tubular microstructures on weathered feldspar surfaces

Mineral surfaces were important during the emergence of life on Earth because the assembly of the necessary complex biomolecules by random collisions in dilute aqueous solutions is implausible. Most silicate mineral surfaces are hydrophilic and organophobic and unsuitable for catalytic reactions, but some silica-rich surfaces of partly dealuminated feldspars and zeolites are organophilic and potentially catalytic. Weathered alkali feldspar crystals from granitic rocks at Shap, north west England, contain abundant tubular etch pits, typically 0.4–0.6 μm wide, forming an orthogonal honeycomb network in a surface zone 50 μm thick, with 2–3 × 106 intersections per mm2 of crystal surface. Surviving metamorphic rocks demonstrate that granites and acidic surface water were present on the Earth’s surface by ∼3.8 Ga. By analogy with Shap granite, honeycombed feldspar has considerable potential as a natural catalytic surface for the start of biochemical evolution. Biomolecules should have become available by catalysis of amino acids, etc. The honeycomb would have provided access to various mineral inclusions in the feldspar, particularly apatite and oxides, which contain phosphorus and transition metals necessary for energetic life. The organized environment would have protected complex molecules from dispersion into dilute solutions, from hydrolysis, and from UV radiation. Sub-micrometer tubes in the honeycomb might have acted as rudimentary cell walls for proto-organisms, which ultimately evolved a lipid lid giving further shelter from the hostile outside environment. A lid would finally have become a complete cell wall permitting detachment and flotation in primordial “soup.” Etch features on weathered alkali feldspar from Shap match the shape of overlying soil bacteria.

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.

Reconstitution of p53-ubiquitinylation reactions from purified components: the role of human ubiquitin-conjugating enzyme UBC4 and E6-associated protein (E6AP).

The E6 protein of the high-risk human papillomaviruses inactivates the tumor suppressor protein p53 by stimulating its ubiquitinylation and subsequent degradation. Ubiquitinylation is a multistep process involving a ubiquitin-activating enzyme, one of many distinct ubiquitin-conjugating enzymes, and in certain cases, a ubiquitin ligase. In human papillomavirus-infected cells, E6 and the E6-associated protein are thought to act as a ubiquitin-protein ligase in the ubiquitinylation of p53. Here we describe the cloning of a human ubiquitin-conjugating enzyme that specifically ubiquitinylates E6-associated protein. Furthermore, we define the biochemical pathway of p53 ubiquitinylation and demonstrate that in vivo inhibition of various components in the pathway leads to an inhibition of E6-stimulated p53 degradation.