Assembly of Biomolecular Gigastructures and Visualization with the Vulkan Graphics API.

Building and displaying all-atom models of biomolecular structures with millions or billions of atoms, like virus particles or cells, remain a challenge due to the sheer size of the data, the required levels of automated building, and the visualization limits of today's graphics hardware. Based on concepts introduced with the CellPack program, we report new algorithms to create such large-scale models using an intermediate coarse-grained "pet representation" of biomolecules with 1/10th the normal size. Pet atoms are placed such that they optimally trace the surface of the original molecule with just ∼1/50th the original atom number and are joined with covalent bonds. Molecular dynamics simulations of pet molecules allow for efficient packing optimization, as well as the generation of realistic DNA/RNA conformations. This pet world can be expanded back to the all-atom representation to be explored and visualized with full details. Essential for the efficient interactive visualization of gigastructures is the use of multiple levels of detail (LODs), where distant molecules are drawn with a heavily reduced polygon count. We present a grid-based algorithm to create such LODs for all common molecular graphics styles (including ball-and-sticks, ribbons, and cartoons) that do not require monochrome molecules to hide LOD transitions. As a practical application, we built all-atom models of SARS-CoV-2, HIV, and an entire presynaptic bouton with 1 µm diameter and 3.6 billion atoms, using modular building blocks to significantly reduce GPU memory requirements through instancing. We employ the Vulkan graphics API to maximize performance on consumer grade hardware and describe how to use the mmCIF format to efficiently store such giant models. An implementation is available as part of the YASARA molecular modeling and simulation program from www.YASARA.org. The free YASARA View program can be used to explore the presented models, which can be downloaded from www.YASARA.org/petworld, a Creative Commons platform for sharing giant biomolecular structures.

A large-scale nanoscopy and biochemistry analysis of postsynaptic dendritic spines.

Dendritic spines, the postsynaptic compartments of excitatory neurotransmission, have different shapes classified from 'stubby' to 'mushroom-like'. Whereas mushroom spines are essential for adult brain function, stubby spines disappear during brain maturation. It is still unclear whether and how they differ in protein composition. To address this, we combined electron microscopy and quantitative biochemistry with super-resolution microscopy to annotate more than 47,000 spines for more than 100 synaptic targets. Surprisingly, mushroom and stubby spines have similar average protein copy numbers and topologies. However, an analysis of the correlation of each protein to the postsynaptic density mass, used as a marker of synaptic strength, showed substantially more significant results for the mushroom spines. Secretion and trafficking proteins correlated particularly poorly to the strength of stubby spines. This suggests that stubby spines are less likely to adequately respond to dynamic changes in synaptic transmission than mushroom spines, which possibly explains their loss during brain maturation.

A comparative analysis of the mobility of 45 proteins in the synaptic bouton.

Many proteins involved in synaptic transmission are well known, and their features, as their abundance or spatial distribution, have been analyzed in systematic studies. This has not been the case, however, for their mobility. To solve this, we analyzed the motion of 45 GFP-tagged synaptic proteins expressed in cultured hippocampal neurons, using fluorescence recovery after photobleaching, particle tracking, and modeling. We compared synaptic vesicle proteins, endo- and exocytosis cofactors, cytoskeleton components, and trafficking proteins. We found that movement was influenced by the protein association with synaptic vesicles, especially for membrane proteins. Surprisingly, protein mobility also correlated significantly with parameters as the protein lifetimes, or the nucleotide composition of their mRNAs. We then analyzed protein movement thoroughly, taking into account the spatial characteristics of the system. This resulted in a first visualization of overall protein motion in the synapse, which should enable future modeling studies of synaptic physiology.

Supersaturated proteins are enriched at synapses and underlie cell and tissue vulnerability in Alzheimer's disease.

Neurodegenerative disorders progress across the brain in characteristic spatio-temporal patterns. A better understanding of the factors underlying the specific cell and tissue vulnerability responsible for such patterns could help identify the molecular origins of these conditions. To investigate these factors, based on the observation that neurodegenerative disorders are closely associated with the presence of aberrant protein deposits, we made the hypothesis that the vulnerability of cells and tissues is associated to the overall levels of supersaturated proteins, which are those most metastable against aggregation. By analyzing single-cell transcriptomic and subcellular proteomics data on healthy brains of ages much younger than those typical of disease onset, we found that the most supersaturated proteins are enriched in cells and tissues that succumb first to neurodegeneration. Then, by focusing the analysis on a metastable subproteome specific to Alzheimer's disease, we show that it is possible to recapitulate the pattern of disease progression using data from healthy brains. We found that this metastable subproteome is significantly enriched for synaptic processes and mitochondrial energy metabolism, thus rendering the synaptic environment dangerous for aggregation. The present identification of protein supersaturation as a signature of cell and tissue vulnerability in neurodegenerative disorders could facilitate the search for effective treatments by providing clearer points of intervention.

The codon sequences predict protein lifetimes and other parameters of the protein life cycle in the mouse brain.

The homeostasis of the proteome depends on the tight regulation of the mRNA and protein abundances, of the translation rates, and of the protein lifetimes. Results from several studies on prokaryotes or eukaryotic cell cultures have suggested that protein homeostasis is connected to, and perhaps regulated by, the protein and the codon sequences. However, this has been little investigated for mammals in vivo. Moreover, the link between the coding sequences and one critical parameter, the protein lifetime, has remained largely unexplored, both in vivo and in vitro. We tested this in the mouse brain, and found that the percentages of amino acids and codons in the sequences could predict all of the homeostasis parameters with a precision approaching experimental measurements. A key predictive element was the wobble nucleotide. G-/C-ending codons correlated with higher protein lifetimes, protein abundances, mRNA abundances and translation rates than A-/U-ending codons. Modifying the proportions of G-/C-ending codons could tune these parameters in cell cultures, in a proof-of-principle experiment. We suggest that the coding sequences are strongly linked to protein homeostasis in vivo, albeit it still remains to be determined whether this relation is causal in nature.

Precisely measured protein lifetimes in the mouse brain reveal differences across tissues and subcellular fractions.

The turnover of brain proteins is critical for organism survival, and its perturbations are linked to pathology. Nevertheless, protein lifetimes have been difficult to obtain in vivo. They are readily measured in vitro by feeding cells with isotopically labeled amino acids, followed by mass spectrometry analyses. In vivo proteins are generated from at least two sources: labeled amino acids from the diet, and non-labeled amino acids from the degradation of pre-existing proteins. This renders measurements difficult. Here we solved this problem rigorously with a workflow that combines mouse in vivo isotopic labeling, mass spectrometry, and mathematical modeling. We also established several independent approaches to test and validate the results. This enabled us to measure the accurate lifetimes of ~3500 brain proteins. The high precision of our data provided a large set of biologically significant observations, including pathway-, organelle-, organ-, or cell-specific effects, along with a comprehensive catalog of extremely long-lived proteins (ELLPs).

Clathrin coat controls synaptic vesicle acidification by blocking vacuolar ATPase activity.

Newly-formed synaptic vesicles (SVs) are rapidly acidified by vacuolar adenosine triphosphatases (vATPases), generating a proton electrochemical gradient that drives neurotransmitter loading. Clathrin-mediated endocytosis is needed for the formation of new SVs, yet it is unclear when endocytosed vesicles acidify and refill at the synapse. Here, we isolated clathrin-coated vesicles (CCVs) from mouse brain to measure their acidification directly at the single vesicle level. We observed that the ATP-induced acidification of CCVs was strikingly reduced in comparison to SVs. Remarkably, when the coat was removed from CCVs, uncoated vesicles regained ATP-dependent acidification, demonstrating that CCVs contain the functional vATPase, yet its function is inhibited by the clathrin coat. Considering the known structures of the vATPase and clathrin coat, we propose a model in which the formation of the coat surrounds the vATPase and blocks its activity. Such inhibition is likely fundamental for the proper timing of SV refilling.

Aptamer Stainings for Super-resolution Microscopy.

Fluorescence microscopy is an invaluable tool to visualize molecules in their biological context with ease and flexibility. However, studies using conventional light microscopy have been limited to the resolution that light diffraction allows (i.e., ~200 nm). This limitation has been recently circumvented by several types of advanced fluorescence microscopy techniques, which have achieved resolutions of up to ~10 nm. The resulting enhanced imaging precision has helped to find important cellular details that were not visible using diffraction-limited instruments. However, it has also revealed that conventional stainings using large affinity tags, such as antibodies, are not accurate enough for these imaging techniques. Since aptamers are substantially smaller than antibodies, they could provide a real advantage in super-resolution imaging. Here we compare the live staining of transferrin receptors (TfnR) obtained with different fluorescently labeled affinity probes: aptamers, specific monoclonal antibodies, or the natural receptor ligand transferrin. We observed negligible differences between these staining strategies when imaging is performed with conventional light microscopy (i.e., laser scanning confocal microscopy). However, a clear superiority of the aptamer tag over antibodies became apparent in super-resolved images obtained with stimulated emission depletion (STED) microscopy.

Composition of isolated synaptic boutons reveals the amounts of vesicle trafficking proteins.

Synaptic vesicle recycling has long served as a model for the general mechanisms of cellular trafficking. We used an integrative approach, combining quantitative immunoblotting and mass spectrometry to determine protein numbers; electron microscopy to measure organelle numbers, sizes, and positions; and super-resolution fluorescence microscopy to localize the proteins. Using these data, we generated a three-dimensional model of an "average" synapse, displaying 300,000 proteins in atomic detail. The copy numbers of proteins involved in the same step of synaptic vesicle recycling correlated closely. In contrast, copy numbers varied over more than three orders of magnitude between steps, from about 150 copies for the endosomal fusion proteins to more than 20,000 for the exocytotic ones.

The highly processive kinesin-8, Kip3, switches microtubule protofilaments with a bias toward the left.

Kinesin-1 motor proteins walk parallel to the protofilament axes of microtubules as they step from one tubulin dimer to the next. Is protofilament tracking an inherent property of processive kinesin motors, like kinesin-1, and what are the structural determinants underlying protofilament tracking? To address these questions, we investigated the tracking properties of the processive kinesin-8, Kip3. Using in vitro gliding motility assays, we found that Kip3 rotates microtubules counterclockwise around their longitudinal axes with periodicities of ∼1 µm. These rotations indicate that the motors switch protofilaments with a bias toward the left. Molecular modeling suggests 1), that the protofilament switching may be due to kinesin-8 having a longer neck linker than kinesin-1, and 2), that the leftward bias is due the asymmetric geometry of the motor neck linker complex.

Quantitative proteomics of the Cav2 channel nano-environments in the mammalian brain.

Local Ca(2+) signaling occurring within nanometers of voltage-gated Ca(2+) (Cav) channels is crucial for CNS function, yet the molecular composition of Cav channel nano-environments is largely unresolved. Here, we used a proteomic strategy combining knockout-controlled multiepitope affinity purifications with high-resolution quantitative MS for comprehensive analysis of the molecular nano-environments of the Cav2 channel family in the whole rodent brain. The analysis shows that Cav2 channels, composed of pore-forming alpha1 and auxiliary beta subunits, are embedded into protein networks that may be assembled from a pool of approximately 200 proteins with distinct abundance, stability of assembly, and preference for the three Cav2 subtypes. The majority of these proteins have not previously been linked to Cav channels; about two-thirds are dedicated to the control of intracellular Ca(2+) concentration, including G protein-coupled receptor-mediated signaling, to activity-dependent cytoskeleton remodeling or Ca(2+)-dependent effector systems that comprise a high portion of the priming and release machinery of synaptic vesicles. The identified protein networks reflect the cellular processes that can be initiated by Cav2 channel activity and define the molecular framework for organization and operation of local Ca(2+) signaling by Cav2 channels in the brain.

Interplay between lipids and the proteinaceous membrane fusion machinery.

For membrane fusion to occur, opposed lipid bilayers initially establish a fusion pore, often followed by complete mixing of the fusing membranes. Contemporary views suggest that during fusion lipid bilayers are continuous passive platforms that are disrupted and remodeled by catalytic proteins. Some models propose that even the architecture and composition of the fusion pore might be dominated by proteins rather than lipids. Hence, lipids have no regulatory contribution to this process; they simply adapt their shape passively for filling space between otherwise autonomous protein machineries. However, an increasing number of experimental findings indicate that membrane fusion critically depends on a variety of lipids and lipid derivatives. Therefore, a purely proteocentric view describes fusion mechanisms insufficiently. Instead, lipids have functions probably at different levels, as (i) a general influence on the propensity of lipid bilayers to fuse, (ii) a role in recruiting exocytotic proteins to the plasma membrane, (iii) a role in organizing membrane domains for fusion and (iv) direct regulatory effects on fusion protein complexes. In this review we have made an attempt to bring together the large body of evidence supporting a major role for lipids in membrane fusion either directly or indirectly.

Anatomy and dynamics of a supramolecular membrane protein cluster.

Most plasmalemmal proteins organize in submicrometer-sized clusters whose architecture and dynamics are still enigmatic. With syntaxin 1 as an example, we applied a combination of far-field optical nanoscopy, biochemistry, fluorescence recovery after photobleaching (FRAP) analysis, and simulations to show that clustering can be explained by self-organization based on simple physical principles. On average, the syntaxin clusters exhibit a diameter of 50 to 60 nanometers and contain 75 densely crowded syntaxins that dynamically exchange with freely diffusing molecules. Self-association depends on weak homophilic protein-protein interactions. Simulations suggest that clustering immobilizes and conformationally constrains the molecules. Moreover, a balance between self-association and crowding-induced steric repulsions is sufficient to explain both the size and dynamics of syntaxin clusters and likely of many oligomerizing membrane proteins that form supramolecular structures.

Molecular anatomy of a trafficking organelle.

Membrane traffic in eukaryotic cells involves transport of vesicles that bud from a donor compartment and fuse with an acceptor compartment. Common principles of budding and fusion have emerged, and many of the proteins involved in these events are now known. However, a detailed picture of an entire trafficking organelle is not yet available. Using synaptic vesicles as a model, we have now determined the protein and lipid composition; measured vesicle size, density, and mass; calculated the average protein and lipid mass per vesicle; and determined the copy number of more than a dozen major constituents. A model has been constructed that integrates all quantitative data and includes structural models of abundant proteins. Synaptic vesicles are dominated by proteins, possess a surprising diversity of trafficking proteins, and, with the exception of the V-ATPase that is present in only one to two copies, contain numerous copies of proteins essential for membrane traffic and neurotransmitter uptake.