A model-based approach for assessing in vivo combination therapy interactions

We present an approach for evaluating the efficacy of combination antitumor agent schedules that accounts for order and timing of drug administration. Our model-based approach compares in vivo tumor volume data over a time course and offers a quantitative definition for additivity of drug effects, relative to which synergism and antagonism are interpreted. We begin by fitting data from individual mice receiving at most one drug to a differential equation tumor growth/drug effect model and combine individual parameter estimates to obtain population statistics. Using two null hypotheses: (i) combination therapy is consistent with additivity or (ii) combination therapy is equivalent to treating with the more effective single agent alone, we compute predicted tumor growth trajectories and their distribution for combination treated animals. We illustrate this approach by comparing entire observed and expected tumor volume trajectories for a data set in which HER-2/neu-overexpressing MCF-7 human breast cancer xenografts are treated with a humanized, anti-HER-2 monoclonal antibody (rhuMAb HER-2), doxorubicin, or one of five proposed combination therapy schedules.

A transgenic mouse model to analyze CD8+ effector T cell differentiation in vivo

Antigen-specific effector T cells are prerequisite to immune protection, but because of the lack of effector cell-specific markers, their generation and differentiation has been difficult to study. We report that effector cells are highly enriched in a T cell subset that can be specifically identified in transgenic (T-GFP) mice expressing green fluorescent protein (GFP) under control of the murine CD4 promoter and proximal enhancer. Consistent with previous studies of these transcriptional control elements, GFP was strongly and specifically expressed in nearly all resting and short-term activated CD4+ and CD8+ T cells. However, when T-GFP mice were challenged with vaccinia virus, allogeneic tumor cells, or staphylococcal enterotoxin A, the cytotoxic and IFN-γ-producing T cells lost GFP expression. Upon T cell receptor (TCR) ligation by αCD3, sorted GFP+ cells fluxed calcium and proliferated vigorously. In contrast, GFP− effector cells showed a diminished calcium flux and did not proliferate. Instead, they underwent apoptosis unless supplied with exogenous IL-2. By reverse transcription–PCR analysis, the GFP− cells up-regulated the pro-apoptotic molecule, Fas-L, and down-regulated gene expression of the proximal TCR signaling molecule, CD3ζ, and c-jun, a component of the AP-1 transcription factor. Thus, differential regulation of TCR signaling may explain the divergent responses of naïve and effector T cells to antigen stimulation.

A model system to study genomic imprinting of human genes

Somatic-cell hybrids have been shown to maintain the correct epigenetic chromatin states to study developmental globin gene expression as well as gene expression on the active and inactive X chromosomes. This suggests the potential use of somatic-cell hybrids containing either a maternal or a paternal human chromosome as a model system to study known imprinted genes and to identify as-yet-unknown imprinted genes. Testing gene expression by using reverse transcription followed by PCR, we show that functional imprints are maintained at four previously characterized 15q11–q13 loci in hybrids containing a single human chromosome 15 and at two chromosome 11p15 loci in hybrids containing a single chromosome 11. In contrast, three γ-aminobutyric acid type A receptor subunit genes in 15q12–q13 are nonimprinted. Furthermore, we have found that differential DNA methylation imprints at the SNRPN promoter and at a CpG island in 11p15 are also maintained in somatic-cell hybrids. Somatic-cell hybrids therefore are a valid and powerful system for studying known imprinted genes as well as for rapidly identifying new imprinted genes.

Differential Distribution of Dynamin Isoforms in Mammalian Cells

Dynamins are 100-kDa GTPases that are essential for clathrin-coated vesicle formation during receptor-mediated endocytosis. To date, three different dynamin genes have been identified, with each gene expressing at least four different alternatively spliced forms. Currently, it is unclear whether these different dynamin gene products perform distinct or redundant cellular functions. Therefore, the focus of this study was to identify additional spliced variants of dynamin from rat tissues and to define the distribution of the dynamin family members in a cultured rat epithelial cell model (Clone 9 cells). After long-distance reverse transcription (RT)-PCR of mRNA from different rat tissues, the full-length cDNAs encoding the different dynamin isoforms were sequenced and revealed four additional spliced variants for dynamin I and nine for dynamin III. Thus, in rat tissues there are a total of at least 25 different mRNAs produced from the three dynamin genes. Subsequently, we generated stably transfected Clone 9 cells expressing full-length cDNAs of six different spliced forms tagged with green fluorescent protein. Confocal or fluorescence microscopy of these transfected cells revealed that many of the dynamin proteins associate with distinct membrane compartments, which include clathrin-coated pits at the plasma membrane and the Golgi apparatus, and several undefined vesicle populations. These results indicate that the dynamin family is more extensive than was originally predicted and suggest that the different dynamin proteins are localized to distinct cytoplasmic or membrane compartments.

Differential quantal release of histamine and 5-hydroxytryptamine from mast cells of vesicular monoamine transporter 2 knockout mice

The recent availability of mice lacking the neuronal form of the vesicular monoamine transporter 2 (VMAT2) affords the opportunity to study its roles in storage and release. Carbon fiber microelectrodes were used to measure individual secretory events of histamine and 5-hydroxytryptamine (5-HT) from VMAT2-expressing mast cells as a model system for quantal release. VMAT2 is indispensable for monoamine storage because mast cells from homozygous (VMAT2−/−) mice, while undergoing granule-cell fusion, do not release monoamines. Cells from heterozygous animals (VMAT2+/−) secrete lower amounts of monoamine per granule than cells from wild-type controls. Investigation of corelease of histamine and 5-HT from granules in VMAT2+/− cells revealed 5-HT quantal size was reduced more than that of histamine. Thus, although vesicular transport is the limiting factor determining quantal size of 5-HT and histamine release, intragranular association with the heparin matrix also plays a significant role.

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.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Regulation of Oleoresinosis in Grand Fir (Abies grandis)1 : Differential Transcriptional Control of Monoterpene, Sesquiterpene, and Diterpene Synthase Genes in Response to Wounding

Grand fir (Abies grandis Lindl.) has been developed as a model system for the study of wound-induced oleoresinosis in conifers as a response to insect attack. Oleoresin is a roughly equal mixture of turpentine (85% monoterpenes [C10] and 15% sesquiterpenes [C15]) and rosin (diterpene [C20] resin acids) that acts to seal wounds and is toxic to both invading insects and their pathogenic fungal symbionts. The dynamic regulation of wound-induced oleoresin formation was studied over 29 d at the enzyme level by in vitro assay of the three classes of synthases directly responsible for the formation of monoterpenes, sesquiterpenes, and diterpenes from the corresponding C10, C15, and C20 prenyl diphosphate precursors, and at the gene level by RNA-blot hybridization using terpene synthase class-directed DNA probes. In overall appearance, the shapes of the time-course curves for all classes of synthase activities are similar, suggesting coordinate formation of all of the terpenoid types. However, closer inspection indicates that the monoterpene synthases arise earlier, as shown by an abbreviated time course over 6 to 48 h. RNA-blot analyses indicated that the genes for all three classes of enzymes are transcriptionally activated in response to wounding, with the monoterpene synthases up-regulated first (transcripts detectable 2 h after wounding), in agreement with the results of cell-free assays of monoterpene synthase activity, followed by the coordinately regulated sesquiterpene synthases and diterpene synthases (transcription beginning on d 3–4). The differential timing in the production of oleoresin components of this defense response is consistent with the immediate formation of monoterpenes to act as insect toxins and their later generation at solvent levels for the mobilization of resin acids responsible for wound sealing.

Differential replication of a single, UV-induced lesion in the leading or lagging strand by a human cell extract: fork uncoupling or gap formation.

We have constructed simian virus 40 minireplicons containing uniquely placed cis,syn-thymine dimers (T <> T) for the analysis of leading- and lagging-strand bypass replication. Assaying for replication in a human cell-free extract through the analysis of full-size labeled product molecules and restriction fragments spanning the T <> T site resulted in the following findings: (i) The primary site of synthesis blockage with T <> T in either the leading or lagging strand was one nucleotide before the lesion. (ii) Replicative bypass of T <> T was detected in both leading and lagging strands. The efficiency of synthesis past T <> T was 22% for leading-strand T <> T and 13% for lagging-strand T <> T. (iii) The lagging-strand T <> T resulted in blocked retrograde synthesis with the replication fork proceeding past the lesion, resulting in daughter molecules containing small gaps (form II' DNA). (iv) With T <> T in the leading-strand template, both the leading and lagging strands were blocked, representing a stalled replication fork. Uncoupling of the concerted synthesis of the two strands of the replication fork was observed, resulting in preferential elongation of the undamaged lagging strand. These data support a model of selective reinitiation downstream from the lesion on lagging strands due to Okazaki synthesis, with no reinitiation close to the damage site on leading strands [Meneghini, R. & Hanawalt, P.C. (1976) Biochim. Biophys. Acta 425, 428-437].

A mechanism for the differential regulation of gonadotropin subunit gene expression by gonadotropin-releasing hormone.

The hypothalamic hormone gonadotropin-releasing hormone (GnRH) is released in a pulsatile fashion, with its frequency varying throughout the reproductive cycle. Varying pulse frequencies and amplitudes differentially regulate the biosynthesis and secretion of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) by pituitary gonadotropes. The mechanism by which this occurs remains a major question in reproductive physiology. Previous studies have been limited by lack of available cell lines that express the LH and FSH subunit genes and respond to GnRH. We have overcome this limitation by transfecting the rat pituitary GH3 cell line with rat GnRH receptor (GnRHR) cDNA driven by a heterologous promoter. These cells, when cotransfected with regulatory regions of the common alpha, LH beta, or FSH beta subunit gene fused to a luciferase reporter gene, respond to GnRH with an increase in luciferase activity. Using this model, we demonstrate that different cell surface densities of the GnRHR result in the differential regulation of LH and FSH subunit gene expression by GnRH. This suggests that the differential regulation of gonadotropin subunit gene expression by GnRH observed in vivo in rats may, in turn, be mediated by varying gonadotrope cell surface GnRHR concentrations. This provides a physiologic mechanism by which a single ligand can act through a single receptor to regulate differentially the production of two hormones in the same cell.

Comparative expressed-sequence-tag analysis of differential gene expression profiles in PC-12 cells before and after nerve growth factor treatment.

Nerve growth factor-induced differentiation of adrenal chromaffin PC-12 cells to a neuronal phenotype involves alterations in gene expression and represents a model system to study neuronal differentiation. We have used the expressed-sequence-tag approach to identify approximately 600 differentially expressed mRNAs in untreated and nerve growth factor-treated PC-12 cells that encode proteins with diverse structural and biochemical functions. Many of these mRNAs encode proteins belonging to cellular pathways not previously known to be regulated by nerve growth factor. Comparative expressed-sequence-tag analysis provides a basis for surveying global changes in gene-expression patterns in response to biological signals at an unprecedented scale, is a powerful tool for identifying potential interactions between different cellular pathways, and allows the gene-expression profiles of individual genes belonging to a particular pathway to be followed.