A new platform to study the molecular mechanisms of exocytosis

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.

Bayesian inference for parameter estimation in lactoferrin-mediated iron transport across blood-brain barrier

BACKGROUND

In neurodegenerative diseases such as Alzheimer's and Parkinson's, excessive irons as well as lactoferrin (Lf), but not transferrin (Tf), have been found in and around the affected regions of the brain. These evidences suggest that lactoferrin plays a critical role during neurodegenerative diseases, although Lf-mediated iron transport across blood-brain barrier (BBB) is negligible compared to that of transferrin in normal condition. However, the kinetics of lactoferrins and lactoferrin-mediated iron transport are still unknown.

METHOD

To determine the kinetic rate constants of lactoferrin-mediated iron transport through BBB, a mass-action based ordinary differential equation model has been presented. A Bayesian framework is developed to estimate the kinetic rate parameters from posterior probability density functions. The iron transport across BBB is studied by considering both Lf- and Tf-mediated pathways for both normal and pathologic conditions.

RESULTS

Using the point estimates of kinetic parameters, our model can effectively reproduce the experimental data of iron transport through BBB endothelial cells. The robustness of the model and parameter estimation process are further verified by perturbation of kinetic parameters. Our results show that surge in high-affinity receptor density increases lactoferrin as well as iron in the brain.

CONCLUSIONS

Due to the lack of a feedback loop such as iron regulatory proteins (IRPs) for lactoferrin, iron can transport to the brain continuously, which might increase brain iron to pathological levels and can contribute to neurodegeneration.

GENERAL SIGNIFICANCE

This study provides an improved understanding of presence of lactoferrin and iron in the brain during neurodegenerative diseases.

Iron transport kinetics through blood-brain barrier endothelial cells

BACKGROUND

Transferrin and its receptors play an important role during the uptake and transcytosis of iron through blood-brain barrier (BBB) endothelial cells (ECs) to maintain iron homeostasis in BBB endothelium and brain. Any disruptions in the cell environment may change the distribution of transferrin receptors on the cell surface, which eventually alter the homeostasis and initiate neurodegenerative disorders. In this paper, we developed a comprehensive mathematical model that considers the necessary kinetics for holo-transferrin internalization and acidification, apo-transferrin recycling, and exocytosis of free iron and transferrin-bound iron through basolateral side of BBB ECs.

METHODS

Ordinary differential equations are formulated based on the first order reaction kinetics to model the iron transport considering their interactions with transferrin and transferrin receptors. Unknown kinetics rate constants are determined from experimental data by applying a non-linear optimization technique.

RESULTS

Using the estimated kinetic rate constants, the presented model can effectively reproduce the experimental data of iron transports through BBB ECs for many in-vitro studies. Model results also suggest that the BBB ECs can regulate the extent of the two possible iron transport pathways (free and transferrin-bound iron) by controlling the receptor expression, internalization of holo-transferrin-receptor complexes and acidification of holo-transferrin inside the cell endosomes.

CONCLUSION

The comprehensive mathematical model described here can predict the iron transport through BBB ECs considering various possible routes from blood side to brain side. The model can also predict the transferrin and iron transport behavior in iron-enriched and iron-depleted cells, which has not been addressed in previous work.

Dynamical Organization of Syntaxin-1A at the Presynaptic Active Zone

Synaptic vesicle fusion is mediated by SNARE proteins forming in between synaptic vesicle (v-SNARE) and plasma membrane (t-SNARE), one of which is Syntaxin-1A. Although exocytosis mainly occurs at active zones, Syntaxin-1A appears to cover the entire neuronal membrane. By using STED super-resolution light microscopy and image analysis of Drosophila neuro-muscular junctions, we show that Syntaxin-1A clusters are more abundant and have an increased size at active zones. A computational particle-based model of syntaxin cluster formation and dynamics is developed. The model is parametrized to reproduce Syntaxin cluster-size distributions found by STED analysis, and successfully reproduces existing FRAP results. The model shows that the neuronal membrane is adjusted in a way to strike a balance between having most syntaxins stored in large clusters, while still keeping a mobile fraction of syntaxins free or in small clusters that can efficiently search the membrane or be traded between clusters. This balance is subtle and can be shifted toward almost no clustering and almost complete clustering by modifying the syntaxin interaction energy on the order of only 1 kBT. This capability appears to be exploited at active zones. The larger active-zone syntaxin clusters are more stable and provide regions of high docking and fusion capability, whereas the smaller clusters outside may serve as flexible reserve pool or sites of spontaneous ectopic release.

A computational analysis framework for molecular cell dynamics: case-study of exocytosis

One difficulty in conducting biologically meaningful dynamic analysis at the systems biology level is that in vivo system regulation is complex. Meanwhile, many kinetic rates are unknown, making global system analysis intractable in practice. In this article, we demonstrate a computational pipeline to help solve this problem, using the exocytotic process as an example. Exocytosis is an essential process in all eukaryotic cells that allows communication in cells through vesicles that contain a wide range of intracellular molecules. During this process a set of proteins called SNAREs acts as an engine in this vesicle-membrane fusion, by forming four-helical bundle complex between (membrane) target-specific and vesicle-specific SNAREs. As expected, the regulatory network for exocytosis is very complex. Based on the current understanding of the protein-protein interaction network related to exocytosis, we mathematically formulated the whole system, by the ordinary differential equations (ODE). We then applied a mathematical approach (called inverse problem) to estimating the kinetic parameters in the fundamental subsystem (without regulation) from limited in vitro experimental data, which fit well with the reports by the conventional assay. These estimates allowed us to conduct an efficient stability analysis under a specified parameter space for the exocytotic process with or without regulation. Finally, we discuss the potential of this approach to explain experimental observations and to make testable hypotheses for further experimentation.

Super-resolution imaging reveals the internal architecture of nano-sized syntaxin clusters

Key synaptic proteins from the soluble SNARE (N-ethylmaleimide-sensitive factor attachment protein receptor) family, among many others, are organized at the plasma membrane of cells as clusters containing dozens to hundreds of protein copies. However, the exact membranal distribution of proteins into clusters or as single molecules, the organization of molecules inside the clusters, and the clustering mechanisms are unclear due to limitations of the imaging and analytical tools. Focusing on syntaxin 1 and SNAP-25, we implemented direct stochastic optical reconstruction microscopy together with quantitative clustering algorithms to demonstrate a novel approach to explore the distribution of clustered and nonclustered molecules at the membrane of PC12 cells with single-molecule precision. Direct stochastic optical reconstruction microscopy images reveal, for the first time, solitary syntaxin/SNAP-25 molecules and small clusters as well as larger clusters. The nonclustered syntaxin or SNAP-25 molecules are mostly concentrated in areas adjacent to their own clusters. In the clusters, the density of the molecules gradually decreases from the dense cluster core to the periphery. We further detected large clusters that contain several density gradients. This suggests that some of the clusters are formed by unification of several clusters that preserve their original organization or reorganize into a single unit. Although syntaxin and SNAP-25 share some common distributional features, their clusters differ markedly from each other. SNAP-25 clusters are significantly larger, more elliptical, and less dense. Finally, this study establishes methodological tools for the analysis of single-molecule-based super-resolution imaging data and paves the way for revealing new levels of membranal protein organization.

Network motif comparison rationalizes Sec1/Munc18-SNARE regulation mechanism in exocytosis

BACKGROUND

Network motifs, recurring subnetwork patterns, provide significant insight into the biological networks which are believed to govern cellular processes.

METHODS

We present a comparative network motif experimental approach, which helps to explain complex biological phenomena and increases the understanding of biological functions at the molecular level by exploring evolutionary design principles of network motifs.

RESULTS

Using this framework to analyze the SM (Sec1/Munc18)-SNARE (N-ethylmaleimide-sensitive factor activating protein receptor) system in exocytic membrane fusion in yeast and neurons, we find that the SM-SNARE network motifs of yeast and neurons show distinct dynamical behaviors. We identify the closed binding mode of neuronal SM (Munc18-1) and SNARE (syntaxin-1) as the key factor leading to mechanistic divergence of membrane fusion systems in yeast and neurons. We also predict that it underlies the conflicting observations in SM overexpression experiments. Furthermore, hypothesis-driven lipid mixing assays validated the prediction.

CONCLUSION

Therefore this study provides a new method to solve the discrepancies and to generalize the functional role of SM proteins.

Friends and foes in synaptic transmission: the role of tomosyn in vesicle priming

Priming is the process by which vesicles become available for fusion at nerve terminals and is modulated by numerous proteins and second messengers. One of the prominent members of this diverse family is tomosyn. Tomosyn has been identified as a syntaxin-binding protein; it inhibits vesicle priming, but its mode of action is not fully understood. The inhibitory activity of tomosyn depends on its N-terminal WD40-repeat domain and is regulated by the binding of its SNARE motif to syntaxin. Here, we describe new physiological information on the function of tomosyn and address possible interpretations of these results in the framework of the recently described crystal structure of the yeast tomosyn homolog Sro7. We also present possible molecular scenarios for vesicle priming and the involvement of tomosyn in these processes.

Vesicle priming and recruitment by ubMunc13-2 are differentially regulated by calcium and calmodulin

Ca2+ regulates multiple processes in nerve terminals, including synaptic vesicle recruitment, priming, and fusion. Munc13s, the mammalian homologs of Caenorhabditis elegans Unc13, are essential vesicle-priming proteins and contain multiple regulatory domains that bind second messengers such as diacylglycerol and Ca2+/calmodulin (Ca2+/CaM). Binding of Ca2+/CaM is necessary for the regulatory effect that allows Munc13-1 and ubMunc13-2 to promote short-term synaptic plasticity. However, the relative contributions of Ca2+ and Ca2+/CaM to vesicle priming and recruitment by Munc13 are not known. Here, we investigated the effect of Ca2+/CaM binding on ubMunc13-2 activity in chromaffin cells via membrane-capacitance measurements and a detailed simulation of the exocytotic machinery. Stimulating secretion under various basal Ca2+ concentrations from cells overexpressing either ubMunc13-2 or a ubMunc13-2 mutant deficient in CaM binding enabled a distinction between the effects of Ca2+ and Ca2+/CaM. We show that vesicle priming by ubMunc13-2 is Ca2+ dependent but independent of CaM binding to ubMunc13-2. However, Ca2+/CaM binding to ubMunc13-2 specifically promotes vesicle recruitment during ongoing stimulation. Based on the experimental data and our simulation, we propose that ubMunc13-2 is activated by two Ca2+-dependent processes: a slow activation mode operating at low Ca2+ concentrations, in which ubMunc13-2 acts as a priming switch, and a fast mode at high Ca2+ concentrations, in which ubMunc13-2 is activated in a Ca2+/CaM-dependent manner and accelerates vesicle recruitment and maturation during stimulation. These different Ca2+ activation steps determine the kinetic properties of exocytosis and vesicle recruitment and can thus alter plasticity and efficacy of transmitter release.

Systematic search for the rate constants that control the exocytotic process from chromaffin cells by a Genetic Algorithm

We have recently created a kinetic model that reproduces the dynamics of exocytosis with high accuracy. The reconstruction necessitated a search, in a multi-dimensional parameter space, for 37 parameters that described the system, with no assurance that the parameters, which reconstructed the observations, are a unique set. In the present study, a Genetic Algorithm (GA) was used for a thorough search in the unknown parameter space, using a strategy of gradual increase of the complexity of the analyzed input data. Upon systematic incorporation of one to four measurable parameters, used as input signals for the analysis, the constraint set on the GA search imposed the convergence of the free parameters into a single narrow range. The mean values for each adjustable parameter represent a minimum for the fitness function in the multi-dimensional parameter space. The GA search demonstrates that the parameters that control the kinetics of exocytosis are the rate constants of the steps downstream to synaptotagmin binding, and that the equilibrium constant of the binding of calcium to Munc13 controls the calcium-dependent priming process. Thus, the systematic use of the GA creates a link between specific reactions in the process of exocytosis and experimental phenotypes.

Somatostatin receptor type 5 modulates somatostatin receptor type 2 regulation of adrenocorticotropin secretion

Somatostatin inhibits adrenocorticotropin (ACTH) secretion from pituitary tumor cells. To assess the contribution of somatostatin receptor subtype 5 (SST5) to somatostatin receptor subtype 2 (SST2) action in these cells, we assessed multipathway responses to novel highly monoreceptor-selective peptide agonists and multireceptor agonists, including octreotide and somatostatin-28. Octreotide and somatostatin-28 cell membrane binding affinities correlated with their respective SST2-selective peptide ligand. Although octreotide had similar inhibiting potency (picomolar) for cAMP accumulation and ACTH secretion as an SST2-selective agonist, somatostatin-28 exhibited a higher potency (femtomolar). Baseline spontaneous calcium oscillations assessed by fluorescent confocal microscopy revealed two distinct effects: SST2 activation reduced oscillations at femtomolar concentrations reflected by high inhibiting potency of averaged normalized oscillation amplitude, whereas SST5 activation induces brief oscillation pauses and increased oscillation amplitude. Octreotide exhibits an integrated effect of both receptors; however, somatostatin-28 exhibited a complex response with two separate inhibitory potencies. SST2 internalization was visualized with SST2-selective agonist at lower concentrations than for octreotide or somatostatin-28, whereas SST5 did not internalize. Using monoreceptor-selective peptide agonists, the results indicate that, in AtT-20 cells, SST5 regulates the dominant SST2 action, attenuating SST2 effects on intracellular calcium oscillation and internalization. This may explain superior somatostatin-28 potency and provides a rationale for somatostatin ligand design to treat ACTH-secreting pituitary tumors.