Water Deficit and Spatial Pattern of Leaf Development. Variability in Responses Can Be Simulated Using a Simple Model of Leaf Development1

We analyzed the effect of short-term water deficits at different periods of sunflower (Helianthus annuus L.) leaf development on the spatial and temporal patterns of tissue expansion and epidermal cell division. Six water-deficit periods were imposed with similar and constant values of soil water content, predawn leaf water potential and [ABA] in the xylem sap, and with negligible reduction of the rate of photosynthesis. Water deficit did not affect the duration of expansion and division. Regardless of their timing, deficits reduced relative expansion rate by 36% and relative cell division rate by 39% (cells blocked at the G0-G1 phase) in all positions within the leaf. However, reductions in final leaf area and cell number in a given zone of the leaf largely differed with the timing of deficit, with a maximum effect for earliest deficits. Individual cell area was only affected during the periods when division slowed down. These behaviors could be simulated in all leaf zones and for all timings by assuming that water deficit affects relative cell division rate and relative expansion rate independently, and that leaf development in each zone follows a stable three-phase pattern in which duration of each phase is stable if expressed in thermal time (C. Granier and F. Tardieu [1998b] Plant Cell Environ 21: 695–703).

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Synergistic inhibition of tumor growth in a murine mammary adenocarcinoma model by combinational gene therapy using IL-12, pro-IL-18, and IL-1β converting enzyme cDNA

We report here that a cancer gene therapy protocol using a combination of IL-12, pro-IL-18, and IL-1β converting enzyme (ICE) cDNA expression vectors simultaneously delivered via gene gun can significantly augment antitumor effects, evidently by generating increased levels of bioactive IL-18 and consequently IFN-γ. First, we compared the levels of IFN-γ secreted by mouse splenocytes stimulated with tumor cells transfected with various test genes, including IL-12 alone; pro-IL-18 alone; pro-IL-18 and ICE; IL-12 and pro-IL-18; and IL-12, pro-IL-18, and ICE. Among these treatments, the combination of IL-12, pro-IL-18, and ICE cDNA resulted in the highest level of IFN-γ production from splenocytes in vitro, and similar results were obtained when these same treatments were delivered to the skin of a mouse by gene gun and IFN-γ levels were measured at the skin transfection site in vivo. Furthermore, the triple gene combinatorial gene therapy protocol was the most effective among all tested groups at suppressing the growth of TS/A (murine mammary adenocarcinoma) tumors previously implanted intradermally at the skin site receiving DNA transfer by gene gun on days 6, 8, 10, and 12 after tumor implantation. Fifty percent of mice treated with the combined three-gene protocol underwent complete tumor regression. In vivo depletion experiments showed that this antitumor effect was CD8+ T cell-mediated and partially IFN-γ-dependent. These results suggest that a combinatorial gene therapy protocol using a mixture of IL-12, pro-IL-18, and ICE cDNAs can confer potent antitumor activities against established TS/A tumors via cytotoxic CD8+ T cells and IFN-γ-dependent pathways.

A human model for multigenic inheritance: Phenotypic expression in Hirschsprung disease requires both the RET gene and a new 9q31 locus

Reduced penetrance in genetic disorders may be either dependent or independent of the genetic background of gene carriers. Hirschsprung disease (HSCR) demonstrates a complex pattern of inheritance with ≈50% of familial cases being heterozygous for mutations in the receptor tyrosine kinase RET. Even when identified, the penetrance of RET mutations is only 50–70%, gender-dependent, and varies with the extent of aganglionosis. We searched for additional susceptibility genes which, in conjunction with RET, lead to phenotypic expression by studying 12 multiplex HSCR families. Haplotype analysis and extensive mutation screening demonstrated three types of families: six families harboring severe RET mutations (group I); and the six remaining families, five of which are RET-linked families with no sequence alterations and one RET-unlinked family (group II). Although the presence of RET mutations in group I families is sufficient to explain HSCR inheritance, a genome scan reveals a new susceptibility locus on 9q31 exclusively in group II families. As such, the gene at 9q31 is a modifier of HSCR penetrance. These observations imply that identification of new susceptibility factors in a complex disease may depend on classification of families by mutational type at known susceptibility genes.

A model for individual and collective cell movement in Dictyostelium discoideum

The cellular slime mold Dictyostelium discoideum is a widely used model system for studying a variety of basic processes in development, including cell–cell signaling, signal transduction, pattern formation, cell motility, and the movement of tissue-like aggregates of cells. Many aspects of cell motion are poorly understood, including how individual cell behavior produces the collective motion of cells observed within the mound and slug. Herein, we describe a biologically realistic model for motile D. discoideum cells that can generate active forces, that interact via surface molecules, and that can detect and respond to chemotactic signals. We model the cells as deformable viscoelastic ellipsoids and incorporate signal transduction and cell–cell signaling by using a previously developed model. The shape constraint restricts the admissible deformations but makes the simulation of a large number of interacting cells feasible. Because the model is based on known processes, the parameters can be estimated or measured experimentally. We show that this model can reproduce the observations on the chemotactic behavior of single cells, streaming during aggregation, and the collective motion of an aggregate of cells driven by a small group of pacemakers. The model predicts that the motion of two-dimensional slugs [Bonner, J. T. (1998) Proc. Natl. Acad. Sci. USA 95, 9355–9359] results from the same behaviors that are exhibited by individual cells; it is not necessary to invoke different mechanisms or behaviors. Our computational experiments also suggest previously uncharacterized phenomena that may be experimentally observable.

Geometric constraints explain much of the species richness pattern in African birds

The world contains boundaries (e.g., continental edge for terrestrial taxa) that impose geometric constraints on the distribution of species ranges. Thus, contrary to traditional thinking, the expected species richness pattern in absence of ecological or physiographical factors is unlikely to be uniform. Species richness has been shown to peak in the middle of a bounded one-dimensional domain, even in the absence of ecological or physiographical factors. Because species ranges are not linear, an extension of the approach to two dimensions is necessary. Here we present a two-dimensional null model accounting for effects of geometric constraints. We use the model to examine the effects of continental edge on the distribution of terrestrial animals in Africa and compare the predictions with the observed pattern of species richness in birds endemic to the continent. Latitudinal, longitudinal, and two-dimensional patterns of species richness are predicted well from the modeled null effects alone. As expected, null effects are of high significance for wide ranging species only. Our results highlight the conceptual significance of an until recently neglected constraint from continental shape alone and support a more cautious analysis of species richness patterns at this scale.

Evidence in support of a four transmembrane-pore-transmembrane topology model for the Arabidopsis thaliana Na+/K+ translocating AtHKT1 protein, a member of the superfamily of K+ transporters

The Arabidopsis thaliana AtHKT1 protein, a Na+/K+ transporter, is capable of mediating inward Na+ currents in Xenopus laevis oocytes and K+ uptake in Escherichia coli. HKT1 proteins are members of a superfamily of K+ transporters. These proteins have been proposed to contain eight transmembrane segments and four pore-forming regions arranged in a mode similar to that of a K+ channel tetramer. However, computer analysis of the AtHKT1 sequence identified eleven potential transmembrane segments. We have investigated the membrane topology of AtHKT1 with three different techniques. First, a gene fusion alkaline phosphatase study in E. coli clearly defined the topology of the N-terminal and middle region of AtHKT1, but the model for membrane folding of the C-terminal region had to be refined. Second, with a reticulocyte-lysate supplemented with dog-pancreas microsomes, we demonstrated that N-glycosylation occurs at position 429 of AtHKT1. An engineered unglycosylated protein variant, N429Q, mediated Na+ currents in X. laevis oocytes with the same characteristics as the wild-type protein, indicating that N-glycosylation is not essential for the functional expression and membrane targeting of AtHKT1. Five potential glycosylation sites were introduced into the N429Q. Their pattern of glycosylation supported the model based on the E. coli-alkaline phosphatase data. Third, immunocytochemical experiments with FLAG-tagged AtHKT1 in HEK293 cells revealed that the N and C termini of AtHKT1, and the regions containing residues 135–142 and 377–384, face the cytosol, whereas the region of residues 55–62 is exposed to the outside. Taken together, our results show that AtHKT1 contains eight transmembrane-spanning segments.

Molecular genetics of pattern formation in the inner ear: Do compartment boundaries play a role?

The membranous labyrinth of the inner ear establishes a precise geometrical topology so that it may subserve the functions of hearing and balance. How this geometry arises from a simple ectodermal placode is under active investigation. The placode invaginates to form the otic cup, which deepens before pinching off to form the otic vesicle. By the vesicle stage many genes expressed in the developing ear have assumed broad, asymmetrical expression domains. We have been exploring the possibility that these domains may reflect developmental compartments that are instrumental in specifying the location and identity of different parts of the ear. The boundaries between compartments are proposed to be the site of inductive interactions required for this specification. Our work has shown that sensory organs and the endolymphatic duct each arise near the boundaries of broader gene expression domains, lending support to this idea. A further prediction of the model, that the compartment boundaries will also represent lineage-restriction compartments, is supported in part by fate mapping the otic cup. Our data suggest that two lineage-restriction boundaries intersect at the dorsal pole of the otocyst, a convergence that may be critical for the specification of endolymphatic duct outgrowth. We speculate that the patterning information necessary to establish these two orthogonal boundaries may emanate, in part, from the hindbrain. The compartment boundary model of ear development now needs to be tested through a variety of experimental perturbations, such as the removal of boundaries, the generation of ectopic boundaries, and/or changes in compartment identity.

Synaptic pattern formation during cellular recognition

Cell–cell recognition often requires the formation of a highly organized pattern of receptor proteins (a synapse) in the intercellular junction. Recent experiments [e.g., Monks, C. R. F., Freiberg, B. A., Kupfer, H., Sciaky, N. & Kupfer, A. (1998) Nature (London) 395, 82–86; Grakoui, A., Bromley, S. K., Sumen, C., Davis, M. M., Shaw, A. S., Allen, P. M. & Dustin, M. L. (1999) Science 285, 221–227; and Davis, D. M., Chiu, I., Fassett, M., Cohen, G. B., Mandelboim, O. & Strominger, J. L. (1999) Proc. Natl. Acad. Sci. USA 96, 15062–15067] vividly demonstrate a complex evolution of cell shape and spatial receptor–ligand patterns (several microns in size) in the intercellular junction during immunological synapse formation. The current view is that this dynamic rearrangement of proteins into organized supramolecular activation clusters is driven primarily by active cytoskeletal processes [e.g., Dustin, M. L. & Cooper, J. A. (2000) Nat. Immunol. 1, 23–29; and Wulfing, C. & Davis, M. M. (1998) Science 282, 2266–2269]. Here, aided by a quantitative analysis of the relevant physico-chemical processes, we demonstrate that the essential characteristics of synaptic patterns observed in living cells can result from spontaneous self-assembly processes. Active cellular interventions are superimposed on these self-organizing tendencies and may also serve to regulate the spontaneous processes. We find that the protein binding/dissociation characteristics, protein mobilities, and membrane constraints measured in the cellular environment are delicately balanced such that the length and time scales of spontaneously evolving patterns are in near-quantitative agreement with observations for synapse formation between T cells and supported membranes [Grakoui, A., Bromley, S. K., Sumen, C., Davis, M. M., Shaw, A. S., Allen, P. M. & Dustin, M. L. (1999) Science 285, 221–227]. The model we present provides a common way of analyzing immunological synapse formation in disparate systems (e.g., T cell/antigen-presenting cell junctions with different MHC-peptides, natural killer cells, etc.).

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

A model of self-thinning through local competition.

Explanations of self-thinning in plant populations have focused on plant shape and packing. A dynamic model based on the structure of local interactions successfully reproduces the pattern and can be approximated to identify key parameters and relationships. The approach generates testable new explanations for differences between species and populations, unifies self-thinning with other patterns in plant population dynamics, and indicates why organisms other than plants can follow the law.

A joint hazard and time scaling model to compare survival curves.

To provide a more general method for comparing survival experience, we propose a model that independently scales both hazard and time dimensions. To test the curve shape similarity of two time-dependent hazards, h1(t) and h2(t), we apply the proposed hazard relationship, h12(tKt)/ h1(t) = Kh, to h1. This relationship doubly scales h1 by the constant hazard and time scale factors, Kh and Kt, producing a transformed hazard, h12, with the same underlying curve shape as h1. We optimize the match of h12 to h2 by adjusting Kh and Kt. The corresponding survival relationship S12(tKt) = [S1(t)]KtKh transforms S1 into a new curve S12 of the same underlying shape that can be matched to the original S2. We apply this model to the curves for regional and local breast cancer contained in the National Cancer Institute's End Results Registry (1950-1973). Scaling the original regional curves, h1 and S1 with Kt = 1.769 and Kh = 0.263 produces transformed curves h12 and S12 that display congruence with the respective local curves, h2 and S2. This similarity of curve shapes suggests the application of the more complete curve shapes for regional disease as templates to predict the long-term survival pattern for local disease. By extension, this similarity raises the possibility of scaling early data for clinical trial curves according to templates of registry or previous trial curves, projecting long-term outcomes and reducing costs. The proposed model includes as special cases the widely used proportional hazards (Kt = 1) and accelerated life (KtKh = 1) models.

A structural model for a homeotic protein-extradenticle-DNA complex accounts for the choice of HOX protein in the heterodimer.

The genes of the homeotic complex (HOX) encode DNA binding homeodomain proteins that control developmental fates by differentially regulating the transcription of downstream target genes. Despite their unique in vivo functions, disparate HOX proteins often bind to very similar DNA sequences in vitro. Thus, a critical question is how HOX proteins select the correct sets of target genes in vivo. The homeodomain proteins encoded by the Drosophila extradenticle gene and its mammalian homologues, the pbx genes, contribute to HOX specificity by cooperatively binding to DNA with HOX proteins. For example, the HOX protein labial cooperatively binds with extradenticle protein to a 20-bp oligonucleotide that is sufficient to direct a labial-like expression pattern in Drosophila embryos. Here we have analyzed the protein-DNA interactions that are important for forming the labial-extradenticle-DNA complex. The data suggest a model in which labial and extradenticle, separated by only 4 bp, bind this DNA as a heterodimer in a head-to-tail orientation. We have confirmed several aspects of this model by characterizing extradenticle-HOX binding to mutant oligonucleotides. Most importantly, mutations in base pairs predicted to contact the HOX N-terminal arm resulted in a change in HOX preference in the heterodimer, from labial to Ultrabithorax. These results demonstrate that extradenticle prefers to bind cooperatively with different HOX proteins depending on subtle differences in the heterodimer binding site.

A general model of invariant chain association with class II major histocompatibility complex proteins.

The binding of invariant chain to major histocompatibility complex (MHC) proteins is an important step in processing of MHC class II proteins and in antigen presentation. The question of how invariant chain can bind to all MHC class II proteins is central to understanding these processes. We have employed molecular modeling to predict the structure of class II-associated invariant chain peptide (CLIP)-MHC protein complexes and to ask whether the predicted mode of association could be general across all MHC class II proteins. CLIP fits identically into the MHC class II alleles HLA-DR3, I-Ak, I-Au, and I-Ad, with a consistent pattern of hydrogen bonds, contacts, and hydrophobic burial and without bad contacts. Our model predicts the burial of CLIP residues Met-91 and Met-99 in the deep P1 and P9 anchor pockets and other detailed interactions, which we have compared with available data. The predicted pattern of I-A allele-specific effects on CLIP binding is very similar to that observed experimentally by alanine-scanning mutations of CLIP. Together, these results indicate that CLIP may bind in a single, general way across products of MHC class II alleles.