Far-red chemigenetic biosensors for multi-dimensional and super-resolved kinase activity imaging

Fluorescent biosensors revolutionized biomedical science by enabling the direct measurement of signaling activities in living cells, yet the current technology is limited in resolution and dimensionality. Here, we introduce highly sensitive chemigenetic kinase activity biosensors that combine the genetically encodable self-labeling protein tag HaloTag7 with bright far-red-emitting synthetic fluorophores. This technology enables five-color biosensor multiplexing, 4D activity imaging, and functional super-resolution imaging via stimulated emission depletion (STED) microscopy.

Harmonizing the Generation and Pre-publication Stewardship of FAIR Image data

Together with the molecular knowledge of genes and proteins, biological images promise to significantly enhance the scientific understanding of complex cellular systems and to advance predictive and personalized therapeutic products for human health. For this potential to be realized, quality-assured image data must be shared among labs at a global scale to be compared, pooled, and reanalyzed, thus unleashing untold potential beyond the original purpose for which the data was generated. There are two broad sets of requirements to enable image data sharing in the life sciences. One set of requirements is articulated in the companion White Paper entitled "Enabling Global Image Data Sharing in the Life Sciences," which is published in parallel and addresses the need to build the cyberinfrastructure for sharing the digital array data (arXiv:2401.13023 [q-bio.OT], https://doi.org/10.48550/arXiv.2401.13023). In this White Paper, we detail a broad set of requirements, which involves collecting, managing, presenting, and propagating contextual information essential to assess the quality, understand the content, interpret the scientific implications, and reuse image data in the context of the experimental details. We start by providing an overview of the main lessons learned to date through international community activities, which have recently made considerable progress toward generating community standard practices for imaging Quality Control (QC) and metadata. We then provide a clear set of recommendations for amplifying this work. The driving goal is to address remaining challenges, and democratize access to common practices and tools for a spectrum of biomedical researchers, regardless of their expertise, access to resources, and geographical location.

Obesity causes mitochondrial fragmentation and dysfunction in white adipocytes due to RalA activation

Mitochondrial dysfunction is a characteristic trait of human and rodent obesity, insulin resistance and fatty liver disease. Here we show that high-fat diet (HFD) feeding causes mitochondrial fragmentation in inguinal white adipocytes from male mice, leading to reduced oxidative capacity by a process dependent on the small GTPase RalA. RalA expression and activity are increased in white adipocytes after HFD. Targeted deletion of RalA in white adipocytes prevents fragmentation of mitochondria and diminishes HFD-induced weight gain by increasing fatty acid oxidation. Mechanistically, RalA increases fission in adipocytes by reversing the inhibitory Ser637 phosphorylation of the fission protein Drp1, leading to more mitochondrial fragmentation. Adipose tissue expression of the human homolog of Drp1, DNM1L, is positively correlated with obesity and insulin resistance. Thus, chronic activation of RalA plays a key role in repressing energy expenditure in obese adipose tissue by shifting the balance of mitochondrial dynamics toward excessive fission, contributing to weight gain and metabolic dysfunction.

In situ tumour arrays reveal early environmental control of cancer immunity

The immune phenotype of a tumour is a key predictor of its response to immunotherapy^1-4. Patients who respond to checkpoint blockade generally present with immune-inflamed^5-7 tumours that are highly infiltrated by T cells. However, not all inflamed tumours respond to therapy, and even lower response rates occur among tumours that lack T cells (immune desert) or that spatially exclude T cells to the periphery of the tumour lesion (immune excluded)⁸. Despite the importance of these tumour immune phenotypes in patients, little is known about their development, heterogeneity or dynamics owing to the technical difficulty of tracking these features in situ. Here we introduce skin tumour array by microporation (STAMP)-a preclinical approach that combines high-throughput time-lapse imaging with next-generation sequencing of tumour arrays. Using STAMP, we followed the development of thousands of arrayed tumours in vivo to show that tumour immune phenotypes and outcomes vary between adjacent tumours and are controlled by local factors within the tumour microenvironment. Particularly, the recruitment of T cells by fibroblasts and monocytes into the tumour core was supportive of T cell cytotoxic activity and tumour rejection. Tumour immune phenotypes were dynamic over time and an early conversion to an immune-inflamed phenotype was predictive of spontaneous or therapy-induced tumour rejection. Thus, STAMP captures the dynamic relationships of the spatial, cellular and molecular components of tumour rejection and has the potential to translate therapeutic concepts into successful clinical strategies.

MitoTNT: Mitochondrial Temporal Network Tracking for 4D live-cell fluorescence microscopy data

Mitochondria form a network in the cell that rapidly changes through fission, fusion, and motility. Dysregulation of this four-dimensional (4D: x,y,z,time) network is implicated in numerous diseases ranging from cancer to neurodegeneration. While lattice light-sheet microscopy has recently made it possible to image mitochondria in 4D, quantitative analysis methods for the resulting datasets have been lacking. Here we present MitoTNT, the first-in-class software for Mitochondrial Temporal Network Tracking in 4D live-cell fluorescence microscopy data. MitoTNT uses spatial proximity and network topology to compute an optimal tracking assignment. To validate the accuracy of tracking, we created a reaction-diffusion simulation to model mitochondrial network motion and remodeling events. We found that our tracking is >90% accurate for ground-truth simulations and agrees well with published motility results for experimental data. We used MitoTNT to quantify 4D mitochondrial networks from human induced pluripotent stem cells. First, we characterized sub-fragment motility and analyzed network branch motion patterns. We revealed that the skeleton node motion is correlated along branch nodes and is uncorrelated in time. Second, we identified fission and fusion events with high spatiotemporal resolution. We found that mitochondrial skeleton nodes near the fission/fusion sites move nearly twice as fast as random skeleton nodes and that microtubules play a role in mediating selective fission/fusion. Finally, we developed graph-based transport simulations that model how material would distribute on experimentally measured mitochondrial temporal networks. We showed that pharmacological perturbations increase network reachability but decrease network resilience through a combination of altered mitochondrial fission/fusion dynamics and motility. MitoTNT's easy-to-use tracking module, interactive 4D visualization capability, and powerful post-tracking analyses aim at making temporal network tracking accessible to the wider mitochondria research community.

Branched actin networks are organized for asymmetric force production during clathrin-mediated endocytosis in mammalian cells

Actin assembly facilitates vesicle formation in several trafficking pathways, including clathrin-mediated endocytosis (CME). Interestingly, actin does not assemble at all CME sites in mammalian cells. How actin networks are organized with respect to mammalian CME sites and how assembly forces are harnessed, are not fully understood. Here, branched actin network geometry at CME sites was analyzed using three different advanced imaging approaches. When endocytic dynamics of unperturbed CME sites are compared, sites with actin assembly show a distinct signature, a delay between completion of coat expansion and vesicle scission, indicating that actin assembly occurs preferentially at stalled CME sites. In addition, N-WASP and the Arp2/3 complex are recruited to one side of CME sites, where they are positioned to stimulate asymmetric actin assembly and force production. We propose that actin assembles preferentially at stalled CME sites where it pulls vesicles into the cell asymmetrically, much as a bottle opener pulls off a bottle cap.

Mechanistic insights into actin force generation during vesicle formation from cryo-electron tomography

Actin assembly provides force for a multitude of cellular processes. Compared to actin-assembly-based force production during cell migration, relatively little is understood about how actin assembly generates pulling forces for vesicle formation. Here, cryo-electron tomography identified actin filament number, organization, and orientation during clathrin-mediated endocytosis in human SK-MEL-2 cells, showing that force generation is robust despite variance in network organization. Actin dynamics simulations incorporating a measured branch angle indicate that sufficient force to drive membrane internalization is generated through polymerization and that assembly is triggered from ∼4 founding "mother" filaments, consistent with tomography data. Hip1R actin filament anchoring points are present along the entire endocytic invagination, where simulations show that it is key to pulling force generation, and along the neck, where it targets filament growth and makes internalization more robust. Actin organization described here allowed direct translation of structure to mechanism with broad implications for other actin-driven processes.

Load adaptation by endocytic actin networks

Clathrin-mediated endocytosis (CME) robustness under elevated membrane tension is maintained by actin assembly-mediated force generation. However, whether more actin assembles at endocytic sites in response to increased load has not previously been investigated. Here actin network ultrastructure at CME sites was examined under low and high membrane tension. Actin and N-WASP spatial organization indicate that actin polymerization initiates at the base of clathrin-coated pits and that the network then grows away from the plasma membrane. Actin network height at individual CME sites was not coupled to coat shape, raising the possibility that local differences in mechanical load feed back on assembly. By manipulating membrane tension and Arp2/3 complex activity we tested the hypothesis that actin assembly at CME sites increases in response to elevated load. Indeed, in response to elevated membrane tension, actin grew higher, resulting in greater coverage of the clathrin coat, and CME slowed. When membrane tension was elevated and the Arp2/3 complex was inhibited, shallow clathrin-coated pits accumulated, indicating that this adaptive mechanism is especially crucial for coat curvature generation. We propose that actin assembly increases in response to increased load to ensure CME robustness over a range of plasma membrane tensions. [Media: see text] [Media: see text] [Media: see text] [Media: see text] [Media: see text].

Direct comparison of clathrin-mediated endocytosis in budding and fission yeast reveals conserved and evolvable features

Conserved proteins drive clathrin-mediated endocytosis (CME), which from yeast to humans involves a burst of actin assembly. To gain mechanistic insights into this process, we performed a side-by-side quantitative comparison of CME in two distantly related yeast species. Though endocytic protein abundance in S. pombe and S. cerevisiae is more similar than previously thought, membrane invagination speed and depth are two-fold greater in fission yeast. In both yeasts, accumulation of ~70 WASp molecules activates the Arp2/3 complex to drive membrane invagination. In contrast to budding yeast, WASp-mediated actin nucleation plays an essential role in fission yeast endocytosis. Genetics and live-cell imaging revealed core CME spatiodynamic similarities between the two yeasts, although the assembly of two zones of actin filaments is specific for fission yeast and not essential for CME. These studies identified conserved CME mechanisms and species-specific adaptations with broad implications that are expected to extend from yeast to humans.

ATP-dependent force generation and membrane scission by ESCRT-III and Vps4

The endosomal sorting complexes required for transport (ESCRTs) catalyze reverse-topology scission from the inner face of membrane necks in HIV budding, multivesicular endosome biogenesis, cytokinesis, and other pathways. We encapsulated ESCRT-III subunits Snf7, Vps24, and Vps2 and the AAA+ ATPase (adenosine triphosphatase) Vps4 in giant vesicles from which membrane nanotubes reflecting the correct topology of scission could be pulled. Upon ATP release by photo-uncaging, this system generated forces within the nanotubes that led to membrane scission in a manner dependent upon Vps4 catalytic activity and Vps4 coupling to the ESCRT-III proteins. Imaging of scission revealed Snf7 and Vps4 puncta within nanotubes whose presence followed ATP release, correlated with force generation and nanotube constriction, and preceded scission. These observations directly verify long-standing predictions that ATP-hydrolyzing assemblies of ESCRT-III and Vps4 sever membranes.

Bidirectional Control of Autophagy by BECN1 BARA Domain Dynamics

Membrane targeting of the BECN1-containing class III PI 3-kinase (PI3KC3) complexes is pivotal to the regulation of autophagy. The interaction of PI3KC3 complex II and its ubiquitously expressed inhibitor, Rubicon, was mapped to the first β sheet of the BECN1 BARA domain and the UVRAG BARA2 domain by hydrogen-deuterium exchange and cryo-EM. These data suggest that the BARA β sheet 1 unfolds to directly engage the membrane. This mechanism was confirmed using protein engineering, giant unilamellar vesicle assays, and molecular simulations. Using this mechanism, a BECN1 β sheet-1 derived peptide activates both PI3KC3 complexes I and II, while HIV-1 Nef inhibits complex II. These data reveal how BECN1 switches on and off PI3KC3 binding to membranes. The observations explain how PI3KC3 inhibition by Rubicon, activation by autophagy-inducing BECN1 peptides, and inhibition by HIV-1 Nef are mediated by the switchable ability of the BECN1 BARA domain to partially unfold and insert into membranes.

4D cell biology: big data image analytics and lattice light-sheet imaging reveal dynamics of clathrin-mediated endocytosis in stem cell-derived intestinal organoids

New methods in stem cell 3D organoid tissue culture, advanced imaging, and big data image analytics now allow tissue-scale 4D cell biology, but currently available analytical pipelines are inadequate for handing and analyzing the resulting gigabytes and terabytes of high-content imaging data. We expressed fluorescent protein fusions of clathrin and dynamin2 at endogenous levels in genome-edited human embryonic stem cells, which were differentiated into hESC-derived intestinal epithelial organoids. Lattice light-sheet imaging with adaptive optics (AO-LLSM) allowed us to image large volumes of these organoids (70 × 60 × 40 µm xyz) at 5.7 s/frame. We developed an open-source data analysis package termed pyLattice to process the resulting large (∼60 Gb) movie data sets and to track clathrin-mediated endocytosis (CME) events. CME tracks could be recorded from ∼35 cells at a time, resulting in ∼4000 processed tracks per movie. On the basis of their localization in the organoid, we classified CME tracks into apical, lateral, and basal events and found that CME dynamics is similar for all three classes, despite reported differences in membrane tension. pyLattice coupled with AO-LLSM makes possible quantitative high temporal and spatial resolution analysis of subcellular events within tissues.

It takes two transducins to activate the cGMP-phosphodiesterase 6 in retinal rods

Among cyclic nucleotide phosphodiesterases (PDEs), PDE6 is unique in serving as an effector enzyme in G protein-coupled signal transduction. In retinal rods and cones, PDE6 is membrane-bound and activated to hydrolyse its substrate, cGMP, by binding of two active G protein α-subunits (Gα*). To investigate the activation mechanism of mammalian rod PDE6, we have collected functional and structural data, and analysed them by reaction-diffusion simulations. Gα* titration of membrane-bound PDE6 reveals a strong functional asymmetry of the enzyme with respect to the affinity of Gα* for its two binding sites on membrane-bound PDE6 and the enzymatic activity of the intermediary 1 : 1 Gα* · PDE6 complex. Employing cGMP and its 8-bromo analogue as substrates, we find that Gα* · PDE6 forms with high affinity but has virtually no cGMP hydrolytic activity. To fully activate PDE6, it takes a second copy of Gα* which binds with lower affinity, forming Gα* · PDE6 · Gα*. Reaction-diffusion simulations show that the functional asymmetry of membrane-bound PDE6 constitutes a coincidence switch and explains the lack of G protein-related noise in visual signal transduction. The high local concentration of Gα* generated by a light-activated rhodopsin molecule efficiently activates PDE6, whereas the low density of spontaneously activated Gα* fails to activate the effector enzyme.

The dynamic Atg13-free conformation of the Atg1 EAT domain is required for phagophore expansion

Yeast macroautophagy begins with the de novo formation of a double-membrane phagophore at the preautophagosomal structure/phagophore assembly site (PAS), followed by its expansion into the autophagosome responsible for cargo engulfment. The kinase Atg1 is recruited to the PAS by Atg13 through interactions between the EAT domain of the former and the tMIM motif of the latter. Mass-spectrometry data have shown that, in the absence of Atg13, the EAT domain structure is strikingly dynamic, but the function of this Atg13-free dynamic state has been unclear. We used structure-based mutational analysis and quantitative and superresolution microscopy to show that Atg1 is present on autophagic puncta at, on average, twice the stoichiometry of Atg13. Moreover, Atg1 colocalizes with the expanding autophagosome in a manner dependent on Atg8 but not Atg13. We used isothermal titration calorimetry and crystal structure information to design an EAT domain mutant allele ATG1^DD that selectively perturbs the function of the Atg13-free state. Atg1^DD shows reduced PAS formation and does not support phagophore expansion, showing that the EAT domain has an essential function that is separate from its Atg13-dependent role in autophagy initiation.

Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations

Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25-109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.

Lipid-mediated PX-BAR domain recruitment couples local membrane constriction to endocytic vesicle fission

Clathrin-mediated endocytosis (CME) involves membrane-associated scaffolds of the bin-amphiphysin-rvs (BAR) domain protein family as well as the GTPase dynamin, and is accompanied and perhaps triggered by changes in local lipid composition. How protein recruitment, scaffold assembly and membrane deformation is spatiotemporally controlled and coupled to fission is poorly understood. We show by computational modelling and super-resolution imaging that phosphatidylinositol 3,4-bisphosphate [PI(3,4)P₂] synthesis within the clathrin-coated area of endocytic intermediates triggers selective recruitment of the PX-BAR domain protein SNX9, as a result of complex interactions of endocytic proteins competing for phospholipids. The specific architecture induces positioning of SNX9 at the invagination neck where its self-assembly regulates membrane constriction, thereby providing a template for dynamin fission. These data explain how lipid conversion at endocytic pits couples local membrane constriction to fission. Our work demonstrates how computational modelling and super-resolution imaging can be combined to unravel function and mechanisms of complex cellular processes.

The spatiotemporal regulation of membrane scaffolds recruitment and coupling between membrane deformation and fission in endocytosis are unclear. Here the authors show that lipid conversion at endocytic pits recruits SNX9, which couples local membrane constriction to fission in endocytosis.

Ballistic Anisotropic Magnetoresistance of Single-Atom Contacts

Anisotropic magnetoresistance, that is, the sensitivity of the electrical resistance of magnetic materials on the magnetization direction, is expected to be strongly enhanced in ballistic transport through nanoscale junctions. However, unambiguous experimental evidence of this effect is difficult to achieve. We utilize single-atom junctions to measure this ballistic anisotropic magnetoresistance (AMR). Single Co and Ir atoms are deposited on domains and domain walls of ferromagnetic Fe layers on W(110) to control their magnetization directions. They are contacted with nonmagnetic tips in a low-temperature scanning tunneling microscope to measure the junction conductances. Large changes of the magnetoresistance occur from the tunneling to the ballistic regime due to the competition of localized and delocalized d-orbitals, which are differently affected by spin-orbit coupling. This work shows that engineering the AMR at the single atom level is feasible.

Explicit spatiotemporal simulation of receptor-G protein coupling in rod cell disk membranes

Dim-light vision is mediated by retinal rod cells. Rhodopsin (R), a G-protein-coupled receptor, switches to its active form (R(∗)) in response to absorbing a single photon and activates multiple copies of the G-protein transducin (G) that trigger further downstream reactions of the phototransduction cascade. The classical assumption is that R and G are uniformly distributed and freely diffusing on disk membranes. Recent experimental findings have challenged this view by showing specific R architectures, including RG precomplexes, nonuniform R density, specific R arrangements, and immobile fractions of R. Here, we derive a physical model that describes the first steps of the photoactivation cascade in spatiotemporal detail and single-molecule resolution. The model was implemented in the ReaDDy software for particle-based reaction-diffusion simulations. Detailed kinetic in vitro experiments are used to parametrize the reaction rates and diffusion constants of R and G. Particle diffusion and G activation are then studied under different conditions of R-R interaction. It is found that the classical free-diffusion model is consistent with the available kinetic data. The existence of precomplexes between inactive R and G is only consistent with the data if these precomplexes are weak, with much larger dissociation rates than suggested elsewhere. Microarchitectures of R, such as dimer racks, would effectively immobilize R but have little impact on the diffusivity of G and on the overall amplification of the cascade at the level of the G protein.

Dynamical Organization of Syntaxin-1A at the Presynaptic Active Zone

Synaptic vesicle fusion is mediated by SNARE proteins forming in between synaptic vesicle (v-SNARE) and plasma membrane (t-SNARE), one of which is Syntaxin-1A. Although exocytosis mainly occurs at active zones, Syntaxin-1A appears to cover the entire neuronal membrane. By using STED super-resolution light microscopy and image analysis of Drosophila neuro-muscular junctions, we show that Syntaxin-1A clusters are more abundant and have an increased size at active zones. A computational particle-based model of syntaxin cluster formation and dynamics is developed. The model is parametrized to reproduce Syntaxin cluster-size distributions found by STED analysis, and successfully reproduces existing FRAP results. The model shows that the neuronal membrane is adjusted in a way to strike a balance between having most syntaxins stored in large clusters, while still keeping a mobile fraction of syntaxins free or in small clusters that can efficiently search the membrane or be traded between clusters. This balance is subtle and can be shifted toward almost no clustering and almost complete clustering by modifying the syntaxin interaction energy on the order of only 1 kBT. This capability appears to be exploited at active zones. The larger active-zone syntaxin clusters are more stable and provide regions of high docking and fusion capability, whereas the smaller clusters outside may serve as flexible reserve pool or sites of spontaneous ectopic release.

Higher-order architecture of rhodopsin in intact photoreceptors and its implication for phototransduction kinetics

The visual pigment rhodopsin belongs to the family of G protein-coupled receptors that can form higher oligomers. It is controversial whether rhodopsin forms oligomers and whether oligomers are functionally relevant. Here, we study rhodopsin organization in cryosections of dark-adapted mouse rod photoreceptors by cryoelectron tomography. We identify four hierarchical levels of organization. Rhodopsin forms dimers; at least ten dimers form a row. Rows form pairs (tracks) that are aligned parallel to the disk incisures. Particle-based simulation shows that the combination of tracks with fast precomplex formation, i.e. rapid association and dissociation between inactive rhodopsin and the G protein transducin, leads to kinetic trapping: rhodopsin first activates transducin from its own track, whereas recruitment of transducin from other tracks proceeds more slowly. The trap mechanism could produce uniform single-photon responses independent of rhodopsin lifetime. In general, tracks might provide a platform that coordinates the spatiotemporal interaction of signaling molecules.

Simulation tools for particle-based reaction-diffusion dynamics in continuous space

Particle-based reaction-diffusion algorithms facilitate the modeling of the diffusional motion of individual molecules and the reactions between them in cellular environments. A physically realistic model, depending on the system at hand and the questions asked, would require different levels of modeling detail such as particle diffusion, geometrical confinement, particle volume exclusion or particle-particle interaction potentials. Higher levels of detail usually correspond to increased number of parameters and higher computational cost. Certain systems however, require these investments to be modeled adequately. Here we present a review on the current field of particle-based reaction-diffusion software packages operating on continuous space. Four nested levels of modeling detail are identified that capture incrementing amount of detail. Their applicability to different biological questions is discussed, arching from straight diffusion simulations to sophisticated and expensive models that bridge towards coarse grained molecular dynamics.

ReaDDy--a software for particle-based reaction-diffusion dynamics in crowded cellular environments

We introduce the software package ReaDDy for simulation of detailed spatiotemporal mechanisms of dynamical processes in the cell, based on reaction-diffusion dynamics with particle resolution. In contrast to other particle-based reaction kinetics programs, ReaDDy supports particle interaction potentials. This permits effects such as space exclusion, molecular crowding and aggregation to be modeled. The biomolecules simulated can be represented as a sphere, or as a more complex geometry such as a domain structure or polymer chain. ReaDDy bridges the gap between small-scale but highly detailed molecular dynamics or Brownian dynamics simulations and large-scale but little-detailed reaction kinetics simulations. ReaDDy has a modular design that enables the exchange of the computing core by efficient platform-specific implementations or dynamical models that are different from Brownian dynamics.