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.

Heterogeneity in molecular recognition by restriction endonucleases: osmotic and hydrostatic pressure effects on BamHI, Pvu II, and EcoRV specificity.

The cleavage specificity of the Pvu II and BamHI restriction endonucleases is found to be dramatically reduced at elevated osmotic pressure. Relaxation in specificity of these otherwise highly accurate and specific enzymes, previously termed "star activity," is uniquely correlated with osmotic pressure between 0 and 100 atmospheres. No other colligative solvent property exhibits a uniform correlation with star activity for all of the compounds tested. Application of hydrostatic pressure counteracts the effects of osmotic pressure and restores the natural selectivity of the enzymes for their canonical recognition sequences. These results indicate that water solvation plays an important role in the site-specific recognition of DNA by many restriction enzymes. Osmotic pressure did not induce an analogous effect on the specificity of the EcoRV endonuclease, implying that selective hydration effects do not participate in DNA recognition in this system. Hydrostatic pressure was found to have little effect on the star activity induced by changes in ionic strength, pH, or divalent cation, suggesting that distinct mechanisms may exist for these observed alterations in specificity. Recent evidence has indicated that BamHI and EcoRI share similar structural motifs, while Pvu II and EcoRV belong to a different structural family. Evidently, the use of hydration water to assist in site-specific recognition is a motif neither limited to nor defined by structural families.

Adeno-associated virus (AAV) site-specific integration: Formation of AAV–AAVS1 junctions in an in vitro system

An in vitro system to study the mechanism of site-specific integration of adeno-associated virus (AAV) was developed. This system is based on two substrates, a linear or circular AAV donor and a circular acceptor containing the preintegration locus AAVS1. In the presence of HeLa extract and the His-Tag-purified Rep68 protein, specific covalent junctions between AAV and AAVS1 were formed and detected by PCR. The majority of the junctions were located within the Rep binding site of both the AAV and the AAVS1 substrates, underlining the involvement of the Rep protein. A limited amount of replication and the presence of nuclear factors promoted the efficiency of the reaction. The process was ATP-dependent, indicating that the helicase activity of Rep may be important in the formation of the junctions. According to current models of integration, the formation of the junctions would represent a first step in the process of AAV integration. This step could be crucial for the site specificity of the recombination event that leads to the integration of AAV into human chromosome 19 in vivo.

An in Vitro System from Maize Seedlings for Tryptophan-Independent Indole-3-Acetic Acid Biosynthesis1

The enzymatic synthesis of indole-3-acetic acid (IAA) from indole by an in vitro preparation from maize (Zea mays L.) that does not use tryptophan (Trp) as an intermediate is described. Light-grown seedlings of normal maize and the maize mutant orange pericarp were shown to contain the necessary enzymes to convert [14C]indole to IAA. The reaction was not inhibited by unlabeled Trp and neither [14C]Trp nor [14C]serine substituted for [14C]indole in this in vitro system. The reaction had a pH optimum greater than 8.0, required a reducing environment, and had an oxidation potential near that of ascorbate. The results obtained with this in vitro enzyme preparation provide strong, additional evidence for the presence of a Trp-independent IAA biosynthesis pathway in plants.

Unexpected signals in a system subject to kinetic proofreading

When multivalent ligands attach to IgEs bound to the receptors with high affinity for IgE on mast cells, the receptors aggregate, tyrosines on the receptors become phosphorylated, and a variety of cellular responses are stimulated. Prior studies, confirmed here, demonstrated that the efficiency with which later events are generated from earlier ones is inversely related to the dissociation rate of the aggregating ligand. This finding suggests that the cellular responses are constrained by a “kinetic proofreading” regimen. We have now observed an apparent exception to this rule. Doses of the rapidly or slowly dissociating ligands that generated equivalent levels of tyrosine-phosphorylated receptors comparably stimulated a putatively distal event: transcription of the gene for monocyte chemoattractant protein 1. Possible explanations of this apparent anomaly were explored.

Fibroblast growth factor (FGF) homologous factors: new members of the FGF family implicated in nervous system development.

Four new members of the fibroblast growth factor (FGF) family, referred to as fibroblast growth factor homologous factors (FHFs), have been identified by a combination of random cDNA sequencing, data base searches, and degenerate PCR. Pairwise comparisons between the four FHFs show between 58% and 71% amino acid sequence identity, but each FHF shows less than 30% identity when compared with other FGFs. Like FGF-1 (acidic FGF) and FGF-2 (basic FGF), the FHFs lack a classical signal sequence and contain clusters of basic residues that can act as nuclear localization signals. In transiently transfected 293 cells FHF-1 accumulates in the nucleus and is not secreted. Each FHF is expressed in the developing and adult nervous systems, suggesting a role for this branch of the FGF family in nervous system development and function.

Functional complementation of xeroderma pigmentosum complementation group E by replication protein A in an in vitro system.

Xeroderma pigmentosum (XP) is caused by a defect in nucleotide excision repair. Patients in the complementation group E (XP-E) have the mildest form of the disease and the highest level of residual repair activity. About 20% of the cell strains derived from XP-E patients lack a damaged DNA-binding protein (DDB) activity that binds to ultraviolet-induced (6-4) photoproducts with high affinity. We report here that cell-free extracts prepared from XP-E cell strains that either lacked or contained DDB activity were severely defective in excising DNA damage including (6-4) photoproducts. However, this excision activity defect was not restored by addition of purified DDB that, in fact, inhibited removal of (6-4) photoproducts by the human excision nuclease reconstituted from purified proteins. Extensive purification of correcting activity from HeLa cells revealed that the correcting activity is inseparable from the human replication/repair protein A [RPA (also known as human single stranded DNA binding protein, HSSB)]. Indeed, supplementing XP-E extracts with recombinant human RPA purified from Escherichia coli restored excision activity. However, no mutation was found in the genes encoding the three subunits of RPA in an XP-E (DDB-) cell line. It is concluded that RPA functionally complements XP-E extracts in vitro, but it is not genetically altered in XP-E patients.

Emergence of order in visual system development.

Neural connections in the adult central nervous system are highly precise. In the visual system, retinal ganglion cells send their axons to target neurons in the lateral geniculate nucleus (LGN) in such a way that axons originating from the two eyes terminate in adjacent but nonoverlapping eye-specific layers. During development, however, inputs from the two eyes are intermixed, and the adult pattern emerges gradually as axons from the two eyes sort out to form the layers. Experiments indicate that the sorting-out process, even though it occurs in utero in higher mammals and always before vision, requires retinal ganglion cell signaling; blocking retinal ganglion cell action potentials with tetrodotoxin prevents the formation of the layers. These action potentials are endogenously generated by the ganglion cells, which fire spontaneously and synchronously with each other, generating "waves" of activity that travel across the retina. Calcium imaging of the retina shows that the ganglion cells undergo correlated calcium bursting to generate the waves and that amacrine cells also participate in the correlated activity patterns. Physiological recordings from LGN neurons in vitro indicate that the quasiperiodic activity generated by the retinal ganglion cells is transmitted across the synapse between ganglion cells to drive target LGN neurons. These observations suggest that (i) a neural circuit within the immature retina is responsible for generating specific spatiotemporal patterns of neural activity; (ii) spontaneous activity generated in the retina is propagated across central synapses; and (iii) even before the photoreceptors are present, nerve cell function is essential for correct wiring of the visual system during early development. Since spontaneously generated activity is known to be present elsewhere in the developing CNS, this process of activity-dependent wiring could be used throughout the nervous system to help refine early sets of neural connections into their highly precise adult patterns.

A sampling problem in molecular dynamics simulations of macromolecules.

Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.

Dominant transformation by mutated human ras genes in vitro requires more than 100 times higher expression than is observed in cancers

The gene-mutation-cancer hypothesis holds that mutated cellular protooncogenes, such as point-mutated proto-ras, “play a dominant part in cancer,” because they are sufficient to transform transfected mouse cell lines in vitro [Alberts, B., Bray, D., Lewis, J., Raff, M., Roberts, K. & Watson, J. D. (1994) Molecular Biology of the Cell (Garland, New York)]. However, in cells transformed in vitro mutated human ras genes are expressed more than 100-fold than in the cancers from which they are isolated. In view of the discrepancy between the very low levels of ras transcription in cancers and the very high levels in cells transformed in vitro, we have investigated the minimal level of human ras expression for transformation in vitro. Using point-mutated human ras genes recombined with different promoters from either human metallothionein-IIA or human fibronectin or from retroviruses we found dominant in vitro transformation of the mouse C3H cell line only with ras genes linked to viral promoters. These ras genes were expressed more than 120-fold higher than are native ras genes of C3H cells. The copy number of transfected ras genes ranged from 2–6 in our system. In addition, nondominant transformation was observed in a small percentage (2–7%) of C3H cells transfected with ras genes that are expressed less than 20 times higher than native C3H ras genes. Because over 90% of cells expressing ras at this moderately enhanced level were untransformed, transformation must follow either a nondominant ras mechanism or a non-ras mechanism. We conclude that the mutated, but normally expressed, ras genes found in human and animal cancers are not likely to “play a dominant part in cancer.” The conclusion that mutated ras genes are not sufficient or dominant for cancer is directly supported by recent discoveries of mutated ras in normal animals, and in benign human tissue, “which has little potential to progress” [Jen, J., Powell, S. M., Papadopoulos, N., Smith, K. J., Hamilton, S. R., Vogelstein, B. & Kinzler, K. W. (1994) Cancer Res. 54, 5523–5526]. Even the view that mutated ras is necessary for cancer is hard to reconcile with (i) otherwise indistinguishable cancers with and without ras mutations, (ii) metastases of the same human cancers with and without ras mutations, (iii) retroviral ras genes that are oncogenic without point mutations, and (iv) human tumor cells having spontaneously lost ras mutation but not tumorigencity.

Development of an in vitro mRNA decay system for Escherichia coli: Poly(A) polymerase I is necessary to trigger degradation

Using a novel Escherichia coli in vitro decay system in which polysomes are the source of both enzymes and mRNA, we demonstrate a requirement for poly(A) polymerase I (PAP I) in mRNA turnover. The in vitro decay of two different mRNAs (trxA and lpp) is triggered by the addition of ATP only when polysomes are prepared from a strain carrying the wild-type gene for PAP I (pcnB+). The relative decay rates of these two messages are similar in vitro and in vivo. Poly(A) tails are formed on both mRNAs, but no poly(A) tails are detected on the 3′ end of mature 23S rRNA. The size distribution of poly(A) tails generated in vitro, averaging 50 nt in length, is comparable to that previously reported in vivo. PAP I activity is associated exclusively with the polysomes. Exogenously added PAP I does not restore mRNA decay to PAP I− polysomes, suggesting that, in vivo, PAP I may be part of a multiprotein complex. The potential of this in vitro system for analyzing mRNA decay in E. coli is discussed.

Analysis of a peptide hormone–receptor interaction in the yeast two-hybrid system

Interaction between a peptide hormone and extracellular domains of its receptor is a crucial step for initiation of hormone action. We have developed a modification of the yeast two-hybrid system to study this interaction and have used it to characterize the interaction of insulin-like growth factor 1 (IGF-1) with its receptor by using GAL4 transcriptional regulation with a β-galactosidase assay as readout. In this system, IGF-1 and proIGF-1 bound to the cysteine-rich domain, extracellular domain, or entire IGF-1 proreceptor. This interaction was specific. Thus, proinsulin showed no significant interaction with the IGF-1 receptor, while a chimeric proinsulin containing the C-peptide of IGF-1 had an intermediate interaction, consistent with its affinity for the IGF-1 receptor. Over 2000 IGF-1 mutants were generated by PCR and screened for interaction with the color assay. About 40% showed a strong interaction, 20% showed an intermediate interaction, and 40% give little or no signal. Of 50 mutants that were sequenced, several (Leu-5 → His, Glu-9 → Val, Arg-37 → Gly, and Met-59 → Leu) appeared to enhance receptor association, others resulted in weaker receptor interaction (Tyr-31 → Phe and Ile-43 → Phe), and two gave no detectable signal (Leu-14 → Arg and Glu-46 → Ala). Using PCR-based mutagenesis with proinsulin, we also identified a gain of function mutant (proinsulin Leu-17 → Pro) that allowed for a strong IGF-1–receptor interaction. These data demonstrate that the specificity of the interaction between a hormone and its receptor can be characterized with high efficiency in the two-hybrid system and that novel hormone analogues may be found by this method.

A Cell-free System to Study Regulation of Focal Adhesions and of the Connected Actin Cytoskeleton

Assembly and modulation of focal adhesions during dynamic adhesive processes are poorly understood. We describe here the use of ventral plasma membranes from adherent fibroblasts to explore mechanisms regulating integrin distribution and function in a system that preserves the integration of these receptors into the plasma membrane. We find that partial disruption of the cellular organization responsible for the maintenance of organized adhesive sites allows modulation of integrin distribution by divalent cations. High Ca2+ concentrations induce quasi-reversible diffusion of β1 integrins out of focal adhesions, whereas low Ca2+ concentrations induce irreversible recruitment of β1 receptors along extracellular matrix fibrils, as shown by immunofluorescence and electron microscopy. Both effects are independent from the presence of actin stress fibers in this system. Experiments with cells expressing truncated β1 receptors show that the cytoplasmic portion of β1 is required for low Ca2+-induced recruitment of the receptors to matrix fibrils. Analysis with function-modulating antibodies indicates that divalent cation-mediated receptor distribution within the membrane correlates with changes in the functional state of the receptors. Moreover, reconstitution experiments show that purified α-actinin colocalizes and redistributes with β1 receptors on ventral plasma membranes depleted of actin, implicating binding of α-actinin to the receptors. Finally, we found that recruitment of exogenous actin is specifically restricted to focal adhesions under conditions in which new actin polymerization is inhibited. Our data show that the described system can be exploited to investigate the mechanisms of integrin function in an experimental setup that permits receptor redistribution. The possibility to uncouple, under cell-free conditions, events involved in focal adhesion and actin cytoskeleton assembly should facilitate the comprehension of the underlying molecular mechanisms.

A ribozyme-mediated, gene “knockdown” strategy for the identification of gene function in zebrafish

The zebrafish system offers many unique opportunities for the study of molecular biology. To date, only random mutagenesis, and not directed gene knockouts, have been demonstrated in this system. To more fully develop the potential of the zebrafish system, an approach to effectively inhibit the expression of any targeted gene in the developing zebrafish embryo has been developed. This approach uses a transient, cytoplasmic, T7 expression system, injected into the fertilized zebrafish egg to rapidly produce high levels of a ribozyme directed against the mRNA encoded by the targeted gene to inhibit its expression. In a demonstration of this strategy, expression of the recessive dominant zebrafish no tail gene was effectively inhibited by using this strategy to yield a phenotype identical to that resulting from a known defective mutation in this same gene. This, ribozyme-mediated, message deletion strategy may have use in determining the function of genetic coding sequences of unknown function.