Markov Random Field Modeling of the Spatial Distribution of Proteins on Cell Membranes

Cell membranes display a range of receptors that bind ligands and activate signaling pathways. Signaling is characterized by dramatic changes in membrane molecular topography, including the co-clustering of receptors with signaling molecules and the segregation of other signaling molecules away from receptors. Electron microscopy of immunogold-labeled membranes is a critical technique to generate topographical information at the 5–10 nm resolution needed to understand how signaling complexes assemble and function. However, due to experimental limitations, only two molecular species can usually be labeled at a time. A formidable challenge is to integrate experimental data across multiple experiments where there are from 10 to 100 different proteins and lipids of interest and only the positions of two species can be observed simultaneously. As a solution, we propose the use of Markov random field (MRF) modeling to reconstruct the distribution of multiple cell membrane constituents from pair-wise data sets. MRFs are a powerful mathematical formalism for modeling correlations between states associated with neighboring sites in spatial lattices. The presence or absence of a protein of a specific type at a point on the cell membrane is a state. Since only two protein types can be observed, i.e., those bound to particles, and the rest cannot be observed, the problem is one of deducing the conditional distribution of a MRF with unobservable (hidden) states. Here, we develop a multiscale MRF model and use mathematical programming techniques to infer the conditional distribution of a MRF for proteins of three types from observations showing the spatial relationships between only two types. Application to synthesized data shows that the spatial distributions of three proteins can be reliably estimated. Application to experimental data provides the first maps of the spatial relationship between groups of three different signaling molecules. The work is an important step toward a more complete understanding of membrane spatial organization and dynamics during signaling.     

Clustering of membrane raft proteins by the actin cytoskeleton

Cell membranes are laterally organized into functionally discrete domains that include the cholesterol-dependent membrane "rafts." However, how membrane domains are established and maintained remains unresolved and controversial but often requires the actin cytoskeleton. In this study, we used fluorescence resonance energy transfer to measure the role of the actin cytoskeleton in the co-clustering of membrane raft-associated fluorescent proteins (FPs) and FPs targeted to the nonraft membrane fraction. By fitting the fluorescence resonance energy transfer data to an isothermal binding equation, we observed a specific co-clustering of raft-associated donor and acceptor probes that was sensitive to latrunculin B (Lat B), which disrupts the actin cytoskeleton. Conversely, treating with jasplakinolide to enhance actin polymerization increased co-clustering of the raft-associated FPs over that of the nonraft probes. We also observed by immunoblotting experiments that the actin-dependent co-clustering coincided with regulation of the raft-associated Src family kinase Lck. Specifically, Lat B decreased the phosphorylation of the C-terminal regulatory tyrosine of Lck (Tyr505), and combining the Lat B with filipin further decreased the Tyr505 phosphorylation. Furthermore, the Lat B-dependent changes in Lck regulation required CD45 because no significant changes occurred in treated T cells lacking CD45 expression. These data define a role for the actin cytoskeleton in promoting co-clustering of raft-associated proteins and show that this property is important toward regulating raft-associated signaling proteins such as Lck.     

Quantitative lifetime unmixing of multiexponentially decaying fluorophores using single-frequency fluorescence lifetime imaging microscopy

Fluorescence lifetime imaging microscopy (FLIM) is a quantitative microscopy technique for imaging nanosecond decay times of fluorophores. In the case of frequency-domain FLIM, several methods have been described to resolve the relative abundance of two fluorescent species with different fluorescence decay times. Thus far, single-frequency FLIM methods generally have been limited to quantifying two species with monoexponential decay. However, multiexponential decays are the norm rather than the exception, especially for fluorescent proteins and biological samples. Here, we describe a novel method for determining the fractional contribution in each pixel of an image of a sample containing two (multiexponentially) decaying species using single-frequency FLIM. We demonstrate that this technique allows the unmixing of binary mixtures of two spectrally identical cyan or green fluorescent proteins, each with multiexponential decay. Furthermore, because of their spectral identity, quantitative images of the relative molecular abundance of these fluorescent proteins can be generated that are independent of the microscope light path. The method is rigorously tested using samples of known composition and applied to live cell microscopy using cells expressing multiple (multiexponentially decaying) fluorescent proteins.     

Correlation functions quantify super-resolution images and estimate apparent clustering due to over-counting

We present an analytical method using correlation functions to quantify clustering in super-resolution fluorescence localization images and electron microscopy images of static surfaces in two dimensions. We use this method to quantify how over-counting of labeled molecules contributes to apparent self-clustering and to calculate the effective lateral resolution of an image. This treatment applies to distributions of proteins and lipids in cell membranes, where there is significant interest in using electron microscopy and super-resolution fluorescence localization techniques to probe membrane heterogeneity. When images are quantified using pair auto-correlation functions, the magnitude of apparent clustering arising from over-counting varies inversely with the surface density of labeled molecules and does not depend on the number of times an average molecule is counted. In contrast, we demonstrate that over-counting does not give rise to apparent co-clustering in double label experiments when pair cross-correlation functions are measured. We apply our analytical method to quantify the distribution of the IgE receptor (FcεRI) on the plasma membranes of chemically fixed RBL-2H3 mast cells from images acquired using stochastic optical reconstruction microscopy (STORM/dSTORM) and scanning electron microscopy (SEM). We find that apparent clustering of FcεRI-bound IgE is dominated by over-counting labels on individual complexes when IgE is directly conjugated to organic fluorophores. We verify this observation by measuring pair cross-correlation functions between two distinguishably labeled pools of IgE-FcεRI on the cell surface using both imaging methods. After correcting for over-counting, we observe weak but significant self-clustering of IgE-FcεRI in fluorescence localization measurements, and no residual self-clustering as detected with SEM. We also apply this method to quantify IgE-FcεRI redistribution after deliberate clustering by crosslinking with two distinct trivalent ligands of defined architectures, and we evaluate contributions from both over-counting of labels and redistribution of proteins.     

Role of sphingolipids in the biogenesis of membrane domains

In recent years, a huge interest in sphingolipid- and cholesterol-enriched membrane domains has risen, after their involvement in fundamental membrane-associated events such as signal transmission, cell adhesion and lipid/protein sorting was postulated. Theoretical considerations and several experimental data suggest that sphingolipids play an important role in the biogenesis and function of domains. In fact, their physicochemical features, different from those of other membrane lipids, allow their interaction either with other sphingolipids or with other membrane components and external ligands. Owing to these features, sphingolipids may undergo segregation and represent a nucleation point for co-clustering with other lipids and proteins in a complex, functional domain. Moreover, sphingolipids confer dynamic properties on domains, a fundamental feature for the modulation of their postulated functions.     

Cytoskeletal modulation of lipid interactions regulates Lck kinase activity

The actin cytoskeleton promotes clustering of proteins associated with cholesterol-dependent rafts, but its effect on lipid interactions that form and maintain rafts is not understood. We addressed this question by determining the effect of disrupting the cytoskeleton on co-clustering of dihexadecyl-(C(16))-anchored DiO and DiI, which co-enrich in ordered lipid environments such as rafts. Co-clustering was assayed by fluorescence resonance energy transfer (FRET) in labeled T cells, where rafts function in the phosphoregulation of the Src family kinase Lck. Our results show that probe co-clustering was sensitive to depolymerization of actin filaments with latrunculin B (Lat B), inhibition of myosin II with blebbistatin, and treatment with neomycin to sequester phosphatidylinositol 4,5-bisphosphate. Cytoskeletal effects on lipid interactions were not restricted to order-preferring label because co-clustering of C(16)-anchored DiO with didodecyl (C(12))-anchored DiI, which favors disordered lipids, was also reduced by Lat B and blebbistatin. Furthermore, conditions that disrupted probe co-clustering resulted in activation of Lck. These data show that the cytoskeleton globally modulates lipid interactions in the plasma membrane, and this property maintains rafts that function in Lck regulation.     

Quantitative analysis of synaptic phosphorylation and protein expression

The postsynaptic density (PSD) signaling machinery contains proteins with diverse functions. Brain region-specific variations in PSD components mediate distinct physiological responses to synaptic activation. We have developed mass spectrometry-based methods to comprehensively compare both relative protein expression and phosphorylation status from proteins present in biochemical preparations of postsynaptic density. Using these methods, we determined the relative expression of 2159 proteins and 1564 phosphorylation sites in PSD preparations from murine cortex, midbrain, cerebellum, and hippocampus. These experiments were conducted twice using independent biological replicates, which allowed us to assess the experimental and biological variability in this system. Concerning protein expression, cluster analysis revealed that known functionally associated proteins display coordinated synaptic expression. Therefore, proteins identified as co-clustering with known protein complexes are prime candidates for assignment as previously unrecognized components. Concerning degree of phosphorylation, we observed more extensive phosphorylation sites on N-methyl-D-aspartate (NMDA) receptors than alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors, consistent with the central role of N-methyl-D-aspartate receptors in processing synaptic transmission patterns. Average kinase and phosphatase levels were highest in the hippocampus, correlating with a higher overall phosphopeptide abundance present in this brain region. These findings suggest that the hippocampus utilizes reversible protein phosphorylation to a greater extent than other brain regions when modifying synaptic strength.     

NSOM/QD-Based Direct Visualization of CD3-Induced and CD28-Enhanced Nanospatial Coclustering of TCR and Coreceptor in Nanodomains in T Cell Activation

Direct molecular imaging of nano-spatial relationship between T cell receptor (TCR)/CD3 and CD4 or CD8 co-receptor before and after activation of a primary T cell has not been reported. We have recently innovated application of near-field scanning optical microscopy (NSOM) and immune-labeling quantum dots (QD) to image Ag-specific TCR response during in vivo clonal expansion, and now up-graded the NSOM/QD-based nanotechnology through dipole-polarization and dual-color imaging. Using this imaging system scanning cell-membrane molecules at a best-optical lateral resolution, we demonstrated that CD3, CD4 or CD8 molecules were distinctly distributed as single QD-bound molecules or nano-clusters equivalent to 2–4 QD fluorescence-intensity/size on cell-membrane of un-stimulated primary T cells, and ∼6–10% of CD3 were co-clustering with CD4 or CD8 as 70–110 nm nano-clusters without forming nano-domains. The ligation of TCR/CD3 on CD4 or CD8 T cells led to CD3 nanoscale co-clustering or interaction with CD4 or CD8 co-receptors forming 200–500 nm nano-domains or >500 nm micro-domains. Such nano-spatial co-clustering of CD3 and CD4 or CD3 and CD8 appeared to be an intrinsic event of TCR/CD3 ligation, not purely limited to MHC engagement, and be driven by Lck phosphorylation. Importantly, CD28 co-stimulation remarkably enhanced TCR/CD3 nanoscale co-clustering or interaction with CD4 co-receptor within nano- or micro-domains on the membrane. In contrast, CD28 co-stimulation did not enhance CD8 clustering or CD3–CD8 co-clustering in nano-domains although it increased molecular number and density of CD3 clustering in the enlarged nano-domains. These nanoscale findings provide new insights into TCR/CD3 interaction with CD4 or CD8 co-receptor in T-cell activation.     

Endocytic Accessory Proteins Are Functionally Distinguished by Their Differential Effects on the Maturation of Clathrin-coated Pits

Diverse cargo molecules (i.e., receptors and ligand/receptor complexes) are taken into the cell by clathrin-mediated endocytosis (CME) utilizing a core machinery consisting of cargo-specific adaptors, clathrin and the GTPase dynamin. Numerous endocytic accessory proteins are also required, but their differential roles and functional hierarchy during CME are not yet understood. Here, we used a combination of quantitative live-cell imaging by total internal reflection fluorescence microscopy (TIR-FM), and decomposition of the lifetime distributions of clathrin-coated pits (CCPs) to measure independent aspects of CCP dynamics, including the turnover of abortive and productive CCP species and their relative contributions. Capitalizing on the sensitivity of this assay, we have examined the effects of specific siRNA-mediated depletion of endocytic accessory proteins on CME progression. Of the 12 endocytic accessory proteins examined, we observed seven qualitatively different phenotypes upon protein depletion. From this data we derive a temporal hierarchy of protein function during early steps of CME. Our results support the idea that a subset of accessory proteins, which mediate coat assembly, membrane curvature, and cargo selection, can provide input into an endocytic restriction point/checkpoint mechanism that monitors CCP maturation.     

Co-clustering of Golgi complex and other cytoplasmic organelles to crescentic region of half-moon nuclei during apoptosis

Early apoptosis is defined by stereotypic morphological changes, especially evident in the nucleus, where chromatin condenses and compacts, and assumes a globular, half-moon or crescent-shaped morphology. Accumulating evidence suggests that cytoplasmic organelles such as mitochondria and the Golgi complex are major sites of integration of pro-apoptotic signaling. In this study, cytoplasmic organelles including Golgi complex, mitochondria, endosomes, lysosomes, and peroxisomes were shown to condense at the same unique region adjacent to the crescentic nucleus during a relatively early stage of apoptosis induced by staurosporine or other agents. The co-clustering phenomenon may be caused by shrinkage of cytoplasm during apoptosis although cytoskeletal markers actin and tubulin were not condensed and appeared excluded. These data suggest the co-clustering of cytoplasmic organelles plays an interesting role during the progression of the apoptotic process. It is possible that modification of pro-apoptotic proteins may arise as a result of the interplay of these cytoplasmic organelles.     

Modular organization of protein interaction networks

MOTIVATION: Accumulating evidence suggests that biological systems are composed of interacting, separable, functional modules. Identifying these modules is essential to understand the organization of biological systems. RESULT: In this paper, we present a framework to identify modules within biological networks. In this approach, the concept of degree is extended from the single vertex to the sub-graph, and a formal definition of module in a network is used. A new agglomerative algorithm was developed to identify modules from the network by combining the new module definition with the relative edge order generated by the Girvan-Newman (G-N) algorithm. A JAVA program, MoNet, was developed to implement the algorithm. Applying MoNet to the yeast core protein interaction network from the database of interacting proteins (DIP) identified 86 simple modules with sizes larger than three proteins. The modules obtained are significantly enriched in proteins with related biological process Gene Ontology terms. A comparison between the MoNet modules and modules defined by Radicchi et al. (2004) indicates that MoNet modules show stronger co-clustering of related genes and are more robust to ties in betweenness values. Further, the MoNet output retains the adjacent relationships between modules and allows the construction of an interaction web of modules providing insight regarding the relationships between different functional modules. Thus, MoNet provides an objective approach to understand the organization and interactions of biological processes in cellular systems. AVAILABILITY: MoNet is available upon request from the authors.     

Imaging Amyloid Fibrils within Cells Using a Se-Labelling Strategy

The process of aggregation leading to amyloid formation by peptides and proteins is associated with diseases ranging from systemic amyloidoses to neurodegenerative disorders such as Alzheimer's disease. A key question in understanding the link between amyloid formation and its pathological consequences is the ultrastructural localisation and morphological form of amyloid species within the cellular environment. The acquisition of such information has proven to be challenging, but we report here a novel approach that enables amyloid fibrils to be visualised directly within a cell. First, fibrils are assembled from selenium analogues of the sulfur-containing cysteine peptides, and then, atomic number contrast transmission electron microscopy is used to detect the selenium doped species selectively within the carbon-rich background of the cell. We demonstrate the power of this approach by imaging human monocyte-derived macrophage cells that have been exposed to fibrils from an amyloidogenic fragment of the disease-associated protein transthyretin. The ready incorporation of seleno-cysteine and methionine instead of their natural sulfur-containing analogues, a feature that is already commonly used in X-ray diffraction studies of proteins, suggests that this method can be used as a general strategy to image specific peptides and proteins within the cellular environment using electron microscopy.     

Stable evolutionary signal in a Yeast protein interaction network

Background  

The recently emerged protein interaction network paradigm can provide novel and important insights into the innerworkings of a cell. Yet, the heavy burden of both false positive and false negative protein-protein interaction data casts doubt on the broader usefulness of these interaction sets. Approaches focusing on one-protein-at-a-time have been powerfully employed to demonstrate the high degree of conservation of proteins participating in numerous interactions; here, we expand his 'node' focused paradigm to investigate the relative persistence of 'link' based evolutionary signals in a protein interaction network of S. cerevisiae and point out the value of this relatively untapped source of information.     

Soluble N-cadherin stimulates fibroblast growth factor receptor dependent neurite outgrowth and N-cadherin and the fibroblast growth factor receptor co-cluster in cells

A chimeric molecule consisting of the extracellular domain of the adhesion molecule, N-cadherin, fused to the Fc region of human IgG (NCAD-Fc) supports calcium-dependent cell adhesion and promotes neurite outgrowth following affinity-capture to a tissue culture substrate. When presented to cerebellar neurons as a soluble molecule, the NCAD-Fc stimulated neurite outgrowth in a manner equivalent to that seen for N-cadherin expressed as a cell surface glycoprotein. Neurons expressing a dominant-negative version of the fibroblast growth factor (FGF) receptor did not respond to soluble NCAD-Fc. In cells transfected with full-length N-cadherin and the FGF receptor, antibody-clustering of N-cadherin resulted in a co-clustering of the FGF receptor to discrete patches in the cell membrane. The data demonstrate that the ability of N-cadherin to stimulate neurite outgrowth can be dissociated from its ability to function as a substrate associated adhesion molecule. The N-cadherin and the FGF receptor co-clustering in cells provides a basis for the neurite outgrowth response stimulated by N-cadherin being dependent on FGF receptor function.     

Sequencing our way to more accurate community abundance

Over the last two decades, there has been a huge increase in our understanding of microbial diversity, structure and composition enabled by high-throughput sequencing technologies. Yet, it is unclear how the number of sequences translates to the number of cells or species within the community. In some cases, additional observational data may be required to ensure relative abundance patterns from sequence reads are biologically meaningful. The goal of DNA-based methods for biodiversity assessments is to obtain robust community abundance data, simultaneously, from environmental samples. In this issue of Molecular Ecology Resources, Pierella Karlusich et al. (2022) describe a new method for quantifying phytoplankton cell abundance. Using Tara Oceans data sets, the authors propose the photosynthetic gene psbO for reporting accurate relative abundance of the entire phytoplankton community from metagenomic data. The authors demonstrate higher correlations with traditional optical methods (including microscopy and flow cytometry), using their new method, improving upon molecular abundance assessments using multicopy marker genes. Furthermore, to facilitate application of their approach, the authors curated a psbO gene database for accessible taxonomic queries. This is an important step towards improving species abundance estimates from molecular data and eventually reporting of absolute species abundance, enhancing our understanding of community dynamics.     

Single-molecule pull-down for studying protein interactions

This protocol describes a single-molecule pull-down (SiMPull) assay for analyzing physiological protein complexes. The assay combines the conventional pull-down assay with single-molecule total internal reflection fluorescence (TIRF) microscopy and allows the probing of single macromolecular complexes directly from cell or tissue extracts. In this method, antibodies against the protein of interest are immobilized on a passivated microscope slide. When cell extracts are applied, the surface-tethered antibody captures the protein together with its physiological interaction partners. After washing away the unbound components, single-molecule fluorescence microscopy is used to probe the pulled-down proteins. Captured proteins are visualized through genetically encoded fluorescent protein tags or through antibody labeling. Compared with western blot analysis, this ultrasensitive assay requires considerably less time and reagents and provides quantitative data. Furthermore, SiMPull can distinguish between multiple association states of the same protein. SiMPull is generally applicable to proteins from a variety of cellular contexts and to endogenous proteins. Starting with the cell extracts and passivated slides, the assay requires 1.5-2.5 h for data acquisition and analysis.     

Some biological features of Candida albicans mutants for genes coding fungal proteins containing the CFEM domain

Several biological features of Candida albicans genes (PGA10, RBT5 and CSA1) coding for putative polypeptide species belonging to a subset of fungal proteins containing an eight-cysteine domain referred as common in several fungal extracellular membrane (CFEM) are described. The deletion of these genes resulted in a cascade of pleiotropic effects. Thus, mutant strains exhibited higher cell surface hydrophobicity levels and an increased ability to bind to inert or biological substrates. Confocal scanning laser microscopy using concanavalin A-Alexafluor 488 (which binds to mannose and glucose residues) and FUN-1 (a cytoplasmic fluorescent probe for cell viability) dyes showed that mutant strains formed thinner and more fragile biofilms. These apparently contained lower quantities of extracellular matrix material and less metabolically active cells than their parental strain counterpart, although the relative percentage of mycelial forms was similar in all cases. The cell surface of C. albicans strains harbouring deletions for genes coding CFEM-domain proteins appeared to be severely altered according to atomic force microscopy observations. Assessment of the relative gene expression within individual C. albicans cells revealed that CFEM-coding genes were upregulated in mycelium, although these genes were shown not to affect virulence in animal models. Overall, this study has demonstrated that CFEM domain protein-encoding genes are pleiotropic, influencing cell surface characteristics and biofilm formation.     

An evaluation of electrophoresis as a taxonomic method using comparative data from the macropodidae (Marsupialia)

Comparative electrophoretic data collected for seven proteins from thirty-four species of marsupial have been examined to define the advantages and shortcomings of electrophoresis as a taxonomic method. The problems of relative rates of change, convergence and parallelism are examined and the range of taxonomic ranks which can usefully be studied with this technique is determined. The relative merits of this type of data, when combined with various methods of taxonomic analysis, are discussed.     

Site localization of membrane-bound proteins on whole cell level using atomic force microscopy

This study presents molecular recognition method, which is based on specific force measurements between modified AFM (atomic force microscopy) tip and mammalian cell. The presented method allows recognition of specific cell surface proteins and receptor sites by nanometer accuracy level. Here we demonstrate specific recognition of membrane-bound Osteopontin (OPN) sites on preosteogenic cell membrane. By merging specific force detection map of the proteins and topography image of the cell, we create a new image (recognition image), which demonstrates the exact locations of the proteins relative to the cell membrane. The recognition results indicate the strong affinity between the modified tip and the target molecules, therefore, it enables the use of an AFM as a remarkable nanoscale tracking tool on the whole cell level.     

Biochemical and immunocytochemical analysis of vitellogenesis in the olive fruit fly Dacus (bactrocera) oleae (Diptera: Tephritidae)

Synthesis and selective accumulation of the major yolk proteins in the developing oocytes of the species Dacus oleae (Diptera: Tephritidae) was studied biochemically and by immunoelectron microscopy. In the hemolymph of adult females, two yolk proteins precursors (or vitellogenins) have been detected. They each exhibit a similar molecular weight and isoelectric point to their respective mature yolk proteins (or vitellins), while electrophoretic analysis of their synthetic profile shows that their levels in the hemolymph increase rapidly during development. Immunogold electron microscopy of ovarian sections, revealed that the hemolymph vitellogenins reach the oocyte through enlarged inter-follicular spaces and demonstrated vitellogenin synthesis by the follicle cells of the vitellogenic follicles. The newly synthesized vitellogenins follow a distinct secretory pathway into these cells as compared to other components being synthesized at the same time (e.g. the vitelline envelope proteins), since they were found in secretory vesicles that appeared to be differentiated from those destined to participate in the vitelline envelope. The vitellogenin-containing vesicles exocytose their contents directionally into the follicle cell/vitelline envelope boundary, and subsequently the vitellogenins diffuse among the gaps of the forming vitelline envelope and reach the oocyte plasma membrane. Their internalization by the oocyte includes the formation of an endocytic complex consisting of coated pits, coated vesicles, endosomes, transitional yolk bodies, and finally mature yolk bodies, in which the storage of the vitellins and other yolk proteins occur. These results are discussed in relation to data obtained from other Dipteran species.