Comparison of electron paramagnetic resonance methods to determine distances between spin labels on human carbonic anhydrase II

Four doubly spin-labeled variants of human carbonic anhydrase II and corresponding singly labeled variants were prepared by site-directed spin labeling. The distances between the spin labels were obtained from continuous-wave electron paramagnetic resonance spectra by analysis of the relative intensity of the half-field transition, Fourier deconvolution of line-shape broadening, and computer simulation of line-shape changes. Distances also were determined by four-pulse double electron-electron resonance. For each variant, at least two methods were applicable and reasonable agreement between methods was obtained. Distances ranged from 7 to 24 A. The doubly spin-labeled samples contained some singly labeled protein due to incomplete labeling. The sensitivity of each of the distance determination methods to the non-interacting component was compared.     

Effect of nucleotides and their analogues on essential light chains in myosin head

The domain movement in myosin head plays a decisive role in the energy transduction process of the muscle contraction. During hydrolysis of ATP, the specific formation of strong binding of myosin head for actin causes conformational changes. As a consequence, the light chain-binding motif generates the powerstroke. In our work maleimide spin labels were covalently attached to Cys-177 residue of ELC in subfragment-1 (S1). Our goal was to study the orientation dependence and the motion of S1, which were incorporated into glycerinated skeletal muscle fibres. The electron paramagnetic resonance spectroscopy (EPR) spectra of the probes depended strongly on the orientation of the fibre axis relative to the magnetic field, indicating that the essential light chain (ELC) and the neck were ordered. The probes were undergoing rapid motion within a cone. The half-width of the cone was estimated to be 65+/-5 degrees (SD, n=8). Addition of ADP affected little the hyperfine splitting and the angular spread of the probe distribution. In the presence of ADP and orthovanadate the intensity of the spectra decreased, which showed the dissociation of S1 and this was accompanied with the disappearance of the orientation dependence.     

Orientation of spin-labeled myosin heads in glycerinated muscle fibers.

We have used electron paramagnetic resonance (EPR) spectra to study spin labels selectively and rigidly attached to myosin heads in glycerinated rabbit psoas muscle fibers. Because the angle between the magnetic field and the principal axis of the probe determines the position of the EPR absorption line, spectra from labeled fibers oriented parallel to the magnetic field yielded directly the distribution of spin label orientations relative to the fiber axis. Two spin labels, having reactivities resembling iodoacetamide (IASL) and maleimide (MSL), were used. In rigor fibers with complete filament overlap, both labels displayed a narrow angular distribution, full width at half maximum approximately 15 degrees, centered at angles of 68 degrees (IASL) and 82 degrees (MSL). Myosin subfragments (heavy meromyosin and subfragment-1) were labeled and allowed to diffuse into fibers. The resulting spectra showed the same sharp angular distribution that was found for the labeled fibers. Thus is appears that virtually all myosin heads in a rigor fiber have the same orientation relative to the fiber axis, and this orientation is determined by the actomyosin bond. Experiments with stretched fibers indicated that the spin labels on the fraction of heads not interacting with actin filaments had a broad angular distribution. Addition of ATP to unstretched fibers under relaxing conditions produced orientational disorder, resulting in a spectrum almost indistinguishable from that of an isotropic distribution of probes. Addition of either an ATP analog (AMPPNP) or pyrophosphate produced partial disorder. That is a fraction of the probes remained sharply oriented as in rigor while a second fraction was in a disordered distribution similar to that of relaxed fibers.     

Simulation vs. reality: a comparison of in silico distance predictions with DEER and FRET measurements

Site specific incorporation of molecular probes such as fluorescent- and nitroxide spin-labels into biomolecules, and subsequent analysis by F?rster resonance energy transfer (FRET) and double electron-electron resonance (DEER) can elucidate the distance and distance-changes between the probes. However, the probes have an intrinsic conformational flexibility due to the linker by which they are conjugated to the biomolecule. This property minimizes the influence of the label side chain on the structure of the target molecule, but complicates the direct correlation of the experimental inter-label distances with the macromolecular structure or changes thereof. Simulation methods that account for the conformational flexibility and orientation of the probe(s) can be helpful in overcoming this problem. We performed distance measurements using FRET and DEER and explored different simulation techniques to predict inter-label distances using the Rpo4/7 stalk module of the M. jannaschii RNA polymerase. This is a suitable model system because it is rigid and a high-resolution X-ray structure is available. The conformations of the fluorescent labels and nitroxide spin labels on Rpo4/7 were modeled using in vacuo molecular dynamics simulations (MD) and a stochastic Monte Carlo sampling approach. For the nitroxide probes we also performed MD simulations with explicit water and carried out a rotamer library analysis. Our results show that the Monte Carlo simulations are in better agreement with experiments than the MD simulations and the rotamer library approach results in plausible distance predictions. Because the latter is the least computationally demanding of the methods we have explored, and is readily available to many researchers, it prevails as the method of choice for the interpretation of DEER distance distributions.     

Effect of ADP on the orientation of spin-labeled myosin heads in muscle fibers: a high-resolution study with deuterated spin labels

We have used electron paramagnetic resonance (EPR) to determine the effects of ADP on the orientational distribution of nitroxide spin labels attached to myosin heads in skinned rabbit psoas muscle fibers. To maximize the specificity of labeling, we spin-labeled isolated myosin heads (subfragment 1) on a single reactive thiol (SH1) and diffused them into unlabeled muscle fibers. To maximize spectral and orientational resolution, we used perdeuterated spin labels, 2H-MSL and 2H-IASL, eliminating superhyperfine broadening and thus narrowing the line widths. Two different spin labels were used, with different orientation relative to the myosin head, to ensure that the results are not affected by unfavorable probe orientation. In rigor, a very narrow three-line spectrum was observed for both spin labels, indicating a narrow orientational distribution, as reported previously (Thomas & Cooke, 1980). ADP induced very slight changes in the spectrum, corresponding to very slight (but significant) changes in the orientational distribution. These changes were quantified by a digital analysis of the spectra, using a two-step simplex fitting procedure (Fajer et al., 1990). First, the magnetic tensor values and line widths were determined by fitting the spectrum of a randomly oriented sample. Then the spectrum of oriented fibers was fit to a model by assuming a Gaussian distribution of the tilt angle (theta) and twist angle (phi) of the nitroxide principal axes relative to the fiber axis. A single-Gaussian distribution resulted in inadequate fits, but a two-component model gave excellent results. ADP induces a small (less than 5 degrees) rotation of the major components for both spin labels, along with a similarly small increase of disorder about the average positions.     

Determination of interspin distances between spin labels attached to insulin: comparison of electron paramagnetic resonance data with the X-ray structure

A method was developed to determine the interspin distances of two or more nitroxide spin labels attached to specific sites in proteins. This method was applied to different conformations of spin-labeled insulins. The electron paramagnetic resonance (EPR) line broadening due to dipolar interaction is determined by fitting simulated EPR powder spectra to experimental data, measured at temperatures below 200 K to freeze the protein motion. The experimental spectra are composed of species with different relative nitroxide orientations and interspin distances because of the flexibility of the spin label side chain and the variety of conformational substates of proteins in frozen solution. Values for the average interspin distance and for the distance distribution width can be determined from the characteristics of the dipolar broadened line shape. The resulting interspin distances determined for crystallized insulins in the R6 and T6 structure agree nicely with structural data obtained by x-ray crystallography and by modeling of the spin-labeled samples. The EPR experiments reveal slight differences between crystal and frozen solution structures of the B-chain amino termini in the R6 and T6 states of hexameric insulins. The study of interspin distances between attached spin labels can be applied to obtain structural information on proteins under conditions where other methods like two-dimensional nuclear magnetic resonance spectroscopy or x-ray crystallography are not applicable.     

Molecular distances from dipolar coupled spin-labels: the global analysis of multifrequency continuous wave electron paramagnetic resonance data.

For immobilized nitroxide spin-labels with a well-defined interprobe geometry, resolved dipolar splittings can be observed in continuous wave electron paramagnetic resonance (CW-EPR) spectra for interelectron distances as large as 30 A using perdeuterated probes. In this work, algorithms are developed for calculating CW-EPR spectra of immobilized, dipolar coupled nitroxides, and then used to define the limits of sensitivity to the interelectron distance as a function of geometry and microwave frequency. Secondly, the CW-EPR spectra of N epsilon-spin-labeled coenzyme NAD+ bound to microcrystalline, tetrameric glyceraldehyde-3-phosphate dehydrogenase (GAPDH) have been collected at 9.8, 34, and 94 GHz. These data have been analyzed, using a combination of simulated annealing and global analysis, to obtain a unique fit to the data. The values of the intermitroxide distance and the five angles defining the relative orientation of the two nitroxides are in reasonable agreement with a molecular model built from the known crystal structure. Finally, the effect of rigid body isotropic rotational diffusion on the CW-EPR spectra of dipolar coupled nitroxides has been investigated using an algorithm based on Brownian dynamics trajectories. These calculations demonstrate the sensitivity of CW-EPR spectra to dipolar coupling in the presence of rigid body rotational diffusion.     

Reappraisal of the electron spin resonance spectra of maleimide and iodoacetamide spin labels in erythrocyte ghosts.

The effects of sulfhydryl inhibitors (iodoacetamide and N-ethylmaleimide) on the electron spin resonance spectra of two maleimide and two iodoacetamide spin labels in erythrocyte ghosts were found to correlate with their relative “lipid”/water partition coefficients. But the spectral characteristics of the maleimide spin labels, and their ghost concentrations after iodoacetamide inhibition, are not consistent with the hypothesis that interprets their spectra solely on the basis of a heterogenous membrane distribution. An alternative hypothesis is suggested which is compatible with relative “lipid solubilities” and the iodoacetamide inhibition spectra.     

De novo high-resolution protein structure determination from sparse spin-labeling EPR data

As many key proteins evade crystallization and remain too large for nuclear magnetic resonance spectroscopy, electron paramagnetic resonance (EPR) spectroscopy combined with site-directed spin labeling offers an alternative approach for obtaining structural information. Such information must be translated into geometric restraints to be used in computer simulations. Here, distances between spin labels are converted into distance ranges between beta carbons by using a "motion-on-a-cone" model, and a linear-correlation model links spin-label accessibility to the number of neighboring residues. This approach was tested on T4-lysozyme and alphaA-crystallin with the de novo structure prediction algorithm Rosetta. The results demonstrate the feasibility of obtaining highly accurate, atomic-detail models from EPR data by yielding 1.0 A and 2.6 A full-atom models, respectively. Distance restraints between amino acids far apart in sequence but close in space are most valuable for structure determination. The approach can be extended to other experimental techniques such as fluorescence spectroscopy, substituted cysteine accessibility method, or mutational studies.     

Calculation of electron paramagnetic resonance spectra from Brownian dynamics trajectories: application to nitroxide side chains in proteins.

We present a method to simulate electron paramagnetic resonance spectra of spin-labeled proteins that explicitly includes the protein structure in the vicinity of the attached spin label. The method is applied to a spin-labeled polyleucine alpha-helix trimer. From short (6 ns) stochastic dynamics simulations of this trimer, an effective potential energy function is calculated. Interaction with secondary and tertiary structures determine the reorientational motion of the spin label side chains. After reduction to a single particle problem, long stochastic dynamic trajectories (700 ns) of the spin label side-chain reorientation are calculated from which the Lamor frequency trajectory and subsequently the electron paramagnetic resonance spectrum is determined. The simulated spectra agree well with experimental electron paramagnetic resonance spectra of bacteriorhodopsin mutants with spin labels in similar secondary and tertiary environments as in the polyleucine.     

A spin label that binds to myosin heads in muscle fibers with its principal axis parallel to the fiber axis.

We have used an indane-dione spin label (2-[-oxyl-2,2,5,5-tetramethyl-3-pyrrolin-3-yl)methenyl]in dane-1,3-dione), designated InVSL, to study the orientation of myosin heads in bundles of chemically skinned rabbit psoas muscle fibers, with electron paramagnetic resonance (EPR) spectroscopy. After reversible preblocking with 5,5'-dithiobis(2-nitro-benzoic acid) (DTNB), we were able to attach most of the spin label covalently and rigidly to either Cys 707 (SH1) or Cys 697 (SH2) on myosin heads. EPR spectra of labeled fibers contained substantial contributions from both oriented and disordered populations of spin labels. Similar spectra were obtained from fibers decorated with InVSL-labeled myosin heads (subfragment 1), indicating that virtually all the spin labels in labeled fibers are on the myosin head. We specifically labeled SH2 with InVSL after reversible preblocking of the SH1 sites with 1-fluoro-2,4-dinitrobenzene (FDNB), resulting in a spectrum that indicated only disordered spin labels. Therefore, the oriented and disordered populations correspond to labels on SH1 and SH2, respectively. The spectrum of SH2-bound labels was subtracted to produce a spectrum corresponding to SH1-bound labels, which was used for further analysis. For this corrected spectrum, the angle between the fiber axis and the principal axis of the spin label was fitted well by a Gaussian distribution centered at theta o = 11 +/- 1 degree, with a full width at half-maximum of delta theta = 15 +/- 2 degrees. The unique orientation of InVSL, with its principal axis almost parallel to the fiber axis, makes it complementary to spin labels previously studied in this system. This label can provide unambiguous information about axial rotations of myosin heads, since any axial rotation of the head must be reflected in the same axial rotation of the principal axis of the probe, thus changing the hyperfine splitting. Therefore, InVSL-labeled fibers have ideal properties needed for further exploration myosin head orientation and rotational motion in muscle.     

Electron spin resonance of TOAC labeled peptides: folding transitions and high frequency spectroscopy

The unnatural, conformationally constrained nitroxide amino acid TOAC (2,2,6,6-tetramethylpiperidine-1-oxyl-4-amino-4-carboxylic acid) stabilizes helical structure and provides a means for studying rigidly spin labeled peptides by electron spin resonance (ESR). Two new directions in TOAC research are described. The first investigates intermediates formed during alpha-helix unfolding. Double TOAC labeled alpha-helical peptides were unfolded at low temperature in aqueous solution with increasing concentrations of guanidine hydrochloride. Comparison of ESR spectra from two doubly labeled peptides suggests that 3(10)-helix emerges as an intermediate. The second research direction involves the use of high frequency ESR (140 GHz) at low temperature to assess dipolar couplings and, hence, distances between TOAC pairs in a series of 3(10)-helical peptides. Preliminary simulations suggest that high frequency ESR is able to extract correct distances between 6 and 11 A. In addition, the spectra appear to be very sensitive to the relative orientation of the TOAC labels.     

The sensitivity of saturation transfer electron paramagnetic resonance spectra to restricted amplitude uniaxial rotational diffusion.

Computational methods have been developed to model the effects of constrained or restricted amplitude uniaxial rotational diffusion (URD) on saturation transfer electron paramagnetic resonance (ST-EPR) signals observed from nitroxide spin labels. These methods, which have been developed to model the global rotational motion of intrinsic membrane proteins that can interact with the cytoskeleton or other peripheral proteins, are an extension of previous work that described computationally efficient algorithms for calculating ST-EPR spectra for unconstrained URD (Hustedt and Beth, 1995, Biophys. J. 69:1409-1423). Calculations are presented that demonstrate the dependence of the ST-EPR signal (V'(2)) on the width (Delta) of a square-well potential as a function of the microwave frequency, the correlation time for URD, and the orientation of the spin-label with respect to the URD axis. At a correlation time of 10 micros, the V'(2) signal is very sensitive to Delta in the range from 0 to 60 degrees, marginally sensitive from 60 degrees to 90 degrees, and insensitive beyond 90 degrees. Sensitivity to Delta depends on the correlation time for URD with higher sensitivity to large values of Delta at the shorter correlation times, on the microwave frequency, and on the orientation of the spin-label relative to the URD axis. The computational algorithm has been incorporated into a global nonlinear least-squares analysis approach, based upon the Marquardt-Levenberg method (Blackman et al., 2001, Biophys. J. 81:3363-3376). This has permitted determination of the correlation time for URD and the width of the square-well potential by automated fitting of experimental ST-EPR data sets obtained from a spin-labeled membrane protein and provided a new automated method for analysis of data obtained from any system that exhibits restricted amplitude URD.     

Myosin regulatory domain orientation in skeletal muscle fibers: application of novel electron paramagnetic resonance spectral decomposition and molecular modeling methods

Reorientation of the regulatory domain of the myosin head is a feature of all current models of force generation in muscle. We have determined the orientation of the myosin regulatory light chain (RLC) using a spin-label bound rigidly and stereospecifically to the single Cys-154 of a mutant skeletal isoform. Labeled RLC was reconstituted into skeletal muscle fibers using a modified method that results in near-stoichiometric levels of RLC and fully functional muscle. Complex electron paramagnetic resonance spectra obtained in rigor necessitated the development of a novel decomposition technique. The strength of this method is that no specific model for a complex orientational distribution was presumed. The global analysis of a series of spectra, from fibers tilted with respect to the magnetic field, revealed two populations: one well-ordered (+/-15 degrees ) with the spin-label z axis parallel to actin, and a second population with a large distribution (+/-60 degrees ). A lack of order in relaxed or nonoverlap fibers demonstrated that regulatory domain ordering was defined by interaction with actin rather than the thick filament surface. No order was observed in the regulatory domain during isometric contraction, consistent with the substantial reorientation that occurs during force generation. For the first time, spin-label orientation has been interpreted in terms of the orientation of a labeled domain. A Monte Carlo conformational search technique was used to determine the orientation of the spin-label with respect to the protein. This in turn allows determination of the absolute orientation of the regulatory domain with respect to the actin axis. The comparison with the electron microscopy reconstructions verified the accuracy of the method; the electron paramagnetic resonance determined that axial orientation was within 10 degrees of the electron microscopy model.     

Selective spin-labeling of the ribosomal proteins of 70S ribosomes from Escherichia coli.

We have used a series of N-(1-oxyl-2,2,5,5-tetramethyl-3-pyrrolidinyl) maleimide spin labels of different length to label, covalently and selectively, the most reactive sulfhydryl groups of 70S ribosomal proteins of Escherichia coli. Under short periods of labeling (1--2 min), less than two spin labels per ribosome are incorporated and were shown to be distributed mainly on five ribosomal proteins in the following order: S18 greater than S21, L27 greater than S17, and S12. With a long period of labeling (3 h) up to 13 spin labels are attached to the ribosome, and protein S1 is the most labeled. The shape of the electron paramagnetic resonance (epr) signal shows two components with a predominance for the strongly immobilized orientation, and the percentage of these components in each spectra has been evaluated. When the distance between the nitroxide group and the maleimide-attaching group exceeds 6 A (1 A = 0.1 nm) the strongly immobilized orientation disappears. The effect of magnesium ions on these selectively spinlabeled ribosomes shows that the dissociation into subunits does not affect the epr signal, but more spin labels are incorporated into the subunits if labeling is performed under conditions of dissociation.     

A study of direct moxibustion using mathematical methods

Direct moxibustion is an important and widely used treatment method in traditional medical science. The use of a mathematical method to analyse direct moxibustion treatment is necessary and helpful in exploring the new direct moxibustion instruments and their standardisation. Thus, this paper aims to use a mathematical method to study direct moxibustion in skin to demonstrate a direct relationship between direct moxibustion and skin stimuli. In this paper, the transient thermal response of skin layers is analysed to study direct moxibustion using the data got from standardised method to measure the temperature of a burning moxa cone. Numerical simulations based on an appropriate finite element model are developed to predict the heat transfer, thermal damage and thermal stress distribution of barley moxa cones and jujube moxa cones in the skin tissue. The results are verified by the ancient literatures of traditional Chinese medicine and clinical application, and showed that mathematical method can be a good interface between moxa cone and skin tissue providing the numerical value basis for moxibustion.     

Orientation of spin-labeled light chain 2 of myosin heads in muscle fibers

Electron paramagnetic resonance (e.p.r.) spectroscopy has been used to monitor the orientation of spin labels attached rigidly to a reactive SH residue on the light chain 2 (LC2) of myosin heads in muscle fibers. e.p.r. spectra from spin-labeled myosin subfragment-1 (S1), allowed to diffuse into unlabeled rigor (ATP-free) fibers, were roughly approximated by a narrow angular distribution of spin labels centered at 66 degrees relative to the fiber axis, indicating a uniform orientation of S1 bound to actin. On the other hand, spectra from spin-labeled heavy meromyosin (HMM) were roughly approximated by two narrow angular distributions centered at 42 degrees and 66 degrees, suggesting that the LC2 domains of the two HMM heads have different orientations. In contrast to S1 or HMM, the spectra from rigor fibers, in which LC2 of endogenous myosin heads was labeled, showed a random orientation which may be due to distortion imposed by the structure of the filament lattice and the mismatch of the helical periodicities of the thick and thin filaments. However, spectra from the fibers in the presence of ATP analog 5'-adenylyl imidodiphosphate (AMPPNP) were approximated by two narrow angular distributions similar to those obtained with HMM. Thus, AMPPNP may cause the LC2 domain to be less flexible and/or the S2 portion to be more flexible, so as to release the distortion of the LC2 domain and make it return to its natural position. At high ionic strength, AMPPNP disoriented the spin labels as ATP did under relaxing conditions, suggesting that the myosin head is detached from and/or weakly (flexibly) attached to a thin filament.     

A study of direct moxibustion using mathematical methods

Direct moxibustion is an important and widely used treatment method in traditional medical science. The use of a mathematical method to analyse direct moxibustion treatment is necessary and helpful in exploring the new direct moxibustion instruments and their standardisation. Thus, this paper aims to use a mathematical method to study direct moxibustion in skin to demonstrate a direct relationship between direct moxibustion and skin stimuli. In this paper, the transient thermal response of skin layers is analysed to study direct moxibustion using the data got from standardised method to measure the temperature of a burning moxa cone. Numerical simulations based on an appropriate finite element model are developed to predict the heat transfer, thermal damage and thermal stress distribution of barley moxa cones and jujube moxa cones in the skin tissue. The results are verified by the ancient literatures of traditional Chinese medicine and clinical application, and showed that mathematical method can be a good interface between moxa cone and skin tissue providing the numerical value basis for moxibustion.     

Rapid Simulation of Unprocessed DEER Decay Data for Protein Fold Prediction

Despite advances in sampling and scoring strategies, Monte Carlo modeling methods still struggle to accurately predict de novo the structures of large proteins, membrane proteins, or proteins of complex topologies. Previous approaches have addressed these shortcomings by leveraging sparse distance data gathered using site-directed spin labeling and electron paramagnetic resonance spectroscopy to improve protein structure prediction and refinement outcomes. However, existing computational implementations entail compromises between coarse-grained models of the spin label that lower the resolution and explicit models that lead to resource-intense simulations. These methods are further limited by their reliance on distance distributions, which are calculated from a primary refocused echo decay signal and contain uncertainties that may require manual refinement. Here, we addressed these challenges by developing RosettaDEER, a scoring method within the Rosetta software suite capable of simulating double electron-electron resonance spectroscopy decay traces and distance distributions between spin labels fast enough to fold proteins de novo. We demonstrate that the accuracy of resulting distance distributions match or exceed those generated by more computationally intensive methods. Moreover, decay traces generated from these distributions recapitulate intermolecular background coupling parameters even when the time window of data collection is truncated. As a result, RosettaDEER can discriminate between poorly folded and native-like models by using decay traces that cannot be accurately converted into distance distributions using regularized fitting approaches. Finally, using two challenging test cases, we demonstrate that RosettaDEER leverages these experimental data for protein fold prediction more effectively than previous methods. These benchmarking results confirm that RosettaDEER can effectively leverage sparse experimental data for a wide array of modeling applications built into the Rosetta software suite.     

Distinct calcium signaling within neuronal growth cones and filopodia

Previous findings indicate that spatial restriction of intracellular calcium levels within growth cones can regulate growth cone behavior at many levels, ranging from filopodial disposition to neurite extension. By combining techniques for focal stimulation of growth cones with those for measurement of filopodia and for capturing low intensity calcium signals, we demonstrate that filopodia on individual growth cones can respond to imposed stimuli independently from one another. Moreover, filopodia and their parent growth cones appear to represent functionally and morphologically distinct domains of calcium regulation, possessing distinct calcium sources and sinks. Both are sensitive to calcium influx; however, application of the calcium ionophore A23187 to cells in calcium-free medium demonstrated the presence of potential intracellular calcium pools in the growth cone proper, but not in isolated filopodia. Thapsigargin significantly reduced the rise in growth cone calcium levels associated with excitatory neurotransmitters, further implicating release from calcium pools as one component of growth cone calcium regulation. The relative contributions of these pools were examined in response to excitatory neurotransmitters by quantitative calcium measurements made in both growth cones and isolated filopodia. Striking differences were observed; filopodia were sensitive to a low concentration of dopamine and serotonin, while growth cones displayed an amplified rise at a higher concentration. The spatial distribution of organelles that could serve as morphological correlates to such calcium amplification was examined using confocal microscopy. While the majority of organelles were located in the central core of the growth cone proper, peripheral organelles were detected at the base of a subset of filopodia. The distinctive distribution of calcium regulation within motile growth cones suggests one mechanism by which growth cones may regulate their complex behavior. © 1996 John Wiley & Sons, Inc.