Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Molecular cross talk between misfolded proteins in animal models of Alzheimer's and prion diseases.

The central event in protein misfolding disorders (PMDs) is the accumulation of a misfolded form of a naturally expressed protein. Despite the diversity of clinical symptoms associated with different PMDs, many similarities in their mechanism suggest that distinct pathologies may cross talk at the molecular level. The main goal of this study was to analyze the interaction of the protein misfolding processes implicated in Alzheimer's and prion diseases. For this purpose, we inoculated prions in an Alzheimer's transgenic mouse model that develop typical amyloid plaques and followed the progression of pathological changes over time. Our findings show a dramatic acceleration and exacerbation of both pathologies. The onset of prion disease symptoms in transgenic mice appeared significantly faster with a concomitant increase on the level of misfolded prion protein in the brain. A striking increase in amyloid plaque deposition was observed in prion-infected mice compared with their noninoculated counterparts. Histological and biochemical studies showed the association of the two misfolded proteins in the brain and in vitro experiments showed that protein misfolding can be enhanced by a cross-seeding mechanism. These results suggest a profound interaction between Alzheimer's and prion pathologies, indicating that one protein misfolding process may be an important risk factor for the development of a second one. Our findings may have important implications to understand the origin and progression of PMDs.

Molecular basis of Corynebacterium diphtheriae virulence and infection in the Caenorhabditis elegans model host

Corynebacterium diphtheriae is the causative agent of cutaneous and pharyngeal diphtheria in humans. While lethality is certainly caused by diphtheria toxin, corynebacterial colonization may primarily require proteinaceous fibers called pili, which mediate adherence to specific tissues. The type strain of C. diphtheriae possesses three distinct pilus structures, namely the SpaA, SpaD, and SpaH-type pili, which are encoded by three distinct pilus gene clusters. The pilus is assembled onto the bacterial peptidoglycan by a specific transpeptidase enzyme called sortase. Although the SpaA pili are shown to be specific for pharyngeal cells in vitro, little is known about functions of the three pili in bacterial pathogenesis. This is mainly due to lack of in vivo models of corynebacterial infection. As an alternative to mouse models as mice do not have functional receptors for diphtheria toxin, in this study I use Caenorhabditis elegans as a model host for C. diphtheriae. A simple C. elegans model would be beneficial in determining the specific role of each pilus-type and the literature suggests that C. elegans infection model can be used to study a variety of bacterial species giving insight into bacterial virulence and host-pathogen interactions. My study examines the hypothesis that pili and toxin are major virulent determinants of C. diphtheriae in the C. elegans model host.

Repeat length and RNA expression level are not primary determinants in CUG expansion toxicity in Drosophila models.

Evidence for an RNA gain-of-function toxicity has now been provided for an increasing number of human pathologies. Myotonic dystrophies (DM) belong to a class of RNA-dominant diseases that result from RNA repeat expansion toxicity. Specifically, DM of type 1 (DM1), is caused by an expansion of CUG repeats in the 3'UTR of the DMPK protein kinase mRNA, while DM of type 2 (DM2) is linked to an expansion of CCUG repeats in an intron of the ZNF9 transcript (ZNF9 encodes a zinc finger protein). In both pathologies the mutant RNA forms nuclear foci. The mechanisms that underlie the RNA pathogenicity seem to be rather complex and not yet completely understood. Here, we describe Drosophila models that might help unravelling the molecular mechanisms of DM1-associated CUG expansion toxicity. We generated transgenic flies that express inducible repeats of different type (CUG or CAG) and length (16, 240, 480 repeats) and then analyzed transgene localization, RNA expression and toxicity as assessed by induced lethality and eye neurodegeneration. The only line that expressed a toxic RNA has a (CTG)(240) insertion. Moreover our analysis shows that its level of expression cannot account for its toxicity. In this line, (CTG)(240.4), the expansion inserted in the first intron of CG9650, a zinc finger protein encoding gene. Interestingly, CG9650 and (CUG)(240.4) expansion RNAs were found in the same nuclear foci. In conclusion, we suggest that the insertion context is the primary determinant for expansion toxicity in Drosophila models. This finding should contribute to the still open debate on the role of the expansions per se in Drosophila and in human pathogenesis of RNA-dominant diseases.

Cellular dynamic simulator: an event driven molecular simulation environment for cellular physiology.

In this paper, we present the Cellular Dynamic Simulator (CDS) for simulating diffusion and chemical reactions within crowded molecular environments. CDS is based on a novel event driven algorithm specifically designed for precise calculation of the timing of collisions, reactions and other events for each individual molecule in the environment. Generic mesh based compartments allow the creation / importation of very simple or detailed cellular structures that exist in a 3D environment. Multiple levels of compartments and static obstacles can be used to create a dense environment to mimic cellular boundaries and the intracellular space. The CDS algorithm takes into account volume exclusion and molecular crowding that may impact signaling cascades in small sub-cellular compartments such as dendritic spines. With the CDS, we can simulate simple enzyme reactions; aggregation, channel transport, as well as highly complicated chemical reaction networks of both freely diffusing and membrane bound multi-protein complexes. Components of the CDS are generally defined such that the simulator can be applied to a wide range of environments in terms of scale and level of detail. Through an initialization GUI, a simple simulation environment can be created and populated within minutes yet is powerful enough to design complex 3D cellular architecture. The initialization tool allows visual confirmation of the environment construction prior to execution by the simulator. This paper describes the CDS algorithm, design implementation, and provides an overview of the types of features available and the utility of those features are highlighted in demonstrations.

Molecular mechanism by which the nucleoid occlusion factor, SlmA, keeps cytokinesis in check.

In Escherichia coli, cytokinesis is orchestrated by FtsZ, which forms a Z-ring to drive septation. Spatial and temporal control of Z-ring formation is achieved by the Min and nucleoid occlusion (NO) systems. Unlike the well-studied Min system, less is known about the anti-DNA guillotining NO process. Here, we describe studies addressing the molecular mechanism of SlmA (synthetic lethal with a defective Min system)-mediated NO. SlmA contains a TetR-like DNA-binding fold, and chromatin immunoprecipitation analyses show that SlmA-binding sites are dispersed on the chromosome except the Ter region, which segregates immediately before septation. SlmA binds DNA and FtsZ simultaneously, and the SlmA-FtsZ structure reveals that two FtsZ molecules sandwich a SlmA dimer. In this complex, FtsZ can still bind GTP and form protofilaments, but the separated protofilaments are forced into an anti-parallel arrangement. This suggests that SlmA may alter FtsZ polymer assembly. Indeed, electron microscopy data, showing that SlmA-DNA disrupts the formation of normal FtsZ polymers and induces distinct spiral structures, supports this. Thus, the combined data reveal how SlmA derails Z-ring formation at the correct place and time to effect NO.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

Molecular mechanisms for the insulin sensitizing activity of rexinoids in skeletal muscle of diabetic (Db/Db) mice

Rexinoids are synthetic agonists for the retinoid X receptors (RXRs), a member of the nuclear receptor family of ligand-activated transcription factors. Rexinoids have been shown to lower serum glucose and insulin levels in animal models of type 2 diabetes. However the mechanisms that are responsible for the insulin-sensitizing action of rexinoids are largely unknown. Skeletal muscle accounts for the majority of insulin-regulated whole-body glucose disposal and impaired insulin action in muscle is an important contributor to the pathophysiology of type 2 diabetes. Glucose transport is a rate-limiting step in glucose utilization. The goal of these studies is to examine the mechanisms of the anti-diabetic activity of rexinoids in skeletal muscle of diabetic db/db mice. The results we have obtained showed that treatment of db/db mice with rexinoids for two weeks resulted in a significant increase in insulin-stimulated glucose transport activity in skeletal muscle. Insulin stimulates glucose transport in muscle via the regulation of both the insulin receptor substrate-1 (IRS-1)/Akt pathway and the Cbl-associated protein (CAP)/Cbl pathway. Rexinoids increased the insulin-stimulated IRS-1 tyrosine phosphorylation and Akt phosphorylation without effects on the activity of the CAP/Cbl pathway. The effects of rexinoids on the IRS-1/Akt pathway were associated with a decrease in the level of IRS-1 Serine 307 phosphorylation as well as qualitative and quantitative alterations in the fatty acyl-CoAs present within the muscle cells. In addition, rexinoids increased the expression of uncoupling protein 3 (UCP3) and activation of AMPK in diabetic muscle. This effect may also enhance the IRS-1/Akt signaling. We believe that it is the concerted activation of the IRS-1/Akt and AMPK signaling systems, a pharmacological mechanism that as far as we know, is unique to rexinoids, that results in the anti-diabetic effects of these drugs. Our results also suggest that the glucose-lowering mechanism of rexinoids is distinct from that of the thiazolidinediones (TZDs), peroxisome proliferator-activated receptor γ (PPARγ) agonists with well-characterized anti-diabetic activity. Rexinoids appear to represent a novel class of insulin sensitizers, with potential applications for the treatment of type 2 diabetes. ^

Molecular genetic analyses of extracellular signaling pathways during murine retinal development

Extracellular signaling pathways initiated by secreted proteins are important in the co-ordination of tissue interactions in multi-cellular organisms, particularly during embryonic development. These signaling cascades direct diverse cellular events, including proliferation, differentiation and migration, in both autocrine and paracrine modes. In adult animals, abnormal function of these proteins often results in degenerative and tumourigenic syndromes. In this study, I have focused on elucidating the role of Bone Morphogenetic Protein (Bmp) signal transduction during neuronal specification and differentiation in the vertebrate embryo, using the mouse retina as a model. Using tissue-specific conditional knock-out approaches, the consequences of genetic loss-of-function of this signaling pathway on retinal physiology were examined. Mutant mice lacking Bmp type I receptor function displayed a range of retinal phenotypes, each of which appeared to be regulated at a different threshold of Bmp receptor activity. Novel essential functions for Bmp signaling were uncovered for retinal neurogenesis, cell survival, and axonal pathfinding at the optic disc. Further, BmprIa and BmprIa exhibited genetic interactions suggestive of functional redundancy. To further characterize the underlying molecular bases for the pleiotropic effects of Bmp receptors, retina-specific loss-of-function mutants of the obligate Bmp-activated transcriptional mediator Smad4 were generated. A comparison of the retina-specific Smad4 mutant phenotypes with those of the Bmp receptor mutant retina revealed that only a subset of retinal phenotypes, namely optic disc axon pathfinding and axial patterning were common for both classes of mutant animals. Thus, these results suggest that, contrary to the classic scheme of Bmp signal transduction, Smad4-independent pathways may be operative downstream of the type I receptors. Indeed, such alternative intracellular signaling cascades may constitute a molecular basis for the multiple cellular responses elicited by Bmp signaling. Finally, I tested whether the potential Bmp pathway targets, the extracellular ligands Fgf9 and Fgf15, mediate essential cellular processes in the retina. The analyses of Fgf9 −/−; Fgf15−/− mutant mice posit a novel shared role for these genes in intra-retinal axon pathfinding. Collectively, these studies have elucidated part of the molecular machinery directing mammalian neuro-retinal development, and provided useful in vivo models to study visual function. ^

Molecular mechanisms of prostacyclin receptor signaling: The intracellular domains in G protein coupling

To better understand the mechanisms of how the human prostacyclin receptor (1P) mediates vasodilation and platelet anti-aggregation through Gs protein coupling, a strategy integrating multiple approaches including high resolution NMR experiments, synthetic peptide, fluorescence spectroscopy, molecular modeling, and recombinant protein was developed and used to characterize the structure/function relationship of important segments and residues of the IP receptor and the α-subunit of the Gs protein (Gαs). The first (iLP1) and third (iLP3) intracellular loops of the IP receptor, as well as the Gαs C-terminal domain, relevant to the Gs-mediated IP receptor signaling, were first identified by observation of the effects of the mini gene-expressed corresponding protein segments in HEK293 cells which co-expressed the receptor and Gαs. Evidence of the IP iLP1 domain interacted with the Gαs C-terminal domain was observed by fluorescence and NMR spectroscopic studies using a constrained synthetic peptide, which mimicked the IP iLP1 domain, and the synthetic peptide, which mimicked Gαs C-terminal domain. The solution structural models and the peptide-peptide interaction of the two synthetic protein segments were determined by high resolution NMR spectroscopy. The important residues in the corresponding domains of the IP receptor and the Gαs predicted by NMR chemical shift mapping were used to guide the identification of their protein-protein interaction in cells. A profile of the residues Arg42 - Ala48 of the IP iLP1 domain and the three residues Glu392 ∼ Leu394 of the Gαs C-terminal domain involved in the IP/Gs protein coupling were confirmed by recombinant proteins. The data revealed an intriguing speculation on the mechanisms of how the signal of the ligand-activated IP receptor is transmitted to the Gs protein in regulating vascular functions and homeostasis, and also provided substantial insights into other prostanoid receptor signaling. ^