Flattening the curve-How to get better results with small deep-mutational-scanning datasets.

Proteins are used in various biotechnological applications, often requiring the optimization of protein properties by introducing specific amino-acid exchanges. Deep mutational scanning (DMS) is an effective high-throughput method for evaluating the effects of these exchanges on protein function. DMS data can then inform the training of a neural network to predict the impact of mutations. Most approaches use some representation of the protein sequence for training and prediction. As proteins are characterized by complex structures and intricate residue interaction networks, directly providing structural information as input reduces the need to learn these features from the data. We introduce a method for encoding protein structures as stacked 2D contact maps, which capture residue interactions, their evolutionary conservation, and mutation-induced interaction changes. Furthermore, we explored techniques to augment neural network training performance on smaller DMS datasets. To validate our approach, we trained three neural network architectures originally used for image analysis on three DMS datasets, and we compared their performances with networks trained solely on protein sequences. The results confirm the effectiveness of the protein structure encoding in machine learning efforts on DMS data. Using structural representations as direct input to the networks, along with data augmentation and pretraining, significantly reduced demands on training data size and improved prediction performance, especially on smaller datasets, while performance on large datasets was on par with state-of-the-art sequence convolutional neural networks. The methods presented here have the potential to provide the same workflow as DMS without the experimental and financial burden of testing thousands of mutants. Additionally, we present an open-source, user-friendly software tool to make these data analysis techniques accessible, particularly to biotechnology and protein engineering researchers who wish to apply them to their mutagenesis data.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2.

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

Minimum information guidelines for experiments structurally characterizing intrinsically disordered protein regions.

An unambiguous description of an experiment, and the subsequent biological observation, is vital for accurate data interpretation. Minimum information guidelines define the fundamental complement of data that can support an unambiguous conclusion based on experimental observations. We present the Minimum Information About Disorder Experiments (MIADE) guidelines to define the parameters required for the wider scientific community to understand the findings of an experiment studying the structural properties of intrinsically disordered regions (IDRs). MIADE guidelines provide recommendations for data producers to describe the results of their experiments at source, for curators to annotate experimental data to community resources and for database developers maintaining community resources to disseminate the data. The MIADE guidelines will improve the interpretability of experimental results for data consumers, facilitate direct data submission, simplify data curation, improve data exchange among repositories and standardize the dissemination of the key metadata on an IDR experiment by IDR data sources.

The sex-specific metabolic signature of C57BL/6NRj mice during aging.

Due to intact reactive oxygen species homeostasis and glucose metabolism, C57BL/6NRj mice are especially suitable to study cellular alterations in metabolism. We applied Nuclear Magnetic resonance spectroscopy to analyze five different tissues of this mouse strain during aging and included female and male mice aged 3, 6, 12, and 24 months. Metabolite signatures allowed separation between the age groups in all tissues, and we identified the most prominently changing metabolites in female and male tissues. A refined analysis of individual metabolite levels during aging revealed an early onset of age-related changes at 6 months, sex-specific differences in the liver, and a biphasic pattern for various metabolites in the brain, heart, liver, and lung. In contrast, a linear decrease of amino acids was apparent in muscle tissues. Based on these results, we assume that age-related metabolic alterations happen at a comparably early aging state and are potentially associated with a metabolic switch. Moreover, identified differences between female and male tissues stress the importance of distinguishing between sexes when studying age-related changes and developing new treatment approaches. Besides, metabolomic features seem to be highly dependent on the genetic background of mouse strains.

Developmental dynamics of RNA translation in the human brain.

The precise regulation of gene expression is fundamental to neurodevelopment, plasticity and cognitive function. Although several studies have profiled transcription in the developing human brain, there is a gap in understanding of accompanying translational regulation. In this study, we performed ribosome profiling on 73 human prenatal and adult cortex samples. We characterized the translational regulation of annotated open reading frames (ORFs) and identified thousands of previously unknown translation events, including small ORFs that give rise to human-specific and/or brain-specific microproteins, many of which we independently verified using proteomics. Ribosome profiling in stem-cell-derived human neuronal cultures corroborated these findings and revealed that several neuronal activity-induced non-coding RNAs encode previously undescribed microproteins. Physicochemical analysis of brain microproteins identified a class of proteins that contain arginine-glycine-glycine (RGG) repeats and, thus, may be regulators of RNA metabolism. This resource expands the known translational landscape of the human brain and illuminates previously unknown brain-specific protein products.

Discovering molecular features of intrinsically disordered regions by using evolution for contrastive learning.

A major challenge to the characterization of intrinsically disordered regions (IDRs), which are widespread in the proteome, but relatively poorly understood, is the identification of molecular features that mediate functions of these regions, such as short motifs, amino acid repeats and physicochemical properties. Here, we introduce a proteome-scale feature discovery approach for IDRs. Our approach, which we call "reverse homology", exploits the principle that important functional features are conserved over evolution. We use this as a contrastive learning signal for deep learning: given a set of homologous IDRs, the neural network has to correctly choose a held-out homolog from another set of IDRs sampled randomly from the proteome. We pair reverse homology with a simple architecture and standard interpretation techniques, and show that the network learns conserved features of IDRs that can be interpreted as motifs, repeats, or bulk features like charge or amino acid propensities. We also show that our model can be used to produce visualizations of what residues and regions are most important to IDR function, generating hypotheses for uncharacterized IDRs. Our results suggest that feature discovery using unsupervised neural networks is a promising avenue to gain systematic insight into poorly understood protein sequences.

FOXO transcription factors differ in their dynamics and intra/intermolecular interactions.

Transcription factors play key roles in orchestrating a plethora of cellular mechanisms and controlling cellular homeostasis. Transcription factors share distinct DNA binding domains, which allows to group them into protein families. Among them, the Forkhead box O (FOXO) family contains transcription factors crucial for cellular homeostasis, longevity and response to stress. The dysregulation of FOXO signaling is linked to drug resistance in cancer therapy or cellular senescence, however, selective drugs targeting FOXOs are limited, thus knowledge about structure and dynamics of FOXO proteins is essential. Here, we provide an extensive study of structure and dynamics of all FOXO family members. We identify residues accounting for different dynamic and structural features. Furthermore, we show that the auto-inhibition of FOXO proteins by their C-terminal trans-activation domain is conserved throughout the family and that these interactions are not only possible intra-, but also inter-molecularly. This indicates a model in which FOXO transcription factors would modulate their activities by interacting mutually.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder.

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

O-Linked-N-Acetylglucosaminylation of the RNA-Binding Protein EWS N-Terminal Low Complexity Region Reduces Phase Separation and Enhances Condensate Dynamics.

Many membraneless organelles are thought to be biomolecular condensates formed by phase separation of proteins and other biopolymers. Post-translational modifications (PTMs) can impact protein phase separation behavior, although for many PTMs this aspect of their function is unknown. O-linked ß-D-N-acetylglucosaminylation (O-GlcNAcylation) is an abundant form of intracellular glycosylation whose roles in regulating biomolecular condensate assembly and dynamics have not been delineated. Using an in vitro approach, we found that O-GlcNAcylation reduces the phase separation propensity of the EWS N-terminal low complexity region (LCRN) under different conditions, including in the presence of the arginine- and glycine-rich RNA-binding domains (RBD). O-GlcNAcylation enhances fluorescence recovery after photobleaching (FRAP) within EWS LCRN condensates and causes the droplets to exhibit more liquid-like relaxation following fusion. Following extended incubation times, EWS LCRN+RBD condensates exhibit diminished FRAP, indicating a loss of fluidity, while condensates containing the O-GlcNAcylated LCRN do not. In HeLa cells, EWS is less O-GlcNAcylated following OGT knockdown, which correlates with its increased accumulation in a filter retardation assay. Relative to the human proteome, O-GlcNAcylated proteins are enriched with regions that are predicted to phase separate, suggesting a general role of O-GlcNAcylation in regulation of biomolecular condensates.

A weakened interface in the P182L variant of HSP27 associated with severe Charcot-Marie-Tooth neuropathy causes aberrant binding to interacting proteins.

HSP27 is a human molecular chaperone that forms large, dynamic oligomers and functions in many aspects of cellular homeostasis. Mutations in HSP27 cause Charcot-Marie-Tooth (CMT) disease, the most common inherited disorder of the peripheral nervous system. A particularly severe form of CMT disease is triggered by the P182L mutation in the highly conserved IxI/V motif of the disordered C-terminal region, which interacts weakly with the structured core domain of HSP27. Here, we observed that the P182L mutation disrupts the chaperone activity and significantly increases the size of HSP27 oligomers formed in vivo, including in motor neurons differentiated from CMT patient-derived stem cells. Using NMR spectroscopy, we determined that the P182L mutation decreases the affinity of the HSP27 IxI/V motif for its own core domain, leaving this binding site more accessible for other IxI/V-containing proteins. We identified multiple IxI/V-bearing proteins that bind with higher affinity to the P182L variant due to the increased availability of the IxI/V-binding site. Our results provide a mechanistic basis for the impact of the P182L mutation on HSP27 and suggest that the IxI/V motif plays an important, regulatory role in modulating protein-protein interactions.

Identifying molecular features that are associated with biological function of intrinsically disordered protein regions.

In previous work, we showed that intrinsically disordered regions (IDRs) of proteins contain sequence-distributed molecular features that are conserved over evolution, despite little sequence similarity that can be detected in alignments (Zarin et al., 2019). Here, we aim to use these molecular features to predict specific biological functions for individual IDRs and identify the molecular features within them that are associated with these functions. We find that the predictable functions are diverse. Examining the associated molecular features, we note some that are consistent with previous reports and identify others that were previously unknown. We experimentally confirm that elevated isoelectric point and hydrophobicity, features that are positively associated with mitochondrial localization, are necessary for mitochondrial targeting function. Remarkably, increasing isoelectric point in a synthetic IDR restores weak mitochondrial targeting. We believe feature analysis represents a new systematic approach to understand how biological functions of IDRs are specified by their protein sequences.

Phase Separation as a Missing Mechanism for Interpretation of Disease Mutations.

It is unclear how disease mutations impact intrinsically disordered protein regions (IDRs), which lack a stable folded structure. These mutations, while prevalent in disease, are frequently neglected or annotated as variants of unknown significance. Biomolecular phase separation, a physical process often mediated by IDRs, has increasingly appreciated roles in cellular organization and regulation. We find that autism spectrum disorder (ASD)- and cancer-associated proteins are enriched for predicted phase separation propensities, suggesting that IDR mutations disrupt phase separation in key cellular processes. More generally, we hypothesize that combinations of small-effect IDR mutations perturb phase separation, potentially contributing to "missing heritability" in complex disease susceptibility.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans.

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5' splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

Whence Blobs? Phylogenetics of functional protein condensates.

What do we know about the molecular evolution of functional protein condensation? The capacity of proteins to form biomolecular condensates (compact, protein-rich states, not bound by membranes, but still separated from the rest of the contents of the cell) appears in many cases to be bestowed by weak, transient interactions within one or between proteins. Natural selection is expected to remove or fix amino acid changes, insertions or deletions that preserve and change this condensation capacity when doing so is beneficial to the cell. A few recent studies have begun to explore this frontier of phylogenetics at the intersection of biophysics and cell biology.

Automated assignment of methyl NMR spectra from large proteins.

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy.

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (Rh) of biomolecules in solution. For unfolded proteins, the Rh provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the Rh of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by 13C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that Rh values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the Rh of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Automatic structure-based NMR methyl resonance assignment in large proteins.

Isotopically labeled methyl groups provide NMR probes in large, otherwise deuterated proteins. However, the resonance assignment constitutes a bottleneck for broader applicability of methyl-based NMR. Here, we present the automated MethylFLYA method for the assignment of methyl groups that is based on methyl-methyl nuclear Overhauser effect spectroscopy (NOESY) peak lists. MethylFLYA is applied to five proteins (28-358 kDa) comprising a total of 708 isotope-labeled methyl groups, of which 612 contribute NOESY cross peaks. MethylFLYA confidently assigns 488 methyl groups, i.e. 80% of those with NOESY data. Of these, 459 agree with the reference, 6 were different, and 23 were without reference assignment. MethylFLYA assigns significantly more methyl groups than alternative algorithms, has an average error rate of 1%, modest runtimes of 0.4-1.2 h, and can handle arbitrary isotope labeling patterns and data from other types of NMR spectra.

Local unfolding of the HSP27 monomer regulates chaperone activity.

The small heat-shock protein HSP27 is a redox-sensitive molecular chaperone that is expressed throughout the human body. Here, we describe redox-induced changes to the structure, dynamics, and function of HSP27 and its conserved α-crystallin domain (ACD). While HSP27 assembles into oligomers, we show that the monomers formed upon reduction are highly active chaperones in vitro, but are susceptible to self-aggregation. By using relaxation dispersion and high-pressure nuclear magnetic resonance (NMR) spectroscopy, we observe that the pair of ß-strands that mediate dimerisation partially unfold in the monomer. We note that numerous HSP27 mutations associated with inherited neuropathies cluster to this dynamic region. High levels of sequence conservation in ACDs from mammalian sHSPs suggest that the exposed, disordered interface present in free monomers or oligomeric subunits may be a general, functional feature of sHSPs.

Entropy and Information within Intrinsically Disordered Protein Regions.

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Automatic Assignment of Methyl-NMR Spectra of Supramolecular Machines Using Graph Theory.

Methyl groups are powerful probes for the analysis of structure, dynamics and function of supramolecular assemblies, using both solution- and solid-state NMR. Widespread application of the methodology has been limited due to the challenges associated with assigning spectral resonances to specific locations within a biomolecule. Here, we present Methyl Assignment by Graph Matching (MAGMA), for the automatic assignment of methyl resonances. A graph matching protocol examines all possibilities for each resonance in order to determine an exact assignment that includes a complete description of any ambiguity. MAGMA gives 100% accuracy in confident assignments when tested against both synthetic data, and 9 cross-validated examples using both solution- and solid-state NMR data. We show that this remarkable accuracy enables a user to distinguish between alternative protein structures. In a drug discovery application on HSP90, we show the method can rapidly and efficiently distinguish between possible ligand binding modes. By providing an exact and robust solution to methyl resonance assignment, MAGMA can facilitate significantly accelerated studies of supramolecular machines using methyl-based NMR spectroscopy.