Models of neocortical layer 5b pyramidal cells capturing a wide range of dendritic and perisomatic active properties

The thick-tufted layer 5b pyramidal cell extends its dendritic tree to all six layers of the mammalian neocortex and serves as a major building block for the cortical column. L5b pyramidal cells have been the subject of extensive experimental and modeling studies, yet conductance-based models of these cells that faithfully reproduce both their perisomatic Na(+)-spiking behavior as well as key dendritic active properties, including Ca(2+) spikes and back-propagating action potentials, are still lacking. Based on a large body of experimental recordings from both the soma and dendrites of L5b pyramidal cells in adult rats, we characterized key features of the somatic and dendritic firing and quantified their statistics. We used these features to constrain the density of a set of ion channels over the soma and dendritic surface via multi-objective optimization with an evolutionary algorithm, thus generating a set of detailed conductance-based models that faithfully replicate the back-propagating action potential activated Ca(2+) spike firing and the perisomatic firing response to current steps, as well as the experimental variability of the properties. Furthermore, we show a useful way to analyze model parameters with our sets of models, which enabled us to identify some of the mechanisms responsible for the dynamic properties of L5b pyramidal cells as well as mechanisms that are sensitive to morphological changes. This automated framework can be used to develop a database of faithful models for other neuron types. The models we present provide several experimentally-testable predictions and can serve as a powerful tool for theoretical investigations of the contribution of single-cell dynamics to network activity and its computational capabilities.     

The generation of phase differences and frequency changes in a network model of inferior olive subthreshold oscillations

It is commonly accepted that the Inferior Olive (IO) provides a timing signal to the cerebellum. Stable subthreshold oscillations in the IO can facilitate accurate timing by phase-locking spikes to the peaks of the oscillation. Several theoretical models accounting for the synchronized subthreshold oscillations have been proposed, however, two experimental observations remain an enigma. The first is the observation of frequent alterations in the frequency of the oscillations. The second is the observation of constant phase differences between simultaneously recorded neurons. In order to account for these two observations we constructed a canonical network model based on anatomical and physiological data from the IO. The constructed network is characterized by clustering of neurons with similar conductance densities, and by electrical coupling between neurons. Neurons inside a cluster are densely connected with weak strengths, while neurons belonging to different clusters are sparsely connected with stronger connections. We found that this type of network can robustly display stable subthreshold oscillations. The overall frequency of the network changes with the strength of the inter-cluster connections, and phase differences occur between neurons of different clusters. Moreover, the phase differences provide a mechanistic explanation for the experimentally observed propagating waves of activity in the IO. We conclude that the architecture of the network of electrically coupled neurons in combination with modulation of the inter-cluster coupling strengths can account for the experimentally observed frequency changes and the phase differences.     

Sustained induction of ERK,protein kinase B,and p70 S6 kinase regulates cell spreading and formation of F-actin microspikes upon ligation of integrins by galectin-8, a mammalian lectin

Galectin-8, a member of the galectin family of mammalian lectins, is a secreted protein that promotes cell adhesion and migration upon binding to a subset of integrins through sugar-protein interactions. Ligation of integrins by galectin-8 triggers a distinct pattern of cytoskeletal organization, including formation of F-actin-containing microspikes. This is associated with activation of integrin-mediated signaling cascades (ERK and phosphatidylinositol 3 kinase (PI3K)) that are much more robust and are of longer duration than those induced upon cell adhesion to fibronectin. Indeed, formation of microspikes is enhanced 40% in cells that overexpress protein kinase B, the downstream effector of PI3K. Inhibition of PI3K activity induced by wortmannin partially inhibits cell adhesion and spreading while largely inhibiting microspike formation in cells adherent to galectin-8. Furthermore, the inhibitory effects of wortmannin are markedly accentuated in cells overexpressing PKB or p70S6K (CHO(PKB) and CHO(p70S6K) cells), whose adhesion and spreading on galectin-8 (but not on fibronectin) is inhibited approximately 25-35% in the presence of wortmannin. The above results suggest that galectin-8 is an extracellular matrix protein that triggers a unique repertoire of integrin-mediated signals, which leads to a distinctive cytoskeletal organization and microspike formation. They further suggest that downstream effectors of PI3K, including PKB and p70 S6 kinase, in part mediate cell adhesion, spreading, and microspike formation induced by galectin-8.     

Replication and/or separation of some human telomeres is delayed beyond S-phase in pre-senescent cells

Cultured primary human cells, which lack telomerase, enter a state of replicative senescence after a characteristic number of population doublings. During this process telomeres shorten to a critical length of approximately 5-7 kb. The mechanistic relationship between advanced cell passage, cellular senescence and telomeric function has yet to be fully elucidated. In the study described here, we investigated the relationship between changes in telomeric replication timing and/or sister chromatid separation at telomeric regions and advanced cell passage. Using fluorescence in situ hybridization, we analyzed the appearance of double hybridization signals (doublets), which indicate that the region of interest has replicated and the replicated products have separated sufficiently to be resolved as two distinct signals. The results showed that the replication and separation of several telomeric regions occurs during the second half of S-phase and that a delay in replication and/or separation of sister chromatids at these regions occurs in pre-senescent human fibroblasts. Surprisingly, in a significant percentage of pre-senescent cells, several telomeric regions did not hybridize as doublets even in metaphase chromosomes. This delay was not associated with extensive changes in methylation levels at subtelomeric regions and was circumvented in human fibroblasts expressing ectopic telomerase. We propose that incomplete replication and/or separation of telomeric regions in metaphase may be associated with proliferative arrest of senescent cells. This cell growth arrest may result from the activation of a mitotic checkpoint, or from chromosomal instability consequent to progression in the cell cycle despite failure to replicate and/or separate these regions completely.     

Modeling the electrical behavior of anatomically complex neurons using a network analysis program: Passive membrane

We describe the application of a popular and widely available electrical circuit simulation program called SPICE to modeling the electrical behavior of neurons with passive membrane properties and arbitrarily complex dendritic trees. Transient responses may be calculated at any location in the cell model following current, voltage or conductance perturbations at any point. A numbering method is described for binary trees which is helpful in transforming complex dendritic structures into a coded list of short cylindrical dendritic segments suitable for input to SPICE. Individual segments are modeled as isopotential compartments comprised of a parallel resistor and capacitor, representing the transmembrane impedance, in series with one or two core resistors. Synaptic current is modeled by a current source controlled by the local membrane potential and an alpha-shaped voltage, thus simulating a conductance change in series with a driving potential. Extensively branched test cell circuits were constructed which satisfied the equivalent cylinder constraints (Rall 1959). These model neurons were perturbed by independent current sources and by synaptic currents. Responses calculated by SPICE are compared with analytical results. With appropriately chosen model parameters, extremely accurate transient calculations may be obtained. Details of the SPICE circuit elements are presented, along with illustrative examples sufficient to allow implementation of passive nerve cell models on a number of common computers. Methods for modeling excitable membrane are presented in the companion paper (Bunow et al. 1985).     

Premature termination codon at the sterol 27-hydroxylase gene causes cerebrotendinous xanthomatosis in a French family

Cerebrotendinous xanthomatosis (CTX) is an autosomal recessive lipid-storage disease caused by mutations in the sterol 27-hydroxylase gene (CYP27). So far several mutations causing CTX have been identified and characterized. A new mutation creating an insertion of cytosine at position 6 in the cDNA, which is expected to result in a frameshift and a premature termination codon at codon 179, has been identified in a French family. The mutation creates a new site for the restriction endonuclease HaeIII.     

Phylogenetic analysis of Sec7-domain-containing Arf nucleotide exchangers

The eukaryotic family of ADP-ribosylation factor (Arf) GTPases plays a key role in the regulation of protein trafficking, and guanine-nucleotide exchange is crucial for Arf function. Exchange is stimulated by members of another family of proteins characterized by a 200-amino acid Sec7 domain, which alone is sufficient to catalyze exchange on Arf. Here, we analyzed the phylogeny of Sec7-domain-containing proteins in seven model organisms, representing fungi, plants, and animals. The phylogenetic tree has seven main groups, of which two include members from all seven model systems. Three groups are specific for animals, whereas two are specific for fungi. Based on this grouping, we propose a phylogenetically consistent set of names for members of the Sec7-domain family. Each group, except for one, contains proteins with known Arf exchange activity, implying that all members of this family have this activity. Contrary to the current convention, the sensitivity of Arf exchange activity to the inhibitor brefeldin A probably cannot be predicted by group membership. Multiple alignment reveals group-specific domains outside the Sec7 domain and a set of highly conserved amino acids within it. Determination of the importance of these conserved elements in Arf exchange activity and other cellular functions is now possible.     

Hemodynamic Response Imaging: A Potential Tool for the Assessment of Angiogenesis in Brain Tumors

Blood oxygenation level dependence (BOLD) imaging under either hypercapnia or hyperoxia has been used to study neuronal activation and for assessment of various brain pathologies. We evaluated the benefit of a combined protocol of BOLD imaging during both hyperoxic and hypercapnic challenges (termed hemodynamic response imaging (HRI)). Nineteen healthy controls and seven patients with primary brain tumors were included: six with glioblastoma (two newly diagnosed and four with recurrent tumors) and one with atypical-meningioma. Maps of percent signal intensity changes (ΔS) during hyperoxia (carbogen; 95%O2+5%CO2) and hypercapnia (95%air+5%CO2) challenges and vascular reactivity mismatch maps (VRM; voxels that responded to carbogen with reduced/absent response to CO2) were calculated. VRM values were measured in white matter (WM) and gray matter (GM) areas of healthy subjects and used as threshold values in patients. Significantly higher response to carbogen was detected in healthy subjects, compared to hypercapnia, with a GM/WM ratio of 3.8 during both challenges. In patients with newly diagnosed/treatment-naive tumors (n = 3), increased response to carbogen was detected with substantially increased VRM response (compared to threshold values) within and around the tumors. In patients with recurrent tumors, reduced/absent response during both challenges was demonstrated. An additional finding in 2 of 4 patients with recurrent glioblastoma was a negative response during carbogen, distant from tumor location, which may indicate steal effect. In conclusion, the HRI method enables the assessment of blood vessel functionality and reactivity. Reference values from healthy subjects are presented and preliminary results demonstrate the potential of this method to complement perfusion imaging for the detection and follow up of angiogenesis in patients with brain tumors.     

Arp2/3 branched actin network mediates filopodia-like bundles formation in vitro

During cellular migration, regulated actin assembly takes place at the cell leading edge, with continuous disassembly deeper in the cell interior. Actin polymerization at the plasma membrane results in the extension of cellular protrusions in the form of lamellipodia and filopodia. To understand how cells regulate the transformation of lamellipodia into filopodia, and to determine the major factors that control their transition, we studied actin self-assembly in the presence of Arp2/3 complex, WASp-VCA and fascin, the major proteins participating in the assembly of lamellipodia and filopodia. We show that in the early stages of actin polymerization fascin is passive while Arp2/3 mediates the formation of dense and highly branched aster-like networks of actin. Once filaments in the periphery of an aster get long enough, fascin becomes active, linking the filaments into bundles which emanate radially from the aster's surface, resulting in the formation of star-like structures. We show that the number of bundles nucleated per star, as well as their thickness and length, is controlled by the initial concentration of Arp2/3 complex ([Arp2/3]). Specifically, we tested several values of [Arp2/3] and found that for given initial concentrations of actin and fascin, the number of bundles per star, as well as their length and thickness are larger when [Arp2/3] is lower. Our experimental findings can be interpreted and explained using a theoretical scheme which combines Kinetic Monte Carlo simulations for aster growth, with a simple mechanistic model for bundles' formation and growth. According to this model, bundles emerge from the aster's (sparsely branched) surface layer. Bundles begin to form when the bending energy associated with bringing two filaments into contact is compensated by the energetic gain resulting from their fascin linking energy. As time evolves the initially thin and short bundles elongate, thus reducing their bending energy and allowing them to further associate and create thicker bundles, until all actin monomers are consumed. This process is essentially irreversible on the time scale of actin polymerization. Two structural parameters, L, which is proportional to the length of filament tips at the aster periphery and b, the spacing between their origins, dictate the onset of bundling; both depending on [Arp2/3]. Cells may use a similar mechanism to regulate filopodia formation along the cell leading edge. Such a mechanism may allow cells to have control over the localization of filopodia by recruiting specific proteins that regulate filaments length (e.g., Dia2) to specific sites along lamellipodia.