3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.Comment: 39 pages, 14 figures. High resolution figures and supplemental movies available upon reques

The macroscopic effects of microscopic heterogeneity

Over the past decade, advances in super-resolution microscopy and particle-based modeling have driven an intense interest in investigating spatial heterogeneity at the level of single molecules in cells. Remarkably, it is becoming clear that spatiotemporal correlations between just a few molecules can have profound effects on the signaling behavior of the entire cell. While such correlations are often explicitly imposed by molecular structures such as rafts, clusters, or scaffolds, they also arise intrinsically, due strictly to the small numbers of molecules involved, the finite speed of diffusion, and the effects of macromolecular crowding. In this chapter we review examples of both explicitly imposed and intrinsic correlations, focusing on the mechanisms by which microscopic heterogeneity is amplified to macroscopic effect.Comment: 20 pages, 5 figures. To appear in Advances in Chemical Physic

Polar Chemoreceptor Clustering by Coupled Trimers of Dimers

Receptors of bacterial chemotaxis form clusters at the cell poles, where clusters act as "antennas" to amplify small changes in ligand concentration. Interestingly, chemoreceptors cluster at multiple length scales. At the smallest scale, receptors form dimers, which assemble into stable timers of dimers. At a large scale, trimers form large polar clusters composed of thousands of receptors. Although much is known about the signaling properties emerging from receptor clusters, it is unknown how receptors localize at the cell poles and what the cluster-size determining factors are. Here, we present a model of polar receptor clustering based on coupled trimers of dimers, where cluster size is determined as a minimum of the cluster-membrane free energy. This energy has contributions from the cluster-membrane elastic energy, penalizing large clusters due to their high intrinsic curvature, and receptor-receptor coupling favoring large clusters. We find that the reduced cluster-membrane curvature mismatch at the curved cell poles leads to large and robust polar clusters in line with experimental observation, while lateral clusters are efficiently suppressed.Comment: 11 pages, 6 figures, and 1 tabl

Bistability in Apoptosis by Receptor Clustering

Apoptosis is a highly regulated cell death mechanism involved in many physiological processes. A key component of extrinsically activated apoptosis is the death receptor Fas, which, on binding to its cognate ligand FasL, oligomerize to form the death-inducing signaling complex. Motivated by recent experimental data, we propose a mathematical model of death ligand-receptor dynamics where FasL acts as a clustering agent for Fas, which form locally stable signaling platforms through proximity-induced receptor interactions. Significantly, the model exhibits hysteresis, providing an upstream mechanism for bistability and robustness. At low receptor concentrations, the bistability is contingent on the trimerism of FasL. Moreover, irreversible bistability, representing a committed cell death decision, emerges at high concentrations, which may be achieved through receptor pre-association or localization onto membrane lipid rafts. Thus, our model provides a novel theory for these observed biological phenomena within the unified context of bistability. Importantly, as Fas interactions initiate the extrinsic apoptotic pathway, our model also suggests a mechanism by which cells may function as bistable life/death switches independently of any such dynamics in their downstream components. Our results highlight the role of death receptors in deciding cell fate and add to the signal processing capabilities attributed to receptor clustering.Comment: Accepted by PLoS Comput Bio

Structure and dynamics of the E. coli chemotaxis core signaling complex by cryo-electron tomography and molecular simulations

To enable the processing of chemical gradients, chemotactic bacteria possess large arrays of transmembrane chemoreceptors, the histidine kinase CheA, and the adaptor protein CheW, organized as coupled core-signaling units (CSU). Despite decades of study, important questions surrounding the molecular mechanisms of sensory signal transduction remain unresolved, owing especially to the lack of a high-resolution CSU structure. Here, we use cryo-electron tomography and sub-tomogram averaging to determine a structure of the Escherichia coli CSU at sub-nanometer resolution. Based on our experimental data, we use molecular simulations to construct an atomistic model of the CSU, enabling a detailed characterization of CheA conformational dynamics in its native structural context. We identify multiple, distinct conformations of the critical P4 domain as well as asymmetries in the localization of the P3 bundle, offering several novel insights into the CheA signaling mechanism

Selective disruption of high sensitivity heat activation but not capsaicin activation of TRPV1 channels by pore turret mutations.

The capsaicin receptor transient receptor potential vanilloid (TRPV)1 is a highly heat-sensitive ion channel. Although chemical activation and heat activation of TRPV1 elicit similar pungent, painful sensation, the molecular mechanism underlying synergistic activation remains mysterious. In particular, where the temperature sensor is located and whether heat and capsaicin share a common activation pathway are debated. To address these fundamental issues, we searched for channel mutations that selectively affected one form of activation. We found that deletion of the first 10 amino acids of the pore turret significantly reduced the heat response amplitude and shifted the heat activation threshold, whereas capsaicin activation remained unchanged. Removing larger portions of the turret disrupted channel function. Introducing an artificial sequence to replace the deleted region restored sensitive capsaicin activation in these nonfunctional channels. The heat activation, however, remained significantly impaired, with the current exhibiting diminishing heat sensitivity to a level indistinguishable from that of a voltage-gated potassium channel, Kv7.4. Our results demonstrate that heat and capsaicin activation of TRPV1 are structurally and mechanistically distinct processes, and the pore turret is an indispensible channel structure involved in the heat activation process but is not part of the capsaicin activation pathway. Synergistic effect of heat and capsaicin on TRPV1 activation may originate from convergence of the two pathways on a common activation gate

Protein connectivity in chemotaxis receptor complexes

The chemotaxis sensory system allows bacteria such as Escherichia coli to swim towards nutrients and away from repellents. The underlying pathway is remarkably sensitive in detecting chemical gradients over a wide range of ambient concentrations. Interactions among receptors, which are predominantly clustered at the cell poles, are crucial to this sensitivity. Although it has been suggested that the kinase CheA and the adapter protein CheW are integral for receptor connectivity, the exact coupling mechanism remains unclear. Here, we present a statistical-mechanics approach to model the receptor linkage mechanism itself, building on nanodisc and electron cryotomography experiments. Specifically, we investigate how the sensing behavior of mixed receptor clusters is affected by variations in the expression levels of CheA and CheW at a constant receptor density in the membrane. Our model compares favorably with dose-response curves from in vivo Förster resonance energy transfer (FRET) measurements, demonstrating that the receptor-methylation level has only minor effects on receptor cooperativity. Importantly, our model provides an explanation for the non-intuitive conclusion that the receptor cooperativity decreases with increasing levels of CheA, a core signaling protein associated with the receptors, whereas the receptor cooperativity increases with increasing levels of CheW, a key adapter protein. Finally, we propose an evolutionary advantage as explanation for the recently suggested CheW-only linker structures

Cellular mechanisms for cargo delivery and polarity maintenance at different polar domains in plant cells

The asymmetric localization of proteins in the plasma membrane domains of eukaryotic cells is a fundamental manifestation of cell polarity that is central to multicellular organization and developmental patterning. In plants, the mechanisms underlying the polar localization of cargo proteins are still largely unknown and appear to be fundamentally distinct from those operating in mammals. Here, we present a systematic, quantitative comparative analysis of the polar delivery and subcellular localization of proteins that characterize distinct polar plasma membrane domains in plant cells. The combination of microscopic analyses and computational modeling revealed a mechanistic framework common to diverse polar cargos and underlying the establishment and maintenance of apical, basal, and lateral polar domains in plant cells. This mechanism depends on the polar secretion, constitutive endocytic recycling, and restricted lateral diffusion of cargos within the plasma membrane. Moreover, our observations suggest that polar cargo distribution involves the individual protein potential to form clusters within the plasma membrane and interact with the extracellular matrix. Our observations provide insights into the shared cellular mechanisms of polar cargo delivery and polarity maintenance in plant cells