Automated segmentation of molecular subunits in electron cryomicroscopy density maps

Electron cryomicroscopy (cryoEM) is capable of imaging large macromolecular machines composed of multiple components. However, it is currently only possible to achieve moderate resolution at which it may be possible to computationally extract the individual components in the machine. In this work, we present application details of an automated method for detecting and segmenting the components of a large machine in an experimentally determined density map. This method is applicable to object with and without symmetry and takes advantage of global and local symmetry axes if present. We have applied this segmentation algorithm to several cryoEM data sets already deposited in EMDB with various complexities, symmetries and resolutions and validated the results using manually segmented density and available structures of the components in the PDB. As such, automated segmentation could become a useful tool for the analysis of the ever-increasing number of structures of macromolecular machines derived from cryoEM.     

Cell biology in neuroscience: Cellular and molecular mechanisms underlying axon formation,growth, and branching

Proper brain wiring during development is pivotal for adult brain function. Neurons display a high degree of polarization both morphologically and functionally, and this polarization requires the segregation of mRNA, proteins, and lipids into the axonal or somatodendritic domains. Recent discoveries have provided insight into many aspects of the cell biology of axonal development including axon specification during neuronal polarization, axon growth, and terminal axon branching during synaptogenesis.

Introduction

Axon development can be divided into three main steps: (1) axon specification during neuronal polarization, (2) axon growth and guidance, and (3) axon branching and presynaptic differentiation (Fig. 1; Barnes and Polleux, 2009; Donahoo and Richards, 2009). These three steps are exemplified during neocortical development in the mouse: upon neurogenesis, newly born neurons engage long-range migration and polarize (Fig. 1, A and B) by adopting a bipolar morphology with a leading and a trailing process (Fig. 1 C). During migration (approximately from embryonic day [E]11 to E18 in the mouse cortex), the trailing process becomes the axon and extends rapidly while being guided to its final destination (lasts until around postnatal day [P]7 in mouse corticofugal axons with distant targets like the spinal cord; Fig. 1, D–F). Finally, upon reaching its target area, extensive axonal branching occurs during the formation of presynaptic contacts with specific postsynaptic partners (during the second and third postnatal week in the mouse cortex; Fig. 1, G–I). Disruption of any of these steps is thought to lead to various neurodevelopmental disorders ranging from mental retardation and infantile epilepsy to autism spectrum disorders (Zoghbi and Bear, 2012). This review will provide an overview of some of the cellular and molecular mechanisms underlying axon specification, growth, and branching.Open in a separate windowFigure 1.Axon specification, growth, and branching during mouse cortical development. Three stages of the development of callosal axons of cortical pyramidal neurons from the superficial layers 2/3 of the somatosensory cortex in the mouse visualized using long-term in utero cortical electroporation. For this class of model axons, development can be divided in three main stages: (1) neurogenesis and axon specification, occurring mostly at embryonic ages (A–C); (2) axon growth/guidance during the first postnatal week (D–F); and (3) axon branching and synapse formation until approximately the end of the third postnatal week (G–I). A, D, and G show coronal sections of mouse cortex at the indicated ages after in utero cortical electroporation of a GFP-coding plasmid at E15.5 in superficial neuron precursors in one brain hemisphere only (GFP signal in inverted color, dotted line indicates the limits of the brain). B, E, and H are a schematic representation of the main morphological changes observed in callosally projecting axons (red) at the corresponding ages. C shows the typical bipolar morphology of a migrating neuron emitting a trailing process (TP) and a leading process (LP) that will ultimately become the axon and dendrite, respectively. F and I show typical axon projections of layer 2/3 neurons located in the primary somatosensory area at P8 and P21, respectively. Neurons and axons in C, F, and I are visualized by GFP expression (inverted color). Image in C is modified from Barnes et al. (2007) with permission from Elsevier. Images in D, F, G, and I are reprinted from Courchet et al. (2013) with permission from Elsevier.

Neuronal polarization and axon specification

Neuronal polarization is the process of breaking symmetry in the newly born cell to create the asymmetry inherent to the formation of the axonal and somatodendritic compartments (Dotti and Banker, 1987). The mechanisms underlying this process have been studied extensively in vitro and more recently in vivo, but the exact sequence of events has remained elusive (Neukirchen and Bradke, 2011) partly because it is studied in various neuronal cell types that might not use the same extrinsic/intrinsic mechanisms to polarize. It is highly likely that at least three factors underlie neuronal polarization: extracellular cues, intracellular signaling cascades, and subcellular organelle localization. The partition-defective proteins (PARs) are a highly conserved family of proteins including two dyads (Par3/Par6 adaptor proteins and the Par4/Par1 serine/threonine kinases) that are required for polarization and axon formation (Shi et al., 2003, 2004; Barnes et al., 2007; Shelly et al., 2007; Chen et al., 2013), while many other intracellular signaling molecules also support axon formation (Oliva et al., 2006; Rašin et al., 2007; Barnes and Polleux, 2009; Shelly et al., 2010; Cheng et al., 2011; Hand and Polleux, 2011; Cheng and Poo, 2012; Gärtner et al., 2012). These intracellular signaling pathways are influenced by localized extracellular cues that instruct which neurite becomes the axon by either directly promoting axon extension or repressing axon growth in favor of dendritic growth (Adler et al., 2006; Yi et al., 2010; Randlett et al., 2011b; Shelly et al., 2011).The role of organelle subcellular localization during neuronal polarization is a more controversial issue. Initially, the orientation of organelles, including the Golgi complex, centrosomes, mitochondria, and endosomes, was shown to correlate with the neurite that becomes the axon in vitro (Bradke and Dotti, 1997; de Anda et al., 2005, 2010) and in vivo (de Anda et al., 2010). However, more recent studies suggest that the positioning of the centrosome is not necessary for neuronal polarization (Distel et al., 2010; Nguyen et al., 2011). Centrosome localization is likely constrained by microtubule organization within the cell, and therefore the centrosome position within the cell changes dynamically during different stages of polarization (Sakakibara et al., 2013). The question of how the interplay between extracellular cues, intracellular signaling, and organelle localization lead to polarization has pushed the field to perform more extensive in vivo imaging studies as in vitro systems/models have a difficult time recapitulating the complex environment and rely on neurons that were previously polarized in vivo.Like other epithelial cells, neural progenitors present a high degree of polarization along the apico-basal axis (Götz and Huttner, 2005). One of the major questions still needing to be addressed is how, or if, newly born mammalian neurons inherit some level of asymmetry from their parent progenitors (Barnes and Polleux, 2009). Recent studies have attempted to answer this question in vivo but have found that the answer might vary in each neuronal subtype. Retinal ganglion cells (RGCs), retinal bipolar neurons, and tegmental hindbrain nuclei neurons seem to inherit the apical/basolateral polarity from their progenitors (Morgan et al., 2006; Zolessi et al., 2006; Distel et al., 2010; Randlett et al., 2011a). In cortical neurons, hippocampal neurons, and cerebellar granule neurons, this relationship is unclear, in part because newly born cortical neurons first exhibit a multipolar morphology with dynamic neurites emerging from the cell body before adopting a bipolar morphology, suggesting they may not retain a predisposed parental polarity (Hand et al., 2005; Barnes et al., 2007). Other factors also suggest that different neuronal subtypes use different mechanisms during polarization. One such factor is the position where neurons specify their axon relative to the original apical/basolateral axis of their progenitors. As an example, cortical neurons in the mouse brain protrude an axon from the membrane facing the original apical surface toward the ventricular zone (Hand et al., 2005; Barnes et al., 2007; Shelly et al., 2007), whereas zebrafish RGCs form their axon from the membrane on the basolateral side (Zolessi et al., 2006; Randlett et al., 2011b). Another significant difference between cortical neurons and RGCs is related to the timing of axogenesis and dendrogenesis. RGCs tend to form their axons and dendrites at the same time during migration (Zolessi et al., 2006; Randlett et al., 2011b). However, cortical neurons form a long axon during migration before significant dendrite arborization takes place. These differences in the regulation of polarization and sequence of axon versus dendrite outgrowth may be linked to the localization of extracellular cues relative to the immature neuron during polarization (Yi et al., 2010).

Neuronal polarization, cytoskeletal dynamics, and polarized transport

What exactly makes the axonal compartment distinct from the somatodendritic domain? This can most easily be illustrated by focusing on the cytoskeleton that forms the framework of the developing axon. The cytoskeleton is composed of microtubules, actin filaments, and intermediate filaments (also called neurofilaments) along with their associated binding partners. Microtubules are composed of α- and β-tubulin subunits that polymerize to form a long filament intrinsically polarized by the addition of tubulin subunits to only one side of the growing filament called the plus end, while on the opposite side depolymerization occurs. It was discovered more than thirty years ago that the axon of a neuron contains a very uniform distribution of microtubules with the plus end facing away from the cell body (Heidemann et al., 1981). Through the years this observation was confirmed in many neuron cell types, and it was determined that dendrites do not have this uniform plus-end out network of microtubules (Fig. 2; Baas et al., 1988). Dendrites appear to have a complex array of microtubule orientations that may vary between species and/or neuronal subtypes. Current research shows that proximal dendrites are composed of mainly minus-end out microtubules, whereas more distal dendrites transition from an equal distribution of minus-end out and plus-end out microtubules to mainly plus-end out microtubules (Stone et al., 2008; Yin et al., 2011; Ori-McKenney et al., 2012). The orientation of microtubules matters greatly because it determines the relative contribution of microtubule-dependent motor proteins (kinesins and dyneins), which are the main motor proteins carrying various cargoes within cells and in particular are responsible for long-range transport in very large cells such as neurons. Dynein (a minus end–directed microtubule motor) is known to be responsible for both the transport of microtubules away from the cell body and for the uniform polarity of microtubules in the axon (Ahmad et al., 1998; Zheng et al., 2008). Recently, it was discovered that kinesin-1 (a plus end–directed microtubule motor) is required for the minus-end out orientation of microtubules in the dendrites of Caenorhabditis elegans DA9 neurons through selective transport of plus-end out microtubule fragments out of the dendrite (Yan et al., 2013). Another hallmark that differentiates the axonal and somatodendritic compartments is the microtubule-associated proteins (MAPs) that decorate microtubules to regulate their bundling and stability (Hirokawa et al., 2010). Microtubules in the axon are mainly decorated by Tau and MAP1B, whereas microtubules in the dendrites are labeled by proteins of the MAP2a-c family. The role of Tau in axon elongation remains controversial because early reports (Harada et al., 1994; Tint et al., 1998; Dawson et al., 2001) of Tau knockout alone suggested that axons were unaffected, but this apparent lack of phenotype might originate from the functional redundancy between MAPs as Tau/MAP1b double knockout mice show clear axon growth defects (Takei et al., 2000).Open in a separate windowFigure 2.Polarity maintenance and trafficking of somatodendritic and axonal proteins. Neurons are polarized into two main compartments: the somatodendritic domain and the axon. These domains are characterized by the underlying cytoskeleton and their unique protein signatures. The axonal cytoskeleton is defined by its uniform microtubule orientation where each microtubule is oriented with its plus end away from the cell body, while the dendrites contain a mixture of microtubules oriented in both directions. The proximal axon is characterized by a structure known as the axon initial segment (AIS, see inset). This highly ordered structure creates a diffusion barrier between the axonal compartment and the rest of the cell. F-actin is responsible for the cytoplasmic barrier, while sodium channels anchored by Ankyrin G form the basis of the plasma membrane barrier. Tau is retained in the axon by a microtubule-based filter at the AIS. Molecular motors (including kinesin, dynein, and myosin) then use the underlying cytoskeleton to restrict cargo transport to either the axon (such as Cntn1 and L1) or the dendrites (such as PSD95, AMPARs, and NMDARs).The dynamics of actin polymerization into actin filaments (F-actin) also play an important role in defining the axonal compartment, and contain an intrinsic polarity based on the polymerization of the free G-actin subunits (Hirokawa et al., 2010). Beyond the well-documented early role of F-actin dynamics in neurite outgrowth, multiple groups have shown that the disruption of actin polymerization allows dendritically localized proteins to incorrectly enter the axonal compartment (Winckler et al., 1999; Lewis et al., 2009; Song et al., 2009). The existence of a “diffusion barrier” in the proximal part of newly formed axons (Song et al., 2009) was long suspected. One of the current hypotheses is that a dense F-actin meshwork creates a cytoplasmic diffusion barrier shortly after polarization, which in part separates the axonal compartment from the neuronal cell body (Fig. 2, inset). Based on functional analysis and electron microscopy analysis, this “F-actin–based filter” is oriented so that the plus ends point toward the cell body while the minus ends point into the axon (Lewis et al., 2009, 2011; Watanabe et al., 2012). Two recent papers show via high resolution imaging techniques that indeed the axon has a unique F-actin network that is not found in dendrites (Watanabe et al., 2012; Xu et al., 2013). The development of this F-actin meshwork appears to directly precede the formation of the axon initial segment (AIS; Song et al., 2009; Galiano et al., 2012). An intra-axonal diffusion barrier, composed of Spectrins and Ankyrin B, defines the eventual position of the AIS. This boundary excludes Ankyrin G, which instead clusters in the most proximal part of the axon close to the cell body, where the AIS will form (Galiano et al., 2012). Ankyrin G is required for AIS formation and maintenance, and its loss causes the axon to start forming protrusions resembling dendritic spines (Hedstrom et al., 2008). Microtubules also play an important role at the AIS, as recent evidence suggests that Tau is retained in the axon through a microtubule-based diffusion barrier independently of the F-actin based filter (X. Li et al., 2011). The AIS is important in the formation of a plasma membrane barrier between the axonal and somatodendritic domains and its disruption affects both neuronal polarity and function because it is critical for clustering of voltage-dependent sodium channels and action potential initiation (Rasband, 2010).One of the critical cellular mechanisms underlying neuronal polarization is the polarized transport of various cargoes in axons and dendrites. Transport of proteins and various organelles is performed by the microtubule-dependent motor proteins kinesin and dynein (Hirokawa et al., 2010). Studies from many laboratories have demonstrated that kinesin motors can carry cargo to both the axonal and dendritic compartments (Burack et al., 2000; Nakata and Hirokawa, 2003). The mechanism for how selection occurs is not completely understood, but it probably incorporates both the affinity of the kinesin head for microtubules and the cargo bound to the motor protein (Nakata and Hirokawa, 2003; Song et al., 2009; Jenkins et al., 2012). In the axon, dynein works to bring cargo and retrograde signals back to the cell body, whereas in the dendrites it is responsible for much of the transport from the soma into the dendrites (Zheng et al., 2008; Kapitein et al., 2010; Harrington and Ginty, 2013). Additionally, the F-actin–dependent myosin motors can affect the polarized transport of cargos by using the F-actin–based cytoplasmic filter at the AIS to deny or facilitate entry of vesicles into the axon. Loss of the actin filter or myosin Va activity (a plus end–directed motor) allows dendritic cargos into the axon, whereas myosin VI (a minus end–directed motor) both removes axonal proteins from the dendritic surface and funnels vesicles containing axonal proteins through the actin filter at the AIS (Lewis et al., 2009, 2011; Al-Bassam et al., 2012). A current working hypothesis is that vesicles composed of multiple cargoes contain binding sites for each of these motors, and that through unknown mechanisms the activity of the motors can be differentially regulated to control the directionality of transport. An interesting example of how the interplay between different motors and cargo adaptors could lead to polarized transport was recently described for mitochondria (van Spronsen et al., 2013).

Axon growth

Microtubule dynamics regulate axon growth.

After axon specification, axon growth constitutes the second step of axonal development and is tightly linked to axon guidance toward the proper postsynaptic targets. Axon elongation by the growth cone is the product of two opposite forces (Fig. 3): slow axonal transport and the polymerization of microtubules providing a pushing force from the axon shaft, and the retrograde flow of actin providing a pulling force at the front of the growth cone (Letourneau et al., 1987; Suter and Miller, 2011). Although coordinated actin and microtubule dynamics are required for the proper function of the growth cone, it was reported that agents disrupting the actin cytoskeleton have limited consequences on axon elongation and are rather involved in axon guidance in vitro (Marsh and Letourneau, 1984; Ruthel and Hollenbeck, 2000) and in vivo (Bentley and Toroian-Raymond, 1986). Furthermore, local disruption of actin organization in the growth cone of minor neurites allows them to turn into axons (Bradke and Dotti, 1999; Kunda et al., 2001), indicating that the dense actin network present at the periphery of an immature neuronal cell body and in immature neurites may prevent microtubule protrusion and elongation necessary for axon specification.Open in a separate windowFigure 3.Cytoskeletal changes during axon elongation and branching. Representation of axon elongation and collateral branch formation in a cultured neuron. Axon growth is a discontinuous process, and collateral branches often originate from sites where the growth cone paused (gray dotted line), after it has resumed its progression. Other modalities of branch formation can occur through the formation of filopodia and lamellipodia. Red box shows a magnification of the main growth cone. Microtubules from the axon shaft spread into the central (C) zone. Some microtubules pass through the transition (T) zone, containing F-actin arcs, to explore filopodia from the peripheral (P) zone. Upon the proper stimulation by extracellular guidance cues or growth-promoting cues, microtubules are stabilized and invade the P-zone where they provide a pushing force, which, combined with the traction force from the actin treadmilling, provides the force required for growth cone extension. Green box shows the cytoskeletal changes occurring during collateral branch formation in the axon. Filopodia and lamellipodia are primarily F-actin–based protrusions that get invaded by microtubules, then elongate upon microtubule bundling. At later developmental stages, axon branches are stabilized or retracted (blue box) by mechanisms relying on the access to extracellular neurotrophins and/or neuronal activity and synapse formation.Contrary to actin, microtubule polymerization is required to sustain axon elongation and branching (Letourneau et al., 1987; Baas and Ahmad, 1993). Axonal proteins and cytoskeletal elements are transported along the axon through slow axonal transport by molecular motors (Yabe et al., 1999; Xia et al., 2003). It is still controversial whether tubulin and other cytoskeletal elements are transported in the axon as monomers and/or as polymers (Roy et al., 2000; Terada et al., 2000; Wang et al., 2000; Brown, 2003; Terada, 2003). Nonetheless, disruption of the slow transport of tubulin impairs the pushing force resulting from microtubule polymerization and impairs axon elongation (Suter and Miller, 2011). Therefore, it is not surprising that axon growth is affected in vitro and in vivo by disruption of plus-end microtubule-binding proteins such as APC (Shi et al., 2004; Zhou et al., 2004; Yokota et al., 2009; Chen et al., 2011) or EB1 and EB3 (Zhou et al., 2004; Jiménez-Mateos et al., 2005; Geraldo et al., 2008), microtubule-associated proteins such as MAP1B (Black et al., 1994; Takei et al., 2000; Dajas-Bailador et al., 2012; Tortosa et al., 2013), or proteins regulating microtubule severing and reorganization such as KIF2A (Homma et al., 2003), katanin, and spastin (Karabay et al., 2004; Yu et al., 2005; Wood et al., 2006; Butler et al., 2010).The contribution of microtubule dynamics to axon growth is not limited to growth cone dynamics but also involves axonal transport and polymerization along the axon shaft. Moreover, changing the balance between microtubule stabilization and depolymerization by local application of microtubule stabilizing agents is sufficient to instruct one neurite to grow and adopt an axon fate (Witte et al., 2008). Many kinase pathways converge on Tau and other axonal MAPs to regulate their function by phosphorylation (Morris et al., 2011). Among them, the MARK kinases regulate microtubule stability and axonal transport through phosphorylation of Tau (Drewes et al., 1997; Mandelkow et al., 2004). Interestingly, MARK-related kinases SAD-A/B control axon specification in part through phosphorylation of Tau (Barnes et al., 2007) and have very recently been linked to the growth and branching of the axons of sensory neurons (Lilley et al., 2013). Our work recently demonstrated that another family member related to MARKs and SAD kinases, called NUAK1, controls axon branching of mouse cortical neurons through the regulation of presynaptic mitochondria capture (Courchet et al., 2013). To what extent the regulation of Tau and other MAPs by the MARKs, SADs, and NUAK1 kinases contributes to axon elongation remains to be explored.

Where does axon elongation take place?

Growth cone progression and guidance constitute the main driver of axonal growth during development, but this process is unlikely to account for the totality of axon elongation. This is especially true after the axon has reached its final target but the axon shaft keeps growing in proportion to the rest of the body. One mechanism that may contribute to this “interstitial” form of axon elongation during brain/body size increase (see Fig. 1 for an example during postnatal cortex growth) is axon stretching, a process that can induce axon elongation in vitro (Smith et al., 2001; Pfister et al., 2004; Loverde et al., 2011) and in vivo (Abe et al., 2004). Aside from extreme stretching performed through the application of external forces, stretching could also contribute to the natural elongation of the axon in response to the tension resulting from growth cone progression (Suter and Miller, 2011).Axon elongation requires the addition of new lipids, proteins, cytoskeleton elements, and organelles along the axon. Where does the synthesis and incorporation of new components take place? Polysaccharides and cholesterol synthesis mostly occur in the cell body; however, lipid synthesis and/or incorporation can occur along the axon as well (Posse De Chaves et al., 2000; Hayashi et al., 2004). The growth cone is also a site of endocytosis, membrane recycling, and exocytosis (Kamiguchi and Yoshihara, 2001; Winckler and Yap, 2011; Nakazawa et al., 2012). One of the best studied examples of endocytosis and its role in axon growth and neuronal survival is the retrograde trafficking of TrkA receptor by target-derived nerve growth factor (NGF) in the peripheral nervous system (Harrington and Ginty, 2013).

Axon branching and presynaptic differentiation

Where do axon branches form?

The last step of axon development is terminal branching, which allows a single axon to connect to a broad set of postsynaptic targets. Collateral branches are formed along the axon through two distinct mechanisms: the first modality of branching is through splitting or bifurcation of the growth cone, which is linked to axon guidance and to the capacity of one single neuron to reach two targets that are far apart with one single axon. Growth cone splitting is observed in vivo in various neuron types including cortical neurons (Sato et al., 1994; Bastmeyer and O’Leary, 1996; Dent et al., 1999; Tang and Kalil, 2005), sympathetic neurons (Letourneau et al., 1986), motorneurons (Matheson and Levine, 1999), sensory neurons (Ma and Tessier-Lavigne, 2007), or mushroom body neurons in Drosophila (Wang et al., 2002). The second modality, known as interstitial branching, occurs through the formation of collateral branches directly along the axon shaft. Contrary to growth cone splitting, interstitial branching serves the purpose of raising axon coverage locally in order to define their “presynaptic territory”, and may contribute to increased network connectivity (Portera-Cailliau et al., 2005). Although both mechanisms can occur simultaneously in the same neuron, the relative importance of splitting versus interstitial branching is highly divergent from one neuron type to the other (Bastmeyer and O’Leary, 1996; Matheson and Levine, 1999; Portera-Cailliau et al., 2005).In culture, the axon grows in a non-continuous fashion with frequent growth cone pausing. Time-lapse imaging of sensorimotor neurons revealed that interstitial branching often occurs at the site where the growth cone paused, shortly after it has continued its course (Szebenyi et al., 1998). Accordingly, treatments with neurotrophins that slow the growth cone correlate with increased axon branching (Szebenyi et al., 1998). This suggests that growth cone pausing leaves a “mark” along the axon shaft that might predetermine future sites of branching (Kalil et al., 2000). However, local applications of neurotrophins shows that aside from growth cone pause sites the axon shaft remains competent to form collateral branches upon stimulation by extracellular factors (Gallo and Letourneau, 1998; Szebenyi et al., 2001), through the formation of filopodia or lamellipodia. Similar observations in vivo revealed that cortical axons are highly dynamic during development and form multiple filopodia that are the origin of collateral branches (Bastmeyer and O’Leary, 1996). Lamellipodia can be observed as motile, F-actin–dependent “waves” along the axon in vitro (Ruthel and Banker, 1998) and in vivo (Flynn et al., 2009). Moreover, disruption of microtubule organization impairs lamellipodia formation along the axon and is correlated with decreased axon branching (Dent and Kalil, 2001; Tint et al., 2009).

Cytoskeleton dynamics and axon branch formation.

Regardless of what type of protrusion gives rise to a branch, cytoskeletal reorganization in the nascent branch generally follows a similar sequence (Fig. 3): initially F-actin filament reorganization gives rise to a protrusion (filopodia, lamellipodia), followed by microtubule invasion of this otherwise transient structure to consolidate it, before the mature branch starts elongating through microtubule bundling (Gallo, 2011). Actin filaments accumulate along the axon and form “patches” that serve as nucleators for axon protrusions such as filopodia and lamellipodia (Korobova and Svitkina, 2008; Mingorance-Le Meur and O’Connor, 2009; Ketschek and Gallo, 2010). The mRNA for β-actin and regulators of actin polymerization such as Wave1 or Cortactin accumulate along the axons of sensory neurons and form hot-spots of local translation that are associated to NGF-dependent branching (Spillane et al., 2012; Donnelly et al., 2013). Subsequently, microtubules in the axon shaft undergo fragmentation at branch points as a prelude to branch invasion by microtubules (Yu et al., 1994, 2008; Gallo and Letourneau, 1998; Dent et al., 1999; Hu et al., 2012), a process that may disrupt transport locally to help trap molecules and organelles at branch points. Moreover, severed microtubules are transported into branches, a process required for branch stabilization (Gallo and Letourneau, 1999; Ahmad et al., 2006; Qiang et al., 2010; Hu et al., 2012). Interestingly, it is clear that, just like growth cone–mediated axon elongation, F-actin and microtubule reorganization events are interconnected to sustain axon branching (Dent and Kalil, 2001). As an example, microtubule-severing enzymes can also contribute to actin nucleation and filopodia formation (Hu et al., 2012).

Is axon branching linked to axon elongation?

Like in the growth cone, cytoskeleton reorganization constitutes the backbone of branch formation. It is therefore not surprising that many manipulations of the cytoskeleton affect both axon elongation and branch formation (Homma et al., 2003; Chen et al., 2011). Moreover, conditions that primarily disrupt axon elongation could secondarily disrupt branching by impairing the ability of the nascent branch to grow. However, axon elongation and axon branching can be considered as two separate phenomena and can be operationally separated because conditions disrupting one do not systematically affect the other. As an example, the microtubule-severing proteins katanin and spastin have differential consequences on axon elongation (primarily dependent upon katanin function) and branching (mostly spastin mediated; Qiang et al., 2010), taxol treatment (which stabilizes microtubules) affects axon elongation but not branching (Gallo and Letourneau, 1999), and disruption of TrkA endocytosis by knock-down of Unc51-like kinase (ULK1/2) proteins has opposite effects on axon elongation and branching (Zhou et al., 2007). In vivo, superficial layer cortical neurons initially go through a phase of elongation through the corpus callosum without branching (see Fig. 1), then stop elongating and form collateral branches in the contralateral cortex (Mizuno et al., 2007; Wang et al., 2007). It is conceivable that even before myelination, axons are actively prevented from branching at places and stages when they elongate (for example in the white matter of the neocortex) where they tend to be highly fasciculated. The identities of the molecules that inhibit interstitial branching along the axon shaft are currently unknown.

Regulation of axon branching by activity.

Immature neurons display spontaneous activity in the form of calcium waves (Gu et al., 1994; Gomez and Spitzer, 1999; Gomez et al., 2001) and spontaneous vesicular release long before they have completed axon development, which suggested a role for early neuronal activity in axon development and guidance (Catalano and Shatz, 1998). Cell-autonomous silencing of neurons in vivo by transfection of the hyperpolarizing inward-rectifying potassium channel Kir2.1 in olfactory neurons (Yu et al., 2004), in RGCs (Hua et al., 2005) or in cortical pyramidal neurons (Mizuno et al., 2007; Wang et al., 2007), or in vitro through infusion of tetrodotoxin (which blocks action potentials generation) in co-cultures of thalamo-cortical projecting neurons (Uesaka et al., 2007) results in a decrease in terminal axon branching, indicating that synaptic activity is required for axons to fully develop their branching pattern. Moreover, inhibition of synaptic release by expression of tetanus toxin light chain (TeTN-LC; Wang et al., 2007) also abolished terminal axon branching, suggesting that the formation of functional presynaptic release sites is required cell autonomously to control terminal axon branching. However, one potential limitation of the experiments involving TeTN-LC is that it blocks most VAMP-mediated vesicular release (with the exception of VAMP7, also called tetanus toxin–independent VAMP, or TI-VAMP). Therefore TeTN-LC action may not be limited to blocking synaptic vesicle release, but could also inhibit peptide release through vesicles containing neurotrophins for example, or other important trophic factors required for axon branching. More recent experiments through silencing of postsynaptic targets revealed that branching of callosal or thalamocortical axons is also dependent upon the activity of the postsynaptic targets (Mizuno et al., 2010; Yamada et al., 2010), albeit activity of the presynaptic neuron is required earlier during the branching process than activity of the postsynaptic targets (Mizuno et al., 2010). Activity is also required in some neurons at the phase of axon elongation through a feedback loop involving the activity-dependent up-regulation of guidance molecules (Mire et al., 2012).How much does spontaneous or evoked neuronal activity contribute to branching? Reduction of neuronal activity through hyperpolarization induced by overexpression of Kir2.1 significantly reduces axon branching without completely eliminating it (Hua et al., 2005; Mizuno et al., 2007; Wang et al., 2007). Activity seems to serve as a competitive regulator of axon branching with regard to its neighbors because silencing of neighboring axons restores normal branching (Hua et al., 2005). Interestingly, neuronal activity induces neurotrophin expression locally, suggesting that activity can contribute to branching partly through activation of activity-independent branching mechanisms (Calinescu et al., 2011).Neuronal activity can regulate branching through modification of the actin cytoskeleton via RhoA activation (Ohnami et al., 2008), and mRNA accumulates at presynaptic sites, indicating a correlation between local translation and synaptic activity (Lyles et al., 2006; Taylor et al., 2013). Neuronal activity is associated with changes in intracellular Ca²⁺ signaling, which has been shown to play a deterministic function in axon growth (Gomez and Spitzer, 1999). Calcium signaling activates the Ca²⁺/calmodulin-dependent kinases (CAMKs) that are known to regulate axon branching in vitro (Wayman et al., 2004; Ageta-Ishihara et al., 2009) and in vivo (Ageta-Ishihara et al., 2009).

Stabilization and refinement of the axonal arborization.

Axon branches are often formed in excess during development, then later refined to select for specific neural circuits (Luo and O’Leary, 2005). Long-range axon branch retraction has long been observed in layer V cortical neurons that initially project to the midbrain, hindbrain, and spinal cord (O’Leary and Terashima, 1988; Bastmeyer and O’Leary, 1996). At later stages, pyramidal neurons from the primary visual cortex will retract their spinal projection through axon pruning, whereas pyramidal neurons from the primary motor cortex will stabilize this projection but retract their axonal branches growing toward visual targets such as the superior colliculus. The molecular mechanisms controlling this area-specific pattern of axon branch pruning are still poorly understood, but seem to involve extracellular cues such as semaphorins and Rac1-dependent signaling (Bagri et al., 2003; Low et al., 2008; Riccomagno et al., 2012). Another example is the well-characterized refinement of retino-geniculate axons during the selective elimination of binocular input of RGC axon synapses onto relay neurons in the dorsal lateral geniculate nucleus (Muir-Robinson et al., 2002). Interestingly, some axons use caspase-dependent pathways locally to induce the selective retraction of axon branches during the process of pruning (Nikolaev et al., 2009; Simon et al., 2012).Circuit refinement and selective branch retraction can be observed in vivo at the level of the neuromuscular junction where individual branches of motor axons are eliminated asynchronously (Keller-Peck et al., 2001). In the developing CNS, neurotrophin-induced branch retraction can also be observed in a context of competition between neighboring axons (Singh et al., 2008). One other way of stabilizing axon branches is through the formation of synapses with postsynaptic targets. In the visual system, the initial axon arbor is refined to establish ocular dominance through activity-dependent retraction of less active branches (Ruthazer et al., 2003). Time-lapse imaging of RGC axons in zebrafish or in Xenopus tadpole revealed that the formation of presynaptic sites occurs concomitantly to axon branching, and branches that form presynaptic structures are less likely to retract (Meyer and Smith, 2006; Ruthazer et al., 2006). The stabilization of axon branches through formation of synaptic contacts parallels with the stabilization of dendritic branches through synapse formation and stabilization (Niell et al., 2004; J. Li et al., 2011). The role of presynaptic bouton formation goes beyond the stabilization of axonal branches because in vivo, new axon branches can emerge from existing presynaptic terminals (Alsina et al., 2001; Javaherian and Cline, 2005; Panzer et al., 2006).In conclusion, axon growth and branching can be regulated by both activity-dependent and activity-independent mechanisms during development. However, for mammalian CNS axons, much more work is needed to define (1) the precise molecular mechanisms underlying axon branching; (2) the cellular and molecular mechanisms regulating the key transition between axon growth and branching when axons start forming presynaptic contacts with their postsynaptic partners; and (3) the mechanisms regulating axon pruning during synapse elimination.     

Virtual screening, molecular interaction field, molecular dynamics, docking, density functional, and ADMET properties of novel AChE inhibitors in Alzheimer's disease

Alzheimer's disease (AD) affects approximately 10% of the world's population with 65 years of age, being the most common form of dementia in adults and is characterized by senile plaquets and cholinergic deficits. Many drugs currently used for the treatment of the AD are based on the improvement of cholinergic neurotransmission achieved by Acetylcholinesterase (AChE) inhibition, the enzyme responsible for acetylcholine hydrolysis. We have focused in this work on the usage of computer-aided molecular design by virtual screening, molecular dynamics with implicit and explicit water solvation, density functional, molecular interaction field studies, docking procedures, ADMET predictions in order to propose novel potential AChE inhibitor for the treatment of Alzheimer's disease.     

Proton energy relaxation in an electron gas in a uniform magnetic field

Energy relaxation of single protons in an ultracold electron gas in the presence and absence of a uniform magnetic field is studied by the method of molecular dynamics. For the case with a magnetic field, the situation is considered in which the electron Larmor radius is much smaller than the classical minimal approach distance. The calculations are performed for electron densities of 10⁸–10⁹ cm^–3, magnetic fields of В = 1–3 T, electron temperatures of 10–50 K, and proton energies of 100–300 000 K.     

Electric field allowed molecular transitions for one and two photon excitation microscopy

We propose an excitation technique for observing single and two photon excitation in those molecules for which such transitions are forbidden by the selection rules. This is possible by the application of an external electric field that perturbs the molecular orbitals, thereby resulting in a significant shift of energy levels. Such a shift of energy levels may bring those levels in resonance with the radiation field which is normally forbidden by selection rules. Further, parity of the these states may significantly improve the emission process. The external electric field results in the mixing of excited (short lifetime) and metastable states (long lifetime), thus reducing the lifetime of metastable (or near metastable) states. This may provide an effective channel for allowing transition from the metastable states. An application of electric field may result in the excitation of poorly excitable biomolecules. This excitation technique may find applications in single- and multi-photon fluorescence microscopy, bioimaging and optical devices.     

Spatial and temporal characteristics of electron and ion beams generated by high-voltage nanosecond discharges in dense molecular gases

Observations of ordered-structure electron beams generated by high-voltage nanosecond discharges in relatively dense (0.04–1 torr) molecular gases at high overvoltages are reported. The beams of accelerated ions generated by high-voltage nanosecond discharges are found to exhibit the same ordered structure. The observed structure of ion beams casts doubt on the mechanisms for the formation of a striped electron-beam structure that assume a regular ecton distribution. Parameters characterizing the temporal behavior of the accelerated electron and ion beams are measured.     

The electron impact, chemical ionization and fast atom bombardment positive ion mass spectra of 1,2-bis(sulfonyl)methylhydrazines.

A series of bis(sulfonyl)-1-methylhydrazines were analyzed by positive ion electron impact (EI), chemical ionization (CI) and fast atom bombardment (FAB) mass spectrometry. Since these compounds showed activity against the L1210 leukemia, an understanding of their mass spectral behavior is important should the structural characterization of metabolites be required. FAB proved to be the most useful technique, generally providing abundant protonated molecule ion peaks, in contrast to the weak peaks observed with CI (ammonia or isobutane) and the total absence of molecular ion peaks in the EI mass spectra. In addition, utilizing FAB eliminated the problem of thermal decomposition, which was very difficult to control under EI and CI experimental conditions. Fragments observed in FAB and CI mass spectra were consistent with protonation at the methyl-bearing nitrogen. One can locate the R1 and R2 moieties relative to the methyl-bearing nitrogen in FAB and CI by assigning that nitrogen as the site of protonation, with subsequent elimination of R2SO2H.     

FireDock: fast interaction refinement in molecular docking

Here, we present FireDock, an efficient method for the refinement and rescoring of rigid-body docking solutions. The refinement process consists of two main steps: (1) rearrangement of the interface side-chains and (2) adjustment of the relative orientation of the molecules. Our method accounts for the observation that most interface residues that are important in recognition and binding do not change their conformation significantly upon complexation. Allowing full side-chain flexibility, a common procedure in refinement methods, often causes excessive conformational changes. These changes may distort preformed structural signatures, which have been shown to be important for binding recognition. Here, we restrict side-chain movements, and thus manage to reduce the false-positive rate noticeably. In the later stages of our procedure (orientation adjustments and scoring), we smooth the atomic radii. This allows for the minor backbone and side-chain movements and increases the sensitivity of our algorithm. FireDock succeeds in ranking a near-native structure within the top 15 predictions for 83% of the 30 enzyme-inhibitor test cases, and for 78% of the 18 semiunbound antibody-antigen complexes. Our refinement procedure significantly improves the ranking of the rigid-body PatchDock algorithm for these cases. The FireDock program is fully automated. In particular, to our knowledge, FireDock's prediction results are comparable to current state-of-the-art refinement methods while its running time is significantly lower. The method is available at http://bioinfo3d.cs.tau.ac.il/FireDock/.     

On the graphical display of molecular electrostatic force-fields and gradients of the electron density

Gradients of the electrostatic potential and electron density of methyl fluoride are displayed as coloured regions on isopotential surfaces. Shading is used to achieve 3D perception and also to encode the direction of the gradient. Colouring the density by the magnitude of the gradient is demonstrated to provide an intrinsic mechanism for colouring the molecular density, so that the underlying atomic structure is vividly displayed in hues commonly used in molecular ball and stick models. The algorithm for contouring a 4D object and drawing coloured, shaded isopotential surfaces is described.     

Using light to induce energy and electron transfer or molecular motions in multicomponent systems.

Light-induced processes are at the basis of fundamental natural phenomena as well as of a variety of applications. Since the functions that can arise from the interaction between light and matter depend on the degree of complexity and organization of the receiving 'matter', the research on these processes has progressively moved from molecular to supramolecular (multicomponent) systems, thereby originating the field of supramolecular photochemistry. In this context, examples of photochemical molecular devices and machines--that is, multicomponent chemical systems capable to perform specific functions under light stimulation--have been developed. Here we report examples of molecular devices in solution, where light is employed (i) to create and transport electronic energy, (ii) to displace electrons, and (iii) to set molecular components in motion, mimicking the operation of mechanical machines and motors. These studies are of interest not only for increasing the basic understanding and testing of current theoretical treatments of photoinduced processes but also eventually for the growth of nanoscience.     

Estimation of the electron density and radiative energy losses in a calcium plasma source based on an electron cyclotron resonance discharge

The parameters of a calcium plasma source based on an electron cyclotron resonance (ECR) discharge were calculated. The analysis was performed as applied to an ion cyclotron resonance system designed for separation of calcium isotopes. The plasma electrons in the source were heated by gyrotron microwave radiation in the zone of the inhomogeneous magnetic field. It was assumed that, in such a combined trap, the energy of the extraordinary microwave propagating from the high-field side was initially transferred to a small group of resonance electrons. As a result, two electron components with different transverse temperatures—the hot resonance component and the cold nonresonance component—were created in the plasma. The longitudinal temperatures of both components were assumed to be equal. The entire discharge space was divided into a narrow ECR zone, where resonance electrons acquired transverse energy, and the region of the discharge itself, where the gas was ionized. The transverse energy of resonance electrons was calculated by solving the equations for electron motion in an inhomogeneous magnetic field. Using the law of energy conservation and the balance condition for the number of hot electrons entering the discharge zone and cooled due to ionization and elastic collisions, the density of hot electrons was estimated and the dependence of the longitudinal temperature T _e∥ of the main (cold) electron component on the energy fraction β lost for radiation was obtained.     

Bsoft: image and molecular processing in electron microscopy

Software for the processing of electron micrographs in structural biology suffers from incompatibility between different packages, poor definition and choice of conventions, and a lack of coherence in software development. The solution lies in adopting a common philosophy of interaction and conventions between the packages. To understand the choices required to have such common interfaces, I am developing a package called "Bsoft." Its foundations lie in the variety of different image file formats used in electron microscopy-a continually frustrating experience to the user and programmer alike. In Bsoft, this problem is greatly diminished by support for many different formats (including MRC, SPIDER, IMAGIC, SUPRIM, and PIF) and by separating algorithmic issues from image format-specific issues. In addition, I implemented a generalized functionality for reading the tag-base STAR (self-defining text archiving and retrieval) parameter file format as a mechanism to exchanging parameters between different packages. Bsoft is written in highly portable code (tested on several Unix systems and under VMS) and offers a continually growing range of image processing functionality, such as Fourier transformation, cross-correlation, and interpolation. Finally, prerequisites for software collaboration are explored, which include agreements on information exchange and conventions, and tests to evaluate compatibility between packages.     

Peptide studies using a fast atom bombardment high field mass spectrometer and data system. 3--Negative ionization: mass calibration, data acquisition and structural characterization

Under negative ionization conditions, nominal mass calibration of the fast bombardment high field mass spectrometer and data system was accomplished using cesium iodide/glycerol as a reference. Mass calibration at --8 kV accelerating potential extends from m/z 387 to m/z 2170 using xenon fast atoms. Negative xenon FAB mass spectra for human angiotensin I and human gastrin I complement their positive fast atom bombardment spectra. Negative xenon fast atom bombardment spectra of underivatized peptides exhibit molecular proton-abstracted ion envelopes and structurally significant fragment ions. Peptide mixture analysis under negative xenon fast atom bombardment reveals peptide molecular ion envelopes of higher relative intensities than under positive xenon fast atom bombardment.     

Experimental and theoretical studies of the spatial distribution of the ionization rate in a non-self-sustained discharge controlled by a fast electron beam

The spatial distribution of the current density of fast electrons and the ionization rate in a gap filled with atmospheric-pressure air under the conditions of a non-self-sustained discharge controlled by a fast electron beam were investigated. The experiments were carried out in a gas-discharge chamber with a grid electrode arranged in parallel to the exit window of the ionization source. Spatial variations in the current density of fast electrons resulting from the grid were measured. The propagation of the electron beam through the discharge system was simulated numerically by the Monte Carlo method in the so-called “effective collision” approximation. The calculated results agree well with the experimental data.     

Solvent density and long-range dipole field around a DNA-binding protein studied by molecular dynamics

The distribution and orientation of solvent around a DNA-binding protein, 434 Cro, were investigated by molecular dynamics simulations with a periodic-boundary condition. The protein was treated in two states: charged and neutral. The computed high-density sites of the solvent around the protein correlated well with the experimentally determined crystal-water sites, in both the charged and neutral states. A local density map, introduced to investigate the solvent density around the highly mobile regions of the protein, showed a hydration shell around hydrophobic sidechains and hydrogen-bondable sites around hydrophilic sidechains, and also showed that the solvent density is sensitive to the slight concaves of the sidechain surface. The long-range solvent-dipole field was observed around the protein, where the pattern of the dipole ordering was considerably different between the charged and neutral states. A local solvent-dipole field was introduced, and the pattern of the dipole ordering was different between the hydrophobic and hydrophilic sidechains. The dipole field from the charged state provided a higher correlation to the electrostatic field obtained from the Poisson-Boltzmann's equation than that from the neutral state, although the correlation become weak quickly for the both states with increasing the protein-solvent distance.     

Change in electron and spin density upon electron transfer to haem

Haems are the cofactors of cytochromes and important catalysts of biological electron transfer. They are composed of a planar porphyrin structure with iron coordinated at the centre. It is known from spectroscopy that ferric low-spin haem has one unpaired electron at the iron, and that this spin is paired as the haem receives an electron upon reduction (I. Bertini, C. Luchinat, NMR of Paramagnetic Molecules in Biological Systems, Benjamin/Cummins Publ. Co., Menlo Park, CA, 1986, pp. 165-170; H.M. Goff, in: A.B.P. Lever, H.B. Gray (Eds.), Iron Porphyrins, Part I, Addison-Wesley Publ. Co., Reading, MA, 1983, pp. 237-281; G. Palmer, in: A.B.P. Lever, H.B. Gray (Eds.), Iron Porphyrins, Part II, Addison-Wesley Publ. Co., Reading, MA, 1983, pp. 43-88). Here we show by quantum chemical calculations on a haem a model that upon reduction the spin pairing at the iron is accompanied by effective delocalisation of electrons from the iron towards the periphery of the porphyrin ring, including its substituents. The change of charge of the iron atom is only approx. 0.1 electrons, despite the unit difference in formal oxidation state. Extensive charge delocalisation on reduction is important in order for the haem to be accommodated in the low dielectric of a protein, and may have impact on the distance dependence of the rates of electron transfer. The lost individuality of the electron added to the haem on reduction is another example of the importance of quantum mechanical effects in biological systems.     

Electric field-directed growth and branching of cultured frog nerves: Effects of aminoglycosides and polycations

The direction and rate of earliest nerve growth are critical determinants of neuronal architecture. One extrinsic cue that influences these parameters is a small direct current electric field, although the underlying mechanisms are unclear. We have studied the orientation, rate of growth, and branching behavior of embryonic Xenopus spinal neurites exposed to aminoglycoside antibiotics, to raised external cations, to applied direct current electric fields, and to combinations of these treatments. Field-induced cathodal turning and cathodal branching of neurites were blocked by the aminoglycosides, by raised extracellular calcium ([Ca²⁺]₀) and by raised extracellular magnesium ([Mg²⁺]₀). Neomycin together with high external Ca²⁺, by contrast, induced a reversal in the polarity of turning and branching, with neurites reorienting and branching more frequently anodally. Aminoglycosides decreased neurite growth rates, and for neomycin this was partially reversed by high external Ca²⁺. Raised [Ca²⁺]₀ alone but not raised [Mg²⁺]₀ altered growth rates in a field-strength dependent manner. Modulation of membrane surface charge may underlie altered galvanotropic orientation and branching. Such an effect is insufficient to explain the changes in growth rates, which may result from additional perturbations to Ca²⁺ influx and inositol phospholipid metabolism. © 1995 John Wiley & Sons, Inc.