Understanding the cell: Future views of structural biology.

Determining the structure and mechanisms of all individual functional modules of cells at high molecular detail has often been seen as equal to understanding how cells work. Recent technical advances have led to a flush of high-resolution structures of various macromolecular machines, but despite this wealth of detailed information, our understanding of cellular function remains incomplete. Here, we discuss present-day limitations of structural biology and highlight novel technologies that may enable us to analyze molecular functions directly inside cells. We predict that the progression toward structural cell biology will involve a shift toward conceptualizing a 4D virtual reality of cells using digital twins. These will capture cellular segments in a highly enriched molecular detail, include dynamic changes, and facilitate simulations of molecular processes, leading to novel and experimentally testable predictions. Transferring biological questions into algorithms that learn from the existing wealth of data and explore novel solutions may ultimately unveil how cells work.

Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony.

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Human fetal brain self-organizes into long-term expanding organoids.

Human brain development involves an orchestrated, massive neural progenitor expansion while a multi-cellular tissue architecture is established. Continuously expanding organoids can be grown directly from multiple somatic tissues, yet to date, brain organoids can solely be established from pluripotent stem cells. Here, we show that healthy human fetal brain in vitro self-organizes into organoids (FeBOs), phenocopying aspects of in vivo cellular heterogeneity and complex organization. FeBOs can be expanded over long time periods. FeBO growth requires maintenance of tissue integrity, which ensures production of a tissue-like extracellular matrix (ECM) niche, ultimately endowing FeBO expansion. FeBO lines derived from different areas of the central nervous system (CNS), including dorsal and ventral forebrain, preserve their regional identity and allow to probe aspects of positional identity. Using CRISPR-Cas9, we showcase the generation of syngeneic mutant FeBO lines for the study of brain cancer. Taken together, FeBOs constitute a complementary CNS organoid platform.

Metabolic trade-offs constrain the cell size ratio in a nitrogen-fixing symbiosis.

Biological dinitrogen (N2) fixation is a key metabolic process exclusively performed by prokaryotes, some of which are symbiotic with eukaryotes. Species of the marine haptophyte algae Braarudosphaera bigelowii harbor the N2-fixing endosymbiotic cyanobacteria UCYN-A, which might be evolving organelle-like characteristics. We found that the size ratio between UCYN-A and their hosts is strikingly conserved across sublineages/species, which is consistent with the size relationships of organelles in this symbiosis and other species. Metabolic modeling showed that this size relationship maximizes the coordinated growth rate based on trade-offs between resource acquisition and exchange. Our findings show that the size relationships of N2-fixing endosymbionts and organelles in unicellular eukaryotes are constrained by predictable metabolic underpinnings and that UCYN-A is, in many regards, functioning like a hypothetical N2-fixing organelle (or nitroplast).

The molecular architecture of the nuclear basket.

The nuclear pore complex (NPC) is the sole mediator of nucleocytoplasmic transport. Despite great advances in understanding its conserved core architecture, the peripheral regions can exhibit considerable variation within and between species. One such structure is the cage-like nuclear basket. Despite its crucial roles in mRNA surveillance and chromatin organization, an architectural understanding has remained elusive. Using in-cell cryo-electron tomography and subtomogram analysis, we explored the NPC's structural variations and the nuclear basket across fungi (yeast; S. cerevisiae), mammals (mouse; M. musculus), and protozoa (T. gondii). Using integrative structural modeling, we computed a model of the basket in yeast and mammals that revealed how a hub of nucleoporins (Nups) in the nuclear ring binds to basket-forming Mlp/Tpr proteins: the coiled-coil domains of Mlp/Tpr form the struts of the basket, while their unstructured termini constitute the basket distal densities, which potentially serve as a docking site for mRNA preprocessing before nucleocytoplasmic transport.

SARS-CoV-2-infection- and vaccine-induced antibody responses are long lasting with an initial waning phase followed by a stabilization phase.

It is thought that mRNA-based vaccine-induced immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) wanes quickly, based mostly on short-term studies. Here, we analyzed the kinetics and durability of the humoral responses to SARS-CoV-2 infection and vaccination using >8,000 longitudinal samples collected over a 3-year period in New York City. Upon primary immunization, participants with pre-existing immunity mounted higher antibody responses faster and achieved higher steady-state antibody titers than naive individuals. Antibody kinetics were characterized by two phases: an initial rapid decay, followed by a stabilization phase with very slow decay. Booster vaccination equalized the differences in antibody concentration between participants with and without hybrid immunity, but the peak antibody titers decreased with each successive antigen exposure. Breakthrough infections increased antibodies to similar titers as an additional vaccine dose in naive individuals. Our study provides strong evidence that SARS-CoV-2 antibody responses are long lasting, with initial waning followed by stabilization.

Reversible, tunable epigenetic silencing of TCF1 generates flexibility in the T cell memory decision.

The immune system encodes information about the severity of a pathogenic threat in the quantity and type of memory cells it forms. This encoding emerges from lymphocyte decisions to maintain or lose self-renewal and memory potential during a challenge. By tracking CD8+ T cells at the single-cell and clonal lineage level using time-resolved transcriptomics, quantitative live imaging, and an acute infection model, we find that T cells will maintain or lose memory potential early after antigen recognition. However, following pathogen clearance, T cells may regain memory potential if initially lost. Mechanistically, this flexibility is implemented by a stochastic cis-epigenetic switch that tunably and reversibly silences the memory regulator, TCF1, in response to stimulation. Mathematical modeling shows how this flexibility allows memory T cell numbers to scale robustly with pathogen virulence and immune response magnitudes. We propose that flexibility and stochasticity in cellular decisions ensure optimal immune responses against diverse threats.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle.

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Proximal termination generates a transcriptional state that determines the rate of establishment of Polycomb silencing.

The mechanisms and timescales controlling de novo establishment of chromatin-mediated transcriptional silencing by Polycomb repressive complex 2 (PRC2) are unclear. Here, we investigate PRC2 silencing at Arabidopsis FLOWERING LOCUS C (FLC), known to involve co-transcriptional RNA processing, histone demethylation activity, and PRC2 function, but so far not mechanistically connected. We develop and test a computational model describing proximal polyadenylation/termination mediated by the RNA-binding protein FCA that induces H3K4me1 removal by the histone demethylase FLD. H3K4me1 removal feeds back to reduce RNA polymerase II (RNA Pol II) processivity and thus enhance early termination, thereby repressing productive transcription. The model predicts that this transcription-coupled repression controls the level of transcriptional antagonism to PRC2 action. Thus, the effectiveness of this repression dictates the timescale for establishment of PRC2/H3K27me3 silencing. We experimentally validate these mechanistic model predictions, revealing that co-transcriptional processing sets the level of productive transcription at the locus, which then determines the rate of the ON-to-OFF switch to PRC2 silencing.

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation.

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

CellVis2: a conference on visualizing the molecular cell.

In January 2024, a targeted conference, 'CellVis2', was held at Scripps Research in La Jolla, USA, the second in a series designed to explore the promise, practices, roadblocks, and prospects of creating, visualizing, sharing, and communicating physical representations of entire biological cells at scales down to the atom.

Visualizing intermediate stages of viral membrane fusion by cryo-electron tomography.

Protein-mediated membrane fusion is the dynamic process where specialized protein machinery undergoes dramatic conformational changes that drive two membrane bilayers together, leading to lipid mixing and opening of a fusion pore between previously separate membrane-bound compartments. Membrane fusion is an essential stage of enveloped virus entry that results in viral genome delivery into host cells. Recent studies applying cryo-electron microscopy techniques in a time-resolved fashion provide unprecedented glimpses into the interaction of viral fusion proteins and membranes, revealing fusion intermediate states from the initiation of fusion to release of the viral genome. In combination with complementary structural, biophysical, and computation modeling approaches, these advances are shedding new light on the mechanics and dynamics of protein-mediated membrane fusion.

Real-time assessment of mitochondrial DNA heteroplasmy dynamics at the single-cell level.

Mitochondrial DNA (mtDNA) is present in multiple copies within cells and is required for mitochondrial ATP generation. Even within individual cells, mtDNA copies can differ in their sequence, a state known as heteroplasmy. The principles underlying dynamic changes in the degree of heteroplasmy remain incompletely understood, due to the inability to monitor this phenomenon in real time. Here, we employ mtDNA-based fluorescent markers, microfluidics, and automated cell tracking, to follow mtDNA variants in live heteroplasmic yeast populations at the single-cell level. This approach, in combination with direct mtDNA tracking and data-driven mathematical modeling reveals asymmetric partitioning of mtDNA copies during cell division, as well as limited mitochondrial fusion and fission frequencies, as critical driving forces for mtDNA variant segregation. Given that our approach also facilitates assessment of segregation between intact and mutant mtDNA, we anticipate that it will be instrumental in elucidating the mechanisms underlying the purifying selection of mtDNA.

Higher-order structure and proteoforms of co-occurring C4b-binding protein assemblies in human serum.

The complement is a conserved cascade that plays a central role in the innate immune system. To maintain a delicate equilibrium preventing excessive complement activation, complement inhibitors are essential. One of the major fluid-phase complement inhibitors is C4b-binding protein (C4BP). Human C4BP is a macromolecular glycoprotein composed of two distinct subunits, C4BPα and C4BPß. These associate with vitamin K-dependent protein S (ProS) forming an ensemble of co-occurring higher-order structures. Here, we characterize these C4BP assemblies. We resolve and quantify isoforms of purified human serum C4BP using distinct single-particle detection techniques: charge detection mass spectrometry, and mass photometry accompanied by high-speed atomic force microscopy. Combining cross-linking mass spectrometry, glycoproteomics, and structural modeling, we report comprehensive glycoproteoform profiles and full-length structural models of the endogenous C4BP assemblies, expanding knowledge of this key complement inhibitor's structure and composition. Finally, we reveal that an increased C4BPα to C4BPß ratio coincides with elevated C-reactive protein levels in patient plasma samples. This observation highlights C4BP isoform variation and affirms a distinct role of co-occurring C4BP assemblies upon acute phase inflammation.

Mechanistic regulation of planarian shape during growth and degrowth.

Adult planarians can grow when fed and degrow (shrink) when starved while maintaining their whole-body shape. It is unknown how the morphogens patterning the planarian axes are coordinated during feeding and starvation or how they modulate the necessary differential tissue growth or degrowth. Here, we investigate the dynamics of planarian shape together with a theoretical study of the mechanisms regulating whole-body proportions and shape. We found that the planarian body proportions scale isometrically following similar linear rates during growth and degrowth, but that fed worms are significantly wider than starved worms. By combining a descriptive model of planarian shape and size with a mechanistic model of anterior-posterior and medio-lateral signaling calibrated with a novel parameter optimization methodology, we theoretically demonstrate that the feedback loop between these positional information signals and the shape they control can regulate the planarian whole-body shape during growth. Furthermore, the computational model produced the correct shape and size dynamics during degrowth as a result of a predicted increase in apoptosis rate and pole signal during starvation. These results offer mechanistic insights into the dynamic regulation of whole-body morphologies.

Differential proliferation regulates multi-tissue morphogenesis during embryonic axial extension: integrating viscous modeling and experimental approaches.

A major challenge in biology is to understand how mechanical interactions and cellular behavior affect the shapes of tissues and embryo morphology. The extension of the neural tube and paraxial mesoderm, which form the spinal cord and musculoskeletal system, respectively, results in the elongated shape of the vertebrate embryonic body. Despite our understanding of how each of these tissues elongates independently of the others, the morphogenetic consequences of their simultaneous growth and mechanical interactions are still unclear. Our study investigates how differential growth, tissue biophysical properties and mechanical interactions affect embryonic morphogenesis during axial extension using a 2D multi-tissue continuum-based mathematical model. Our model captures the dynamics observed in vivo by time-lapse imaging of bird embryos, and reveals the underestimated influence of differential tissue proliferation rates. We confirmed this prediction in quail embryos by showing that decreasing the rate of cell proliferation in the paraxial mesoderm affects long-term tissue dynamics, and shaping of both the paraxial mesoderm and the neighboring neural tube. Overall, our work provides a new theoretical platform upon which to consider the long-term consequences of tissue differential growth and mechanical interactions on morphogenesis.

A transit-amplifying progenitor with biphasic behavior contributes to epidermal renewal.

Transit-amplifying (TA) cells are progenitors that undergo an amplification phase followed by transition into an extinction phase. A long postulated epidermal TA progenitor with biphasic behavior has not yet been experimentally observed in vivo. Here, we identify such a TA population using clonal analysis of Aspm-CreER genetic cell-marking in mice, which uncovers contribution to both homeostasis and injury repair of adult skin. This TA population is more frequently dividing than a Dlx1-CreER-marked long-term self-renewing (e.g. stem cell) population. Newly developed generalized birth-death modeling of long-term lineage tracing data shows that both TA progenitors and stem cells display neutral competition, but only the stem cells display neutral drift. The quantitative evolution of a nascent TA cell and its direct descendants shows that TA progenitors indeed amplify the basal layer before transition and that the homeostatic TA population is mostly in extinction phase. This model will be broadly useful for analyzing progenitors whose behavior changes with their clone age. This work identifies a long-missing class of non-self-renewing biphasic epidermal TA progenitors and has broad implications for understanding tissue renewal mechanisms.

Developing long bones respond to surrounding tissues by trans-pairing of periosteal osteoclasts and endocortical osteoblasts.

Developing long bones alter their shape while maintaining uniform cortical thickness via coordinated activity of bone-forming osteoblasts and bone-resorbing osteoclasts at periosteal and endosteal surfaces, a process we designate trans-pairing. Two types of trans-pairing shift cortical bone in opposite orientations: peri-forming trans-pairing (peri-t-p) increases bone marrow space and endo-forming trans-pairing (endo-t-p) decreases it, via paired activity of bone resorption and formation across the cortex. Here, we focused on endo-t-p in growing bones. Analysis of endo-t-p activity in the cortex of mouse fibulae revealed osteoclasts under the periosteum compressed by muscles and expression of RANKL in periosteal cells of the cambium layer. Furthermore, mature osteoblasts were localized on the endosteum, while preosteoblasts were at the periosteum and within cortical canals. X-ray tomographic microscopy revealed the presence of cortical canals more closely associated with endo- than with peri-t-p. Sciatic nerve transection followed by muscle atrophy and unloading induced circumferential endo-t-p with concomitant spread of cortical canals. Such canals likely supply the endosteum with preosteoblasts from the periosteum under endo-t-p, allowing bone shape to change in response to mechanical stress or nerve injury.

Historic obstacles and emerging opportunities in the field of developmental metabolism - lessons from Heidelberg.

The field of developmental metabolism is experiencing a technological revolution that is opening entirely new fields of inquiry. Advances in metabolomics, small-molecule sensors, single-cell RNA sequencing and computational modeling present new opportunities for exploring cell-specific and tissue-specific metabolic networks, interorgan metabolic communication, and gene-by-metabolite interactions in time and space. Together, these advances not only present a means by which developmental biologists can tackle questions that have challenged the field for centuries, but also present young scientists with opportunities to define new areas of inquiry. These emerging frontiers of developmental metabolism were at the center of a highly interactive 2023 EMBO workshop 'Developmental metabolism: flows of energy, matter, and information'. Here, we summarize key discussions from this forum, emphasizing modern developmental biology's challenges and opportunities.

A mutation in F-actin polymerization factor suppresses the distal arthrogryposis type 5 PIEZO2 pathogenic variant in Caenorhabditis elegans.

The mechanosensitive PIEZO channel family has been linked to over 26 disorders and diseases. Although progress has been made in understanding these channels at the structural and functional levels, the underlying mechanisms of PIEZO-associated diseases remain elusive. In this study, we engineered four PIEZO-based disease models using CRISPR/Cas9 gene editing. We performed an unbiased chemical mutagen-based genetic suppressor screen to identify putative suppressors of a conserved gain-of-function variant pezo-1[R2405P] that in human PIEZO2 causes distal arthrogryposis type 5 (DA5; p. R2718P). Electrophysiological analyses indicate that pezo-1(R2405P) is a gain-of-function allele. Using genomic mapping and whole-genome sequencing approaches, we identified a candidate suppressor allele in the C. elegans gene gex-3. This gene is an ortholog of human NCKAP1 (NCK-associated protein 1), a subunit of the Wiskott-Aldrich syndrome protein (WASP)-verprolin homologous protein (WAVE/SCAR) complex, which regulates F-actin polymerization. Depletion of gex-3 by RNAi, or with the suppressor allele gex-3(av259[L353F]), significantly increased brood size and ovulation rate, as well as alleviating the crushed oocyte phenotype of the pezo-1(R2405P) mutant. Expression of GEX-3 in the soma is required to rescue the brood size defects in pezo-1(R2405P) animals. Actin organization and orientation were disrupted and distorted in the pezo-1 mutants. Mutation of gex-3(L353F) partially alleviated these defects. The identification of gex-3 as a suppressor of the pathogenic variant pezo-1(R2405P) suggests that the PIEZO coordinates with the cytoskeleton regulator to maintain the F-actin network and provides insight into the molecular mechanisms of DA5 and other PIEZO-associated diseases.