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.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Spatial modulation of primate inferotemporal responses by eye position.

BACKGROUND: A key aspect of representations for object recognition and scene analysis in the ventral visual stream is the spatial frame of reference, be it a viewer-centered, object-centered, or scene-based coordinate system. Coordinate transforms from retinocentric space to other reference frames involve combining neural visual responses with extraretinal postural information. METHODOLOGY/PRINCIPAL FINDINGS: We examined whether such spatial information is available to anterior inferotemporal (AIT) neurons in the macaque monkey by measuring the effect of eye position on responses to a set of simple 2D shapes. We report, for the first time, a significant eye position effect in over 40% of recorded neurons with small gaze angle shifts from central fixation. Although eye position modulates responses, it does not change shape selectivity. CONCLUSIONS/SIGNIFICANCE: These data demonstrate that spatial information is available in AIT for the representation of objects and scenes within a non-retinocentric frame of reference. More generally, the availability of spatial information in AIT calls into questions the classic dichotomy in visual processing that associates object shape processing with ventral structures such as AIT but places spatial processing in a separate anatomical stream projecting to dorsal structures.

Dynamics of a minimal model of interlocked positive and negative feedback loops of transcriptional regulation by cAMP-response element binding proteins.

cAMP-response element binding (CREB) proteins are involved in transcriptional regulation in a number of cellular processes (e.g., neural plasticity and circadian rhythms). The CREB family contains activators and repressors that may interact through positive and negative feedback loops. These loops can be generated by auto- and cross-regulation of expression of CREB proteins, via CRE elements in or near their genes. Experiments suggest that such feedback loops may operate in several systems (e.g., Aplysia and rat). To understand the functional implications of such feedback loops, which are interlocked via cross-regulation of transcription, a minimal model with a positive and negative loop was developed and investigated using bifurcation analysis. Bifurcation analysis revealed diverse nonlinear dynamics (e.g., bistability and oscillations). The stability of steady states or oscillations could be changed by time delays in the synthesis of the activator (CREB1) or the repressor (CREB2). Investigation of stochastic fluctuations due to small numbers of molecules of CREB1 and CREB2 revealed a bimodal distribution of CREB molecules in the bistability region. The robustness of the stable HIGH and LOW states of CREB expression to stochastic noise differs, and a critical number of molecules was required to sustain the HIGH state for days or longer. Increasing positive feedback or decreasing negative feedback also increased the lifetime of the HIGH state, and persistence of this state may correlate with long-term memory formation. A critical number of molecules was also required to sustain robust oscillations of CREB expression. If a steady state was near a deterministic Hopf bifurcation point, stochastic resonance could induce oscillations. This comparative analysis of deterministic and stochastic dynamics not only provides insights into the possible dynamics of CREB regulatory motifs, but also demonstrates a framework for understanding other regulatory processes with similar network architecture.

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.

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

A reduced model clarifies the role of feedback loops and time delays in the Drosophila circadian oscillator.

Although several detailed models of molecular processes essential for circadian oscillations have been developed, their complexity makes intuitive understanding of the oscillation mechanism difficult. The goal of the present study was to reduce a previously developed, detailed model to a minimal representation of the transcriptional regulation essential for circadian rhythmicity in Drosophila. The reduced model contains only two differential equations, each with time delays. A negative feedback loop is included, in which PER protein represses per transcription by binding the dCLOCK transcription factor. A positive feedback loop is also included, in which dCLOCK indirectly enhances its own formation. The model simulated circadian oscillations, light entrainment, and a phase-response curve with qualitative similarities to experiment. Time delays were found to be essential for simulation of circadian oscillations with this model. To examine the robustness of the simplified model to fluctuations in molecule numbers, a stochastic variant was constructed. Robust circadian oscillations and entrainment to light pulses were simulated with fewer than 80 molecules of each gene product present on average. Circadian oscillations persisted when the positive feedback loop was removed. Moreover, elimination of positive feedback did not decrease the robustness of oscillations to stochastic fluctuations or to variations in parameter values. Such reduced models can aid understanding of the oscillation mechanisms in Drosophila and in other organisms in which feedback regulation of transcription may play an important role.

Epstein-Barr virus infection of Langerhans cell precursors as a mechanism of oral epithelial entry, persistence, and reactivation.

Epstein-Barr virus (EBV) is a ubiquitous human herpesvirus associated with many malignant and nonmalignant human diseases. Life-long latent EBV persistence occurs in blood-borne B lymphocytes, while EBV intermittently productively replicates in mucosal epithelia. Although several models have previously been proposed, the mechanism of EBV transition between these two reservoirs of infection has not been determined. In this study, we present the first evidence demonstrating that EBV latently infects a unique subset of blood-borne mononuclear cells that are direct precursors to Langerhans cells and that EBV both latently and productively infects oral epithelium-resident cells that are likely Langerhans cells. These data form the basis of a proposed new model of EBV transition from blood to oral epithelium in which EBV-infected Langerhans cell precursors serve to transport EBV to the oral epithelium as they migrate and differentiate into oral Langerhans cells. This new model contributes fresh insight into the natural history of EBV infection and the pathogenesis of EBV-associated epithelial disease.

Direct intrathecal implantation of mesenchymal stromal cells leads to enhanced neuroprotection via an NFkappaB-mediated increase in interleukin-6 production.

Mesenchymal stromal cell (MSC) therapy has shown promise for the treatment of traumatic brain injury (TBI). Although the mechanism(s) by which MSCs offer protection is unclear, initial in vivo work has suggested that modulation of the locoregional inflammatory response could explain the observed benefit. We hypothesize that the direct implantation of MSCs into the injured brain activates resident neuronal stem cell (NSC) niches altering the intracerebral milieu. To test our hypothesis, we conducted initial in vivo studies, followed by a sequence of in vitro studies. In vivo: Sprague-Dawley rats received a controlled cortical impact (CCI) injury with implantation of 1 million MSCs 6 h after injury. Brain tissue supernatant was harvested for analysis of the proinflammatory cytokine profile. In vitro: NSCs were transfected with a firefly luciferase reporter for NFkappaB and placed in contact culture and transwell culture. Additionally, multiplex, quantitative PCR, caspase 3, and EDU assays were completed to evaluate NSC cytokine production, apoptosis, and proliferation, respectively. In vivo: Brain supernatant analysis showed an increase in the proinflammatory cytokines IL-1alpha, IL-1beta, and IL-6. In vitro: NSC NFkappaB activity increased only when in contact culture with MSCs. When in contact with MSCs, NSCs show an increase in IL-6 production as well as a decrease in apoptosis. Direct implantation of MSCs enhances neuroprotection via activation of resident NSC NFkappaB activity (independent of PI3 kinase/AKT pathway) leading to an increase in IL-6 production and decrease in apoptosis. In addition, the observed NFkappaB activity depends on direct cell contact.

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.

Lobe specific Ca2+-calmodulin nano-domain in neuronal spines: a single molecule level analysis.

Calmodulin (CaM) is a ubiquitous Ca(2+) buffer and second messenger that affects cellular function as diverse as cardiac excitability, synaptic plasticity, and gene transcription. In CA1 pyramidal neurons, CaM regulates two opposing Ca(2+)-dependent processes that underlie memory formation: long-term potentiation (LTP) and long-term depression (LTD). Induction of LTP and LTD require activation of Ca(2+)-CaM-dependent enzymes: Ca(2+)/CaM-dependent kinase II (CaMKII) and calcineurin, respectively. Yet, it remains unclear as to how Ca(2+) and CaM produce these two opposing effects, LTP and LTD. CaM binds 4 Ca(2+) ions: two in its N-terminal lobe and two in its C-terminal lobe. Experimental studies have shown that the N- and C-terminal lobes of CaM have different binding kinetics toward Ca(2+) and its downstream targets. This may suggest that each lobe of CaM differentially responds to Ca(2+) signal patterns. Here, we use a novel event-driven particle-based Monte Carlo simulation and statistical point pattern analysis to explore the spatial and temporal dynamics of lobe-specific Ca(2+)-CaM interaction at the single molecule level. We show that the N-lobe of CaM, but not the C-lobe, exhibits a nano-scale domain of activation that is highly sensitive to the location of Ca(2+) channels, and to the microscopic injection rate of Ca(2+) ions. We also demonstrate that Ca(2+) saturation takes place via two different pathways depending on the Ca(2+) injection rate, one dominated by the N-terminal lobe, and the other one by the C-terminal lobe. Taken together, these results suggest that the two lobes of CaM function as distinct Ca(2+) sensors that can differentially transduce Ca(2+) influx to downstream targets. We discuss a possible role of the N-terminal lobe-specific Ca(2+)-CaM nano-domain in CaMKII activation required for the induction of synaptic plasticity.

CART: an Hrs/actinin-4/BERP/myosin V protein complex required for efficient receptor recycling.

Altering the number of surface receptors can rapidly modulate cellular responses to extracellular signals. Some receptors, like the transferrin receptor (TfR), are constitutively internalized and recycled to the plasma membrane. Other receptors, like the epidermal growth factor receptor (EGFR), are internalized after ligand binding and then ultimately degraded in the lysosome. Routing internalized receptors to different destinations suggests that distinct molecular mechanisms may direct their movement. Here, we report that the endosome-associated protein hrs is a subunit of a protein complex containing actinin-4, BERP, and myosin V that is necessary for efficient TfR recycling but not for EGFR degradation. The hrs/actinin-4/BERP/myosin V (CART [cytoskeleton-associated recycling or transport]) complex assembles in a linear manner and interrupting binding of any member to its neighbor produces an inhibition of transferrin recycling rate. Disrupting the CART complex results in shunting receptors to a slower recycling pathway that involves the recycling endosome. The novel CART complex may provide a molecular mechanism for the actin-dependence of rapid recycling of constitutively recycled plasma membrane receptors.

Poly(A)-binding protein-interacting protein 1 binds to eukaryotic translation initiation factor 3 to stimulate translation.

Poly(A)-binding protein (PABP) stimulates translation initiation by binding simultaneously to the mRNA poly(A) tail and eukaryotic translation initiation factor 4G (eIF4G). PABP activity is regulated by PABP-interacting (Paip) proteins. Paip1 binds PABP and stimulates translation by an unknown mechanism. Here, we describe the interaction between Paip1 and eIF3, which is direct, RNA independent, and mediated via the eIF3g (p44) subunit. Stimulation of translation by Paip1 in vivo was decreased upon deletion of the N-terminal sequence containing the eIF3-binding domain and upon silencing of PABP or several eIF3 subunits. We also show the formation of ternary complexes composed of Paip1-PABP-eIF4G and Paip1-eIF3-eIF4G. Taken together, these data demonstrate that the eIF3-Paip1 interaction promotes translation. We propose that eIF3-Paip1 stabilizes the interaction between PABP and eIF4G, which brings about the circularization of the mRNA.

Lipid-protein interactions drive membrane protein topogenesis in accordance with the positive inside rule.

Transmembrane domain orientation within some membrane proteins is dependent on membrane lipid composition. Initial orientation occurs within the translocon, but final orientation is determined after membrane insertion by interactions within the protein and between lipid headgroups and protein extramembrane domains. Positively and negatively charged amino acids in extramembrane domains represent cytoplasmic retention and membrane translocation forces, respectively, which are determinants of protein orientation. Lipids with no net charge dampen the translocation potential of negative residues working in opposition to cytoplasmic retention of positive residues, thus allowing the functional presence of negative residues in cytoplasmic domains without affecting protein topology.

Activation of heat shock and antioxidant responses by the natural product celastrol: transcriptional signatures of a thiol-targeted molecule.

Stress response pathways allow cells to sense and respond to environmental changes and adverse pathophysiological states. Pharmacological modulation of cellular stress pathways has implications in the treatment of human diseases, including neurodegenerative disorders, cardiovascular disease, and cancer. The quinone methide triterpene celastrol, derived from a traditional Chinese medicinal herb, has numerous pharmacological properties, and it is a potent activator of the mammalian heat shock transcription factor HSF1. However, its mode of action and spectrum of cellular targets are poorly understood. We show here that celastrol activates Hsf1 in Saccharomyces cerevisiae at a similar effective concentration seen in mammalian cells. Transcriptional profiling revealed that celastrol treatment induces a battery of oxidant defense genes in addition to heat shock genes. Celastrol activated the yeast Yap1 oxidant defense transcription factor via the carboxy-terminal redox center that responds to electrophilic compounds. Antioxidant response genes were likewise induced in mammalian cells, demonstrating that the activation of two major cell stress pathways by celastrol is conserved. We report that celastrol's biological effects, including inhibition of glucocorticoid receptor activity, can be blocked by the addition of excess free thiol, suggesting a chemical mechanism for biological activity based on modification of key reactive thiols by this natural product.