Mathematical simulation of membrane protein clustering for efficient signal transduction

Engineering a Novel Probiotic Toolkit in Escherichia coli Nissle 1917 for Sensing and Mitigating Gut Inflammatory Diseases

Inflammatory bowel disease (IBD) is characterized by chronic intestinal inflammation with no cure and limited treatment options that often have systemic side effects. In this study, we developed a target-specific system to potentially treat IBD by engineering the probiotic bacterium Escherichia coli Nissle 1917 (EcN). Our modular system comprises three components: a transcription factor-based sensor (NorR) capable of detecting the inflammation biomarker nitric oxide (NO), a type 1 hemolysin secretion system, and a therapeutic cargo consisting of a library of humanized anti-TNFα nanobodies. Despite a reduction in sensitivity, our system demonstrated a concentration-dependent response to NO, successfully secreting functional nanobodies with binding affinities comparable to the commonly used drug Adalimumab, as confirmed by enzyme-linked immunosorbent assay and in vitro assays. This newly validated nanobody library expands EcN therapeutic capabilities. The adopted secretion system, also characterized for the first time in EcN, can be further adapted as a platform for screening and purifying proteins of interest. Additionally, we provided a mathematical framework to assess critical parameters in engineering probiotic systems, including the production and diffusion of relevant molecules, bacterial colonization rates, and particle interactions. This integrated approach expands the synthetic biology toolbox for EcN-based therapies, providing novel parts, circuits, and a model for tunable responses at inflammatory hotspots.

Insights gained from computational modeling of YAP/TAZ signaling for cellular mechanotransduction

YAP/TAZ signaling pathway is regulated by a multiplicity of feedback loops, crosstalk with other pathways, and both mechanical and biochemical stimuli. Computational modeling serves as a powerful tool to unravel how these different factors can regulate YAP/TAZ, emphasizing biophysical modeling as an indispensable tool for deciphering mechanotransduction and its regulation of cell fate. We provide a critical review of the current state-of-the-art of computational models focused on YAP/TAZ signaling.

Diversity in Notch ligand-receptor signaling interactions

The Notch signaling pathway uses families of ligands and receptors to transmit signals to nearby cells. These components are expressed in diverse combinations in different cell types, interact in a many-to-many fashion, both within the same cell (in cis) and between cells (in trans), and their interactions are modulated by Fringe glycosyltransferases. A fundamental question is how the strength of Notch signaling depends on which pathway components are expressed, at what levels, and in which cells. Here, we used a quantitative, bottom-up, cell-based approach to systematically characterize trans-activation, cis-inhibition, and cis-activation signaling efficiencies across a range of ligand and Fringe expression levels in two mammalian cell types. Each ligand (Dll1, Dll4, Jag1, and Jag2) and receptor variant (Notch1 and Notch2) analyzed here exhibited a unique profile of interactions, Fringe-dependence, and signaling outcomes. All four ligands were able to bind receptors in cis and in trans, and all ligands trans-activated both receptors, although Jag1-Notch1 signaling was substantially weaker than other ligand-receptor combinations. Cis-interactions were predominantly inhibitory, with the exception of the Dll1- and Dll4-Notch2 pairs, which exhibited cis-activation stronger than trans-activation. Lfng strengthened Delta-mediated trans-activation and weakened Jagged-mediated trans-activation for both receptors. Finally, cis-ligands showed diverse cis-inhibition strengths, which depended on the identity of the trans-ligand as well as the receptor. The map of receptor-ligand-Fringe interaction outcomes revealed here should help guide rational perturbation and control of the Notch pathway.

Modeling the Cluster Size Distribution of Vascular Endothelial Growth Factor (VEGF) Receptors

We previously developed a method of defining receptor clusters in the membrane based on mutual distance and applied it to a set of transmission microscopy images of vascular endothelial growth factor receptors. An optimal length parameter was identified, resulting in cluster identification and a procedure that assigned a geometric shape to each cluster. We showed that the observed particle distribution results were consistent with the random placement of receptors within the clusters and, to a lesser extent, the random placement of the clusters on the cell membrane. Here, we develop and validate a stochastic model of clustering, based on a hypothesis of preexisting domains that have a high affinity for receptors. The proximate objective is to clarify the mechanism behind cluster formation and to estimate the effect on signaling. Receptor-enriched domains may significantly impact signaling pathways that rely on ligand-induced dimerization of receptors. We define a simple statistical model, based on the preexisting domain hypothesis, to predict the probability distribution of cluster sizes. The process yielded sets of parameter values that can readily be used in dynamical calculations as the estimates of the quantitative characteristics of the clustering domains.

Exploring the influence of cytosolic and membrane FAK activation on YAP/TAZ nuclear translocation

Membrane binding and unbinding dynamics play a crucial role in the biological activity of several nonintegral membrane proteins, which have to be recruited to the membrane to perform their functions. By localizing to the membrane, these proteins are able to induce downstream signal amplification in their respective signaling pathways. Here, we present a 3D computational approach using reaction-diffusion equations to investigate the relation between membrane localization of focal adhesion kinase (FAK), Ras homolog family member A (RhoA), and signal amplification of the YAP/TAZ signaling pathway. Our results show that the theoretical scenarios in which FAK is membrane bound yield robust and amplified YAP/TAZ nuclear translocation signals. Moreover, we predict that the amount of YAP/TAZ nuclear translocation increases with cell spreading, confirming the experimental findings in the literature. In summary, our in silico predictions show that when the cell membrane interaction area with the underlying substrate increases, for example, through cell spreading, this leads to more encounters between membrane-bound signaling partners and downstream signal amplification. Because membrane activation is a motif common to many signaling pathways, this study has important implications for understanding the design principles of signaling networks.

Rhodopsin Oligomerization and Aggregation

Rhodopsin is the light receptor in photoreceptor cells of the retina and a prototypical G protein-coupled receptor. Two types of quaternary structures can be adopted by rhodopsin. If rhodopsin folds and attains a proper tertiary structure, it can then form oligomers and nanodomains within the photoreceptor cell membrane. In contrast, if rhodopsin misfolds, it cannot progress through the biosynthetic pathway and instead will form aggregates that can cause retinal degenerative disease. In this review, emerging views are highlighted on the supramolecular organization of rhodopsin within the membrane of photoreceptor cells and the aggregation of rhodopsin that can lead to retinal degeneration.

Surface reaction-diffusion kinetics on lattice at the microscopic scale

Microscopic models of reaction-diffusion processes on the cell membrane can link local spatiotemporal effects to macroscopic self-organized patterns often observed on the membrane. Simulation schemes based on the microscopic lattice method (MLM) can model these processes at the microscopic scale by tracking individual molecules, represented as hard spheres, on fine lattice voxels. Although MLM is simple to implement and is generally less computationally demanding than off-lattice approaches, its accuracy and consistency in modeling surface reactions have not been fully verified. Using the Spatiocyte scheme, we study the accuracy of MLM in diffusion-influenced surface reactions. We derive the lattice-based bimolecular association rates for two-dimensional (2D) surface-surface reaction and one-dimensional (1D) volume-surface adsorption according to the Smoluchowski-Collins-Kimball model and random walk theory. We match the time-dependent rates on lattice with off-lattice counterparts to obtain the correct expressions for MLM parameters in terms of physical constants. The expressions indicate that the voxel size needs to be at least 0.6% larger than the molecule to accurately simulate surface reactions on triangular lattice. On square lattice, the minimum voxel size should be even larger, at 5%. We also demonstrate the ability of MLM-based schemes such as Spatiocyte to simulate a reaction-diffusion model that involves all dimensions: three-dimensional (3D) diffusion in the cytoplasm, 2D diffusion on the cell membrane, and 1D cytoplasm-membrane adsorption. With the model, we examine the contribution of the 2D reaction pathway to the overall reaction rate at different reactant diffusivity, reactivity, and concentrations.

TGF-β family co-receptor function and signaling

Transforming growth factor-β (TGF-β) family members, which include TGF-βs, activins and bone morphogenetic proteins, are pleiotropic cytokines that elicit cell type-specific effects in a highly context-dependent manner in many different tissues. These secreted protein ligands signal via single-transmembrane Type I and Type II serine/threonine kinase receptors and intracellular SMAD transcription factors. Deregulation in signaling has been implicated in a broad array of diseases, and implicate the need for intricate fine tuning in cellular signaling responses. One important emerging mechanism by which TGF-β family receptor signaling intensity, duration, specificity and diversity are regulated and/or mediated is through cell surface co-receptors. Here, we provide an overview of the co-receptors that have been identified for TGF-β family members. While some appear to be specific to TGF-β family members, others are shared with other pathways and provide possible ways for signal integration. This review focuses on novel functions of TGF-β family co-receptors, which continue to be discovered.

Cellular mechanosensing of the biophysical microenvironment: A review of mathematical models of biophysical regulation of cell responses

Cells in vivo reside within complex microenvironments composed of both biochemical and biophysical cues. The dynamic feedback between cells and their microenvironments hinges upon biophysical cues that regulate critical cellular behaviors. Understanding this regulation from sensing to reaction to feedback is therefore critical, and a large effort is afoot to identify and mathematically model the fundamental mechanobiological mechanisms underlying this regulation. This review provides a critical perspective on recent progress in mathematical models for the responses of cells to the biophysical cues in their microenvironments, including dynamic strain, osmotic shock, fluid shear stress, mechanical force, matrix rigidity, porosity, and matrix shape. The review highlights key successes and failings of existing models, and discusses future opportunities and challenges in the field.

Differential mast cell outcomes are sensitive to FcεRI-Syk binding kinetics

Cross-linking of immunoglobulin E-bound FcεRI triggers multiple cellular responses, including degranulation and cytokine production. Signaling is dependent on recruitment of Syk via docking of its dual SH2 domains to phosphorylated tyrosines within the FcεRI immunoreceptor tyrosine-based activation motifs. Using single-molecule imaging in live cells, we directly visualized and quantified the binding of individual mNeonGreen-tagged Syk molecules as they associated with the plasma membrane after FcεRI activation. We found that Syk colocalizes transiently to FcεRI and that Syk-FcεRI binding dynamics are independent of receptor aggregate size. Substitution of glutamic acid for tyrosine between the Syk SH2 domains (Syk-Y130E) led to an increased Syk-FcεRI off-rate, loss of site-specific Syk autophosphorylation, and impaired downstream signaling. Genome edited cells expressing only Syk-Y130E were deficient in antigen-stimulated calcium release, degranulation, and production of some cytokines (TNF-a, IL-3) but not others (MCP-1, IL-4). We propose that kinetic discrimination along the FcεRI signaling pathway occurs at the level of Syk-FcεRI interactions, with key outcomes dependent upon sufficiently long-lived Syk binding events.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Effect of Spatial Inhomogeneities on the Membrane Surface on Receptor Dimerization and Signal Initiation

Important signal transduction pathways originate on the plasma membrane, where microdomains may transiently entrap diffusing receptors. This results in a non-random distribution of receptors even in the resting state, which can be visualized as “clusters” by high resolution imaging methods. Here, we explore how spatial in-homogeneities in the plasma membrane might influence the dimerization and phosphorylation status of ErbB2 and ErbB3, two receptor tyrosine kinases that preferentially heterodimerize and are often co-expressed in cancer. This theoretical study is based upon spatial stochastic simulations of the two-dimensional membrane landscape, where variables include differential distributions and overlap of transient confinement zones (“domains”) for the two receptor species. The in silico model is parameterized and validated using data from single particle tracking experiments. We report key differences in signaling output based on the degree of overlap between domains and the relative retention of receptors in such domains, expressed as escape probability. Results predict that a high overlap of domains, which favors transient co-confinement of both receptor species, will enhance the rate of hetero-interactions. Where domains do not overlap, simulations confirm expectations that homo-interactions are favored. Since ErbB3 is uniquely dependent on ErbB2 interactions for activation of its catalytic activity, variations in domain overlap or escape probability markedly alter the predicted patterns and time course of ErbB3 and ErbB2 phosphorylation. Taken together, these results implicate membrane domain organization as an important modulator of signal initiation, motivating the design of novel experimental approaches to measure these important parameters across a wider range of receptor systems.

Classification of dynamical diffusion states in single molecule tracking microscopy

Single molecule tracking of membrane proteins by fluorescence microscopy is a promising method to investigate dynamic processes in live cells. Translating the trajectories of proteins to biological implications, such as protein interactions, requires the classification of protein motion within the trajectories. Spatial information of protein motion may reveal where the protein interacts with cellular structures, because binding of proteins to such structures often alters their diffusion speed. For dynamic diffusion systems, we provide an analytical framework to determine in which diffusion state a molecule is residing during the course of its trajectory. We compare different methods for the quantification of motion to utilize this framework for the classification of two diffusion states (two populations with different diffusion speed). We found that a gyration quantification method and a Bayesian statistics-based method are the most accurate in diffusion-state classification for realistic experimentally obtained datasets, of which the gyration method is much less computationally demanding. After classification of the diffusion, the lifetime of the states can be determined, and images of the diffusion states can be reconstructed at high resolution. Simulations validate these applications. We apply the classification and its applications to experimental data to demonstrate the potential of this approach to obtain further insights into the dynamics of cell membrane proteins.

Evaluation of fluorophores to label SNAP-tag fused proteins for multicolor single-molecule tracking microscopy in live cells

Single-molecule tracking has become a widely used technique for studying protein dynamics and their organization in the complex environment of the cell. In particular, the spatiotemporal distribution of membrane receptors is an active field of study due to its putative role in the regulation of signal transduction. The SNAP-tag is an intrinsically monovalent and highly specific genetic tag for attaching a fluorescent label to a protein of interest. Little information is currently available on the choice of optimal fluorescent dyes for single-molecule microscopy utilizing the SNAP-tag labeling system. We surveyed 6 green and 16 red excitable dyes for their suitability in single-molecule microscopy of SNAP-tag fusion proteins in live cells. We determined the nonspecific binding levels and photostability of these dye conjugates when bound to a SNAP-tag fused membrane protein in live cells. We found that only a limited subset of the dyes tested is suitable for single-molecule tracking microscopy. The results show that a careful choice of the dye to conjugate to the SNAP-substrate to label SNAP-tag fusion proteins is very important, as many dyes suffer from either rapid photobleaching or high nonspecific staining. These characteristics appear to be unpredictable, which motivated the need to perform the systematic survey presented here. We have developed a protocol for evaluating the best dyes, and for the conditions that we evaluated, we find that Dy 549 and CF 640 are the best choices tested for single-molecule tracking. Using an optimal dye pair, we also demonstrate the possibility of dual-color single-molecule imaging of SNAP-tag fusion proteins. This survey provides an overview of the photophysical and imaging properties of a range of SNAP-tag fluorescent substrates, enabling the selection of optimal dyes and conditions for single-molecule imaging of SNAP-tagged fusion proteins in eukaryotic cell lines.

Kinetics of enzymatic reactions in lipid membranes containing domains

An appreciable part of enzymes operating in vivo is associated with lipid membranes. The function of such enzymes can be influenced by the presence of domains containing proteins and/or composed of different lipids. The corresponding experimental model-system studies can be performed under well controlled conditions, e.g., on a planar supported lipid bilayer or surface-immobilized vesicles. To clarify what may happen in such systems, we propose general kinetic equations describing the enzyme-catalyzed substrate conversion occurring via the Michaelis-Menten (MM) mechanism on a membrane with domains which do not directly participate in reaction. For two generic situations when a relatively slow reaction takes place primarily in or outside domains, we take substrate saturation and lateral substrate-substrate interactions at domains into account and scrutinize the dependence of the reaction rate on the average substrate coverage. With increasing coverage, depending on the details, the reaction rate reaches saturation via an inflection point or monotonously as in the conventional MM case. In addition, we show analytically the types of reaction kinetics occurring primarily at domain boundaries. In the physically interesting situation when the domain growth is fast on the reaction time scale, the latter kinetics are far from conventional. The opposite situation when the reaction is fast and controlled by diffusion has been studied by using the Monte Carlo technique. The corresponding results indicate that the dependence of the reaction kinetics on the domain size may be weak.

Nanodomain organization of rhodopsin in native human and murine rod outer segment disc membranes

Biological membranes display distinct domains that organize membrane proteins and signaling molecules to facilitate efficient and reliable signaling. The organization of rhodopsin, a G protein-coupled receptor, in native rod outer segment disc membranes was investigated by atomic force microscopy. Atomic force microscopy revealed that rhodopsin is arranged into domains of variable size, which we refer to herein as nanodomains, in native membranes. Quantitative analysis of 150 disc membranes revealed that the physical properties of nanodomains are conserved in humans and mice and that the properties of individual disc membranes can be variable. Examining the variable properties of disc membranes revealed some of the factors contributing to the size of rod outer segment discs and the formation of nanodomains in the membrane. The diameter of rod outer segment discs was dependent on the number of rhodopsin molecules incorporated into the membrane but independent of the spatial density of rhodopsin. The number of nanodomains present in a single disc was also dependent on the number of rhodopsin molecules incorporated into the membrane. The size of the nanodomains was largely independent of the number or spatial density of rhodopsin in the membrane.