Intrinsically Disordered Proteins Research Articles

AlphaFold2 modeling and molecular dynamics simulations of an intrinsically disordered protein.

We use AlphaFold2 (AF2) to model the monomer and dimer structures of an intrinsically disordered protein (IDP), Nvjp-1, assisted by molecular dynamics (MD) simulations. We observe relatively rigid dimeric structures of Nvjp-1 when compared with the monomer structures. We suggest that protein conformations from multiple AF2 models and those from MD trajectories exhibit a coherent trend: the conformations of an IDP are deviated from each other and the conformations of a well-folded protein are consistent with each other. We use a residue-residue interaction network (RIN) derived from the contact map which show that the residue-residue interactions in Nvjp-1 are mainly transient; however, those in a well-folded protein are mainly persistent. Despite the variation in 3D shapes, we show that the AF2 models of both disordered and ordered proteins exhibit highly consistent profiles of the pLDDT (predicted local distance difference test) scores. These results indicate a potential protocol to justify the IDPs based on multiple AF2 models and MD simulations.

Identification of Essential Components of RNA Binding Domain of TLS/FUS

TLS/FUS is RNA-binding protein having multiple functions of regulations of genes, homeostasis, and cellular growth. Recent studies show that TLS is involved in phase separation and occasionally forms precipitation related to neurodegenerative diseases like amyotrophic lateral sclerosis (ALS). RNA has been reported to suppress phase separation, droplet formation, and concomitant precipitation of TLS, suggesting that RNA is a possible candidate for ALS drug discovery. Our experiments demonstrated that a long noncoding RNA, promoter-associated noncoding RNA (pncRNA-D), specifically binds TLS and represses its phase separation and precipitation. To obtain competent drug seeds, it is essential to reveal mechanism of action of lncRNAs with specificity to TLS and inhibitory activity on phase separation and related precipitation. For this purpose, several lncRNAs (lncRNAs 1 to 6) were selected upon assays with GST-TLS binding and inhibition on the precipitation. With criteria of binding specificity for TLS, lncRNA3 has been selected for further analysis for RNA-binding ability. Initially, RNA-binding region at TLS amino acid sequence was identified from four fragments of TLS. RNA binding assay with biotinylated lncRNA3 precipitated with avidin magnetic beads indicated clearly that TLS binds the fragment 4 (373-526 aa), C-terminus end of TLS. Then, dissecting fragment 4 presents four regions, RGG2, zinc finger, RGG3, and the nuclear localization signal (NLS) region in this order. Experiments with extensive deletion mutants indicated that just one deletion out of the four regions irs not enough to delete the TLS binding, although combinatorial deletion of zinc finger with other three regions almost wiped off the lncRNA3 binding. Remarkably, each of four regions alone has no binding to TLS, either. Collectively, RGG2, zinc finger, RGG3, and NLS all are essential for binding to lncRNA3, but are required to work synergistically for full binding. These data indicate that dynamic assembly of RNA-binding domain works for action of lncRNAs and possibly has allosteric effect on intrinsically disordered region (IDR) of N-terminus of TLS, implying relation of RNA-binding with phase separation and the resultant precipitation.

Understanding the Energy Landscape of Intrinsically Disordered Protein Ensembles.

A substantial portion of various organisms' proteomes comprises intrinsically disordered proteins (IDPs) that lack a defined three-dimensional structure. These IDPs exhibit a diverse array of conformations, displaying remarkable spatiotemporal heterogeneity and exceptional conformational flexibility. Characterizing the structure or structural ensemble of IDPs presents significant conceptual and methodological challenges owing to the absence of a well-defined native structure. While databases such as the Protein Ensemble Database (PED) provide IDP ensembles obtained through a combination of experimental data and molecular modeling, the absence of reaction coordinates poses challenges in comprehensively understanding pertinent aspects of the system. In this study, we leverage the energy landscape visualization method (JCTC, 6482, 2019) to scrutinize four IDP ensembles sourced from PED. ELViM, a methodology that circumvents the need for a priori reaction coordinates, aids in analyzing the ensembles. The specific IDP ensembles investigated are as follows: two fragments of nucleoporin (NUL: 884-993 and NUS: 1313-1390), yeast sic 1 N-terminal (1-90), and the N-terminal SH3 domain of Drk (1-59). Utilizing ELViM enables the comprehensive validation of ensembles, facilitating the detection of potential inconsistencies in the sampling process. Additionally, it allows for identifying and characterizing the most prevalent conformations within an ensemble. Moreover, ELViM facilitates the comparative analysis of ensembles obtained under diverse conditions, thereby providing a powerful tool for investigating the functional mechanisms of IDPs.

Dynamics of Single-Chain Nanoparticles under Crowding: A Neutron Spin Echo Study.

We present a neutron spin echo (NSE) investigation to examine the impact of macromolecular crowding on the dynamics of single-chain nanoparticles (SCNPs), serving as synthetic models for biomacromolecules with flexibility and internal degrees of freedom, such as intrinsically disordered proteins (IDPs). In particular, we studied the dynamics of a medium-size poly(methyl methacrylate) (PMMA)-based SCNP (33 kDa) in solutions with low- (10 kDa) and high- (100 kDa) molecular weight analogous deuterated PMMA linear crowders. The dynamic structure factors of the SCNPs in dilute solution show certain degrees of freedom, yet the analysis in terms of the Zimm model reveals high internal friction that effectively stiffens the chain-a phenomenon also observed for IDPs. Under crowding conditions, the internal dynamics remains essentially unchanged, but the center-of-mass diffusion slows down. The effective viscosity felt by the SCNPs at the timescales probed by NSE is lower than the macroscopic viscosity of the crowder solution, and it does not depend significantly on the molecular weight.

Unveiling the Binding between the Armadillo-Repeat Domain of Plakophilin 1 and the Intrinsically Disordered Transcriptional Repressor RYBP.

Plakophilin 1 (PKP1), a member of the p120ctn subfamily of the armadillo (ARM)-repeat-containing proteins, is an important structural component of cell-cell adhesion scaffolds although it can also be ubiquitously found in the cytoplasm and the nucleus. RYBP (RING 1A and YY1 binding protein) is a multifunctional intrinsically disordered protein (IDP) best described as a transcriptional regulator. Both proteins are involved in the development and metastasis of several types of tumors. We studied the binding of the armadillo domain of PKP1 (ARM-PKP1) with RYBP by using in cellulo methods, namely immunofluorescence (IF) and proximity ligation assay (PLA), and in vitro biophysical techniques, namely fluorescence, far-ultraviolet (far-UV) circular dichroism (CD), and isothermal titration calorimetry (ITC). We also characterized the binding of the two proteins by using in silico experiments. Our results showed that there was binding in tumor and non-tumoral cell lines. Binding in vitro between the two proteins was also monitored and found to occur with a dissociation constant in the low micromolar range (~10 μM). Finally, in silico experiments provided additional information on the possible structure of the binding complex, especially on the binding ARM-PKP1 hot-spot. Our findings suggest that RYBP might be a rescuer of the high expression of PKP1 in tumors, where it could decrease the epithelial-mesenchymal transition in some cancer cells.

Isolation of nucleic acids using liquid–liquid phase separation of pH-sensitive elastin-like polypeptides

Extraction of nucleic acids (NAs) is critical for many methods in molecular biology and bioanalytical chemistry. NA extraction has been extensively studied and optimized for a wide range of applications and its importance to society has significantly increased. The COVID-19 pandemic highlighted the importance of early and efficient NA testing, for which NA extraction is a critical analytical step prior to the detection by methods like polymerase chain reaction. This study explores simple, new approaches to extraction using engineered smart nanomaterials, namely NA-binding, intrinsically disordered proteins (IDPs), that undergo triggered liquid–liquid phase separation (LLPS). Two types of NA-binding IDPs are studied, both based on genetically engineered elastin-like polypeptides (ELPs), model IDPs that exhibit a lower critical solution temperature in water and can be designed to exhibit LLPS at desired temperatures in a variety of biological solutions. We show that ELP fusion proteins with natural NA-binding domains can be used to extract DNA and RNA from physiologically relevant solutions. We further show that LLPS of pH responsive ELPs that incorporate histidine in their sequences can be used for both binding, extraction and release of NAs from biological solutions, and can be used to detect SARS-CoV-2 RNA in samples from COVID-positive patients.

Hydrodynamic Radii of Intrinsically Disordered Proteins: Fast Prediction by Minimum Dissipation Approximation and Experimental Validation.

The diffusion coefficients of globular and fully unfolded proteins can be predicted with high accuracy solely from their mass or chain length. However, this approach fails for intrinsically disordered proteins (IDPs) containing structural domains. We propose a rapid predictive methodology for estimating the diffusion coefficients of IDPs. The methodology uses accelerated conformational sampling based on self-avoiding random walks and includes hydrodynamic interactions between coarse-grained protein subunits, modeled using the generalized Rotne-Prager-Yamakawa approximation. To estimate the hydrodynamic radius, we rely on the minimum dissipation approximation recently introduced by Cichocki et al. Using a large set of experimentally measured hydrodynamic radii of IDPs over a wide range of chain lengths and domain contributions, we demonstrate that our predictions are more accurate than the Kirkwood approximation and phenomenological approaches. Our technique may prove to be valuable in predicting the hydrodynamic properties of both fully unstructured and multidomain disordered proteins.

DR-BERT: A protein language model to annotate disordered regions

Despite their lack of a rigid structure, intrinsically disordered regions (IDRs) in proteins play important roles in cellular functions, including mediating protein-protein interactions. Therefore, it is important to computationally annotate IDRs with high accuracy. In this study, we present Disordered Region prediction using Bidirectional Encoder Representations from Transformers (DR-BERT), a compact protein language model. Unlike most popular tools, DR-BERT is pretrained on unannotated proteins and trained to predict IDRs without relying on explicit evolutionary or biophysical data. Despite this, DR-BERT demonstrates significant improvement over existing methods on the Critical Assessment of protein Intrinsic Disorder (CAID) evaluation datasetand outperforms competitors on two out of four test cases in the CAID 2 dataset, while maintaining competitiveness in the others. This performance is due to the information learned during pretraining and DR-BERT's ability to use contextual information.

Taxonomy-specific assessment of intrinsic disorder predictions at residue and region levels in higher eukaryotes, protists, archaea, bacteria and viruses

Intrinsic disorder predictors were evaluated in several studies including the two large CAID experiments. However, these studies are biased towards eukaryotic proteins and focus primarily on the residue-level predictions. We provide first-of-its-kind assessment that comprehensively covers the taxonomy and evaluates predictions at the residue and disordered region levels. We curate a benchmark dataset that uniformly covers eukaryotic, archaeal, bacterial, and viral proteins. We find that predictive performance differs substantially across taxonomy, where viruses are predicted most accurately, followed by protists and higher eukaryotes, while bacterial and archaeal proteins suffer lower levels of accuracy. These trends are consistent across predictors. We also find that current tools, except for flDPnn, struggle with reproducing native distributions of the numbers and sizes of the disordered regions. Moreover, analysis of two variants of disorder predictions derived from the AlphaFold2 predicted structures reveals that they produce accurate residue-level propensities for archaea, bacteria and protists. However, they underperform for higher eukaryotes and generally struggle to accurately identify disordered regions. Our results motivate development of new predictors that target bacteria and archaea and which produce accurate results at both residue and region levels. We also stress the need to include the region-level assessments in future assessments.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation.

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

HYPK: A marginally disordered protein sensitive to charge decoration

Intrinsically disordered proteins (IDPs) that lie close to the empirical boundary separating IDPs and folded proteins in Uversky's charge-hydropathy plot may behave as "marginal IDPs" and sensitively switch conformation upon changes in environment (temperature, crowding, and charge screening), sequence, or both. In our search for such a marginal IDP, we selected Huntingtin-interacting protein K (HYPK) near that boundary as a candidate; PKIα, also near that boundary, has lower secondary structure propensity; and Crk1, just across the boundary on the folded side, has higher secondary structure propensity. We used a qualitative Förster resonance energy transfer-based assay together with circular dichroism to simultaneously probe global and local conformation. HYPK shows several unique features indicating marginality: a cooperative transition in end-to-end distance with temperature, like Crk1 and folded proteins, but unlike PKIα; enhanced secondary structure upon crowding, in contrast to Crk1 and PKIα; and a cross-over from salt-induced expansion to compaction at high temperature, likely due to a structure-to-disorder transition not seen in Crk1 and PKIα. We then tested HYPK's sensitivity to charge patterning by designing charge-flipped variants including two specific sequences with identical amino acid composition that markedly differ in their predicted size and response to salt. The experimentally observed trends, also including mutants of PKIα, verify the predictions from sequence charge decoration metrics. Marginal proteins like HYPK show features of both folded and disordered proteins that make them sensitive to physicochemical perturbations and structural control by charge patterning.

GENESIS CGDYN: large-scale coarse-grained MD simulation with dynamic load balancing for heterogeneous biomolecular systems.

Residue-level coarse-grained (CG) molecular dynamics (MD) simulation is widely used to investigate slow biological processes that involve multiple proteins, nucleic acids, and their complexes. Biomolecules in a large simulation system are distributed non-uniformly, limiting computational efficiency with conventional methods. Here, we develop a hierarchical domain decomposition scheme with dynamic load balancing for heterogeneous biomolecular systems to keep computational efficiency even after drastic changes in particle distribution. These schemes are applied to the dynamics of intrinsically disordered protein (IDP) droplets. During the fusion of two droplets, we find that the changes in droplet shape correlate with the mixing of IDP chains. Additionally, we simulate large systems with multiple IDP droplets, achieving simulation sizes comparable to those observed in microscopy. In our MD simulations, we directly observe Ostwald ripening, a phenomenon where small droplets dissolve and their molecules redeposit into larger droplets. These methods have been implemented in CGDYN of the GENESIS software, offering a tool for investigating mesoscopic biological processes using the residue-level CG models.

Adenosine Triphosphate: The Primordial Molecule That Controls Protein Homeostasis and Shapes the Genome-Proteome Interface.

Adenosine triphosphate (ATP) acts as the universal energy currency that drives various biological processes, while nucleic acids function to store and transmit genetic information for all living organisms. Liquid-liquid phase separation (LLPS) represents the common principle for the formation of membrane-less organelles (MLOs) composed of proteins rich in intrinsically disordered regions (IDRs) and nucleic acids. Currently, while IDRs are well recognized to facilitate LLPS through dynamic and multivalent interactions, the precise mechanisms by which ATP and nucleic acids affect LLPS still remain elusive. This review summarizes recent NMR results on the LLPS of human FUS, TDP-43, and the viral nucleocapsid (N) protein of SARS-CoV-2, as modulated by ATP and nucleic acids, revealing the following: (1) ATP binds to folded domains overlapping with nucleic-acid-binding interfaces; (2) ATP and nucleic acids interplay to biphasically modulate LLPS by competitively binding to overlapping pockets of folded domains and Arg/Lys within IDRs; (3) ATP energy-independently induces protein folding with the highest efficiency known so far. As ATP likely emerged in the prebiotic monomeric world, while LLPS represents a pivotal mechanism to concentrate and compartmentalize rare molecules for forming primordial cells, ATP appears to control protein homeostasis and shape genome-proteome interfaces throughout the evolutionary trajectory, from prebiotic origins to modern cells.

Massive ER protein disposal by reticulophagy receptors and selective disposal by TOLLIP

ABSTRACT Proteostasis of the endoplasmic reticulum (ER) is maintained by coordinated action of two major catabolic pathways: proteasome-dependent ER-associated degradation (ERAD) and less characterized lysosomal pathways. Recent studies on ER-specific autophagy (termed “reticulophagy”) have highlighted the importance of lysosomes for ER proteostasis. Key to this process are proteins termed reticulophagy receptors that connect ER fragments and Atg8-family proteins, facilitating the lysosomal degradation of both native and aberrant ER proteins in a relatively nonselective manner. In contrast, our recent work identified TOLLIP as a novel type of cargo receptor specifically dedicated to the lysosomal degradation of aberrant ER membrane proteins. The clients of TOLLIP include an engineered model substrate, which mimics an ER-retained aberrant membrane protein, and motor neuron disease-linked misfolded mutants of VAPB and BSCL2/Seipin. TOLLIP acts as a receptor to connect these aberrant ER membrane proteins and phosphatidylinositol-3-phosphate (PtdIns3P) by recognizing the former through its misfolding-sensing intrinsically disordered region (IDR) and ubiquitin-binding CUE domain, and the latter through its C2 domain. These interactions enable PtdIns3P-dependent vesicular trafficking of aberrant membrane proteins to lysosomes without promoting reticulophagic turnover of bulk ER.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins.

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

Pathogenic mutations of human phosphorylation sites affect protein–protein interactions

Despite their lack of a defined 3D structure, intrinsically disordered regions (IDRs) of proteins play important biological roles. Many IDRs contain short linear motifs (SLiMs) that mediate protein-protein interactions (PPIs), which can be regulated by post-translational modifications like phosphorylation. 20% of pathogenic missense mutations are found in IDRs, and understanding how such mutations affect PPIs is essential for unraveling disease mechanisms. Here, we employ peptide-based interaction proteomics to investigate 36 disease-associated mutations affecting phosphorylation sites. Our results unveil significant differences in interactomes between phosphorylated and non-phosphorylated peptides, often due to disrupted phosphorylation-dependent SLiMs. We focused on a mutation of a serine phosphorylation site in the transcription factor GATAD1, which causes dilated cardiomyopathy. We find that this phosphorylation site mediates interaction with 14-3-3 family proteins. Follow-up experiments reveal the structural basis of this interaction and suggest that 14-3-3 binding affects GATAD1 nucleocytoplasmic transport by masking a nuclear localisation signal. Our results demonstrate that pathogenic mutations of human phosphorylation sites can significantly impact protein-protein interactions, offering insights into potential molecular mechanisms underlying pathogenesis.

Genetically Fusing Order-Promoting and Thermoresponsive Building Blocks to Design Hybrid Biomaterials.

The unique biophysical and biochemical properties of intrinsically disordered proteins (IDPs) and their recombinant derivatives, intrinsically disordered protein polymers (IDPPs) offer opportunities for producing multistimuli-responsive materials; their sequence-encoded disorder and tendency for phase separation facilitate the development of multifunctional materials. This review highlights the strategies for enhancing the structural diversity of elastin-like polypeptides (ELPs) and resilin-like polypeptides (RLPs), and their self-assembled structures via genetic fusion to ordered motifs such as helical or beta sheet domains. In particular, this review describes approaches that harness the synergistic interplay between order-promoting and thermoresponsive building blocks to design hybrid biomaterials, resulting in well-structured, stimuli-responsive supramolecular materials ordered on the nanoscale.

Neurofilament Biophysics: From Structure to Biomechanics.

Neurofilaments (NFs) are multisubunit, neuron-specific intermediate filaments consisting of a 10-nm diameter filament "core" surrounded by a layer of long intrinsically disordered protein (IDP) "tails." NFs are thought to regulate axonal caliber during development and then stabilize the mature axon, with NF subunit misregulation, mutation, and aggregation featuring prominently in multiple neurological diseases. The field's understanding of NF structure, mechanics, and function has been deeply informed by a rich variety of biochemical, cell biological, and mouse genetic studies spanning more than four decades. These studies have contributed much to our collective understanding of NF function in axonal physiology and disease. In recent years, however, there has been a resurgence of interest in NF subunit proteins in two new contexts: as potential blood- and cerebrospinal fluid-based biomarkers of neuronal damage, and as model IDPs with intriguing properties. Here, we review established principles and more recent discoveries in NF structure and function. Where possible, we place these findings in the context of biophysics of NF assembly, interaction, and contributions to axonal mechanics.

Intact Transition Epitope Mapping-Force Interferences by Variable Extensions (ITEM-FIVE).

Investigations on binding strength differences of non-covalent protein complex components were performed by mass spectrometry. T4 fibritin foldon (T4Ff) is a well-studied miniprotein, which together with its biotinylated version served as model system to represent a compactly folded protein to which an Intrinsically Disordered Region (IDR) was attached. The apparent enthalpies of the gas phase dissociation reactions of the homo-trimeric foldon F-F-F and of the homo-trimeric triply biotinylated foldon bF-bF-bF have been determined to be rather similar (3.32 kJ/mol and 3.85 kJ/mol) but quite distinct from those of the singly and doubly biotinylated hetero-trimers F-F-bF and F-bF-bF (1.86 kJ/mol and 1.08 kJ/mol). Molecular dynamics simulations suggest that the ground states of the (biotinylated) T4Ff trimers are highly symmetric and well comparable to each other, indicating that the energy levels of all four (biotinylated) T4Ff trimer ground states are nearly indistinguishable. The experimentally determined differences and/or similarities in enthalpies of the complex dissociation reactions are explained by entropic spring effects, which are noticeable in the T4Ff hetero-trimers but not in the T4Ff homo-trimers. A lowering of the transition state energy levels of the T4Ff hetero-trimers seems likely because the biotin moieties, mimicking intrinsically disordered regions (IDRs), induced asymmetries in the transition states of the biotinylated T4Ff hetero-trimers. This transition state energy level lowering effect is absent in the T4Ff homo-trimer, as well as in the triply biotinylated T4Ff homo-trimer. In the latter, the IDR-associated entropic spring effects on complex stability cancel each other out. ITEM-FIVE enabled semi-quantitative determination of energy differences of complex dissociation reactions, whose differences were modulated by IDRs attached to compactly folded proteins.

Contribution of the ELRs to the development of advanced in vitro models.

Developing in vitro models that accurately mimic the microenvironment of biological structures or processes holds substantial promise for gaining insights into specific biological functions. In the field of tissue engineering and regenerative medicine, in vitro models able to capture the precise structural, topographical, and functional complexity of living tissues, prove to be valuable tools for comprehending disease mechanisms, assessing drug responses, and serving as alternatives or complements to animal testing. The choice of the right biomaterial and fabrication technique for the development of these in vitro models plays an important role in their functionality. In this sense, elastin-like recombinamers (ELRs) have emerged as an important tool for the fabrication of in vitro models overcoming the challenges encountered in natural and synthetic materials due to their intrinsic properties, such as phase transition behavior, tunable biological properties, viscoelasticity, and easy processability. In this review article, we will delve into the use of ELRs for molecular models of intrinsically disordered proteins (IDPs), as well as for the development of in vitro 3D models for regenerative medicine. The easy processability of the ELRs and their rational design has allowed their use for the development of spheroids and organoids, or bioinks for 3D bioprinting. Thus, incorporating ELRs into the toolkit of biomaterials used for the fabrication of in vitro models, represents a transformative step forward in improving the accuracy, efficiency, and functionality of these models, and opening up a wide range of possibilities in combination with advanced biofabrication techniques that remains to be explored.