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.

Neurogranin controls the spatiotemporal pattern of postsynaptic Ca2+/CaM signaling.

Neurogranin (Ng) is a postsynaptic IQ-motif containing protein that accelerates Ca(2+) dissociation from calmodulin (CaM), a key regulator of long-term potentiation and long-term depression in CA1 pyramidal neurons. The exact physiological role of Ng, however, remains controversial. Two genetic knockout studies of Ng showed opposite outcomes in terms of the induction of synaptic plasticity. To understand its function, we test the hypothesis that Ng could regulate the spatial range of action of Ca(2+)/CaM based on its ability to accelerate the dissociation of Ca(2+) from CaM. Using a mathematical model constructed on the known biochemistry of Ng, we calculate the cycle time that CaM molecules alternate between the fully Ca(2+) saturated state and the Ca(2+) unbound state. We then use these results and include diffusion of CaM to illustrate the impact that Ng has on modulating the spatial profile of Ca(2+)-saturated CaM within a model spine compartment. Finally, the first-passage time of CaM to transition from the Ca(2+)-free state to the Ca(2+)-saturated state was calculated with or without Ng present. These analyses suggest that Ng regulates the encounter rate between Ca(2+) saturated CaM and its downstream targets during postsynaptic Ca(2+) transients.

Role of the N- and C-lobes of calmodulin in the activation of Ca(2+)/calmodulin-dependent protein kinase II.

Understanding the principles of calmodulin (CaM) activation of target enzymes will help delineate how this seemingly simple molecule can play such a complex role in transducing Ca (2+)-signals to a variety of downstream pathways. In the work reported here, we use biochemical and biophysical tools and a panel of CaM constructs to examine the lobe specific interactions between CaM and CaMKII necessary for the activation and autophosphorylation of the enzyme. Interestingly, the N-terminal lobe of CaM by itself was able to partially activate and allow autophosphorylation of CaMKII while the C-terminal lobe was inactive. When used together, CaMN and CaMC produced maximal CaMKII activation and autophosphorylation. Moreover, CaMNN and CaMCC (chimeras of the two N- or C-terminal lobes) both activated the kinase but with greater K act than for wtCaM. Isothermal titration calorimetry experiments showed the same rank order of affinities of wtCaM > CaMNN > CaMCC as those determined in the activity assay and that the CaM to CaMKII subunit binding ratio was 1:1. Together, our results lead to a proposed sequential mechanism to describe the activation pathway of CaMKII led by binding of the N-lobe followed by the C-lobe. This mechanism contrasts the typical sequential binding mode of CaM with other CaM-dependent enzymes, where the C-lobe of CaM binds first. The consequence of such lobe specific binding mechanisms is discussed in relation to the differential rates of Ca (2+)-binding to each lobe of CaM during intracellular Ca (2+) oscillations.

Nonviolent factors: Variables associated with ethnic violence, and factors which are associated with nonviolent response to grievances toward members of other ethnic groups

Ethnic violence appears to be the major source of violence in the world. Ethnic hostilities are potentially all-pervasive because most countries in the world are multi-ethnic. Public health's focus on violence documents its increasing role in this issue.^ The present study is based on a secondary analysis of a dataset of responses by 272 individuals from four ethnic groups (Anglo, African, Mexican, and Vietnamese Americans) who answered questions regarding variables related to ethnic violence from a general questionnaire which was distributed to ethnically diverse purposive, nonprobability, self-selected groups of individuals in Houston, Texas, in 1993.^ One goal was psychometric: learning about issues in analysis of datasets with modest numbers, comparison of two approaches to dealing with missing observations not missing at random (conducting analysis on two datasets), transformation analysis of continuous variables for logistic regression, and logistic regression diagnostics.^ Regarding the psychometric goal, it was concluded that measurement model analysis was not possible with a relatively small dataset with nonnormal variables, such as Likert-scaled variables; therefore, exploratory factor analysis was used. The two approaches to dealing with missing values resulted in comparable findings. Transformation analysis suggested that the continuous variables were in the correct scale, and diagnostics that the model fit was adequate.^ The substantive portion of the analysis included the testing of four hypotheses. Hypothesis One proposed that attitudes/efficacy regarding alternative approaches to resolving grievances from the general questionnaire represented underlying factors: nonpunitive social norms and strategies for addressing grievances--using the political system, organizing protests, using the system to punish offenders, and personal mediation. Evidence was found to support all but one factor, nonpunitive social norms.^ Hypothesis Two proposed that the factor variables and the other independent variables--jail, grievance, male, young, and membership in a particular ethnic group--were associated with (non)violence. Jail, grievance, and not using the political system to address grievances were associated with a greater likelihood of intergroup violence.^ No evidence was found to support Hypotheses Three and Four, which proposed that grievance and ethnic group membership would interact with other variables (i.e., age, gender, etc.) to produce variant levels of subgroup (non)violence.^ The generalizability of the results of this study are constrained by the purposive self-selected nature of the sample and small sample size (n = 272).^ Suggestions for future research include incorporating other possible variables or factors predictive of intergroup violence in models of the kind tested here, and the development and evaluation of interventions that promote electoral and nonelectoral political participation as means of reducing interethnic conflict. ^

MYELINATION IN UNDERNOURISHED RAT BRAIN

Several interactive parameters of protein-calorie malnutrition imposed during postnatal ontogeny on the myelination of rat brain wre investigated. Postnatal starvation depresses the rate of myelin protein synthesis to approximately the same extent in all major brain regions examined (cerebral cortex, cerebellum, striatum, hippocampus, hypothalamus, midbrain and medulla), indicating a relatively uniform reduction in myelination throughout the brain. Early starvation from birth through 8 days, as well as starvation occurring late, from 14 to 30 days, produced no lasting deficit in myelin accumulation. Starvation from birth through 14 days or from birth through 20 days produces lasting, significant myelin deficits in all brain regions when examined following ad libitum feeding to 60 days of age. These data, in combination with the metabolic studies of myelin synthesis, show that severe starvation occurring during the 2nd and 3rd weeks of postnatal life produces an immediate reduction in myelin synthesis, and that the subsequent deficit in myelin accumulation is irreversible by nutritional rehabilitation. With respect to the relative severity of nutritional restriction occurring during this "critical" interval of brain ontogeny, additional studies showed that mild undernourishment (producing less than 20 percent growth lag) produces no myelin deficit. There appears to be a threshold effect such that undernutrition producing a growth lag of between 20 to 30 percent first produces a measurable deficit. Increasingly severe regimens of nutritional restriction which produce approximately 30, 40 and 50 percent body weight lags result in initial myelin deficits of 25, 55 and 60 percent, respectively. Initial myelin deficits do not recover following nutritional rehabilitation, although myelin continues to increase in both normal and all undernourished populations. At the cellular level, severe postnatal nutritional restriction appears to depress both the initial synthesis of myelin precursor proteins (as demonstrated for proteolipid protein) as well as their subsequent assembly into myelin membrane. All of the findings of the present studies are consistent with a hypothetical model of undernutrition-induced brain hypomyelination in which the primary defect consists of a failure of oligodendroglia to myelinate a substantial percentage of axons, resulting in a greatly decreased ratio of myelinated to unmyelinated axons. ^

14-3-3zeta overexpression in multiple human cancers plays a critical role in cancer progression

14-3-3 is a family of highly conserved and ubiquitously expressed proteins in eukaryotic organisms. 14-3-3 isoforms bind in a phospho-serine/threonine-dependent manner to a host of proteins involved in essential cellular processes including cell cycle, signal transduction and apoptosis. We fortuitously discovered 14-3-3 zeta overexpression in many human primary cancers, such as breast, lung, and sarcoma, and in a majority of cancer cell lines. To determine 14-3-3 zeta involvement in breast cancer progression, we used immunohistochemical analysis to examine 14-3-3 zeta expression in human primary invasive breast carcinomas. High 14-3-3 zeta expression was significantly correlated with poor prognosis of breast cancer patients. Increased expression of 14-3-3 zeta was also significantly correlated with elevated PKB/Akt activation in patient samples. Thus, 14-3-3 zeta is a marker of poor prognosis in breast cancers. Furthermore, up-regulation of 14-3-3 zeta enhanced malignant transformation of cancer cells in vitro. ^ To determine the biological significance of 14-3-3 zeta in human cancers, small interfering RNAs (siRNA) were used to specifically block 14-3-3 zeta expression in cancer cells. 14-3-3 zeta siRNA inhibited cellular proliferation by inducing a G1 arrest associated with up-regulation of p27 KIP1 and p21CIP1 cyclin dependent kinase inhibitors. Reduced 14-3-3 zeta inhibited PKB/Akt activation while stimulating the p38 signaling pathway. Silencing 14-3-3 zeta expression also increased stress-induced apoptosis by caspase activation. Notably, 14-3-3 zeta siRNA inhibited transformation related properties of breast cancer cells in vitro and inhibited tumor progression of breast cancer cells in vivo. 14-3-3 zeta may be a key regulatory factor controlling multiple signaling pathways leading to tumor progression. ^ The data indicate 14-3-3 zeta is a major regulator of cell growth and apoptosis and may play a critical role in the development of multiple cancer types. Hence, blocking 14-3-3 zeta may be a promising therapeutic approach for numerous cancers. ^

A proposed application of the attributable risk technique in the assessment of unmet need among newly organized populations

In the midst of health care reform, and as health care organizations reorganize to provide more cost-effective healthcare, the population is being shifted into new healthcare delivery systems such as health insurance purchasing alliances, and health maintenance organizations. These new models of delivery are usually organized within resource restricted and data limited environments. Health care planners are faced with the challenge of identifying priorities for preventive and primary care services within these newly organized populations (Medicare HMO, Medicaid HMO, etc.). The author proposes a technique usually employed in epidemiology--attributable risk estimation--as a planning methodology to establish preventive health priorities within newly organized populations. Illustrations of the methodology are provided utilizing the Texas 1992 population. ^

Platelets and Anti-Angiogenic Resistance in Ovarian Carcinoma

Background: Resistance to targeted anti-angiogenic therapy is a growing clinical concern given the disappointing clinical impact of anti-angiogenic. Platelets represent a component of the tumor microenvironment that are implicated in metastasis and represent a significant reservoir of angiogenic regulators. Thrombocytosis has been shown to be caused by malignancy and associated with adverse clinical outcomes, however the causal connections between these associations remain to be identified. Materials and Methods: Following IRB approval, patient data were collected on patients from four U.S. centers and platelet levels through and after therapy were considered as indicators of recurrence of disease. In vitro effects of platelets on cancer cell proliferation, apoptosis, and migration were examined. RNA interference was used to query signaling pathways mediating these effects. The necessity of platelet activation for in vitro effect was analyzed. In vivo orthotopic models were used to query the impact of thrombocytosis and thrombocytopenia on the efficacy of cytotoxic chemotherapy, the effect of aspirin on thrombocytosis and cancer, and platelet effect on anti-angiogenic therapy. Results: Platelets were found to increase at the time of diagnosis of ovarian cancer recurrence in a pattern comparable to CA-125. Platelet co-culture increased proliferation, increased migration, and decreased apoptosis in all cell lines tested. RNA interference implicated platelet derived growth factor alpha (PDGFRA) and transforming growth factor beta-receptor 1 (TGFBR1) signaling. Biodistribution studies suggested minimal platelet sequestration of taxanes. Blockade of platelet activation blocked in vitro effects. In vivo, thrombocytosis blocked chemotherapeutic efficacy, thrombocytopenia increased chemotherapeutic efficacy, and aspirin therapy partially blocked the effects of thrombocytosis. In vivo, withdrawal of anti-angiogenic therapy caused loss of therapeutic benefit with evidence of accelerated disease growth. This effect was blocked by use of a small-molecule inhibitor of Focal Adhesion Kinase. Anti-angiogenic therapy was also associated with increased platelet infiltration into tumor that was not seen to the same degree in the control or FAK-inhibitor-treated mice. Conclusions: Platelets are active participants in the growth and metastasis of tumor, both directly and via facilitation of angiogenesis. Blocking platelets, blocking platelet activation, and blocking platelet trafficking into tumor are novel therapeutic avenues supported by this data. Copyright © 2012 Justin Neal Bottsford-Miller, all rights reserved.

Role of Neurogranin in the regulation of calcium binding to Calmodulin

Role of Neurogranin in the regulation of calcium binding to Calmodulin Anuja Chandrasekar, B.S Advisor: M. Neal Waxham, Ph.D The overall goal of my project was to gain a quantitative understanding of how the interaction between two proteins neurogranin (RC3) and calmodulin (CaM) alters a fundamental property of CaM. CaM, has been extensively studied for more than four decades due to its seminal role in almost all biological functions as a calcium signal transducer. Calcium signals in cardiac and neuronal cells are exquisitely precise and enable activation of some processes while down-regulating others. CaM, with its four calcium binding sites, serves as a central component of calcium signaling in these cells. It is aided in this role as a regulatory hub that differentially activates targets in response to a calcium flux by proteins that alter its calcium binding properties. Neurogranin, also known as RC3, is a member of a family of small neuronal IQ (SNIQ) domain proteins that was originally thought to play a ‘capacitive’ role by sequestering CaM until a calcium influx of sufficient intensity arrived. However, based on earlier work in our lab on neurogranin, we believe that this protein plays a more nuanced role in neurons than simply acting as a CaM buffer. We believe that neurogranin is one of the proteins which, by altering the kinetics of calcium binding allow CaM to decode a variety of signals with fine precision. To quantify the interaction between CaM, neurogranin and calcium, I used biophysical techniques and computational simulations. From my results, I conclude that neurogranin finely regulates the proportion of calcium-saturated CaM and thereby directs CaM’s target specificity.