Cryo-EM images of phase-separated lipid bilayer vesicles analyzed with a machine-learning approach.

Lateral lipid heterogeneity (i.e., raft formation) in biomembranes plays a functional role in living cells. Three-component mixtures of low- and high-melting lipids plus cholesterol offer a simplified experimental model for raft domains in which a liquid-disordered (Ld) phase coexists with a liquid-ordered (Lo) phase. Using such models, we recently showed that cryogenic electron microscopy (cryo-EM) can detect phase separation in lipid vesicles based on differences in bilayer thickness. However, the considerable noise within cryo-EM data poses a significant challenge for accurately determining the membrane phase state at high spatial resolution. To this end, we have developed an image-processing pipeline that utilizes machine learning (ML) to predict the bilayer phase in projection images of lipid vesicles. Importantly, the ML method exploits differences in both the thickness and molecular density of Lo compared to Ld, which leads to improved phase identification. To assess accuracy, we used artificial images of phase-separated lipid vesicles generated from all-atom molecular dynamics simulations of Lo and Ld phases. Synthetic ground-truth data sets mimicking a series of compositions along a tieline of Ld + Lo coexistence were created and then analyzed with various ML models. For all tieline compositions, we find that the ML approach can correctly identify the bilayer phase with >90% accuracy, thus providing a means to isolate the intensity profiles of coexisting Ld and Lo phases, as well as accurately determine domain-size distributions, number of domains, and phase-area fractions. The method described here provides a framework for characterizing nanoscopic lateral heterogeneities in membranes and paves the way for a more detailed understanding of raft properties in biological contexts.

Direct label-free imaging of nanodomains in biomimetic and biological membranes by cryogenic electron microscopy.

The nanoscale organization of biological membranes into structurally and compositionally distinct lateral domains is believed to be central to membrane function. The nature of this organization has remained elusive due to a lack of methods to directly probe nanoscopic membrane features. We show here that cryogenic electron microscopy (cryo-EM) can be used to directly image coexisting nanoscopic domains in synthetic and bioderived membranes without extrinsic probes. Analyzing a series of single-component liposomes composed of synthetic lipids of varying chain lengths, we demonstrate that cryo-EM can distinguish bilayer thickness differences as small as 0.5 Å, comparable to the resolution of small-angle scattering methods. Simulated images from computational models reveal that features in cryo-EM images result from a complex interplay between the atomic distribution normal to the plane of the bilayer and imaging parameters. Simulations of phase-separated bilayers were used to predict two sources of contrast between coexisting ordered and disordered phases within a single liposome, namely differences in membrane thickness and molecular density. We observe both sources of contrast in biomimetic membranes composed of saturated lipids, unsaturated lipids, and cholesterol. When extended to isolated mammalian plasma membranes, cryo-EM reveals similar nanoscale lateral heterogeneities. The methods reported here for direct, probe-free imaging of nanodomains in unperturbed membranes open new avenues for investigation of nanoscopic membrane organization.

Exploring the F-actin/CPEB3 interaction and its possible role in the molecular mechanism of long-term memory.

Dendritic spines are tiny membranous protrusions on the dendrites of neurons. Dendritic spines change shape in response to input signals, thereby strengthening the connections between neurons. The growth and stabilization of dendritic spines is thought to be essential for maintaining long-term memory. Actin cytoskeleton remodeling in spines is a key element of their formation and growth. More speculatively, the aggregation of CPEB3, a functional prion that binds RNA, has been reported to be involved in the maintenance of long-term memory. Here we study the interaction between actin and CPEB3 and propose a molecular model for the complex structure of CPEB3 and an actin filament (F-actin). The results of our computational modeling, including both energetic and structural analyses, are compared with novel data from peptide array experiments. Our model of the CPEB3/F-actin interaction suggests that F-actin potentially triggers the aggregation-prone structural transition of a short CPEB3 sequence by zipping it into a beta-hairpin form. We also propose that the CPEB3/F-actin interaction might be regulated by the SUMOylation of CPEB3, based on bioinformatic searches for potential SUMOylation sites as well as SUMO interacting motifs in CPEB3. On the basis of these results and the existing literature, we put forward a possible molecular mechanism underlying long-term memory that involves CPEB3's binding to actin, its aggregation, and its regulation by SUMOylation.

The role of the Arp2/3 complex in shaping the dynamics and structures of branched actomyosin networks.

Actomyosin networks give cells the ability to move and divide. These networks contract and expand while being driven by active energy-consuming processes such as motor protein walking and actin polymerization. Actin dynamics is also regulated by actin-binding proteins, such as the actin-related protein 2/3 (Arp2/3) complex. This complex generates branched filaments, thereby changing the overall organization of the network. In this work, the spatiotemporal patterns of dynamical actin assembly accompanying the branching-induced reorganization caused by Arp2/3 were studied using a computational model (mechanochemical dynamics of active networks [MEDYAN]); this model simulates actomyosin network dynamics as a result of chemical reactions whose rates are modulated by rapid mechanical equilibration. We show that branched actomyosin networks relax significantly more slowly than do unbranched networks. Also, branched networks undergo rare convulsive movements, "avalanches," that release strain in the network. These avalanches are associated with the more heterogeneous distribution of mechanically linked filaments displayed by branched networks. These far-from-equilibrium events arising from the marginal stability of growing actomyosin networks provide a possible mechanism of the "cytoquakes" recently seen in experiments.

Assemblies of calcium/calmodulin-dependent kinase II with actin and their dynamic regulation by calmodulin in dendritic spines.

Calcium/calmodulin-dependent kinase II (CaMKII) plays a key role in the plasticity of dendritic spines. Calcium signals cause calcium-calmodulin to activate CaMKII, which leads to remodeling of the actin filament (F-actin) network in the spine. We elucidate the mechanism of the remodeling by combining computer simulations with protein array experiments and electron microscopic imaging, to arrive at a structural model for the dodecameric complex of CaMKII with F-actin. The binding interface involves multiple domains of CaMKII. This structure explains the architecture of the micrometer-scale CaMKII/F-actin bundles arising from the multivalence of CaMKII. We also show that the regulatory domain of CaMKII may bind either calmodulin or F-actin, but not both. This frustration, along with the multipartite nature of the binding interface, allows calmodulin transiently to strip CaMKII from actin assemblies so that they can reorganize. This observation therefore provides a simple mechanism by which the structural dynamics of CaMKII establishes the link between calcium signaling and the morphological plasticity of dendritic spines.

Morphology of mitochondria in spatially restricted axons revealed by cryo-electron tomography.

Neurons project axons to local and distal sites and can display heterogeneous morphologies with limited physical dimensions that may influence the structure of large organelles such as mitochondria. Using cryo-electron tomography (cryo-ET), we characterized native environments within axons and presynaptic varicosities to examine whether spatial restrictions within these compartments influence the morphology of mitochondria. Segmented tomographic reconstructions revealed distinctive morphological characteristics of mitochondria residing at the narrowed boundary between presynaptic varicosities and axons with limited physical dimensions (approximately 80 nm), compared to mitochondria in nonspatially restricted environments. Furthermore, segmentation of the tomograms revealed discrete organizations between the inner and outer membranes, suggesting possible independent remodeling of each membrane in mitochondria at spatially restricted axonal/varicosity boundaries. Thus, cryo-ET of mitochondria within axonal subcompartments reveals that spatial restrictions do not obstruct mitochondria from residing within them, but limited available space can influence their gross morphology and the organization of the inner and outer membranes. These findings offer new perspectives on the influence of physical and spatial characteristics of cellular environments on mitochondrial morphology and highlight the potential for remarkable structural plasticity of mitochondria to adapt to spatial restrictions within axons.

Lipidomic atlas of mammalian cell membranes reveals hierarchical variation induced by culture conditions, subcellular membranes, and cell lineages.

Lipid membranes are ubiquitous biological organizers, required for structural and functional compartmentalization of the cell and sub-cellular organelles. Membranes in living cells are compositionally complex, comprising hundreds of dynamically regulated, distinct lipid species. Cellular physiology requires tight regulation of these lipidomic profiles to achieve proper membrane functionality. While some general features of tissue- and organelle-specific lipid complements have been identified, less is known about detailed lipidomic variations caused by cell-intrinsic or extrinsic factors. Here, we use shotgun lipidomics to report detailed, comprehensive lipidomes of a variety of cultured and primary mammalian membrane preparations to identify trends and sources of variation. Unbiased principle component analysis (PCA) shows clear separation between cultured and primary cells, with primary erythrocytes, synaptic membranes, and other mammalian tissue lipidomes sharply diverging from all cultured cell lines and also from one other. Most broadly, cultured cell membrane preparations were distinguished by their paucity of polyunsaturated lipids. Cultured mammalian cell lines were comparatively similar to one another, although we detected clear, highly reproducible lipidomic signatures of individual cell lines and plasma membrane (PM) isolations thereof. These measurements begin to establish a comprehensive lipidomic atlas of mammalian cells and tissues, identifying some major sources of variation. These observations will allow investigation of the regulation and functional significance of mammalian lipidomes in various contexts.

The ubiquitin ligase UBE4B regulates amyloid precursor protein ubiquitination, endosomal trafficking, and amyloid ß42 generation and secretion.

The extracellular accumulation of amyloid ß (Aß) fragments of amyloid precursor protein (APP) in brain parenchyma is a pathological hallmark of Alzheimer's disease (AD). APP can be cleaved into Aß on late endosomes/multivesicular bodies (MVBs). E3 ubiquitin ligases have been linked to Aß production, but specific E3 ligases associated with APP ubiquitination that may affect targeting of APP to endosomes have not yet been described. Using cultured cortical neurons isolated from rat pups, we reconstituted APP movement into the internal vesicles (ILVs) of MVBs. Loss of endosomal sorting complexes required for transport (ESCRT) components inhibited APP movement into ILVs and increased endosomal Aß42 generation, implying a requirement for APP ubiquitination. We identified an ESCRT-binding and APP-interacting endosomal E3 ubiquitin ligase, ubiquitination factor E4B (UBE4B) that regulates APP ubiquitination. Depleting UBE4B in neurons inhibited APP ubiquitination and internalization into MVBs, resulting in increased endosomal Aß42 levels and increased neuronal secretion of Aß42. When we examined AD brains, we found levels of the UBE4B-interacting ESCRT component, hepatocyte growth factor-regulated tyrosine kinase substrate (Hrs), were significantly decreased in AD brains. These data suggest that ESCRT components critical for membrane protein sorting in the endocytic pathway are altered in AD. These results indicate that the molecular machinery underlying endosomal trafficking of APP, including the ubiquitin ligase UBE4B, regulates Aß levels and may play an essential role in AD progression.

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism.

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.

On the Mechanism of Bilayer Separation by Extrusion, or Why Your LUVs Are Not Really Unilamellar.

Extrusion through porous filters is a widely used method for preparing biomimetic model membranes. Of primary importance in this approach is the efficient production of single bilayer (unilamellar) vesicles that eliminate the influence of interlamellar interactions and strictly define the bilayer surface area available to external reagents such as proteins. Submicroscopic vesicles produced using extrusion are widely assumed to be unilamellar, and large deviations from this assumption would impact interpretations from many model membrane experiments. Using three probe-free methods-small angle X-ray and neutron scattering and cryogenic electron microscopy-we report unambiguous evidence of extensive multilamellarity in extruded vesicles composed of neutral phosphatidylcholine lipids, including for the common case of neutral lipids dispersed in physiological buffer and extruded through 100-nm diameter pores. In such preparations, only ∼35% of lipids are externally accessible and this fraction is highly dependent on preparation conditions. Charged lipids promote unilamellarity as does decreasing solvent ionic strength, indicating the importance of electrostatic interactions in determining the lamellarity of extruded vesicles. Smaller extrusion pore sizes also robustly increase the fraction of unilamellar vesicles, suggesting a role for membrane bending. Taken together, these observations suggest a mechanistic model for extrusion, wherein the formation of unilamellar vesicles involves competition between bilayer bending and adhesion energies. The findings presented here have wide-ranging implications for the design and interpretation of model membrane studies, especially ensemble-averaged observations relying on the assumption of unilamellarity.

Distinct mechanisms enable inward or outward budding from late endosomes/multivesicular bodies.

Regulating the residence time of membrane proteins on the cell surface can modify their response to extracellular cues and allow for cellular adaptation in response to changing environmental conditions. The fate of membrane proteins that are internalized from the plasma membrane and arrive at the limiting membrane of the late endosome/multivesicular body (MVB) is dictated by whether they remain on the limiting membrane, bud into internal MVB vesicles, or bud outwardly from the membrane. The molecular details underlying the disposition of membrane proteins that transit this pathway and the mechanisms regulating these trafficking events are unclear. We established a cell-free system that reconstitutes budding of membrane protein cargo into internal MVB vesicles and onto vesicles that bud outwardly from the MVB membrane. Both budding reactions are cytosol-dependent and supported by Saccharomyces cerevisiae (yeast) cytosol. We observed that inward and outward budding from the MVB membrane are mechanistically distinct but may be linked, such that inhibition of inward budding triggers a re-routing of cargo from inward to outward budding vesicles, without affecting the number of vesicles that bud outwardly from MVBs.

Cytoskeletal-like Filaments of Ca2+-Calmodulin-Dependent Protein Kinase II Are Formed in a Regulated and Zn2+-Dependent Manner.

Ca2+-calmodulin-dependent protein kinase II (CaMKII) is highly abundant in neurons, where its concentration reaches that typically found for cytoskeletal proteins. Functional reasons for such a high concentration are not known, but given the multitude of known binding partners for CaMKII, a role as a scaffolding molecule has been proposed. In this report, we provide experimental evidence that demonstrates a novel structural role for CaMKII. We discovered that CaMKII forms filaments that can extend for several micrometers in the presence of certain divalent cations (Zn2+, Cd2+, and Cu2+) but not with others (Ca2+, Mg2+, Co2+, and Ni2+). Once formed, depleting the divalent ion concentration with chelators completely dissociated the filaments, and this process could be repeated by cyclic addition and removal of divalent ions. Using the crystal structure of the CaMKII holoenzyme, we computed an electrostatic potential map of the dodecameric complex to predict divalent ion binding sites. This analysis revealed a potential surface-exposed divalent ion binding site involving amino acids that also participate in calmodulin (CaM) binding and suggested CaM binding might inhibit formation of the filaments. As predicted, Ca2+/CaM binding both inhibited divalent ion-induced filament formation and could disassemble preformed filaments. Interestingly, CaMKII within the filaments retains the capacity to autophosphorylate; however, activity toward exogenous substrates is significantly decreased. Activity is restored upon filament disassembly. We compile our results with structural and mechanistic data from the literature to propose a model of Zn2+-mediated CaMKII filament formation, in which assembly and activity are further regulated by Ca2+/CaM.

Spatiotemporal Analysis of K-Ras Plasma Membrane Interactions Reveals Multiple High Order Homo-oligomeric Complexes.

Self-assembly of plasma membrane-associated Ras GTPases has major implications to the regulation of cell signaling. However, the structural basis of homo-oligomerization and the fractional distribution of oligomeric states remained undetermined. We have addressed these issues by deciphering the distribution of dimers and higher-order oligomers of K-Ras4B, the most frequently mutated Ras isoform in human cancers. We focused on the constitutively active G12V K-Ras and two of its variants, K101E and K101C/E107C, which respectively destabilize and stabilize oligomers. Using raster image correlation spectroscopy and number and brightness analysis combined with fluorescence recovery after photobleaching, fluorescence correlation spectroscopy and electron microscopy in live cells, we show that G12V K-Ras exists as a mixture of monomers, dimers and larger oligomers, while the K101E mutant is predominantly monomeric and K101C/E107C is dominated by oligomers. This observation demonstrates the ability of K-Ras to exist in multiple oligomeric states whose population can be altered by interfacial mutations. Using molecular modeling and simulations we further show that K-Ras uses two partially overlapping interfaces to form compositionally and topologically diverse oligomers. Our results thus provide the first detailed insight into the multiplicity, structure, and membrane organization of K-Ras homomers.

Protein recognition and selection through conformational and mutually induced fit.

Protein-protein interactions drive most every biological process, but in many instances the domains mediating recognition are disordered. How specificity in binding is attained in the absence of defined structure contrasts with well-established experimental and theoretical work describing ligand binding to protein. The signaling protein calmodulin presents a unique opportunity to investigate mechanisms for target recognition given that it interacts with several hundred different targets. By advancing coarse-grained computer simulations and experimental techniques, mechanistic insights were gained in defining the pathways leading to recognition and in how target selectivity can be achieved at the molecular level. A model requiring mutually induced conformational changes in both calmodulin and target proteins was necessary and broadly informs how proteins can achieve both high affinity and high specificity.

Relative Cosolute Size Influences the Kinetics of Protein-Protein Interactions.

Protein signaling occurs in crowded intracellular environments, and while high concentrations of macromolecules are postulated to modulate protein-protein interactions, analysis of their impact at each step of the reaction pathway has not been systematically addressed. Potential cosolute-induced alterations in target association are particularly important for a signaling molecule like calmodulin (CaM), where competition among >300 targets governs which pathways are selectively activated. To explore how high concentrations of cosolutes influence CaM-target affinity and kinetics, we methodically investigated each step of the CaM-target binding mechanism under crowded or osmolyte-rich environments mimicked by ficoll-70, dextran-10, and sucrose. All cosolutes stabilized compact conformers of CaM and modulated association kinetics by affecting diffusion and rates of conformational change; however, the results showed that differently sized molecules had variable effects to enhance or impede unique steps of the association pathway. On- and off-rates were modulated by all cosolutes in a compensatory fashion, producing little change in steady-state affinity. From this work insights were gained on how high concentrations of inert crowding agents and osmolytes fit into a kinetic framework to describe protein-protein interactions relevant for cellular signaling.

Neurogranin alters the structure and calcium binding properties of calmodulin.

Neurogranin (Ng) is a member of the IQ motif class of calmodulin (CaM)-binding proteins, and interactions with CaM are its only known biological function. In this report we demonstrate that the binding affinity of Ng for CaM is weakened by Ca(2+) but to a lesser extent (2-3-fold) than that previously suggested from qualitative observations. We also show that Ng induced a >10-fold decrease in the affinity of Ca(2+) binding to the C-terminal domain of CaM with an associated increase in the Ca(2+) dissociation rate. We also discovered a modest, but potentially important, increase in the cooperativity in Ca(2+) binding to the C-lobe of CaM in the presence of Ng, thus sharpening the threshold for the C-domain to become Ca(2+)-saturated. Domain mapping using synthetic peptides indicated that the IQ motif of Ng is a poor mimetic of the intact protein and that the acidic sequence just N-terminal to the IQ motif plays an important role in reproducing Ng-mediated decreases in the Ca(2+) binding affinity of CaM. Using NMR, full-length Ng was shown to make contacts largely with residues in the C-domain of CaM, although contacts were also detected in residues in the N-terminal domain. Together, our results can be consolidated into a model where Ng contacts residues in the N- and C-lobes of both apo- and Ca(2+)-bound CaM and that although Ca(2+) binding weakens Ng interactions with CaM, the most dramatic biochemical effect is the impact of Ng on Ca(2+) binding to the C-terminal lobe of CaM.

Domain contributions to signaling specificity differences between Ras-guanine nucleotide releasing factor (Ras-GRF) 1 and Ras-GRF2.

Ras-GRF1 (GRF1) and Ras-GRF2 (GRF2) constitute a family of similar calcium sensors that regulate synaptic plasticity. They are both guanine exchange factors that contain a very similar set of functional domains, including N-terminal pleckstrin homology, coiled-coil, and calmodulin-binding IQ domains and C-terminal Dbl homology Rac-activating domains, Ras-exchange motifs, and CDC25 Ras-activating domains. Nevertheless, they regulate different forms of synaptic plasticity. Although both GRF proteins transduce calcium signals emanating from NMDA-type glutamate receptors in the CA1 region of the hippocampus, GRF1 promotes LTD, whereas GRF2 promotes θ-burst stimulation-induced LTP (TBS-LTP). GRF1 can also mediate high frequency stimulation-induced LTP (HFS-LTP) in mice over 2-months of age, which involves calcium-permeable AMPA-type glutamate receptors. To add to our understanding of how proteins with similar domains can have different functions, WT and various chimeras between GRF1 and GRF2 proteins were tested for their abilities to reconstitute defective LTP and/or LTD in the CA1 hippocampus of Grf1/Grf2 double knock-out mice. These studies revealed a critical role for the GRF2 CDC25 domain in the induction of TBS-LTP by GRF proteins. In contrast, the N-terminal pleckstrin homology and/or coiled-coil domains of GRF1 are key to the induction of HFS-LTP by GRF proteins. Finally, the IQ motif of GRF1 determines whether a GRF protein can induce LTD. Overall, these findings show that for the three forms of synaptic plasticity that are regulated by GRF proteins in the CA1 hippocampus, specificity is encoded in only one or two domains, and a different set of domains for each form of synaptic plasticity.

Conformational frustration in calmodulin-target recognition.

Calmodulin (CaM) is a primary calcium (Ca(2+) )-signaling protein that specifically recognizes and activates highly diverse target proteins. We explored the molecular basis of target recognition of CaM with peptides representing the CaM-binding domains from two Ca(2+) -CaM-dependent kinases, CaMKI and CaMKII, by employing experimentally constrained molecular simulations. Detailed binding route analysis revealed that the two CaM target peptides, although similar in length and net charge, follow distinct routes that lead to a higher binding frustration in the CaM-CaMKII complex than in the CaM-CaMKI complex. We discovered that the molecular origin of the binding frustration is caused by intermolecular contacts formed with the C-domain of CaM that need to be broken before the formation of intermolecular contacts with the N-domain of CaM. We argue that the binding frustration is important for determining the kinetics of the recognition process of proteins involving large structural fluctuations.

Precisely tunable engineering of sub-30 nm monodisperse oligonucleotide nanoparticles.

Advancement of RNAi therapies is mainly hindered by the development of efficient delivery vehicles. The ability to create small size (<30 nm) oligonucleotide nanoparticles is essential for many aspects of the delivery process but is often overlooked. In this report, we describe diblock star polymers that can reproducibly complex double-stranded oligonucleotides into monodisperse nanoparticles with 15, 23, or 30 nm in diameter. The polymer-nucleic acid nanoparticles have a core-shell architecture with dense PEG brush coating. We characterized these nanoparticles using ITC, DLS, FRET, FCS, TIRF, and TEM. In addition to small size, these nanoparticles have neutral zeta-potentials, making the presented polymer architecture a very attractive platform for investigation of yet poorly studied polyplex size range for siRNA and antisense oligonucleotide delivery applications.

Phase separation in model lipid membranes investigated with cryogenic electron microscopy.

We describe a method for investigating lateral membrane heterogeneity using cryogenic electron microscopy (cryo-EM) images of liposomes. The method takes advantage of differences in the thickness and molecular density of ordered and disordered phases that are resolvable in phase contrast cryo-EM. Compared to biophysical techniques like FRET or neutron scattering that yield ensemble-averaged information, cryo-EM provides direct visualization of individual vesicles and can therefore reveal variability that would otherwise be obscured by averaging. Moreover, because the contrast mechanism involves inherent properties of the lipid phases themselves, no extrinsic probes are required. We explain and discuss various complementary analyses of spatially resolved thickness and intensity measurements that enable an assessment of the membrane's phase state. The method opens a window to nanodomain structure in synthetic and biological membranes that should lead to an improved understanding of lipid raft phenomena.