Search

Neurotransmitter (NT) release is accomplished by a machinery that unclamps fusion in response to calcium and then fuses the synaptic vesicle and plasma membranes. These are often thought of as distinct tasks assigned to non-overlapping components. Vesicle release rates have a power law dependence on [Ca2+] with an exponent of 3-5, long taken to indicate that 3-5 Ca2+ions bind the calcium sensor Synaptotagmin to trigger release. However, dependencies at low [Ca] are inconsistent with simple sequential binding to a single Ca2+sensor followed by a final fusion step. Here we developed coarse-grained molecular dynamics simulations of the NT release machinery accounting for Synaptotagmin-mediated unclamping and SNARE-mediated fusion. Calcium-triggered unclamping and SNARE-mediated fusion emerged from simulations as contemporaneous, coupled processes. Increasing cytosolic [Ca2+], the instantaneous fusion rate increased as SNAREpins were progressively and reversibly released by dissociation of Synaptotagmin-SNAREpin complexes. Simulations reproduced the observed dependence of release rates on [Ca2+], but the power law was unrelated to the number of Ca2+ions required. Action potential-evoked vesicle release probabilities depended on the number of transiently unclamped SNAREpins, explaining experimental dependencies of release probabilities on both unclamping and membrane-fusing machinery components. These results describe a highly cooperative NT release machinery with intrinsically inseparable unclamping and membrane-fusing functionalities.

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ(3), we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ(3) bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.

SNARE proteins constitute the core of the exocytotic membrane fusion machinery. Fusion occurs when vesicle-associated and target membrane-associated SNAREs zipper into trans-SNARE complexes ('SNAREpins'), but the number required is controversial and the mechanism of cooperative fusion is poorly understood. We developed a highly coarse-grained molecular dynamics simulation to access the long fusion timescales, which revealed a two-stage process. First, zippering energy was dissipated and cooperative entropic forces assembled the SNAREpins into a ring; second, entropic forces expanded the ring, pressing membranes together and catalyzing fusion. We predict that any number of SNAREs fuses membranes, but fusion is faster with more SNAREs.

SummaryUnderstanding how functional lipid domains in live cell membranes are generated has posed a challenge. Here, we show that transbilayer interactions are necessary for the generation of cholesterol-dependent nanoclusters of GPI-anchored proteins mediated by membrane-adjacent dynamic actin filaments. We find that long saturated acyl-chains are required for forming GPI-anchor nanoclusters. Simultaneously, at the inner leaflet, long acyl-chain-containing phosphatidylserine (PS) is necessary for transbilayer coupling. All-atom molecular dynamics simulations of asymmetric multicomponent-membrane bilayers in a mixed phase provide evidence that immobilization of long saturated acyl-chain lipids at either leaflet stabilizes cholesterol-dependent transbilayer interactions forming local domains with characteristics similar to a liquid-ordered (lo) phase. This is verified by experiments wherein immobilization of long acyl-chain lipids at one leaflet effects transbilayer interactions of corresponding lipids at the opposite leaflet. This suggests a general mechanism for the generation and stabilization of nanoscale cholesterol-dependent and actin-mediated lipid clusters in live cell membranes.

Lipid/cholesterol mixtures derived from cell membranes as well as their synthetic reconstitutions exhibit well-defined miscibility phase transitions and critical phenomena near physiological temperatures. This suggests that lipid/cholesterol-mediated phase separation plays a role in the organization of live cell membranes. However, macroscopic lipid-phase separation is not generally observed in cell membranes, and the degree to which properties of isolated lipid mixtures are preserved in the cell membrane remain unknown. A fundamental property of phase transitions is that the variation of tagged particle diffusion with temperature exhibits an abrupt change as the system passes through the transition, even when the two phases are distributed in a nanometer-scale emulsion. We support this using a variety of Monte Carlo and atomistic simulations on model lipid membrane systems. However, temperature-dependent fluorescence correlation spectroscopy of labeled lipids and membrane-anchored proteins in live cell membranes shows a consistently smooth increase in the diffusion coefficient as a function of temperature. We find no evidence of a discrete miscibility phase transition throughout a wide range of temperatures: 14-37 °C. This contrasts the behavior of giant plasma membrane vesicles (GPMVs) blebbed from the same cells, which do exhibit phase transitions and macroscopic phase separation. Fluorescence lifetime analysis of a DiI probe in both cases reveals a significant environmental difference between the live cell and the GPMV. Taken together, these data suggest the live cell membrane may avoid the miscibility phase transition inherent to its lipid constituents by actively regulating physical parameters, such as tension, in the membrane.

Integrin binding to extracellular matrix proteins is regulated by conformational transitions from closed, low affinity states to open, high affinity states. However, the pathways of integrin conformational activation remain incompletely understood. Here, by combining all-atom molecular dynamics simulation, coarse-graining, heterogeneous elastic network modeling, and experimental ligand binding measurements, we test the effect of integrin β mutations that destabilize the closed conformation. Our results support a “deadbolt” model of integrin activation, where extension of the headpiece is not coupled to leg separation, consistent with recent cryo-EM reconstructions of integrin intermediates. Moreover, our results are inconsistent with a “switchblade-like” mechanism. The data show that locally correlated atomistic motions are likely responsible for extension of integrin headpiece before separation of transmembrane legs, without persistence of these correlations across the entire protein. By combining modeling and simulation with experiment, this study provides new insight into the structural basis of full-length integrin activation.

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling. publishedVersion

The molecular mechanism of ethanol and its effects on neurological function is far from clear. In this study, we investigate the effects of ethanol on various structural and dynamical properties of mixed bilayers consisting of different ratios of dipalmitoylphosphatidylcholine (DPPC), sphingomyelin (SM) and cholesterol that are typical constituents of neural cell membranes (Calderon et al., 1995) using molecular dynamics (MD) simulations. The bilayer properties such as thickness, hydrophobic chain order, and diffusive motion of individual lipids as well collective properties like lateral pressure profiles are affected by the presence of ethanol molecules. The simulations show that the percentage of cholesterol present in the bilayers significantly affects the depth of penetration of ethanol molecules. In particular, presence of very high concentration of cholesterol molecules enhances the rigidity of the bilayer and renders them resistant to the penetration of the ethanol molecules, consistent with experiments. Ethanol molecules compete with cholesterol molecules for hydrogen bonding and disrupt cholesterol-lipid interactions, especially those between SM and cholesterol. Ethanol molecules also affect the lateral pressure profiles in the bilayer systems. These results may have implications in understanding the general anesthetic mechanism and role played by cholesterol on partitioning of such anesthetic/alcohol molecules into cell membranes.

The cell membrane is inherently asymmetric and heterogeneous in its composition, a feature that is crucial for its function. Using atomistic molecular dynamics simulations, the physical properties of a 3-component asymmetric mixed lipid bilayer system comprising an unsaturated POPC (palmitoyloleoylphosphatidylcholine), a saturated PSM (palmitoylsphingomyelin), and cholesterol are investigated. Our simulations explore both the dynamics of coarsening following a quench from the mixed phase and the final phase-segregated regime obtained by equilibrating a fully segregated configuration. Following a quench, the membrane quickly enters a coarsening regime, where the initial stages of liquid ordered, l(o), domain formation are observed. These growing domains are found to be highly enriched in cholesterol and PSM. Consistent with this, the final phase-segregated regime contains large l(o) domains at equilibrium, enriched in cholesterol and PSM. Our simulations suggest that the cholesterol molecules may partition into these PSM-dominated regions in the ratio of 3:1 when compared to POPC-dominated regions. PSM molecules exhibit a measurable tilt and long-range tilt correlations within the l(o) domain as a consequence of the asymmetry of the bilayer, with implications to local membrane deformation and budding. Tagged particle diffusion for PSM and cholesterol molecules, which reflects spatial variations in the physical environment encountered by the tagged particle, is computed and compared with recent experimental results obtained from high-resolution microscopy.

Despite the importance of integrin proteins in cell signal transduction and force generation, the mechanism of action is not understood at the atomic level. To date, no experimental method has been able to probe the structure of integrin proteins in their open, force-activated state. To gain a better understanding of the transformation of the inactive, closed, state to an active state of the integrin, we are employing advanced computational techniques that will allow us to more quickly sample the conformational change and probe the pathway of the transformation. To make connection with the experimental studies, we have first studied the effect of single and double amino acid substitutions to the wild type integrin. Our simulations reveal certain mutants which destabilize the closed state by increasing the kink angle of the extracellular region and distance between the transmembrane region of the alpha and beta subunits of the integrin. We have also developed a coarse-grained model based on our atomistic simulation data, and are using that model to probe the transition of integrin to the open state.

Localized contractile configurations or asters spontaneously appear and disappear as emergent structures in the collective stochastic dynamics of active polar actomyosin filaments. Passive parti- cles which (un)bind to the active filaments get advected into the asters, forming transient clusters. We study the phase segregation of such passive advective scalars in a medium of dynamic asters, as a function of the aster density and the ratio of the rates of aster remodeling to particle diffusion. The dynamics of coarsening shows strong violation of Porod behaviour, suggesting diffuse interfaces. The phase segregated steady state shows strongly fluctuations characterized by multiscaling and in- termittency. We expect these unique nonequilibrium features to manifest in the actin-dependent molecular clustering at the cell surface., 5 pages, 4 figures

Membrane proteins that couple to cortical actin show temperature-independent diffusion. The loss of this coupling and perturbation of cortical actomyosin dynamics render the diffusion temperature dependent. The findings suggest that active fluctuations arising from dynamic actin filaments at the cortex can drive diffusion on the cell membrane., Molecular diffusion at the surface of living cells is believed to be predominantly driven by thermal kicks. However, there is growing evidence that certain cell surface molecules are driven by the fluctuating dynamics of cortical cytoskeleton. Using fluorescence correlation spectroscopy, we measure the diffusion coefficient of a variety of cell surface molecules over a temperature range of 24–37°C. Exogenously incorporated fluorescent lipids with short acyl chains exhibit the expected increase of diffusion coefficient over this temperature range. In contrast, we find that GPI-anchored proteins exhibit temperature-independent diffusion over this range and revert to temperature-dependent diffusion on cell membrane blebs, in cells depleted of cholesterol, and upon acute perturbation of actin dynamics and myosin activity. A model transmembrane protein with a cytosolic actin-binding domain also exhibits the temperature-independent behavior, directly implicating the role of cortical actin. We show that diffusion of GPI-anchored proteins also becomes temperature dependent when the filamentous dynamic actin nucleator formin is inhibited. However, changes in cortical actin mesh size or perturbation of branched actin nucleator Arp2/3 do not affect this behavior. Thus cell surface diffusion of GPI-anchored proteins and transmembrane proteins that associate with actin is driven by active fluctuations of dynamic cortical actin filaments in addition to thermal fluctuations, consistent with expectations from an “active actin-membrane composite” cell surface.

We present a study of the bend angle distribution of semiflexible polymers of short and intermediate lengths within the wormlike chain model. This enables us to calculate the elastic response of a stiff molecule to a bending moment. Our results go beyond the Hookean regime and explore the nonlinear elastic behaviour of a single molecule. We present analytical formulae for the bend angle distribution and for the moment-angle relation. Our analytical study is compared against numerical Monte Carlo simulations. The functional forms derived here can be applied to fluorescence microscopic studies on actin and DNA. Our results are relevant to recent studies in "kinks" and cyclization in short and intermediate length DNA strands., 4 pages

A question of considerable interest to cell membrane biology is whether phase segregated domains across an asymmetric bilayer are strongly correlated with each other and whether phase segregation in one leaflet can induce segregation in the other. We answer both these questions in the affirmative, using an atomistic molecular dynamics simulation to study the equilibrium statistical properties of a 3-component {\em asymmetric} lipid bilayer comprising an unsaturated POPC (palmitoyl-oleoyl-phosphatidyl-choline), a saturated SM (sphingomyelin) and cholesterol with different composition ratios. Our simulations are done by fixing the composition of the upper leaflet to be at the coexistence of the liquid ordered ($l_o$) - liquid disordered ($l_d$) phases, while the composition of the lower leaflet is varied from the phase coexistence regime to the mixed $l_d$ phase, across a first-order phase boundary. In the regime of phase coexistence in each leaflet, we find strong transbilayer correlations of the $l_o$ domains across the two leaflets, resulting in {\it bilayer registry}. This transbilayer correlation depends sensitively upon the chain length of the participating lipids and possibly other features of lipid chemistry, such as degree of saturation. We find that the $l_o$ domains in the upper leaflet can {\em induce} phase segregation in the lower leaflet, when the latter is nominally in the mixed ($l_d$) phase., Comment: 6 figures

The spatial distribution and remodeling dynamics of nanoscale clusters of outer-leaflet GPI-anchor proteins (GPI-APs), glycolipids and inner-leaflet Ras proteins are regulated by the activity of the underlying cortical acto-myosin machinery. Recently, we proposed a theoretical framework based on active hydrodynamics of short polar filaments interacting with myosin motors to explain the nanoclustering of the GPI-APs. This framework involves the coupling of cell-surface molecules to these dynamic filaments of actin at the inner leaflet, this in turn leads to their transient clustering. We predicted and experimentally confirmed that transmembrane proteins with cytoplasmic domains capable of directly binding filamentous actin (TM-ABD), would be organized both at nanoscale (