Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Thermodynamics of Antimicrobial Lipopeptide Binding to Membranes: Origins of Affinity and Selectivity

The development of novel antibiotic drugs is one of the most pressing biomedical problems due to the increasing number of antibiotic-resistant pathogens. Antimicrobial peptides and lipopeptides are a promising category of candidates, but the molecular origins of their antimembrane activity is unclear. Here we explore a series of recently developed antimicrobial lipopeptides, using coarse-grained molecular-dynamics simulations and free energy methods to uncover the thermodynamics governing their binding to membranes. Specifically, we quantify C16-KGGK’s binding affinity to the two types of membrane by umbrella sampling. We also examined the origin of C16-KGGK’s selectivity for bacterial versus mammalian membranes by systematically varying the peptide sequence and salt concentration. Our data showed that the C16 hydrophobic tail is the main contributor to its affinity to lipid membrane, whereas the peptide portion is mainly responsible for its selectivity. Furthermore, the electrostatic interaction between the cationic peptide and anionic bacterial membrane plays a significant role in the selectivity.     

Driving force of water entry into hydrophobic channels of carbon nanotubes: entropy or energy?

Spontaneous entry of water molecules inside single-wall carbon nanotubes (SWCNTs) has been confirmed by both simulations and experiments. Using molecular dynamics simulations, we have studied the thermodynamics of filling of a (6,6) carbon nanotube in a temperature range from 273 to 353 K and with different strengths of the nanotube–water interaction. From explicit energy and entropy calculations using the two-phase thermodynamics method, we have presented a thermodynamic understanding of the filling behaviour of a nanotube. We show that both the energy and the entropy of transfer decrease with increasing temperature. On the other hand, scaling down the attractive part of the carbon–oxygen interaction results in increased energy of transfer while the entropy of transfer increases slowly with decreasing the interaction strength. Our results indicate that both energy and entropy favour water entry into (6,6) SWCNTs. Our results are compared with those of several recent studies of water entry into carbon nanotubes.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Quantifying the Entropy of Binding for Water Molecules in Protein Cavities by Computing Correlations

Protein structural analysis demonstrates that water molecules are commonly found in the internal cavities of proteins. Analysis of experimental data on the entropies of inorganic crystals suggests that the entropic cost of transferring such a water molecule to a protein cavity will not typically be greater than 7.0 cal/mol/K per water molecule, corresponding to a contribution of approximately +2.0 kcal/mol to the free energy. In this study, we employ the statistical mechanical method of inhomogeneous fluid solvation theory to quantify the enthalpic and entropic contributions of individual water molecules in 19 protein cavities across five different proteins. We utilize information theory to develop a rigorous estimate of the total two-particle entropy, yielding a complete framework to calculate hydration free energies. We show that predictions from inhomogeneous fluid solvation theory are in excellent agreement with predictions from free energy perturbation (FEP) and that these predictions are consistent with experimental estimates. However, the results suggest that water molecules in protein cavities containing charged residues may be subject to entropy changes that contribute more than +2.0 kcal/mol to the free energy. In all cases, these unfavorable entropy changes are predicted to be dominated by highly favorable enthalpy changes. These findings are relevant to the study of bridging water molecules at protein-protein interfaces as well as in complexes with cognate ligands and small-molecule inhibitors.     

High affinity antigen recognition of the dual specific variants of herceptin is entropy-driven in spite of structural plasticity

The antigen-binding site of Herceptin, an anti-human Epidermal Growth Factor Receptor 2 (HER2) antibody, was engineered to add a second specificity toward Vascular Endothelial Growth Factor (VEGF) to create a high affinity two-in-one antibody bH1. Crystal structures of bH1 in complex with either antigen showed that, in comparison to Herceptin, this antibody exhibited greater conformational variability, also called "structural plasticity". Here, we analyzed the biophysical and thermodynamic properties of the dual specific variants of Herceptin to understand how a single antibody binds two unrelated protein antigens. We showed that while bH1 and the affinity-improved bH1-44, in particular, maintained many properties of Herceptin including binding affinity, kinetics and the use of residues for antigen recognition, they differed in the binding thermodynamics. The interactions of bH1 and its variants with both antigens were characterized by large favorable entropy changes whereas the Herceptin/HER2 interaction involved a large favorable enthalpy change. By dissecting the total entropy change and the energy barrier for dual interaction, we determined that the significant structural plasticity of the bH1 antibodies demanded by the dual specificity did not translate into the expected increase of entropic penalty relative to Herceptin. Clearly, dual antigen recognition of the Herceptin variants involves divergent antibody conformations of nearly equivalent energetic states. Hence, increasing the structural plasticity of an antigen-binding site without increasing the entropic cost may play a role for antibodies to evolve multi-specificity. Our report represents the first comprehensive biophysical analysis of a high affinity dual specific antibody binding two unrelated protein antigens, furthering our understanding of the thermodynamics that drive the vast antigen recognition capacity of the antibody repertoire.     

Entropy rectifies the Brownian steps of kinesin

Kinesin is a stepping motor that successively produces forward and backward 8-nm steps along microtubules. Under physiological conditions, the steps powering kinesin's motility are biased in one direction and drive various biological motile processes. The physical mechanism underlying the unidirectional bias of the kinesin steps is not fully understood. Here we explored the mechanical kinetics and thermodynamics of forward and backward kinesin steps by analyzing their temperature and load dependence. Results show that the frequency asymmetry between forward and backward steps is produced by entropy. Furthermore, the magnitude of the entropic asymmetry is 6 k(B)T, more than three times greater than expected from a current model, in which a mechanical conformational change within the kinesin molecular structure directly biases the kinesin steps forward. We propose that the stepping direction of kinesin is preferably caused by an entropy asymmetry resulting from the compatibility between the kinesin and microtubule interaction based on their polar structures.     

Enhanced Human Tissue Microdialysis Using Hydroxypropyl-?-Cyclodextrin as Molecular Carrier

Microdialysis sampling of lipophilic molecules in human tissues is challenging because protein binding and adhesion to the membrane limit recovery. Hydroxypropyl-ß-cyclodextrin (HP-ß-CD) forms complexes with hydrophobic molecules thereby improving microdialysis recovery of lipophilic molecules in vitro and in rodents. We tested the approach in human subjects. First, we determined HP-ß-CD influences on metabolite stability, delivery, and recovery in vitro. Then, we evaluated HP-ß-CD as microdialysis perfusion fluid supplement in 20 healthy volunteers. We placed 20 kDa microdialysis catheters in subcutaneous abdominal adipose tissue and in the vastus lateralis muscle. We perfused catheters with lactate free Ringer solution with or without 10% HP-ß-CD at flow rates of 0.3–2.0 µl/min. We assessed tissue metabolites, ultrafiltration effects, and blood flow. In both tissues, metabolite concentrations with Ringer+HP-ß-CD perfusate were equal or higher compared to Ringer alone. Addition of HP-ß-CD increased dialysate volume by 10%. Adverse local or systemic reactions to HP-ß-CD did not occur and analytical methods were not disturbed. HP-ß-CD addition allowed to measure interstitial anandamide concentrations, a highly lipophilic endogenous molecule. Our findings suggest that HP-ß-CD is a suitable supplement in clinical microdialysis to enhance recovery of lipophilic molecules from human interstitial fluid.     

Dynamics of intracellular neonatal Fc receptor–ligand interactions in primary macrophages using biophysical fluorescence techniques

The neonatal Fc receptor (FcRn) is responsible for the recycling of endocytosed albumin and IgG, and contributes to their long plasma half-life. We recently identified an FcRn-dependent recycling pathway from macropinosomes in macrophages; however, little is known about the dynamics of intracellular FcRn–ligand interactions to promote recycling. Here we demonstrate a multiplexed biophysical fluorescent microscopy approach to resolve the spatiotemporal dynamics of albumin–FcRn interactions in living bone marrow–derived macrophages (BMDMs). We used the phasor approach to fluorescence lifetime imaging microscopy (FLIM) of Förster resonance energy transfer (FRET) to detect the interaction of a FcRn–mCherry fusion protein with endocytosed Alexa Fluor 488–labeled human serum albumin (HSA–AF488) in BMDMs, and raster image correlation spectroscopy (RICS) analysis of single fluorescent-labeled albumin molecules to monitor the diffusion kinetics of internalized albumin. Our data identified a major fraction of immobile HSA–AF488 molecules in endosomal structures of human FcRn-positive mouse macrophages and an increase in FLIM-FRET following endocytosis, including detection of FRET in tubular-like structures. A nonbinding mutant of albumin showed minimum FLIM-FRET and high mobility. These data reveal the kinetics of FcRn–ligand binding within endosomal structures for recruitment into transport carriers for recycling. These approaches have wide applicability for analyses of intracellular ligand–receptor interactions.     

Solution and Crystal Molecular Dynamics Simulation Study of m4-Cyanovirin-N Mutants Complexed with Di-Mannose

Cyanovirin-N (CVN) is a highly potent anti-HIV carbohydrate-binding agent that establishes its microbicide activity through interaction with mannose-rich glycoprotein gp120 on the virion surface. The m4-CVN and P51G-m4-CVN mutants represent simple models for studying the high-affinity binding site, B^M. A recently determined 1.35 Å high-resolution structure of P51G-m4-CVN provided details on the di-mannose binding mechanism, and suggested that the Arg-76 and Glu-41 residues are critical components of high mannose specificity and affinity. We performed molecular-dynamics simulations in solution and a crystal environment to study the role of Arg-76. Network analysis and clustering were used to characterize the dynamics of Arg-76. The results of our explicit solvent solution and crystal simulations showed a significant correlation with conformations of Arg-76 proposed from x-ray crystallographic studies. However, the crystal simulation showed that the crystal environment strongly biases conformational sampling of the Arg-76 residue. The solution simulations demonstrated no conformational preferences for Arg-76, which would support its critical role as the residue that locks the ligand in the bound state. Instead, a comparative analysis of trajectories from >50 ns of simulation for two mutants revealed the existence of a very stable eight-hydrogen-bond network between the di-mannose ligand and predominantly main-chain atoms. This network may play a key role in the specific recognition and strong binding of mannose oligomers in CVN and its homologs.     

Thermodynamic characterization of the multivalent binding of chartreusin to DNA

Characterization of the thermodynamics of DNA– drug interactions is a very useful part in rational drug design. Isothermal titration calorimetry (ITC), differential scanning calorimetry (DSC) and UV melting experiments have been used to analyze the multivalent (intercalation plus minor groove) binding of the antitumor antibiotic chartreusin to DNA. Using DNA UV melting studies in the presence of the ligand and the binding enthalpy determined by ITC, we determined that the binding constant for the interaction was 3.6 × 10⁵ M^–1 at 20°C, in a solution containing 18 mM Na⁺. The DNA–drug interaction was enthalpy driven, with a ΔH_b of –7.07 kcal/mol at 20°C. Binding enthalpies were determined by ITC in the 20–35°C range and used to calculate a binding-induced change in heat capacity (ΔCp) of –391 cal/mol K. We have obtained a detailed thermodynamic profile for the interaction of this multivalent drug, which makes possible a dissection of ΔG_obs into the component free energy terms. The hydrophobic transfer of the chartreusin chromophore from the solution to the DNA intercalating site is the main contributor to the free energy of binding.     

Hydrophobe-water interactions: methane as a model

Recent molecular-dynamics simulations have demonstrated that the use of an empirical hydrophobic potential displaying two minima, i.e., one for hydrophobes in close contact and one for hydrophobes separated by a hydration layer, leads to a marked improvement in protein structure prediction. This potential is supported by experimental data and simulations, but its physical origin and mathematical formulation have not been derived as yet. Here we show that water-mediated attraction (the “wetting regime”) between two hydrophobic molecules originates in the interaction between the dipoles induced at the surface of the hydrophobes by the surrounding structured water. As an example, we derive the effective hydrophobic potential that describes the interaction between two methane molecules, a classical model of a double-well energy function. We found an excellent agreement with published results from all-atom, explicit solvent molecular-dynamics simulations of this interaction. The approach presented here provides the theoretical basis for implementing an adequate representation of the wetting regime of the hydrophobic interactions in force fields used for structure prediction. The results are useful for modeling both protein folding and binding.     

Atomistic Molecular-Dynamics Simulations Enable Prediction of the Arginine Permeation Pathway through OccD1/OprD from Pseudomonas?aeruginosa

Pseudomonas aeruginosa is a Gram-negative bacterium that does not contain large, nonspecific porins in its outer membrane. Consequently, the outer membrane is highly impermeable to polar solutes and serves as a barrier against the penetration of antimicrobial agents. This is one of the reasons why such bacteria are intrinsically resistant to antibiotics. Polar molecules that permeate across the outer membrane do so through substrate-specific channels proteins. To design antibiotics that target substrate-channel proteins, it is essential to first identify the permeation pathways of their natural substrates. In P. aeruginosa, the largest family of substrate-specific proteins is the OccD (previously reported under the name OprD) family. Here, we employ equilibrium and steered molecular-dynamics simulations to study OccD1/OprD, the archetypical member of the OccD family. We study the permeation of arginine, one of the natural substrates of OccD1, through the protein. The combination of simulation methods allows us to predict the pathway taken by the amino acid, which is enabled by conformational rearrangements of the extracellular loops of the protein. Furthermore, we show that arginine adopts a specific orientation to form the molecular interactions that facilitate its passage through part of the protein. We predict a three-stage permeation process for arginine.     

Thermodynamics of Antimicrobial Lipopeptide Binding to Membranes: Origins of Affinity and Selectivity

The development of novel antibiotic drugs is one of the most pressing biomedical problems due to the increasing number of antibiotic-resistant pathogens. Antimicrobial peptides and lipopeptides are a promising category of candidates, but the molecular origins of their antimembrane activity is unclear. Here we explore a series of recently developed antimicrobial lipopeptides, using coarse-grained molecular-dynamics simulations and free energy methods to uncover the thermodynamics governing their binding to membranes. Specifically, we quantify C16-KGGK’s binding affinity to the two types of membrane by umbrella sampling. We also examined the origin of C16-KGGK’s selectivity for bacterial versus mammalian membranes by systematically varying the peptide sequence and salt concentration. Our data showed that the C16 hydrophobic tail is the main contributor to its affinity to lipid membrane, whereas the peptide portion is mainly responsible for its selectivity. Furthermore, the electrostatic interaction between the cationic peptide and anionic bacterial membrane plays a significant role in the selectivity.     

Competitive inhibition of the classical complement pathway using exogenous single-chain C1q recognition proteins

Complement component C1q is a protein complex of the innate immune system with well-characterized binding partners that constitutes part of the classical complement pathway. In addition, C1q was recently described in the central nervous system as having a role in synapse elimination both in the healthy brain and in neurodegenerative diseases. However, the molecular mechanism of C1q-associated synapse phagocytosis is still unclear. Here, we designed monomer and multimer protein constructs, which comprised the globular interaction recognition parts of mouse C1q (globular part of C1q [gC1q]) as single-chain molecules (sc-gC1q proteins) lacking the collagen-like effector region. These molecules, which can competitively inhibit the function of C1q, were expressed in an Escherichia coli expression system, and their structure and capabilities to bind known complement pathway activators were validated by mass spectrometry, analytical size-exclusion chromatography, analytical ultracentrifugation, CD spectroscopy, and ELISA. We further characterized the interactions between these molecules and immunoglobulins and neuronal pentraxins using surface plasmon resonance spectroscopy. We demonstrated that sc-gC1qs potently inhibited the function of C1q. Furthermore, these sc-gC1qs competed with C1q in binding to the embryonal neuronal cell membrane. We conclude that the application of sc-gC1qs can reveal neuronal localization and functions of C1q in assays in vivo and might serve as a basis for engineering inhibitors for therapeutic purposes.     

Molecular Simulations of a Dynamic Protein Complex: Role of Salt-Bridges and Polar Interactions in Configurational Transitions

Ion charge pairs and hydrogen bonds have been extensively studied for their roles in stabilizing protein complexes and in steering the process of protein association. Recently, it has become clear that some protein complexes are dynamic in that they interconvert between several alternate configurations. We have previously characterized one such system: the EphA2:SHIP2 SAM-SAM heterodimer by solution NMR. Here we carried out extensive all-atom molecular-dynamics simulations on a microsecond time-scale starting with different NMR-derived structures for the complex. Transitions are observed between several discernible configurations at average time intervals of 50–100 ns. The domains reorient relative to one another by substantial rotation and a slight shifting of the interfaces. Bifurcated and intermediary salt-bridge and hydrogen-bond interactions play a role in the transitions in a process that can be described as moving along a “monkey-bar”. We notice an increased density of salt bridges near protein interaction surfaces that appear to enable these transitions, also suggesting why the trajectories can become kinetically hindered in regions where fewer of such interactions are possible. In this context, even microsecond molecular-dynamics simulations are not sufficient to sample the energy landscape unless the structures remain close to their experimentally derived low-energy configurations.     

Sequence requirements for the recognition of tyrosine-based endocytic signals by clathrin AP-2 complexes.

We recently determined that fusion proteins containing tyrosine-based endocytic signals bind to the mu 2 subunit of AP-2, the complex that drives clathrin coat formation and mediates endocytosis from the plasma membrane. Here we analyze the selectivity of peptide recognition by mu 2 and by AP-2 using combinatorial selection methods and surface plasmon resonance. Both mu 2 and AP-2 are shown to interact with various sequences of the form tyrosine-polar-polar-hydrophobic (Yppø) found on receptors that follow the clathrin pathway. The optimal sequence for interaction with mu 2 and with AP-2 has tyrosine as an anchor and prefers arginine at position Y + 2 and leucine at position Y + 3. In contrast, no preferred sequence is detected surrounding the Yppø signal, indicating that recognition of the Yppø endocytic signal does not require a prefolded structure. We conclude that sorting into the endocytic pathway is governed by a surprisingly simple interaction between the mu 2 chain and a tyrosine-containing tetrapeptide sequence.     

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k_2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.     

Dipole-dipole interaction stabilizes the transition state of apurinic/apyrimidinic endonuclease--abasic site interaction

Human apurinic/apyrimidinic (AP) endonuclease (hAPE) initiates the repair of an abasic site (AP site). To gain insight into the mechanisms of damage recognition of hAPE, we conducted surface plasmon resonance spectroscopy to study the thermodynamics and kinetics of its interaction with substrate DNA containing an abasic site (AP DNA). The affinity of hAPE binding toward DNA increased as much as 6-fold after replacing a single adenine (equilibrium dissociation constant, K(D), 5.3 nm) with an AP site (K(D), 0.87 nm). The enzyme-substrate complex formation appears to be thermodynamically stabilized and favored by a large change in Gibbs free energy, DeltaG degrees (-50 kJ/mol). The latter is supported by a high negative change in enthalpy, DeltaH degrees (-43 kJ/mol) and also positive change in entropy, DeltaS degrees (24 J/(K mol)), and thus the binding process is spontaneous at all temperatures. Analysis of kinetic parameters reveals small enthalpy of activation for association, DeltaH degrees++(ass) (-17 kJ/mol), and activation energy for association (E(a), -14 kJ/mol) when compared with the enthalpy of activation for dissociation, DeltaH degrees++(diss) (26 kJ/mol), and activation energy in the reverse direction (E(d), 28 kJ/mol). Furthermore, varying concentration of KCl showed an increase in binding affinity at low concentration but complete abrogation of the binding at higher concentration, implying the importance of hydrophobic, but predominantly ionic, forces in the Michaelis-Menten complex formation. Thus, low activation energy and the enthalpy of activation, which are perhaps a result of dipole-dipole interactions, play critical roles in AP site binding of APE.     

Biophysics of Active Vesicle Transport,an Intermediate Step That Couples Excitation and Exocytosis of Serotonin in the Neuronal Soma

Transmitter exocytosis from the neuronal soma is evoked by brief trains of high frequency electrical activity and continues for several minutes. Here we studied how active vesicle transport towards the plasma membrane contributes to this slow phenomenon in serotonergic leech Retzius neurons, by combining electron microscopy, the kinetics of exocytosis obtained from FM1-43 dye fluorescence as vesicles fuse with the plasma membrane, and a diffusion equation incorporating the forces of local confinement and molecular motors. Electron micrographs of neurons at rest or after stimulation with 1 Hz trains showed cytoplasmic clusters of dense core vesicles at 1.5±0.2 and 3.7±0.3 µm distances from the plasma membrane, to which they were bound through microtubule bundles. By contrast, after 20 Hz stimulation vesicle clusters were apposed to the plasma membrane, suggesting that transport was induced by electrical stimulation. Consistently, 20 Hz stimulation of cultured neurons induced spotted FM1-43 fluorescence increases with one or two slow sigmoidal kinetics, suggesting exocytosis from an equal number of vesicle clusters. These fluorescence increases were prevented by colchicine, which suggested microtubule-dependent vesicle transport. Model fitting to the fluorescence kinetics predicted that 52–951 vesicles/cluster were transported along 0.60–6.18 µm distances at average 11–95 nms⁻¹ velocities. The ATP cost per vesicle fused (0.4–72.0), calculated from the ratio of the ΔG_process/ΔG_ATP, depended on the ratio of the traveling velocity and the number of vesicles in the cluster. Interestingly, the distance-dependence of the ATP cost per vesicle was bistable, with low energy values at 1.4 and 3.3 µm, similar to the average resting distances of the vesicle clusters, and a high energy barrier at 1.6–2.0 µm. Our study confirms that active vesicle transport is an intermediate step for somatic serotonin exocytosis by Retzius neurons and provides a quantitative method for analyzing similar phenomena in other cell types.