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.

Adjustment of conformational flexibility is a key event in the thermal adaptation of proteins

3-Isopropylmalate dehydrogenase (IPMDH, E.C. 1.1.1.85) from the thermophilic bacterium Thermus thermophilus HB8 is homologous to IPMDH from the mesophilic Escherichia coli, but has an approximately 17°C higher melting temperature. Its temperature optimum is 22–25°C higher than that of the E. coli enzyme; however, it is hardly active at room temperature. The increased conformational rigidity required to stabilize the thermophilic enzyme against heat denaturation might explain its different temperature-activity profile. Hydrogen/deuterium exchange studies were performed on this thermophilic-mesophilic enzyme pair to compare their conformational flexibilities. It was found that Th. thermophilus IPMDH is significantly more rigid at room temperature than E. coli IPMDH, whereas the enzymes have nearly identical flexibilities under their respective optimal working conditions, suggesting that evolutionary adaptation tends to maintain a “corresponding state” regarding conformational flexibility. These observations confirm that conformational fluctuations necessary for catalytic function are restricted at room temperature in the thermophilic enzyme, suggesting a close relationship between conformational flexibility and enzyme function.

Purification of ribonucleotide reductase subunits Y1, Y2, Y3, and Y4 from yeast: Y4 plays a key role in diiron cluster assembly

Ribonucleotide reductases (RNRs) catalyze the conversion of nucleotides to deoxynucleotides. Class I RNRs are composed of two types of subunits: RNR1 contains the active site for reduction and the binding sites for the nucleotide allosteric effectors. RNR2 contains the diiron-tyrosyl radical (Y⋅) cofactor essential for the reduction process. Studies in yeast have recently identified four RNR subunits: Y1 and Y3, Y2 and Y4. These proteins have been expressed in Saccharomyces cerevisiae and in Escherichia coli and purified to ≈90% homogeneity. The specific activity of Y1 isolated from yeast and E. coli is 0.03 μmol⋅min−1⋅mg−1 and of (His)6-Y2 [(His)6-Y2-K387N] from yeast is 0.037 μmol⋅min−1⋅mg−1 (0.125 μmol⋅min−1⋅mg−1). Y2, Y3, and Y4 isolated from E. coli have no measurable activity. Efforts to generate Y⋅ in Y2 or Y4 using Fe2+, O2, and reductant have been unsuccessful. However, preliminary studies show that incubation of Y4 and Fe2+ with inactive E. coli Y2 followed by addition of O2 generates Y2 with a specific activity of 0.069 μmol⋅min−1⋅mg−1 and a Y⋅. A similar experiment with (His)6-Y2-K387N, Y4, O2, and Fe2+ results in an increase in its specific activity to 0.30 μmol⋅min−1⋅mg−1. Studies with antibodies to Y4 and Y2 reveal that they can form a complex in vivo. Y4 appears to play an important role in diiron-Y⋅ assembly of Y2.

Identification of a gene encoding an acyl CoA:diacylglycerol acyltransferase, a key enzyme in triacylglycerol synthesis

Triacylglycerols are quantitatively the most important storage form of energy for eukaryotic cells. Acyl CoA:diacylglycerol acyltransferase (DGAT, EC 2.3.1.20) catalyzes the terminal and only committed step in triacylglycerol synthesis, by using diacylglycerol and fatty acyl CoA as substrates. DGAT plays a fundamental role in the metabolism of cellular diacylglycerol and is important in higher eukaryotes for physiologic processes involving triacylglycerol metabolism such as intestinal fat absorption, lipoprotein assembly, adipose tissue formation, and lactation. DGAT is an integral membrane protein that has never been purified to homogeneity, nor has its gene been cloned. We identified an expressed sequence tag clone that shared regions of similarity with acyl CoA:cholesterol acyltransferase, an enzyme that also uses fatty acyl CoA as a substrate. Expression of a mouse cDNA for this expressed sequence tag in insect cells resulted in high levels of DGAT activity in cell membranes. No other acyltransferase activity was detected when a variety of substrates, including cholesterol, were used as acyl acceptors. The gene was expressed in all tissues examined; during differentiation of NIH 3T3-L1 cells into adipocytes, its expression increased markedly in parallel with increases in DGAT activity. The identification of this cDNA encoding a DGAT will greatly facilitate studies of cellular glycerolipid metabolism and its regulation.

Nitric oxide plays a key role in adaptive control of locomotion in cat

Diverse roles in cellular functions have been ascribed to nitric oxide (NO), and its involvement in induction of long-term depression in cerebellar Purkinje cells has been demonstrated. Manipulations of NO concentration or its synthesis in cerebellar tissues therefore provide a means for investigating roles of NO in cerebellar functions at both cellular and behavioral levels. We tested adaptive control of locomotion to perturbation in cats, and found that this form of motor learning was abolished by application of either an inhibitor of NO synthase or a scavenger of NO to the cerebellar cortical locomotion area. This finding supports the view that NO in the cerebellum plays a key role in motor learning.

Reciprocal interactions between β1-integrin and epidermal growth factor receptor in three-dimensional basement membrane breast cultures: A different perspective in epithelial biology

Anchorage and growth factor independence are cardinal features of the transformed phenotype. Although it is logical that the two pathways must be coregulated in normal tissues to maintain homeostasis, this has not been demonstrated directly. We showed previously that down-modulation of β1-integrin signaling reverted the malignant behavior of a human breast tumor cell line (T4–2) derived from phenotypically normal cells (HMT-3522) and led to growth arrest in a three-dimensional (3D) basement membrane assay in which the cells formed tissue-like acini (14). Here, we show that there is a bidirectional cross-modulation of β1-integrin and epidermal growth factor receptor (EGFR) signaling via the mitogen-activated protein kinase (MAPK) pathway. The reciprocal modulation does not occur in monolayer (2D) cultures. Antibody-mediated inhibition of either of these receptors in the tumor cells, or inhibition of MAPK kinase, induced a concomitant down-regulation of both receptors, followed by growth-arrest and restoration of normal breast tissue morphogenesis. Cross-modulation and tissue morphogenesis were associated with attenuation of EGF-induced transient MAPK activation. To specifically test EGFR and β1-integrin interdependency, EGFR was overexpressed in nonmalignant cells, leading to disruption of morphogenesis and a compensatory up-regulation of β1-integrin expression, again only in 3D. Our results indicate that when breast cells are spatially organized as a result of contact with basement membrane, the signaling pathways become coupled and bidirectional. They further explain why breast cells fail to differentiate in monolayer cultures in which these events are mostly uncoupled. Moreover, in a subset of tumor cells in which these pathways are misregulated but functional, the cells could be “normalized” by manipulating either pathway.

Cloning of chlorophyllase, the key enzyme in chlorophyll degradation: Finding of a lipase motif and the induction by methyl jasmonate

Chlorophyllase (Chlase) is the first enzyme involved in chlorophyll (Chl) degradation and catalyzes the hydrolysis of ester bond to yield chlorophyllide and phytol. In the present study, we isolated the Chlase cDNA. We synthesized degenerate oligo DNA probes based on the internal amino acid sequences of purified Chlase from Chenopodium album, screened the C. album cDNA library, and cloned a cDNA (CaCLH, C. album chlorophyll-chlorophyllido hydrolase). The deduced amino acid sequence (347 aa residues) had a lipase motif overlapping with an ATP/GTP-binding motif (P-loop). CaCLH possibly was localized in the extraplastidic part of the cell, because a putative signal sequence for endoplasmic reticulum is at the N terminus. The amino acid sequence shared 37% identity with a function-unknown gene whose mRNA is inducible by coronatine and methyl jasmonate (MeJA) in Arabidopsis thaliana (AtCLH1). We expressed the gene products of AtCLH1 and of CaCLH in Escherichia coli, and they similarly exhibited Chlase activity. Moreover, we isolated another full-length cDNA based on an Arabidopsis genomic fragment and expressed it in E. coli, demonstrating the presence of the second Arabidopsis CLH gene (AtCLH2). No typical feature of signal sequence was identified in AtCLH1, whereas AtCLH2 had a typical signal sequence for chloroplast. AtCLH1 mRNA was induced rapidly by a treatment of MeJA, which is known to promote senescence and Chl degradation in plants, and a high mRNA level was maintained up to 9 h. AtCLH2, however, did not respond to MeJA.

A targeted mutation in the murine gene encoding the high density lipoprotein (HDL) receptor scavenger receptor class B type I reveals its key role in HDL metabolism

Plasma high density lipoprotein (HDL), which protects against atherosclerosis, is thought to remove cholesterol from peripheral tissues and to deliver cholesteryl esters via a selective uptake pathway to the liver (reverse cholesterol transport) and steroidogenic tissues (e.g., adrenal gland for storage and hormone synthesis). Despite its physiologic and pathophysiologic importance, the cellular metabolism of HDL has not been well defined. The class B, type I scavenger receptor (SR-BI) has been proposed to play an important role in HDL metabolism because (i) it is a cell surface HDL receptor which mediates selective cholesterol uptake in cultured cells, (ii) its physiologically regulated expression is most abundant in the liver and steroidogenic tissues, and (iii) hepatic overexpression dramatically lowers plasma HDL. To test directly the normal role of SR-BI in HDL metabolism, we generated mice with a targeted null mutation in the SR-BI gene. In heterozygous and homozygous mutants relative to wild-type controls, plasma cholesterol concentrations were increased by ≈31% and 125%, respectively, because of the formation of large, apolipoprotein A-I (apoA-I)-containing particles, and adrenal gland cholesterol content decreased by 42% and 72%, respectively. The plasma concentration of apoA-I, the major protein in HDL, was unchanged in the mutants. This, in conjunction with the increased lipoprotein size, suggests that the increased plasma cholesterol in the mutants was due to decreased selective cholesterol uptake. These results provide strong support for the proposal that in mice the gene encoding SR-BI plays a key role in determining the levels of plasma lipoprotein cholesterol (primarily HDL) and the accumulation of cholesterol stores in the adrenal gland. If it has a similar role in controlling plasma HDL in humans, SR-BI may influence the development and progression of atherosclerosis and may be an attractive candidate for therapeutic intervention in this disease.

Signaling through fibroblast growth factor receptor 2b plays a key role in the development of the exocrine pancreas

The development of the pancreas depends on epithelial-mesenchymal interactions. Fibroblast growth factors (FGFs) and their receptors (FGFRs 1–4) have been identified as mediators of epithelial-mesenchymal interactions in different organs. We show here that FGFR-2 IIIb and its ligands FGF-1, FGF-7, and FGF-10 are expressed throughout pancreatic development. We also show that in mesenchyme-free cultures of embryonic pancreatic epithelium FGF-1, FGF-7, and FGF-10 stimulate the growth, morphogenesis, and cytodifferentiation of the exocrine cells of the pancreas. The role of FGFs signaling through FGFR-2 IIIb was further investigated by inhibiting FGFR-2 IIIb signaling in organocultures of pancreatic explants (epithelium + mesenchyme) by using either antisense FGFR-2 IIIb oligonucleotides or a soluble recombinant FGFR-2 IIIb protein. Abrogation of FGFR-2 IIIb signaling resulted in a considerable reduction in the size of the explants and in a 2-fold reduction of the development of the exocrine cells. These results demonstrate that FGFs signaling through FGFR-2 IIIb play an important role in the development of the exocrine pancreas.

Pax8 has a key role in thyroid cell differentiation

Transformation of rat thyroid cells with polyoma virus middle T antigen results in loss of the thyroid-differentiated phenotype, measured as the expression of the thyroglobulin (Tg), thyroperoxidase (TPO), and sodium/iodide symporter (NIS) genes. Among the transcription factors involved in the regulation of these genes, TTF-1 and TTF-2 were still detected at nearly wild-type levels, while a specific loss of the paired domain transcription factor Pax8 was observed. In this study, we used the PCPy cell line as a model system to study the role of Pax8 in thyroid differentiation. We demonstrate that the reintroduction of Pax8 in PCPy cells is sufficient to activate expression of the endogenous genes encoding thyroglobulin, thyroperoxidase, and sodium/iodide symporter. Thus, this cell system provides direct evidence for the ability of Pax8 to activate transcription of thyroid-specific genes at their chromosomal locus and strongly suggests a fundamental role of this transcription factor in the maintenance of functional differentiation in thyroid cells. Moreover, we show that Pax8 and TTF-1 cooperate in the activation of the thyroglobulin promoter and that additional thyroid-specific mechanism(s) are involved in such a cooperation. To identify the Pax8 domain able to mediate the specific activation of the thyroglobulin promoter, we transfected in PCPy cells three different Pax8 isoforms. The results of such experiments indicate that for the transcriptional activation of thyroid-specific genes, Pax8 uses an as yet unidentified functional domain.

The hypocone as a key innovation in mammalian evolution.

The hypocone, a cusp added to the primitively triangular upper molar teeth of therian mammals, has evolved convergently > 20 times among mammals during the Cenozoic. Acquisition of the hypocone itself involves little phenotypic change, but subsequent diversification of groups possessing the hypocone may be greatly enhanced. Our analysis of the Cenozoic mammalian radiations, including the Recent fauna, shows that high species diversity of mammals with hypocones and association of the hypocone with herbivory strongly support recognition of the hypocone as a key innovation that has allowed invasion of, and diversification within, herbivorous adaptive zones. In contrast, mammals lacking hypocones show no marked increase in species diversity during the Cenozoic.

Conservation of fibroblast growth factor function in lens regeneration

In urodele amphibians, lens induction during development and regeneration occurs through different pathways. During development, the lens is induced from the mutual interaction of the ectoderm and the optic vesicle, whereas after lentectomy the lens is regenerated through the transdifferentiation of the iris-pigmented epithelial cells. Given the known role of fibroblast growth factors (FGFs) during lens development, we examined whether or not the expression and the effects of exogenous FGF during urodele lens regeneration were conserved. In this paper, we describe expression of FGF-1 and its receptors, FGFR-2 (KGFR and bek variants) and FGFR-3, in newts during lens regeneration. Expression of these genes was readily observed in the dedifferentiating pigmented epithelial cells, and the levels of expression were high in the lens epithelium and the differentiating fibers and lower in the retina. These patterns of expression implied involvement of FGFs in lens regeneration. To further elucidate this function, we examined the effects of exogenous FGF-1 and FGF-4 during lens regeneration. FGF-1 or FGF-4 treatment in lentectomized eyes resulted in the induction of abnormalities reminiscent to the ones induced during lens development in transgenic mice. Effects included transformation of epithelial cells to fiber cells, double lens regeneration, and lenses with abnormal polarity. These results establish that FGF molecules are key factors in fiber differentiation, polarity, and morphogenesis of the lens during regeneration even though the regenerating lens is induced by a different mechanism than in lens development. In this sense, FGF function in lens regeneration and development should be regarded as conserved. Such conservation should help elucidate the mechanisms of lens regeneration in urodeles and its absence in higher vertebrates.

Plant biology in the future

In the beginning of modern plant biology, plant biologists followed a simple model for their science. This model included important branches of plant biology known then. Of course, plants had to be identified and classified first. Thus, there was much work on taxonomy, genetics, and physiology. Ecology and evolution were approached implicitly, rather than explicitly, through paleobotany, taxonomy, morphology, and historical geography. However, the burgeoning explosion of knowledge and great advances in molecular biology, e.g., to the extent that genes for specific traits can be added (or deleted) at will, have created a revolution in the study of plants. Genomics in agriculture has made it possible to address many important issues in crop production by the identification and manipulation of genes in crop plants. The current model of plant study differs from the previous one in that it places greater emphasis on developmental controls and on evolution by differential fitness. In a rapidly changing environment, the current model also explicitly considers the phenotypic variation among individuals on which selection operates. These are calls for the unity of science. In fact, the proponents of “Complexity Theory” think there are common algorithms describing all levels of organization, from atoms all the way to the structure of the universe, and that when these are discovered, the issue of scaling will be greatly simplified! Plant biology must seriously contribute to, among other things, meeting the nutritional needs of the human population. This challenge constitutes a key part of the backdrop against which future evolution will occur. Genetic engineering technologies are and will continue to be an important component of agriculture; however, we must consider the evolutionary implications of these new technologies. Meeting these demands requires drastic changes in the undergraduate curriculum. Students of biology should be trained in molecular, cellular, organismal, and ecosystem biology, including all living organisms.