Hebbian analysis of the transformation of medial entorhinal grid-cell inputs to hippocampal place fields.

The discovery of grid cells in the medial entorhinal cortex (MEC) permits the characterization of hippocampal computation in much greater detail than previously possible. The present study addresses how an integrate-and-fire unit driven by grid-cell spike trains may transform the multipeaked, spatial firing pattern of grid cells into the single-peaked activity that is typical of hippocampal place cells. Previous studies have shown that in the absence of network interactions, this transformation can succeed only if the place cell receives inputs from grids with overlapping vertices at the location of the place cell's firing field. In our simulations, the selection of these inputs was accomplished by fast Hebbian plasticity alone. The resulting nonlinear process was acutely sensitive to small input variations. Simulations differing only in the exact spike timing of grid cells produced different field locations for the same place cells. Place fields became concentrated in areas that correlated with the initial trajectory of the animal; the introduction of feedback inhibitory cells reduced this bias. These results suggest distinct roles for plasticity of the perforant path synapses and for competition via feedback inhibition in the formation of place fields in a novel environment. Furthermore, they imply that variability in MEC spiking patterns or in the rat's trajectory is sufficient for generating a distinct population code in a novel environment and suggest that recalling this code in a familiar environment involves additional inputs and/or a different mode of operation of the network.

Dominance of the proximal coordinate frame in determining the locations of hippocampal place cell activity during navigation.

The place-specific activity of hippocampal cells provides downstream structures with information regarding an animal's position within an environment and, perhaps, the location of goals within that environment. In rodents, recent research has suggested that distal cues primarily set the orientation of the spatial representation, whereas the boundaries of the behavioral apparatus determine the locations of place activity. The current study was designed to address possible biases in some previous research that may have minimized the likelihood of observing place activity bound to distal cues. Hippocampal single-unit activity was recorded from six freely moving rats as they were trained to perform a tone-initiated place-preference task on an open-field platform. To investigate whether place activity was bound to the room- or platform-based coordinate frame (or both), the platform was translated within the room at an "early" and at a "late" phase of task acquisition (Shift 1 and Shift 2). At both time points, CA1 and CA3 place cells demonstrated room-associated and/or platform-associated activity, or remapped in response to the platform shift. Shift 1 revealed place activity that reflected an interaction between a dominant platform-based (proximal) coordinate frame and a weaker room-based (distal) frame because many CA1 and CA3 place fields shifted to a location intermediate to the two reference frames. Shift 2 resulted in place activity that became more strongly bound to either the platform- or room-based coordinate frame, suggesting the emergence of two independent spatial frames of reference (with many more cells participating in platform-based than in room-based representations).

drp, A novel cell density regulated protein

Growing cells are continuously processing signals of all varieties and responding to these signals by changes in cellular gene expression. One signal that cells in close proximity relay to each other is cell-cell contact. Non-transformed cells respond to cell-cell contact by arrest of growth and entry into G$\sb0,$ a process known as contact inhibition. Transformed cells do not respond to contact inhibition and continue to grow to high cell density, forming foci when in cell culture and tumors in the living organism. The events surrounding the generation, transduction, and response to cellular contact are poorly understood. In the present study, a novel gene product, drp, is shown to be expressed at high levels in cultured cells at high cell density. This density regulated protein, drp, has an apparent molecular weight of 70 kDa. Northern analysis shows drp to be highly expressed in cardiac and skeletal muscle and least abundant in lung and kidney tissues. By homology to two independently derived sequence tagged sites (STSs) used in the human genome project, drp or a closely related sequence maps to human chromosome 12. Density-dependent increases in drp expression have been demonstrated in six different cell lines including NIH 3T3, Hela and a human teratocarcinoma cell line, PA-1. Cells exhibit increased drp expression both when they are plated at increasing concentrations per unit area, or plated at low density and allowed to grow naturally to higher cell density. Cells at high density can exhibit several phenotypes including growth arrest, accumulation of soluble factors in the media, and increased numbers of cell contacts. Growth arrest by serum starvation or TGF-$\beta$ treatment fails to produce an increase in drp expression. Similarly, treatment of low density cells with conditioned media from high density cells fails to elicit drp expression. These results argue that neither soluble factors accumulated or expressed at high density nor simple exit from the cell cycle is sufficient to produce an increase in drp expression. The expression of drp appears to be uniquely regulated by cell density alone. ^

4HPR-X radiation mechanisms of interactions in non-small cell lung cancer (NSCLC) cells in vitro

Lung cancer is the leading cause of cancer death. However, poor survival using conventional therapies fuel the search for more rational interventions. The objective of this study was to design and implement a 4HPR-radiation interaction model in NSCLC, employing a traditional clinical modality (radiation), a relatively new, therapeutically unexplored agent (4HPR) and rationally combining them based on molecular mechanistic findings pertaining to their interactions. To test the hypothesis that 4HPR sensitizes cells to radiation-induced cell death via G2+M accumulation, we designed a working model consisting of H522 adenocarcinoma cells (p53, K-ras mutated) derived from an NSCLC patient; 4HPR at concentrations up to 10 μM; and X radiation up to 6 Gy generated by a patient-dedicated Phillips RT-250 X ray unit at 250 KV, 15 mA, 1.85 Gy/min. We found that 4HPR produced time- and dose-dependent morphological changes, growth inhibition, and DNA damage-inducing enhancement of reactive oxygen species. A transient G2+M accumulation of cells maximal at 24 h of continuous 4HPR exposure was used for irradiation time scheduling. Our data demonstrated enhanced cell death (both apoptotic and necrotic) in irradiated cells pre-treated with 4HPR versus those with either stressor alone. 4HPR's effect of increased NSCLC cells' radioresponse was confirmed by clonogenic assay. To explore these practical findings from a molecular mechanistic perspective, we further investigated and showed that levels of cyclin B1 and p34cdc2 kinase—both components of the mitosis promoting factor (MPF) regulating the G2/M transition—did not change following 4HPR treatment. Likewise, cdc25C phosphatase was not altered. However, enhanced p34cdc2 phosphorylation on its Thr14Tyr15 residues—indicative of its inactivation and increased expression of MPF negative regulators chk1 and wee1 kinases—were supportive of explaining 4HPR-treated cells' accumulation. Hence, p34cdc2 phosphorylation, chk1, and wee1 warrant further evaluation as potential molecular targets for 4HPR-X radiation combination. In summary, we (1) demonstrated that 4HPR not only induces cell death by itself, but also increases NSCLC cells' subsequent radioresponse, indicative of potential clinical applicability, and (2) for the first time, shed light on deciphering 4HPR-X radiation molecular mechanisms of interaction, including the finding of 4HPR's role as a p34cdc2 inactivator via Thr14Tyr15 phosphorylation. ^

Mechanism of damage sensing and cell killing following the inhibition of nucleotide excision repair by fludarabine in quiescent cells

Inhibition of DNA repair by the nucleoside of fludarabine (F-ara-A) induces toxicity in quiescent human cells. The sensing and signaling mechanisms following DNA repair inhibition by F-ara-A are unknown. The central hypothesis of this project was that the mechanistic interaction of a DNA repair initiating agent and a nucleoside analog initiates an apoptotic signal in quiescent cells. The purpose of this research was to identify the sensing and signaling mechanism(s) that respond to DNA repair inhibition by F-ara-A. Lymphocytes were treated with F-ara-A, to accumulate the active triphosphate metabolite and subsequently DNA repair was activated by UV irradiation. Pre-incubation of lymphocytes with 3 μM F-ara-A inhibited DNA repair initiated by 2 J/m2 UV and induced greater than additive apoptosis after 24 h. Blocking the incorporation of F-ara-A nucleotide into repairing DNA using 30 μM aphidicolin considerably lowered the apoptotic response. ^ Wild-type quiescent cells showed a significant loss in viability than did cells lacking functional sensor kinase DNA-PKcs or p53 as measured by colony formation assays. The functional status of ATM did not appear to affect the apoptotic outcome. Immunoprecipitation studies showed an interaction between the catalytic sub-unit of DNA-PK and p53 following DNA repair inhibition. Confocal fluorescence microscopy studies have indicated the localization pattern of p53, DNA-PK and γ-H2AX in the nucleus following DNA damage. Foci formation by γ-H2AX was seen as an early event that is followed by interaction with DNA-PKcs. p53 serine-15 phosphorylation and accumulation were detected 2 h after treatment. Fas/Fas ligand expression increased significantly after repair inhibition and was dependent on the functional status of p53. Blocking the interaction between Fas and Fas ligand by neutralizing antibodies significantly rescued the apoptotic fraction of cells. ^ Collectively, these results suggest that incorporation of the nucleoside analog into repair patches is critical for cytotoxicity and that the DNA damage, while being sensed by DNA-PK, may induce apoptosis by a p53-mediated signaling mechanism. Based on the results, a model is proposed for the sensing of F-ara-A-induced DNA damage that includes γ-H2AX, DNA-PKcs, and p53. Targeting the cellular DNA repair mechanism can be a potential means of producing cytotoxicity in a quiescent population of neoplastic cells. These results also provide mechanistic support for the success of nucleoside analogs with cyclophosphamide or other agents that initiate excision repair processes, in the clinic. ^

αvβ8 Integrin Interacts with RhoGDI1 to Regulate Rac1 and Cdc42 Activation and Drive Glioblastoma Cell Invasion

As an interface between the circulatory and central nervous systems, the neurovascular unit is vital to the development and survival of tumors. The malignant brain cancer glioblastoma multiforme (GBM) displays invasive growth behaviors that are major impediments to surgical resection and targeted therapies. Adhesion and signaling pathways that drive GBM cell invasion remain largely uncharacterized. Here we have utilized human GBM cell lines, primary patient samples, and pre-clinical mouse models to demonstrate that integrin αvβ8 is a major driver of GBM cell invasion. β8 integrin is overexpressed in many human GBM cells, with higher integrin expression correlating with increased invasion and diminished patient survival. Silencing β8 integrin in human GBM cells leads to impaired tumor cell invasion due to hyperactivation of the Rho GTPases Rac1 and Cdc42. β8 integrin associates with Rho GDP Dissociation Inhibitor 1 (RhoGDI1), an intracellular signaling effector that sequesters Rho GTPases in their inactive GDP-bound states. Silencing RhoGDI1 expression or uncoupling αvβ8 integrin-RhoGDI1 protein interactions blocks GBM cell invasion due to Rho GTPase hyperactivation. These data reveal for the first time that αvβ8 integrin, via interactions with RhoGDI1, suppresses activation of Rho proteins to promote GBM cell invasiveness. Hence, targeting the αvβ8 integrin-RhoGDI1 signaling axis may be an effective strategy for blocking GBM cell invasion.

TIME SERIES ANALYSIS AS INPUT FOR PREDICTIVE MODELING: PREDICTING CARDIAC ARREST IN A PEDIATRIC INTENSIVE CARE UNIT

The first manuscript, entitled "Time-Series Analysis as Input for Clinical Predictive Modeling: Modeling Cardiac Arrest in a Pediatric ICU" lays out the theoretical background for the project. There are several core concepts presented in this paper. First, traditional multivariate models (where each variable is represented by only one value) provide single point-in-time snapshots of patient status: they are incapable of characterizing deterioration. Since deterioration is consistently identified as a precursor to cardiac arrests, we maintain that the traditional multivariate paradigm is insufficient for predicting arrests. We identify time series analysis as a method capable of characterizing deterioration in an objective, mathematical fashion, and describe how to build a general foundation for predictive modeling using time series analysis results as latent variables. Building a solid foundation for any given modeling task involves addressing a number of issues during the design phase. These include selecting the proper candidate features on which to base the model, and selecting the most appropriate tool to measure them. We also identified several unique design issues that are introduced when time series data elements are added to the set of candidate features. One such issue is in defining the duration and resolution of time series elements required to sufficiently characterize the time series phenomena being considered as candidate features for the predictive model. Once the duration and resolution are established, there must also be explicit mathematical or statistical operations that produce the time series analysis result to be used as a latent candidate feature. In synthesizing the comprehensive framework for building a predictive model based on time series data elements, we identified at least four classes of data that can be used in the model design. The first two classes are shared with traditional multivariate models: multivariate data and clinical latent features. Multivariate data is represented by the standard one value per variable paradigm and is widely employed in a host of clinical models and tools. These are often represented by a number present in a given cell of a table. Clinical latent features derived, rather than directly measured, data elements that more accurately represent a particular clinical phenomenon than any of the directly measured data elements in isolation. The second two classes are unique to the time series data elements. The first of these is the raw data elements. These are represented by multiple values per variable, and constitute the measured observations that are typically available to end users when they review time series data. These are often represented as dots on a graph. The final class of data results from performing time series analysis. This class of data represents the fundamental concept on which our hypothesis is based. The specific statistical or mathematical operations are up to the modeler to determine, but we generally recommend that a variety of analyses be performed in order to maximize the likelihood that a representation of the time series data elements is produced that is able to distinguish between two or more classes of outcomes. The second manuscript, entitled "Building Clinical Prediction Models Using Time Series Data: Modeling Cardiac Arrest in a Pediatric ICU" provides a detailed description, start to finish, of the methods required to prepare the data, build, and validate a predictive model that uses the time series data elements determined in the first paper. One of the fundamental tenets of the second paper is that manual implementations of time series based models are unfeasible due to the relatively large number of data elements and the complexity of preprocessing that must occur before data can be presented to the model. Each of the seventeen steps is analyzed from the perspective of how it may be automated, when necessary. We identify the general objectives and available strategies of each of the steps, and we present our rationale for choosing a specific strategy for each step in the case of predicting cardiac arrest in a pediatric intensive care unit. Another issue brought to light by the second paper is that the individual steps required to use time series data for predictive modeling are more numerous and more complex than those used for modeling with traditional multivariate data. Even after complexities attributable to the design phase (addressed in our first paper) have been accounted for, the management and manipulation of the time series elements (the preprocessing steps in particular) are issues that are not present in a traditional multivariate modeling paradigm. In our methods, we present the issues that arise from the time series data elements: defining a reference time; imputing and reducing time series data in order to conform to a predefined structure that was specified during the design phase; and normalizing variable families rather than individual variable instances. The final manuscript, entitled: "Using Time-Series Analysis to Predict Cardiac Arrest in a Pediatric Intensive Care Unit" presents the results that were obtained by applying the theoretical construct and its associated methods (detailed in the first two papers) to the case of cardiac arrest prediction in a pediatric intensive care unit. Our results showed that utilizing the trend analysis from the time series data elements reduced the number of classification errors by 73%. The area under the Receiver Operating Characteristic curve increased from a baseline of 87% to 98% by including the trend analysis. In addition to the performance measures, we were also able to demonstrate that adding raw time series data elements without their associated trend analyses improved classification accuracy as compared to the baseline multivariate model, but diminished classification accuracy as compared to when just the trend analysis features were added (ie, without adding the raw time series data elements). We believe this phenomenon was largely attributable to overfitting, which is known to increase as the ratio of candidate features to class examples rises. Furthermore, although we employed several feature reduction strategies to counteract the overfitting problem, they failed to improve the performance beyond that which was achieved by exclusion of the raw time series elements. Finally, our data demonstrated that pulse oximetry and systolic blood pressure readings tend to start diminishing about 10-20 minutes before an arrest, whereas heart rates tend to diminish rapidly less than 5 minutes before an arrest.