Frequent transitions in self-assembly across the evolution of a central metabolic enzyme.

Many enzymes assemble into homomeric protein complexes comprising multiple copies of one protein. Because structural form is usually assumed to follow function in biochemistry, these assemblies are thought to evolve because they provide some functional advantage. In many cases, however, no specific advantage is known and, in some cases, quaternary structure varies among orthologs. This has led to the proposition that self-assembly may instead vary neutrally within protein families. The extent of such variation has been difficult to ascertain because quaternary structure has until recently been difficult to measure on large scales. Here, we employ mass photometry, phylogenetics, and structural biology to interrogate the evolution of homo-oligomeric assembly across the entire phylogeny of prokaryotic citrate synthases - an enzyme with a highly conserved function. We discover a menagerie of different assembly types that come and go over the course of evolution, including cases of parallel evolution and reversions from complex to simple assemblies. Functional experiments in vitro and in vivo indicate that evolutionary transitions between different assemblies do not strongly influence enzyme catalysis. Our work suggests that enzymes can wander relatively freely through a large space of possible assemblies and demonstrates the power of characterizing structure-function relationships across entire phylogenies.

Emergence of fractal geometries in the evolution of a metabolic enzyme.

Fractals are patterns that are self-similar across multiple length-scales1. Macroscopic fractals are common in nature2-4; however, so far, molecular assembly into fractals is restricted to synthetic systems5-12. Here we report the discovery of a natural protein, citrate synthase from the cyanobacterium Synechococcus elongatus, which self-assembles into Sierpinski triangles. Using cryo-electron microscopy, we reveal how the fractal assembles from a hexameric building block. Although different stimuli modulate the formation of fractal complexes and these complexes can regulate the enzymatic activity of citrate synthase in vitro, the fractal may not serve a physiological function in vivo. We use ancestral sequence reconstruction to retrace how the citrate synthase fractal evolved from non-fractal precursors, and the results suggest it may have emerged as a harmless evolutionary accident. Our findings expand the space of possible protein complexes and demonstrate that intricate and regulatable assemblies can evolve in a single substitution.

State- and stimulus-specific dynamics of SMAD signaling determine fate decisions in individual cells.

SMAD-mediated signaling regulates apoptosis, cell cycle arrest, and epithelial-to-mesenchymal transition to safeguard tissue homeostasis. However, it remains elusive how the relatively simple pathway can determine such a broad range of cell fate decisions and how it differentiates between varying ligands. Here, we systematically investigate how SMAD-mediated responses are modulated by various ligands of the transforming growth factor ß (TGFß) family and compare these ligand responses in quiescent and proliferating MCF10A cells. We find that the nature of the phenotypic response is mainly determined by the proliferation status, with migration and cell cycle arrest being dominant in proliferating cells for all tested TGFß family ligands, whereas cell death is the major outcome in quiescent cells. In both quiescent and proliferating cells, the identity of the ligand modulates the strength of the phenotypic response proportional to the dynamics of induced SMAD nuclear-to-cytoplasmic translocation and, as a consequence, the corresponding gene expression changes. Interestingly, the proliferation state of a cell has little impact on the set of genes induced by SMAD signaling; instead, it modulates the relative cellular sensitivity to TGFß superfamily members. Taken together, diversity of SMAD-mediated responses is mediated by differing cellular states, which determine ligand sensitivity and phenotypic effects, while the pathway itself merely serves as a quantitative relay from the cell membrane to the nucleus.

Protocol for RNA-seq Expression Analysis in Yeast.

Genome-wide sequencing of RNA (RNA-seq) has become an inexpensive tool to gain key insights into cellular and disease mechanisms. Sample preparation and sequencing are streamlined and allow the acquisition of hundreds of gene expression profiles in a few days; however, in particular, data processing, curation, and analysis involve numerous steps that can be overwhelming to non-experts. Here, the sample preparation, sequencing, and data processing workflow for RNA-seq expression analysis in yeast is described. While this protocol covers only a small portion of the RNA-seq landscape, the principal workflow common to such experiments is described, allowing the reader to adapt the protocol where necessary. Graphic abstract: Basic workflow of RNA-seq expression analysis.

Cryo-EM structure of the cetacean morbillivirus nucleoprotein-RNA complex.

Cetacean morbillivirus (CeMV) is an emerging and highly infectious paramyxovirus that causes outbreaks in cetaceans and occasionally in pinnipeds, representing a major threat to biodiversity and conservation of endangered marine mammal populations in both hemispheres. As for all non-segmented, negative-sense, single-stranded RNA (ssRNA) viruses, the morbilliviral genome is enwrapped by thousands of nucleoprotein (N) protomers. Each bound to six ribonucleotides, N protomers assemble to form a helical ribonucleoprotein (RNP) complex that serves as scaffold for nucleocapsid formation and as template for viral replication and transcription. While the molecular details on RNP complexes elucidated in human measles virus (MeV) served as paradigm model for these processes in all members of the Morbillivirus genus, no structural information has been obtained from other morbilliviruses, nor has any CeMV structure been solved so far. We report the structure of the CeMV RNP complex, reconstituted in vitro upon binding of recombinant CeMV N to poly-adenine ssRNA hexamers and solved to 4.0 Å resolution by cryo-electron microscopy. In spite of the amino acid sequence similarity and consequently similar folding of the N protomer, the CeMV RNP complex exhibits different helical parameters as compared to previously reported MeV orthologs. The CeMV structure reveals exclusive interactions leading to more extensive protomer-RNA and protomer-protomer interfaces. We identified twelve residues, among those varying between CeMV strains, as putatively important for the stabilization of the RNP complex, which highlights the need to study the potential of CeMV N mutations that modulate nucleocapsid assembly to also affect viral phenotype and host adaptation.

Structural insights into photosystem II assembly.

Biogenesis of photosystem II (PSII), nature's water-splitting catalyst, is assisted by auxiliary proteins that form transient complexes with PSII components to facilitate stepwise assembly events. Using cryo-electron microscopy, we solved the structure of such a PSII assembly intermediate from Thermosynechococcus elongatus at 2.94 Å resolution. It contains three assembly factors (Psb27, Psb28 and Psb34) and provides detailed insights into their molecular function. Binding of Psb28 induces large conformational changes at the PSII acceptor side, which distort the binding pocket of the mobile quinone (QB) and replace the bicarbonate ligand of non-haem iron with glutamate, a structural motif found in reaction centres of non-oxygenic photosynthetic bacteria. These results reveal mechanisms that protect PSII from damage during biogenesis until water splitting is activated. Our structure further demonstrates how the PSII active site is prepared for the incorporation of the Mn4CaO5 cluster, which performs the unique water-splitting reaction.

High-resolution structure and biophysical characterization of the nucleocapsid phosphoprotein dimerization domain from the Covid-19 severe acute respiratory syndrome coronavirus 2.

Unprecedented by number of casualties and socio-economic burden occurring worldwide, the coronavirus disease 2019 (Covid-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the worst health crisis of this century. In order to develop adequate countermeasures against Covid-19, identification and structural characterization of suitable antiviral targets within the SARS-CoV-2 protein repertoire is urgently needed. The nucleocapsid phosphoprotein (N) is a multifunctional and highly immunogenic determinant of virulence and pathogenicity, whose main functions consist in oligomerizing and packaging the single-stranded RNA (ssRNA) viral genome. Here we report the structural and biophysical characterization of the SARS-CoV-2 N C-terminal domain (CTD), on which both N homo-oligomerization and ssRNA binding depend. Crystal structures solved at 1.44 Å and 1.36 Å resolution describe a rhombus-shape N CTD dimer, which stably exists in solution as validated by size-exclusion chromatography coupled to multi-angle light scattering and analytical ultracentrifugation. Differential scanning fluorimetry revealed moderate thermal stability and a tendency towards conformational change. Microscale thermophoresis demonstrated binding to a 7-bp SARS-CoV-2 genomic ssRNA fragment at micromolar affinity. Furthermore, a low-resolution preliminary model of the full-length SARS-CoV N in complex with ssRNA, obtained by cryo-electron microscopy, provides an initial understanding of self-associating and RNA binding functions exerted by the SARS-CoV-2 N.

Genetic interaction mapping informs integrative structure determination of protein complexes.

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on ~500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations.

Disentangling Pro-mitotic Signaling during Cell Cycle Progression using Time-Resolved Single-Cell Imaging.

Cells rely on input from extracellular growth factors to control their proliferation during development and adult homeostasis. Such mitogenic inputs are transmitted through multiple signaling pathways that synergize to precisely regulate cell cycle entry and progression. Although the architecture of these signaling networks has been characterized in molecular detail, their relative contribution, especially at later cell cycle stages, remains largely unexplored. By combining quantitative time-resolved measurements of fluorescent reporters in untransformed human cells with targeted pharmacological inhibitors and statistical analysis, we quantify epidermal growth factor (EGF)-induced signal processing in individual cells over time and dissect the dynamic contribution of downstream pathways. We define signaling features that encode information about extracellular ligand concentrations and critical time windows for inducing cell cycle transitions. We show that both extracellular signal-regulated kinase (ERK) and phosphatidylinositol 3-kinase (PI3K) activity are necessary for initial cell cycle entry, whereas only PI3K affects the duration of S phase at later stages of mitogenic signaling.

Cell-specific responses to the cytokine TGFß are determined by variability in protein levels.

The cytokine TGFß provides important information during embryonic development, adult tissue homeostasis, and regeneration. Alterations in the cellular response to TGFß are involved in severe human diseases. To understand how cells encode the extracellular input and transmit its information to elicit appropriate responses, we acquired quantitative time-resolved measurements of pathway activation at the single-cell level. We established dynamic time warping to quantitatively compare signaling dynamics of thousands of individual cells and described heterogeneous single-cell responses by mathematical modeling. Our combined experimental and theoretical study revealed that the response to a given dose of TGFß is determined cell specifically by the levels of defined signaling proteins. This heterogeneity in signaling protein expression leads to decomposition of cells into classes with qualitatively distinct signaling dynamics and phenotypic outcome. Negative feedback regulators promote heterogeneous signaling, as a SMAD7 knock-out specifically affected the signal duration in a subpopulation of cells. Taken together, we propose a quantitative framework that allows predicting and testing sources of cellular signaling heterogeneity.

Design and evaluation of a multimedia electronic patient record "oncoflow" with clinical workflow assistance for head and neck tumor therapy.

OBJECTIVES: The management of patient-specific information is a challenging task for surgeons and physicians because existing clinical information systems are insufficiently integrated into daily clinical routine and contained information entities are distributed across different proprietary databases. Thus, existing information is hardly usable for further electronic processing, workflow support or clinical studies. METHODS: A Web-based clinical information system has been developed that automatically imports patient-specific information from different information systems. The system is tailored to the existing workflow for the treatment of patients with head and neck cancer. In this paper, the clinical assistance functions and a quantitative as well as a qualitative system evaluation are presented. RESULTS: The information system has been deployed at a clinical site and is in use in daily clinical routine. Two evaluation studies show that the information integration, the structured information presentation in the Web browser and the assistance functions improve the physician's workflow. The studies also show that the usage of the new information system does not impair the time physicians need for a process step compared with the usage of the existing information system. CONCLUSIONS: Information integration is crucial for efficient workflow support in the clinic. The central access to information within a modern and structured user interface saves valuable time for the physician. The comprehensive database allows an instant usage of the existing information clinical workflow support or the conduction of trial studies.

Localization of the regulatory particle subunit Sem1 in the 26S proteasome.

The ubiquitin-proteasome system is responsible for regulated protein degradation in the cell with the 26S proteasome acting as its executive arm. The molecular architecture of this 2.5 MDa complex has been established recently, with the notable exception of the small acidic subunit Sem1. Here, we localize the C-terminal helix of Sem1 binding to the PCI domain of the subunit Rpn7 using cryo-electron microscopy single particle reconstruction of proteasomes purified from yeast cells with sem1 deletion. The approximate position of the N-terminal region of Sem1 bridging the cleft between Rpn7 and Rpn3 was inferred based on site-specific cross-linking data of the 26S proteasome. Our structural studies indicate that Sem1 can assume different conformations in different contexts, which supports the idea that Sem1 functions as a molecular glue stabilizing the Rpn3/Rpn7 heterodimer.

Near-atomic resolution structural model of the yeast 26S proteasome.

The 26S proteasome operates at the executive end of the ubiquitin-proteasome pathway. Here, we present a cryo-EM structure of the Saccharomyces cerevisiae 26S proteasome at a resolution of 7.4 Å or 6.7 Å (Fourier-Shell Correlation of 0.5 or 0.3, respectively). We used this map in conjunction with molecular dynamics-based flexible fitting to build a near-atomic resolution model of the holocomplex. The quality of the map allowed us to assign α-helices, the predominant secondary structure element of the regulatory particle subunits, throughout the entire map. We were able to determine the architecture of the Rpn8/Rpn11 heterodimer, which had hitherto remained elusive. The MPN domain of Rpn11 is positioned directly above the AAA-ATPase N-ring suggesting that Rpn11 deubiquitylates substrates immediately following commitment and prior to their unfolding by the AAA-ATPase module. The MPN domain of Rpn11 dimerizes with that of Rpn8 and the C-termini of both subunits form long helices, which are integral parts of a coiled-coil module. Together with the C-terminal helices of the six PCI-domain subunits they form a very large coiled-coil bundle, which appears to serve as a flexible anchoring device for all the lid subunits.

False discovery rate estimation for cross-linked peptides identified by mass spectrometry.

The mass spectrometric identification of chemically cross-linked peptides (CXMS) specifies spatial restraints of protein complexes; these values complement data obtained from common structure-determination techniques. Generic methods for determining false discovery rates of cross-linked peptide assignments are currently lacking, thus making data sets from CXMS studies inherently incomparable. Here we describe an automated target-decoy strategy and the software tool xProphet, which solve this problem for large multicomponent protein complexes.

Molecular architecture of the 26S proteasome holocomplex determined by an integrative approach.

The 26S proteasome is at the executive end of the ubiquitin-proteasome pathway for the controlled degradation of intracellular proteins. While the structure of its 20S core particle (CP) has been determined by X-ray crystallography, the structure of the 19S regulatory particle (RP), which recruits substrates, unfolds them, and translocates them to the CP for degradation, has remained elusive. Here, we describe the molecular architecture of the 26S holocomplex determined by an integrative approach based on data from cryoelectron microscopy, X-ray crystallography, residue-specific chemical cross-linking, and several proteomics techniques. The "lid" of the RP (consisting of Rpn3/5/6/7/8/9/11/12) is organized in a modular fashion. Rpn3/5/6/7/9/12 form a horseshoe-shaped heterohexamer, which connects to the CP and roofs the AAA-ATPase module, positioning the Rpn8/Rpn11 heterodimer close to its mouth. Rpn2 is rigid, supporting the lid, while Rpn1 is conformationally variable, positioned at the periphery of the ATPase ring. The ubiquitin receptors Rpn10 and Rpn13 are located in the distal part of the RP, indicating that they were recruited to the complex late in its evolution. The modular structure of the 26S proteasome provides insights into the sequence of events prior to the degradation of ubiquitylated substrates.

Expanding the chemical cross-linking toolbox by the use of multiple proteases and enrichment by size exclusion chromatography.

Chemical cross-linking in combination with mass spectrometric analysis offers the potential to obtain low-resolution structural information from proteins and protein complexes. Identification of peptides connected by a cross-link provides direct evidence for the physical interaction of amino acid side chains, information that can be used for computational modeling purposes. Despite impressive advances that were made in recent years, the number of experimentally observed cross-links still falls below the number of possible contacts of cross-linkable side chains within the span of the cross-linker. Here, we propose two complementary experimental strategies to expand cross-linking data sets. First, enrichment of cross-linked peptides by size exclusion chromatography selects cross-linked peptides based on their higher molecular mass, thereby depleting the majority of unmodified peptides present in proteolytic digests of cross-linked samples. Second, we demonstrate that the use of proteases in addition to trypsin, such as Asp-N, can additionally boost the number of observable cross-linking sites. The benefits of both SEC enrichment and multiprotease digests are demonstrated on a set of model proteins and the improved workflow is applied to the characterization of the 20S proteasome from rabbit and Schizosaccharomyces pombe.

Localization of the proteasomal ubiquitin receptors Rpn10 and Rpn13 by electron cryomicroscopy.

Two canonical subunits of the 26S proteasome, Rpn10 and Rpn13, function as ubiquitin (Ub) receptors. The mutual arrangement of these subunits--and all other non-ATPase subunits--in the regulatory particle is unknown. Using electron cryomicroscopy, we calculated difference maps between wild-type 26S proteasome from Saccharomyces cerevisiae and deletion mutants (rpn10Δ, rpn13Δ, and rpn10Δrpn13Δ). These maps allowed us to localize the two Ub receptors unambiguously. Rpn10 and Rpn13 mapped to the apical part of the 26S proteasome, above the N-terminal coiled coils of the AAA-ATPase heterodimers Rpt4/Rpt5 and Rpt1/Rpt2, respectively. On the basis of the mutual positions of Rpn10 and Rpn13, we propose a model for polyubiquitin binding to the 26S proteasome.

The proteasomal subunit Rpn6 is a molecular clamp holding the core and regulatory subcomplexes together.

Proteasomes execute the degradation of most cellular proteins. Although the 20S core particle (CP) has been studied in great detail, the structure of the 19S regulatory particle (RP), which prepares ubiquitylated substrates for degradation, has remained elusive. Here, we report the crystal structure of one of the RP subunits, Rpn6, and we describe its integration into the cryo-EM density map of the 26S holocomplex at 9.1 Å resolution. Rpn6 consists of an α-solenoid-like fold and a proteasome COP9/signalosome eIF3 (PCI) module in a right-handed suprahelical configuration. Highly conserved surface areas of Rpn6 interact with the conserved surfaces of the Pre8 (alpha2) and Rpt6 subunits from the alpha and ATPase rings, respectively. The structure suggests that Rpn6 has a pivotal role in stabilizing the otherwise weak interaction between the CP and the RP.

Structure of the 26S proteasome from Schizosaccharomyces pombe at subnanometer resolution.

The structure of the 26S proteasome from Schizosaccharomyces pombe has been determined to a resolution of 9.1 Å by cryoelectron microscopy and single particle analysis. In addition, chemical cross-linking in conjunction with mass spectrometry has been used to identify numerous residue pairs in close proximity to each other, providing an array of spatial restraints. Taken together these data clarify the topology of the AAA-ATPase module in the 19S regulatory particle and its spatial relationship to the α-ring of the 20S core particle. Image classification and variance analysis reveal a belt of high "activity" surrounding the AAA-ATPase module which is tentatively assigned to the reversible association of proteasome interacting proteins and the conformational heterogeneity among the particles. An integrated model is presented which sheds light on the early steps of protein degradation by the 26S complex.

Image restoration and spatial resolution in 7-tesla magnetic resonance imaging.

A good spatial resolution is essential for high precision segmentations of small structures in magnetic resonance images. However, any increase in the spatial resolution results in a decrease of the signal-to-noise ratio (SNR). In this article, this problem is addressed by a new image restoration technique that is used to partly compensate for the loss in SNR. Specifically, a two-stage hybrid image restoration procedure is proposed where the first stage is a Wiener wavelet filter for an initial denoising. The artifacts that will inevitably be produced by this step are subsequently reduced using a recent variant of anisotropic diffusion. The method is applied to magnetic resonance imaging data acquired on a 7-T magnetic resonance imaging scanner and compared with averaged multiple measurements of the same subject. It was found that the effect of image restoration procedure roughly corresponds to averaging across three repeated measurements.