Sequence-Based Prediction of Fuzzy Protein Interactions

Predicting Conformational Ensembles of Intrinsically Disordered Proteins: From Molecular Dynamics to Machine Learning

Intrinsically disordered proteins and regions (IDP/IDRs) are ubiquitous across all domains of life. Characterized by a lack of a stable tertiary structure, IDP/IDRs populate a diverse set of transiently formed structural states that can promiscuously adapt upon binding with specific interaction partners and/or certain alterations in environmental conditions. This malleability is foundational for their role as tunable interaction hubs in core cellular processes such as signaling, transcription, and translation. Tracing the conformational ensemble of an IDP/IDR and its perturbation in response to regulatory cues is thus paramount for illuminating its function. However, the conformational heterogeneity of IDP/IDRs poses several challenges. Here, we review experimental and computational methods devised to disentangle the conformational landscape of IDP/IDRs, highlighting recent computational advances that permit proteome-wide scans of IDP/IDRs conformations. We briefly evaluate selected computational methods using the disordered N-terminal of the human copper transporter 1 as a test case and outline further challenges in IDP/IDRs ensemble prediction.

Context-dependent, fuzzy protein interactions: Towards sequence-based insights

Predicting protein interactions in the cellular environment still remains a challenge in the AlphaFold era. Protein interactions, similarly to their structures, sample a continuum from ordered to disordered states, with specific partners in many bound configurations. A multiplicity of binding modes (MBM) enables transition between these states under different cellular conditions. This review focuses on how the cellular environment affects protein interactions, highlighting the molecular mechanisms, biophysical origin, and sequence-based principles of context-dependent, fuzzy interactions. It summarises experimental and computational approaches to address the challenge of interaction heterogeneity and its contribution to a wide range of biological functions. These insights will help in understanding complex cellular processes, involving conversions between protein assembly states, such as from liquid-like droplet state to the amyloid state.

Intrinsic disorder in the human vitreous proteome

The vitreous is a vital component of the eye, occupying a substantial portion of its volume and maintaining its structure. This study delves into the presence and significance of intrinsically disordered proteins (IDPs) within the vitreous, utilizing a dataset of 1240 vitreous proteins previously discovered in the vitreous proteome by Murthy et al.in five healthy subjects. The results indicate that 26.9 % of vitreous proteins are highly disordered, 68.8 % possess moderate disorder, and only 4.3 % are highly ordered. A complex interaction network among these proteins suggests their biological importance, and approximately 25 % may undergo liquid-liquid phase separation (LLPS). These findings offer new perspectives on the vitreous' molecular composition and behavior, potentially impacting our understanding of eye-related diseases, physiological changes such as vitreous syneresis. Further research is needed to translate these insights into clinical applications, although the intrinsic protein disorder and its association with LLPS appears to play a role in vitreous proteome function.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Deep mutational scanning of the RNase III-like domain in Trypanosoma brucei RNA editing protein KREPB4

Kinetoplastid pathogens including Trypanosoma brucei, T. cruzi, and Leishmania species, are early diverged, eukaryotic, unicellular parasites. Functional understanding of many proteins from these pathogens has been hampered by limited sequence homology to proteins from other model organisms. Here we describe the development of a high-throughput deep mutational scanning approach in T. brucei that facilitates rapid and unbiased assessment of the impacts of many possible amino acid substitutions within a protein on cell fitness, as measured by relative cell growth. The approach leverages several molecular technologies: cells with conditional expression of a wild-type gene of interest and constitutive expression of a library of mutant variants, degron-controlled stabilization of I-SceI meganuclease to mediate highly efficient transfection of a mutant allele library, and a high-throughput sequencing readout for cell growth upon conditional knockdown of wild-type gene expression and exclusive expression of mutant variants. Using this method, we queried the effects of amino acid substitutions in the apparently non-catalytic RNase III-like domain of KREPB4 (B4), which is an essential component of the RNA Editing Catalytic Complexes (RECCs) that carry out mitochondrial RNA editing in T. brucei. We measured the impacts of thousands of B4 variants on bloodstream form cell growth and validated the most deleterious variants containing single amino acid substitutions. Crucially, there was no correlation between phenotypes and amino acid conservation, demonstrating the greater power of this method over traditional sequence homology searching to identify functional residues. The bloodstream form cell growth phenotypes were combined with structural modeling, RECC protein proximity data, and analysis of selected substitutions in procyclic form T. brucei. These analyses revealed that the B4 RNaseIII-like domain is essential for maintenance of RECC integrity and RECC protein abundances and is also involved in changes in RECCs that occur between bloodstream and procyclic form life cycle stages.

Protein structure-function continuum model: Emerging nexuses between specificity, evolution, and structure

The rationale for replacing the old binary of structure-function with the trinity of structure, disorder, and function has gained considerable ground in recent years. A continuum model based on the expanded form of the existing paradigm can now subsume importance of both conformational flexibility and intrinsic disorder in protein function. The disorder is actually critical for understanding the protein-protein interactions in many regulatory processes, formation of membrane-less organelles, and our revised notions of specificity as amply illustrated by moonlighting proteins. While its importance in formation of amyloids and function of prions is often discussed, the roles of intrinsic disorder in infectious diseases and protein function under extreme conditions are also becoming clear. This review is an attempt to discuss how our current understanding of protein function, specificity, and evolution fit better with the continuum model. This integration of structure and disorder under a single model may bring greater clarity in our continuing quest for understanding proteins and molecular mechanisms of their functionality.

Searching for Structure: Characterizing the Protein Conformational Landscape with Clustering-Based Algorithms

The identification and characterization of the main conformations from a protein population are a challenging and inherently high-dimensional problem. Here, we evaluate the performance of the Secondary sTructural Ensembles with machine LeArning (StELa) double-clustering method, which clusters protein structures based on the relationship between the φ and ψ dihedral angles in a protein backbone and the secondary structure of the protein, thus focusing on the local properties of protein structures. The classification of states as vectors composed of the clusters' indices arising naturally from the Ramachandran plot is followed by the hierarchical clustering of the vectors to allow for the identification of the main features of the corresponding free energy landscape (FEL). We compare the performance of StELa with the established root-mean-squared-deviation (RMSD)-based clustering algorithm, which focuses on global properties of protein structures and with Combinatorial Averaged Transient Structure (CATS), the combinatorial averaged transient structure clustering method based on distributions of the φ and ψ dihedral angle coordinates. Using ensembles of conformations from molecular dynamics simulations of intrinsically disordered proteins (IDPs) of various lengths (tau protein fragments) or short fragments from a globular protein, we show that StELa is the clustering method that identifies many of the minima and relevant energy states around the minima from the corresponding FELs. In contrast, the RMSD-based algorithm yields a large number of clusters that usually cover most of the FEL, thus being unable to distinguish between states, while CATS does not sample well the FELs for long IDPs and fragments from globular proteins.

Inherently disordered regions of axonemal dynein assembly factors

The dynein-driven beating of cilia is required to move individual cells and to generate fluid flow across surfaces and within cavities. These motor enzymes are highly complex and can contain upwards of 20 different protein components with a total mass approaching 2 MDa. The dynein heavy chains are enormous proteins consisting of ~4500 residues and ribosomes take approximately 15 min to synthesize one. Studies in a broad array of organisms ranging from the green alga Chlamydomonas to humans has identified 19 cytosolic factors (DNAAFs) that are needed to specifically build axonemal dyneins; defects in many of these proteins lead to primary ciliary dyskinesia in mammals which can result in infertility, severe bronchial problems, and situs inversus. How all these factors cooperate in a spatially and temporally regulated manner to promote dynein assembly in cytoplasm remains very uncertain. These DNAAFs contain a variety of well-folded domains many of which provide protein interaction surfaces. However, many also exhibit large regions that are predicted to be inherently disordered. Here I discuss the nature of these unstructured segments, their predicted propensity for driving protein phase separation, and their potential for adopting more defined conformations during the dynein assembly process.

Protein nanocondensates: the next frontier

Over the past decade, myriads of studies have highlighted the central role of protein condensation in subcellular compartmentalization and spatiotemporal organization of biological processes. Conceptually, protein condensation stands at the highest level in protein structure hierarchy, accounting for the assembly of bodies ranging from thousands to billions of molecules and for densities ranging from dense liquids to solid materials. In size, protein condensates range from nanocondensates of hundreds of nanometers (mesoscopic clusters) to phase-separated micron-sized condensates. In this review, we focus on protein nanocondensation, a process that can occur in subsaturated solutions and can nucleate dense liquid phases, crystals, amorphous aggregates, and fibers. We discuss the nanocondensation of proteins in the light of general physical principles and examine the biophysical properties of several outstanding examples of nanocondensation. We conclude that protein nanocondensation cannot be fully explained by the conceptual framework of micron-scale biomolecular condensation. The evolution of nanocondensates through changes in density and order is currently under intense investigation, and this should lead to the development of a general theoretical framework, capable of encompassing the full range of sizes and densities found in protein condensates.

Towards sequence-based principles for protein phase separation predictions

The phenomenon of protein phase separation, which underlies the formation of biomolecular condensates, has been associated with numerous cellular functions. Recent studies indicate that the amino acid sequences of most proteins may harbour not only the code for folding into the native state but also for condensing into the liquid-like droplet state and the solid-like amyloid state. Here we review the current understanding of the principles for sequence-based methods for predicting the propensity of proteins for phase separation. A guiding concept is that entropic contributions are generally more important to stabilise the droplet state than they are for the native and amyloid states. Although estimating these entropic contributions has proven difficult, we describe some progress that has been recently made in this direction. To conclude, we discuss the challenges ahead to extend sequence-based prediction methods of protein phase separation to include quantitative in vivo characterisations of this process.

A sequence-based foldability score combined with AlphaFold2 predictions to disentangle the protein order/disorder continuum

Order and disorder govern protein functions, but there is a great diversity in disorder, from regions that are-and stay-fully disordered to conditional order. This diversity is still difficult to decipher even though it is encoded in the amino acid sequences. Here, we developed an analytic Python package, named pyHCA, to estimate the foldability of a protein segment from the only information of its amino acid sequence and based on a measure of its density in regular secondary structures associated with hydrophobic clusters, as defined by the hydrophobic cluster analysis (HCA) approach. The tool was designed by optimizing the separation between foldable segments from databases of disorder (DisProt) and order (SCOPe [soluble domains] and OPM [transmembrane domains]). It allows to specify the ratio between order, embodied by regular secondary structures (either participating in the hydrophobic core of well-folded 3D structures or conditionally formed in intrinsically disordered regions) and disorder. We illustrated the relevance of pyHCA with several examples and applied it to the sequences of the proteomes of 21 species ranging from prokaryotes and archaea to unicellular and multicellular eukaryotes, for which structure models are provided in the AlphaFold protein structure database. Cases of low-confidence scores related to disorder were distinguished from those of sequences that we identified as foldable but are still excluded from accurate modeling by AlphaFold2 due to a lack of sequence homologs or to compositional biases. Overall, our approach is complementary to AlphaFold2, providing guides to map structural innovations through evolutionary processes, at proteome and gene scales.

AlphaFold2 and its applications in the fields of biology and medicine

AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.

Alternatively spliced exon regulates context-dependent MEF2D higher-order assembly during myogenesis

During muscle cell differentiation, the alternatively spliced, acidic β-domain potentiates transcription of Myocyte-specific Enhancer Factor 2 (Mef2D). Sequence analysis by the FuzDrop method indicates that the β-domain can serve as an interaction element for Mef2D higher-order assembly. In accord, we observed Mef2D mobile nuclear condensates in C2C12 cells, similar to those formed through liquid-liquid phase separation. In addition, we found Mef2D solid-like aggregates in the cytosol, the presence of which correlated with higher transcriptional activity. In parallel, we observed a progress in the early phase of myotube development, and higher MyoD and desmin expression. In accord with our predictions, the formation of aggregates was promoted by rigid β-domain variants, as well as by a disordered β-domain variant, capable of switching between liquid-like and solid-like higher-order states. Along these lines, NMR and molecular dynamics simulations corroborated that the β-domain can sample both ordered and disordered interactions leading to compact and extended conformations. These results suggest that β-domain fine-tunes Mef2D higher-order assembly to the cellular context, which provides a platform for myogenic regulatory factors and the transcriptional apparatus during the developmental process.

Computational Prediction of Protein Intrinsically Disordered Region Related Interactions and Functions

Intrinsically Disordered Proteins (IDPs) and Regions (IDRs) exist widely. Although without well-defined structures, they participate in many important biological processes. In addition, they are also widely related to human diseases and have become potential targets in drug discovery. However, there is a big gap between the experimental annotations related to IDPs/IDRs and their actual number. In recent decades, the computational methods related to IDPs/IDRs have been developed vigorously, including predicting IDPs/IDRs, the binding modes of IDPs/IDRs, the binding sites of IDPs/IDRs, and the molecular functions of IDPs/IDRs according to different tasks. In view of the correlation between these predictors, we have reviewed these prediction methods uniformly for the first time, summarized their computational methods and predictive performance, and discussed some problems and perspectives.

Protein interactions: anything new?

How do proteins interact in the cellular environment? Which interactions stabilize liquid-liquid phase separated condensates? Are the concepts, which have been developed for specific protein complexes also applicable to higher-order assemblies? Recent discoveries prompt for a universal framework for protein interactions, which can be applied across the scales of protein communities. Here, we discuss how our views on protein interactions have evolved from rigid structures to conformational ensembles of proteins and discuss the open problems, in particular related to biomolecular condensates. Protein interactions have evolved to follow changes in the cellular environment, which manifests in multiple modes of interactions between the same partners. Such cellular context-dependence requires multiplicity of binding modes (MBM) by sampling multiple minima of the interaction energy landscape. We demonstrate that the energy landscape framework of protein folding can be applied to explain this phenomenon, opening a perspective toward a physics-based, universal model for cellular protein behaviors.

Mapping the deformability of natural and designed cellulosomes in solution

Background

Natural cellulosome multi-enzyme complexes, their components, and engineered ‘designer cellulosomes’ (DCs) promise an efficient means of breaking down cellulosic substrates into valuable biofuel products. Their broad uptake in biotechnology relies on boosting proximity-based synergy among the resident enzymes, but the modular architecture challenges structure determination and rational design.

Results

We used small angle X-ray scattering combined with molecular modeling to study the solution structure of cellulosomal components. These include three dockerin-bearing cellulases with distinct substrate specificities, original scaffoldins from the human gut bacterium Ruminococcus champanellensis (ScaA, ScaH and ScaK) and a trivalent cohesin-bearing designer scaffoldin (Scaf20L), followed by cellulosomal complexes comprising these components, and the nonavalent fully loaded Clostridium thermocellum CipA in complex with Cel8A from the same bacterium. The size analysis of R_g and D_max values deduced from the scattering curves and corresponding molecular models highlight their variable aspects, depending on composition, size and spatial organization of the objects in solution.

Conclusions

Our data quantifies variability of form and compactness of cellulosomal components in solution and confirms that this native plasticity may well be related to speciation with respect to the substrate that is targeted. By showing that scaffoldins or components display enhanced compactness compared to the free objects, we provide new routes to rationally enhance their stability and performance in their environment of action.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02165-3.

NMR insights into dynamic, multivalent interactions of intrinsically disordered regions: from discrete complexes to condensates

The spatial and temporal organization of interactions between proteins underlie the regulation of most cellular processes. The requirement for such interactions to be specific predisposes a view that protein-protein interactions are relatively static and are formed through the stable complementarity of the interacting partners. A growing body of reports indicate, however, that many interactions lead to fuzzy complexes with an ensemble of conformations in dynamic exchange accounting for the observed binding. Here, we discuss how NMR has facilitated the characterization of these discrete, dynamic complexes and how such characterization has aided the understanding of dynamic, condensed phases of phase-separating proteins with exchanging multivalent interactions.

The MDMX Acidic Domain Uses Allovalency to Bind Both p53 and MDMX

Autoinhibition of p53 binding to MDMX requires two short-linear motifs (SLiMs) containing adjacent tryptophan (WW) and tryptophan-phenylalanine (WF) residues. NMR spectroscopy was used to show the WW and WF motifs directly compete for the p53 binding site on MDMX and circular dichroism spectroscopy was used to show the WW motif becomes helical when it is bound to the p53 binding domain (p53BD) of MDMX. Binding studies using isothermal titration calorimetry showed the WW motif is a stronger inhibitor of p53 binding than the WF motif when they are both tethered to p53BD by the natural disordered linker. We also investigated how the WW and WF motifs interact with the DNA binding domain (DBD) of p53. Both motifs bind independently to similar sites on DBD that overlap the DNA binding site. Taken together our work defines a model for complex formation between MDMX and p53 where a pair of disordered SLiMs bind overlapping sites on both proteins.

MobiDB: 10 years of intrinsically disordered proteins

The MobiDB database (URL: https://mobidb.org/) is a knowledge base of intrinsically disordered proteins. MobiDB aggregates disorder annotations derived from the literature and from experimental evidence along with predictions for all known protein sequences. MobiDB generates new knowledge and captures the functional significance of disordered regions by processing and combining complementary sources of information. Since its first release 10 years ago, the MobiDB database has evolved in order to improve the quality and coverage of protein disorder annotations and its accessibility. MobiDB has now reached its maturity in terms of data standardization and visualization. Here, we present a new release which focuses on the optimization of user experience and database content. The major advances compared to the previous version are the integration of AlphaFoldDB predictions and the re-implementation of the homology transfer pipeline, which expands manually curated annotations by two orders of magnitude. Finally, the entry page has been restyled in order to provide an overview of the available annotations along with two separate views that highlight structural disorder evidence and functions associated with different binding modes.

New insights into disordered proteins and regions according to the FOD-M model

A collection of intrinsically disordered proteins (IDPs) having regions with the status of intrinsically disordered (IDR) according to the Disprot database was analyzed from the point of view of the structure of hydrophobic core in the structural unit (chain / domain). The analysis includes all the Homo Sapiens as well as Mus Musculus proteins present in the DisProt database for which the structure is available. In the analysis, the fuzzy oil drop modified model (FOD-M) was used, taking into account the external force field, modified by the presence of other factors apart from polar water, influencing protein structuring. The paper presents an alternative to secondary-structure-based classification of intrinsically disordered regions (IDR). The basis of our classification is the ordering of hydrophobic core as calculated by the FOD-M model resulting in FOD-ordered or FOD-unordered IDRs.

Nanomedicines Meet Disordered Proteins: A Shift from Traditional Materials and Concepts to Innovative Polymers

Water-soluble nanomedicines have been widely studied for the targeted delivery of drugs for a very long time. As a notable example, biomaterials based on N-(2-hydroxypropyl) methacrylamide (HPMA) copolymers have been under investigation for nearly half a century. In particular, anticancer drug carriers have been developed under the assumption that the leading mechanism with a therapeutic impact on solid tumors is the enhanced permeability and retention (EPR) effect, which dates back more than three decades. Nevertheless, these (and other) materials and concepts have encountered several barriers in their successful translation into clinical practice, and future nanomedicines need improvements in both passive and active targeting to their site of action. Notions borrowed from recent studies on intrinsically disordered proteins (IDPs) seem promising for enhancing the self-assembly, stimuli-responsiveness, and recognition properties of protein/peptide-based copolymers. Accordingly, IDP-based nanomedicines are ready to give new impetus to more traditional research in this field.

Protein conformation and biomolecular condensates

Protein conformation and cell compartmentalization are fundamental concepts and subjects of vast scientific endeavors. In the last two decades, we have witnessed exciting advances that unveiled the conjunction of these concepts. An avalanche of studies highlighted the central role of biomolecular condensates in membraneless subcellular compartmentalization that permits the spatiotemporal organization and regulation of myriads of simultaneous biochemical reactions and macromolecular interactions. These studies have also shown that biomolecular condensation, driven by multivalent intermolecular interactions, is mediated by order-disorder transitions of protein conformation and by protein domain architecture. Conceptually, protein condensation is a distinct level in protein conformational landscape in which collective folding of large collections of molecules takes place. Biomolecular condensates arise by the physical process of phase separation and comprise a variety of bodies ranging from membraneless organelles to liquid condensates to solid-like conglomerates, spanning lengths from mesoscopic clusters (nanometers) to micrometer-sized objects. In this review, we summarize and discuss recent work on the assembly, composition, conformation, material properties, thermodynamics, regulation, and functions of these bodies. We also review the conceptual framework for future studies on the conformational dynamics of condensed proteins in the regulation of cellular processes.

Challenges in describing the conformation and dynamics of proteins with ambiguous behavior

Traditionally, our understanding of how proteins operate and how evolution shapes them is based on two main data sources: the overall protein fold and the protein amino acid sequence. However, a significant part of the proteome shows highly dynamic and/or structurally ambiguous behavior, which cannot be correctly represented by the traditional fixed set of static coordinates. Representing such protein behaviors remains challenging and necessarily involves a complex interpretation of conformational states, including probabilistic descriptions. Relating protein dynamics and multiple conformations to their function as well as their physiological context (e.g., post-translational modifications and subcellular localization), therefore, remains elusive for much of the proteome, with studies to investigate the effect of protein dynamics relying heavily on computational models. We here investigate the possibility of delineating three classes of protein conformational behavior: order, disorder, and ambiguity. These definitions are explored based on three different datasets, using interpretable machine learning from a set of features, from AlphaFold2 to sequence-based predictions, to understand the overlap and differences between these datasets. This forms the basis for a discussion on the current limitations in describing the behavior of dynamic and ambiguous proteins.

Are casein micelles extracellular condensates formed by liquid-liquid phase separation?

Casein micelles are extracellular polydisperse assemblies of unstructured casein proteins. Caseins are the major component of milk. Within casein micelles, casein molecules are stabilised by binding to calcium phosphate nanoclusters and, by acting as molecular chaperones, through multivalent interactions. In light of such interactions, we discuss whether casein micelles can be considered as extracellular condensates formed by liquid-liquid phase separation. We analyse the sequence, structure and interactions of caseins in comparison to proteins forming intracellular condensates. Furthermore, we review the similarities between caseins and small heat-shock proteins whose chaperone activity is linked to phase separation of proteins. By bringing these observations together, we describe a regulatory mechanism for protein condensates, as exemplified by casein micelles.

FuzDrop on AlphaFold: visualizing the sequence-dependent propensity of liquid-liquid phase separation and aggregation of proteins

Many proteins perform their functions within membraneless organelles, where they form a liquid-like condensed state, also known as droplet state. The FuzDrop method predicts the probability of spontaneous liquid-liquid phase separation of proteins and provides a sequence-based score to identify the regions that promote this process. Furthermore, the FuzDrop method estimates the propensity of conversion of proteins to the amyloid state, and identifies aggregation hot-spots, which can drive the irreversible maturation of the liquid-like droplet state. These predictions can also identify mutations that can induce formation of amyloid aggregates, including those implicated in human diseases. To facilitate the interpretation of the predictions, the droplet-promoting and aggregation-promoting regions can be visualized on protein structures generated by AlphaFold. The FuzDrop server (https://fuzdrop.bio.unipd.it) thus offers insights into the complex behavior of proteins in their condensed states and facilitates the understanding of the functional relationships of proteins.

Fuzzy RNA recognition by the Trypanosoma brucei editosome

The assembly of high molecular mass ribonucleoprotein complexes typically relies on the binary interaction of defined RNA sequences or precisely folded RNA motifs with dedicated RNA-binding domains on the protein side. Here we describe a new molecular recognition principle of RNA molecules by a high molecular mass protein complex. By chemically probing the solvent accessibility of mitochondrial pre-mRNAs when bound to the Trypanosoma brucei editosome, we identified multiple similar but non-identical RNA motifs as editosome contact sites. However, by treating the different motifs as mathematical graph objects we demonstrate that they fit a consensus 2D-graph consisting of 4 vertices (V) and 3 edges (E) with a Laplacian eigenvalue of 0.5477 (λ₂). We establish that synthetic 4V(3E)-RNAs are sufficient to compete for the editosomal pre-mRNA binding site and that they inhibit RNA editing in vitro. Furthermore, we demonstrate that only two topological indices are necessary to predict the binding of any RNA motif to the editosome with a high level of confidence. Our analysis corroborates that the editosome has adapted to the structural multiplicity of the mitochondrial mRNA folding space by recognizing a fuzzy continuum of RNA folds that fit a consensus graph descriptor.

SRC homology 3 domains: multifaceted binding modules

The assembly of complexes following the detection of extracellular signals is often controlled by signaling proteins comprising multiple peptide binding modules. The SRC homology (SH)3 family represents the archetypical modular protein interaction module, with ~300 annotated SH3 domains in humans that regulate an impressive array of signaling processes. We review recent findings regarding the allosteric contributions of SH3 domains host protein context, their phosphoregulation, and their roles in phase separation that challenge the simple model in which SH3s are considered to be portable domains binding to specific proline-rich peptide motifs.

AlphaFold2: A Role for Disordered Protein/Region Prediction?

The development of AlphaFold2 marked a paradigm-shift in the structural biology community. Herein, we assess the ability of AlphaFold2 to predict disordered regions against traditional sequence-based disorder predictors. We find that AlphaFold2 performs well at discriminating disordered regions, but also note that the disorder predictor one constructs from an AlphaFold2 structure determines accuracy. In particular, a naïve, but non-trivial assumption that residues assigned to helices, strands, and H-bond stabilized turns are likely ordered and all other residues are disordered results in a dramatic overestimation in disorder; conversely, the predicted local distance difference test (pLDDT) provides an excellent measure of residue-wise disorder. Furthermore, by employing molecular dynamics (MD) simulations, we note an interesting relationship between the pLDDT and secondary structure, that may explain our observations and suggests a broader application of the pLDDT for characterizing the local dynamics of intrinsically disordered proteins and regions (IDPs/IDRs).

Targeting key proteins involved in transcriptional regulation for cancer therapy: Current strategies and future prospective

The key proteins involved in transcriptional regulation play convergent roles in cellular homeostasis, and their dysfunction mediates aberrant gene expressions that underline the hallmarks of tumorigenesis. As tumor progression is dependent on such abnormal regulation of transcription, it is important to discover novel chemical entities as antitumor drugs that target key tumor-associated proteins involved in transcriptional regulation. Despite most key proteins (especially transcription factors) involved in transcriptional regulation are historically recognized as undruggable targets, multiple targeting approaches at diverse levels of transcriptional regulation, such as epigenetic intervention, inhibition of DNA-binding of transcriptional factors, and inhibition of the protein-protein interactions (PPIs), have been established in preclinically or clinically studies. In addition, several new approaches have recently been described, such as targeting proteasomal degradation and eliciting synthetic lethality. This review will emphasize on accentuating these developing therapeutic approaches and provide a thorough conspectus of the drug development to target key proteins involved in transcriptional regulation and their impact on future oncotherapy.

Intrinsically disordered signaling proteins: Essential hub players in the control of stress responses in Saccharomyces cerevisiae

Cells have developed diverse mechanisms to monitor changes in their surroundings. This allows them to establish effective responses to cope with adverse environments. Some of these mechanisms have been well characterized in the budding yeast Saccharomyces cerevisiae, an excellent experimental model to explore and elucidate some of the strategies selected in eukaryotic organisms to adjust their growth and development in stressful conditions. The relevance of structural disorder in proteins and the impact on their functions has been uncovered for proteins participating in different processes. This is the case of some transcription factors (TFs) and other signaling hub proteins, where intrinsically disordered regions (IDRs) play a critical role in their function. In this work, we present a comprehensive bioinformatic analysis to evaluate the significance of structural disorder in those TFs (170) recognized in S. cerevisiae. Our findings show that 85.2% of these TFs contain at least one IDR, whereas ~30% exhibit a higher disorder level and thus were considered as intrinsically disordered proteins (IDPs). We also found that TFs contain a higher number of IDRs compared to the rest of the yeast proteins, and that intrinsically disordered TFs (IDTFs) have a higher number of protein-protein interactions than those with low structural disorder. The analysis of different stress response pathways showed a high content of structural disorder not only in TFs but also in other signaling proteins. The propensity of yeast proteome to undergo a liquid-liquid phase separation (LLPS) was also analyzed, showing that a significant proportion of IDTFs may undergo this phenomenon. Our analysis is a starting point for future research on the importance of structural disorder in yeast stress responses.

Topological frustration leading to backtracking in a coupled folding-binding process

Intrinsically disordered proteins (IDPs) are abundant in all species. Their discovery challenges the traditional "sequence-structure-function" paradigm of protein science because IDPs play important roles in various biological processes without preformed folded structures. Bioinformatic analysis reveals that the intrinsically conformational disorder of IDPs as well as their conformational transition upon binding to their targets is encoded by their amino acid sequences. The rRNase domain of colicin E3 and the immunity protein Im3 are a pair of proteins involved in bacterial survival. While the N-terminal segment and the central segment of E3 make comparable intermolecular contacts with Im3 in the bound state, binding of E3 with Im3 is dominantly triggered by the central segment of E3. In this work, to further investigate the binding mechanism of disordered E3 with Im3, we performed systematic free energy and transition path analysis through coarse-grained molecular dynamics simulations. We observed backtracking of the N-terminal segment of E3 in the binding process, whose occurrence depends on salt concentration. Conformational analysis revealed that initial binding of the N-terminal segment of E3 to Im3 usually leads to misorientation of a central hairpin of E3 on Im3, which generates topological frustration and results in backtracking of the N-terminal segment. Our results not only provide deeper mechanistic insights into the coupled folding-binding process of the E3/Im3 complex, but also suggest that topological frustration could be present in the coupled folding-binding process of IDPs and play an important role in regulating the binding transition pathways.

Molecular Determinants of Selectivity in Disordered Complexes May Shed Light on Specificity in Protein Condensates

Biomolecular condensates challenge the classical concepts of molecular recognition. The variable composition and heterogeneous conformations of liquid-like protein droplets are bottlenecks for high-resolution structural studies. To obtain atomistic insights into the organization of these assemblies, here we have characterized the conformational ensembles of specific disordered complexes, including those of droplet-driving proteins. First, we found that these specific complexes exhibit a high degree of conformational heterogeneity. Second, we found that residues forming contacts at the interface also sample many conformations. Third, we found that different patterns of contacting residues form the specific interface. In addition, we observed a wide range of sequence motifs mediating disordered interactions, including charged, hydrophobic and polar contacts. These results demonstrate that selective recognition can be realized by variable patterns of weakly defined interaction motifs in many different binding configurations. We propose that these principles also play roles in determining the selectivity of biomolecular condensates.

Characterization of the High-Affinity Fuzzy Complex between the Disordered Domain of the E7 Oncoprotein from High-Risk HPV and the TAZ2 Domain of CBP

The intrinsically disordered N-terminal region of the E7 protein from high-risk human papillomavirus (HPV) strains is responsible for oncogenic transformation of host cells through its interaction with a number of cellular factors, including the TAZ2 domain of the transcriptional coactivator CREB-binding protein. Using a variety of spectroscopic and biochemical tools, we find that despite its nanomolar affinity, the HPV16 E7 complex with TAZ2 is disordered and highly dynamic. The disordered domain of HPV16 E7 protein does not adopt a single conformation on the surface of TAZ2 but engages promiscuously with its target through multiple interactions involving two conserved motifs, termed CR1 and CR2, that occupy an extensive binding surface on TAZ2. The fuzzy nature of the complex is a reflection of the promiscuous binding repertoire of viral proteins, which must efficiently dysregulate host cell processes by binding to a variety of host factors in the cellular environment.

Unveiling induced folding of intrinsically disordered proteins - Protein engineering, frustration and emerging themes

Intrinsically disordered proteins (IDPs) can be generally described as a class of proteins that lack a well-defined ordered structure in isolation at physiological conditions. Upon binding to their physiological ligands, IDPs typically undergo a disorder-to-order transition, which may or may not lead to the complete folding of the IDP. In this short review, we focus on some of the key findings pertaining to the mechanisms of such induced folding. In particular, first we describe the general features of the reaction; then, we discuss some of the most remarkable findings obtained from applying protein engineering in synergy with kinetic studies to induced folding; and finally, we offer a critical view on some of the emerging themes when considering the structural heterogeneity of IDPs vis-à-vis to their inherent frustration.

FuzDB: a new phase in understanding fuzzy interactions

Fuzzy interactions are specific, variable contacts between proteins and other biomolecules (proteins, DNA, RNA, small molecules) formed in accord to the cellular context. Fuzzy interactions have recently been demonstrated to regulate biomolecular condensates generated by liquid-liquid phase separation. The FuzDB v4.0 database (https://fuzdb.org) assembles experimentally identified examples of fuzzy interactions, where disordered regions mediate functionally important, context-dependent contacts between the partners in stoichiometric and higher-order assemblies. The new version of FuzDB establishes cross-links with databases on structure (PDB, BMRB, PED), function (ELM, UniProt) and biomolecular condensates (PhaSepDB, PhaSePro, LLPSDB). FuzDB v4.0 is a source to decipher molecular basis of complex cellular interaction behaviors, including those in protein droplets.

Brain Immunoinformatics: A Symmetrical Link between Informatics, Wet Lab and the Clinic

Breakthrough advances in informatics over the last decade have thoroughly influenced the field of immunology. The intermingling of machine learning with wet lab applications and clinical results has hatched the newly defined immunoinformatics society. Immunoinformatics of the central neural system, referred to as neuroimmunoinformatics (NII), investigates symmetrical and asymmetrical interactions of the brain-immune interface. This interdisciplinary overview on NII is addressed to bioscientists and computer scientists. We delineate the dominating trajectories and field-shaping achievements and elaborate on future directions using bridging language and terminology. Computation, varying from linear modeling to complex deep learning approaches, fuels neuroimmunology through three core directions. Firstly, by providing big-data analysis software for high-throughput methods such as next-generation sequencing and genome-wide association studies. Secondly, by designing models for the prediction of protein morphology, functions, and symmetrical and asymmetrical protein–protein interactions. Finally, NII boosts the output of quantitative pathology by enabling the automatization of tedious processes such as cell counting, tracing, and arbor analysis. The new classification of microglia, the brain’s innate immune cells, was an NII achievement. Deep sequencing classifies microglia in “sensotypes” to accurately describe the versatility of immune responses to physiological and pathological challenges, as well as to experimental conditions such as xenografting and organoids. NII approaches complex tasks in the brain-immune interface, recognizes patterns and allows for hypothesis-free predictions with ultimate targeted individualized treatment strategies, and personalizes disease prognosis and treatment response.

An intrinsically disordered pathological prion variant Y145Stop converts into self-seeding amyloids via liquid-liquid phase separation

Biomolecular condensation via liquid-liquid phase separation of intrinsically disordered proteins/regions (IDPs/IDRs) along with other biomolecules is proposed to control critical cellular functions, whereas aberrant phase transitions are associated with a range of neurodegenerative diseases. Here, we show that a disease-associated stop codon mutation of the prion protein (PrP) at tyrosine 145 (Y145Stop), resulting in a truncated, highly disordered, N-terminal IDR, spontaneously phase-separates into dynamic liquid-like droplets. Phase separation of this highly positively charged N-terminal segment is promoted by the electrostatic screening and a multitude of weak, transient, multivalent, intermolecular interactions. Single-droplet Raman measurements, in conjunction with an array of bioinformatic, spectroscopic, microscopic, and mutagenesis studies, revealed a highly mobile internal organization within the liquid-like condensates. The phase behavior of Y145Stop is modulated by RNA. Lower RNA:protein ratios promote condensation at a low micromolar protein concentration under physiological conditions. At higher concentrations of RNA, phase separation is abolished. Upon aging, these highly dynamic liquid-like droplets gradually transform into ordered, β-rich, amyloid-like aggregates. These aggregates formed via phase transitions display an autocatalytic self-templating characteristic involving the recruitment and binding-induced conformational conversion of monomeric Y145Stop into amyloid fibrils. In contrast to this intrinsically disordered truncated variant, the wild-type full-length PrP exhibits a much lower propensity for both condensation and maturation into amyloids, hinting at a possible protective role of the C-terminal domain. Such an interplay of molecular factors in modulating the protein phase behavior might have much broader implications in cell physiology and disease.

Studying the Structures of Relaxed and Fuzzy Interactions: The Diverse World of S100 Complexes

S100 proteins are small, dimeric, Ca²⁺-binding proteins of considerable interest due to their associations with cancer and rheumatic and neurodegenerative diseases. They control the functions of numerous proteins by forming protein–protein complexes with them. Several of these complexes were found to display “fuzzy” properties. Examining these highly flexible interactions, however, is a difficult task, especially from a structural biology point of view. Here, we summarize the available in vitro techniques that can be deployed to obtain structural information about these dynamic complexes. We also review the current state of knowledge about the structures of S100 complexes, focusing on their often-asymmetric nature.

When Order Meets Disorder: Modeling and Function of the Protein Interface in Fuzzy Complexes

The degree of proteins structural organization ranges from highly structured, compact folding to intrinsic disorder, where each degree of self-organization corresponds to specific functions: well-organized structural motifs in enzymes offer a proper environment for precisely positioned functional groups to participate in catalytic reactions; at the other end of the self-organization spectrum, intrinsically disordered proteins act as binding hubs via the formation of multiple, transient and often non-specific interactions. This review focusses on cases where structurally organized proteins or domains associate with highly disordered protein chains, leading to the formation of interfaces with varying degrees of fuzziness. We present a review of the computational methods developed to provide us with information on such fuzzy interfaces, and how they integrate experimental information. The discussion focusses on two specific cases, microtubules and homologous recombination nucleoprotein filaments, where a network of intrinsically disordered tails exerts regulatory function in recruiting partner macromolecules, proteins or DNA and tuning the atomic level association. Notably, we show how computational approaches such as molecular dynamics simulations can bring new knowledge to help bridging the gap between experimental analysis, that mostly concerns ensemble properties, and the behavior of individual disordered protein chains that contribute to regulation functions.

AlphaFold and Implications for Intrinsically Disordered Proteins

Accurate predictions of the three-dimensional structures of proteins from their amino acid sequences have come of age. AlphaFold, a deep learning-based approach to protein structure prediction, shows remarkable success in independent assessments of prediction accuracy. A significant epoch in structural bioinformatics was the structural annotation of over 98% of protein sequences in the human proteome. Interestingly, many predictions feature regions of very low confidence, and these regions largely overlap with intrinsically disordered regions (IDRs). That over 30% of regions within the proteome are disordered is congruent with estimates that have been made over the past two decades, as intense efforts have been undertaken to generalize the structure-function paradigm to include the importance of conformational heterogeneity and dynamics. With structural annotations from AlphaFold in hand, there is the temptation to draw inferences regarding the "structures" of IDRs and their interactomes. Here, we offer a cautionary note regarding the misinterpretations that might ensue and highlight efforts that provide concrete understanding of sequence-ensemble-function relationships of IDRs. This perspective is intended to emphasize the importance of IDRs in sequence-function relationships (SERs) and to highlight how one might go about extracting quantitative SERs to make sense of how IDRs function.

Intrinsically disordered proteins: modes of binding with emphasis on disordered domains

Our notions of protein function have long been determined by the protein structure-function paradigm. However, the idea that protein function is dictated by a prerequisite complementarity of shapes at the binding interface is becoming increasingly challenged. Interactions involving intrinsically disordered proteins (IDPs) have indicated a significant degree of disorder present in the bound state, ranging from static disorder to complete disorder, termed 'random fuzziness'. This review assesses the anatomy of an IDP and relates how its intrinsic properties permit promiscuity and allow for the various modes of interaction. Furthermore, a mechanistic overview of the types of disordered domains is detailed, while also relating to a recent example and the kinetic and thermodynamic principles governing its formation.

The sequence-ensemble relationship in fuzzy protein complexes

Intrinsically disordered proteins (IDPs) interact with globular proteins through a variety of mechanisms, resulting in the structurally heterogeneous ensembles known as fuzzy complexes. While there exists a reasonable comprehension on how IDP sequence determines the unbound IDP ensemble, little is known about what shapes the structural characteristics of IDPs bound to their targets. Using a statistical thermodynamic model, we show that the target-bound ensembles are determined by a simple code that combines the IDP sequence and the distribution of IDP-target interaction hotspots. These two parameters define the conformational space of target-bound IDPs and rationalize the observed structural heterogeneity of fuzzy complexes. The presented model successfully reproduces the dynamical signatures of target-bound IDPs from the NMR relaxation experiments as well as the changes of interaction affinity and the IDP helicity induced by mutations. The model explains how the target-bound IDP ensemble adapts to mutations in order to achieve an optimal balance between conformational freedom and interaction energy. Taken together, the presented sequence-ensemble relationship of fuzzy complexes explains the different manifestations of IDP disorder in folding-upon-binding processes.

Sequence Determinants of the Aggregation of Proteins Within Condensates Generated by Liquid-liquid Phase Separation

The transition between the native and amyloid states of proteins can proceed via a deposition pathway via oligomeric intermediates or via a condensation pathway involving liquid droplet intermediates generated through liquid-liquid phase separation. While several computational methods are available to perform sequence-based predictions of the propensity of proteins to aggregate via the deposition pathway, much less is known about the physico-chemical principles that underlie aggregation within condensates. Here we investigate the sequence determinants of aggregation via the condensation pathway, and identify three relevant features: droplet-promoting propensity, aggregation-promoting propensity and multimodal interactions quantified by the binding mode entropy. By using this approach, we show that it is possible to predict aggregation-promoting mutations in droplet-forming proteins associated with amyotrophic lateral sclerosis (ALS). This analysis provides insights into the amino acid code for the conversion of proteins between liquid-like and solid-like condensates.

Fuzzy protein theory for disordered proteins

Why proteins are fuzzy? Constant adaptation to the cellular environment requires a wide range of changes in protein structure and interactions. Conformational ensembles of disordered proteins in particular exhibit large shifts to activate or inhibit alternative pathways. Fuzziness is critical for liquid–liquid phase separation and conversion of biomolecular condensates into fibrils. Interpretation of these phenomena presents a challenge for the classical structure-function paradigm. Here I discuss a multi-valued formalism, based on fuzzy logic, which can be applied to describe complex cellular behavior of proteins.