Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

Induction of solid tumor differentiation by the peroxisome proliferator-activated receptor-γ ligand troglitazone in patients with liposarcoma

Agonist ligands for the nuclear receptor peroxisome proliferator-activated receptor-γ have been shown to induce terminal differentiation of normal preadipocytes and human liposarcoma cells in vitro. Because the differentiation status of liposarcoma is predictive of clinical outcomes, modulation of the differentiation status of a tumor may favorably impact clinical behavior. We have conducted a clinical trial for treatment of patients with advanced liposarcoma by using the peroxisome proliferator-activated receptor-γ ligand troglitazone, in which extensive correlative laboratory studies of tumor differentiation were performed. We report here the results of three patients with intermediate to high-grade liposarcomas in whom troglitazone administration induced histologic and biochemical differentiation in vivo. Biopsies of tumors from each of these patients while on troglitazone demonstrated histologic evidence of extensive lipid accumulation by tumor cells and substantial increases in NMR-detectable tumor triglycerides compared with pretreatment biopsies. In addition, expression of several mRNA transcripts characteristic of differentiation in the adipocyte lineage was induced. There was also a marked reduction in immunohistochemical expression of Ki-67, a marker of cell proliferation. Together, these data indicate that terminal adipocytic differentiation was induced in these malignant tumors by troglitazone. These results indicate that lineage-appropriate differentiation can be induced pharmacologically in a human solid tumor.

Induction of ARF tumor suppressor gene expression and cell cycle arrest by transcription factor DMP1

Expression of the DMP1 transcription factor, a cyclin D-binding Myb-like protein, induces growth arrest in mouse embryo fibroblast strains but is devoid of antiproliferative activity in primary diploid fibroblasts that lack the ARF tumor suppressor gene. DMP1 binds to a single canonical recognition site in the ARF promoter to activate gene expression, and in turn, p19ARF synthesis causes p53-dependent cell cycle arrest. Unlike genes such as Myc, adenovirus E1A, and E2F-1, which, when overexpressed, activate the ARF-p53 pathway and trigger apoptosis, DMP1, like ARF itself, does not induce programmed cell death. Therefore, apart from its recently recognized role in protecting cells from potentially oncogenic signals, ARF can be induced in response to antiproliferative stimuli that do not obligatorily lead to apoptosis.

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.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

A subset of skin tumor modifier loci determines survival time of tumor-bearing mice

Studies of mouse models of human cancer have established the existence of multiple tumor modifiers that influence parameters of cancer susceptibility such as tumor multiplicity, tumor size, or the probability of malignant progression. We have carried out an analysis of skin tumor susceptibility in interspecific Mus musculus/Mus spretus hybrid mice and have identified another seven loci showing either significant (six loci) or suggestive (one locus) linkage to tumor susceptibility or resistance. A specific search was carried out for skin tumor modifier loci associated with time of survival after development of a malignant tumor. A combination of resistance alleles at three markers [D6Mit15 (Skts12), D7Mit12 (Skts2), and D17Mit7 (Skts10)], all of which are close to or the same as loci associated with carcinoma incidence and/or papilloma multiplicity, is significantly associated with increased survival of mice with carcinomas, whereas the reverse combination of susceptibility alleles is significantly linked to early mortality caused by rapid carcinoma growth (χ2 = 25.22; P = 5.1 × 10−8). These data indicate that host genetic factors may be used to predict carcinoma growth rate and/or survival of individual backcross mice exposed to the same carcinogenic stimulus and suggest that mouse models may provide an approach to the identification of genetic modifiers of cancer survival in humans.

Growth suppression of glioma cells by PTEN requires a functional phosphatase catalytic domain

Deletions of all or part of chromosome 10 are the most common genetic alterations in high-grade gliomas. The PTEN gene (also called MMAC1 and TEP1) maps to chromosome region 10q23 and has been implicated as a target of alteration in gliomas and also in other cancers such as those of the breast, prostate, and kidney. Here we sought to provide a functional test of its candidacy as a growth suppressor in glioma cells. We used a combination of Northern blot analysis, protein truncation assays, and sequence analysis to determine the types and frequency of PTEN mutations in glioma cell lines so that we could define appropriate recipients to assess the growth suppressive function of PTEN by gene transfer. Introduction of wild-type PTEN into glioma cells containing endogenous mutant alleles caused growth suppression, but was without effect in cells containing endogenous wild-type PTEN. The ectopic expression of PTEN alleles, which carried mutations found in primary tumors and have been shown or are expected to inactivate its phosphatase activity, caused little growth suppression. These data strongly suggest that PTEN is a protein phosphatase that exhibits functional and specific growth-suppressing activity.

Three-dimensional modeling of and ligand docking to vitamin D receptor ligand binding domain

The ligand binding domain of the human vitamin D receptor (VDR) was modeled based on the crystal structure of the retinoic acid receptor. The ligand binding pocket of our VDR model is spacious at the helix 11 site and confined at the β-turn site. The ligand 1α,25-dihydroxyvitamin D3 was assumed to be anchored in the ligand binding pocket with its side chain heading to helix 11 (site 2) and the A-ring toward the β-turn (site 1). Three residues forming hydrogen bonds with the functionally important 1α- and 25-hydroxyl groups of 1α,25-dihydroxyvitamin D3 were identified and confirmed by mutational analysis: the 1α-hydroxyl group is forming pincer-type hydrogen bonds with S237 and R274 and the 25-hydroxyl group is interacting with H397. Docking potential for various ligands to the VDR model was examined, and the results are in good agreement with our previous three-dimensional structure-function theory.

Highly efficient electro-gene therapy of solid tumor by using an expression plasmid for the herpes simplex virus thymidine kinase gene

We report successful electro-gene therapy (EGT) by using plasmid DNA for tumor-bearing mice. Subcutaneously inoculated CT26 tumor was subjected to EGT, which consists of intratumoral injection of a naked plasmid encoding a marker gene or a therapeutic gene, followed by in vivo electroporation (EP). When this treatment modality is carried out with the plasmid DNA for the green fluorescent protein gene, followed by in vivo EP with the optimized pulse parameters, numerous intensely bright green fluorescent signals appeared within the tumor. EGT, by using the “A” fragment of the diphtheria toxin gene significantly inhibited the growth of tumors, by about 30%, on the flank of mice. With the herpes simplex virus thymidine kinase gene, followed by systemic injection of ganciclovir, EGT was far more effective in retarding tumor growth, varying between 50% and 90%, compared with the other controls. Based on these results, it appears that EGT can be used successfully for treating murine solid tumors.

Diverse roles of the tumor necrosis factor family member TRANCE in skeletal physiology revealed by TRANCE deficiency and partial rescue by a lymphocyte-expressed TRANCE transgene

Tumor necrosis factor-related, activation-induced cytokine (TRANCE), a tumor necrosis factor family member, mediates survival of dendritic cells in the immune system and is required for osteoclast differentiation and activation in the skeleton. We report the skeletal phenotype of TRANCE-deficient mice and its rescue by the TRANCE transgene specifically expressed in lymphocytes. TRANCE-deficient mice showed severe osteopetrosis, with no osteoclasts, marrow spaces, or tooth eruption, and exhibited profound growth retardation at several skeletal sites, including the limbs, skull, and vertebrae. These mice had marked chondrodysplasia, with thick, irregular growth plates and a relative increase in hypertrophic chondrocytes. Transgenic overexpression of TRANCE in lymphocytes of TRANCE-deficient mice rescued osteoclast development in two locations in growing long bones: excavation of marrow cavities permitting hematopoiesis in the marrow spaces, and remodeling of osteopetrotic woven bone in the shafts of long bones into histologically normal lamellar bone. However, osteoclasts in these mice failed to appear at the chondroosseous junction and the metaphyseal periosteum of long bones, nor were they present in tooth eruption pathways. These defects resulted in sclerotic metaphyses with persistence of club-shaped long bones and unerupted teeth, and the growth plate defects were largely unimproved by the TRANCE transgene. Thus, TRANCE-mediated regulation of the skeleton is complex, and impacts chondrocyte differentiation and osteoclast formation in a manner that likely requires local delivery of TRANCE.

A meanfield approach to the thermodynamics of a protein-solvent system with application to the oligomerization of the tumor suppressor p53

The thermodynamic stability and oligomerization status of the tumor suppressor p53 tetramerization domain have been studied experimentally and theoretically. A series of hydrophilic mutations at Met-340 and Leu-344 of human p53 were designed to disrupt the hydrophobic dimer–dimer interface of the tetrameric oligomerization domain of p53 (residues 325–355). Meanfield calculations of the free energy of the solvated mutants as a function of interdimer distance were compared with experimental data on the thermal stability and oligomeric state (tetramer, dimer, or equilibrium mixture of both) of each mutant. The calculations predicted a decreasing stability and oligomeric state for the following amino acids at residue 340: Met (tetramer) > Ser Asp, His, Gln, > Glu, Lys (dimer), whereas the experimental results showed the following order: Met (tetramer) > Ser > Gln > His, Lys > Asp, Glu (dimers). For residue 344, the calculated trend was Leu (tetramer) > Ala > Arg, Gln, Lys (dimer), and the experimental trend was Leu (tetramer) > Ala, Arg, Gln, Lys (dimer). The discrepancy for the lysine side chain at residue 340 is attributed to the dual nature of lysine, both hydrophobic and charged. The incorrect prediction of stability of the mutant with Asp at residue 340 is attributed to the fact that within the meanfield approach, we use the wild-type backbone configuration for all mutants, but low melting temperatures suggest a softening of the α-helices at the dimer–dimer interface. Overall, this initial application of meanfield theory toward a protein-solvent system is encouraging for the application of the theoretical model to more complex systems.

A numerical strategy for efficient modeling of the equatorial wave guide

Convection in the tropics is observed to involve a wide-ranging hierarchy of scales from a few kilometers to the planetary scales and also has a profound impact on short-term climate. The mechanisms responsible for this behavior present a major unsolved problem. A promising emerging approach to address these issues is cloud-resolving modeling. Here a family of numerical models is introduced specifically to model the feedback of small-scale deep convection on tropical planetary waves and tropical circulation in a highly efficient manner compatible with the approach through cloud-resolving modeling. Such a procedure is also useful for theoretical purposes. The basic idea in the approach is to use low-order truncation in the meriodonal direction through Gauss–Hermite quadrature projected onto a simple discrete radiation condition. In this fashion, the cloud-resolving modeling of equatorially trapped planetary waves reduces to the solution of a small number of purely zonal two-dimensional wave systems along a few judiciously chosen meriodonal layers that are coupled only by some additional source terms. The approach is analyzed in detail with full mathematical rigor for linearized equatorial primitive equations with source terms.

Embryonic expression of the tumor-associated PAX3-FKHR fusion protein interferes with the developmental functions of Pax3

A unique chromosomal translocation involving the genes PAX3 and FKHR is characteristic of most human alveolar rhabdomyosarcomas. The resultant chimeric protein fuses the PAX3 DNA-binding domains to the transactivation domain of FKHR, suggesting that PAX3-FKHR exerts its role in alveolar rhabdomyosarcomas through dysregulation of PAX3-specific target genes. Here, we have produced transgenic mice in which PAX3-FKHR expression was driven by mouse Pax3 promoter/enhancer sequences. Five independent lines expressed PAX3-FKHR in the dorsal neural tube and lateral dermomyotome. Each line exhibited phenotypes that correlated with PAX3-FKHR expression levels and predominantly involved pigmentary disturbances of the abdomen, hindpaws, and tail, with additional neurological related alterations. Phenotypic severity could be increased by reducing Pax3 levels through matings with Pax3-defective Splotch mice, and interference between PAX3 and PAX3-FKHR was apparent in transcription reporter assays. These data suggest that the tumor-associated PAX3-FKHR fusion protein interferes with normal Pax3 developmental functions as a prelude to transformation.

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

Statistical modeling of large microarray data sets to identify stimulus-response profiles

A statistical modeling approach is proposed for use in searching large microarray data sets for genes that have a transcriptional response to a stimulus. The approach is unrestricted with respect to the timing, magnitude or duration of the response, or the overall abundance of the transcript. The statistical model makes an accommodation for systematic heterogeneity in expression levels. Corresponding data analyses provide gene-specific information, and the approach provides a means for evaluating the statistical significance of such information. To illustrate this strategy we have derived a model to depict the profile expected for a periodically transcribed gene and used it to look for budding yeast transcripts that adhere to this profile. Using objective criteria, this method identifies 81% of the known periodic transcripts and 1,088 genes, which show significant periodicity in at least one of the three data sets analyzed. However, only one-quarter of these genes show significant oscillations in at least two data sets and can be classified as periodic with high confidence. The method provides estimates of the mean activation and deactivation times, induced and basal expression levels, and statistical measures of the precision of these estimates for each periodic transcript.