Disease mutations in disordered regions--exception to the rule?

Biomolecular condensates and disease pathogenesis

Biomolecular condensates or membraneless organelles (MLOs) formed by liquid-liquid phase separation (LLPS) divide intracellular spaces into discrete compartments for specific functions. Dysregulation of LLPS or aberrant phase transition that disturbs the formation or material states of MLOs is closely correlated with neurodegeneration, tumorigenesis, and many other pathological processes. Herein, we summarize the recent progress in development of methods to monitor phase separation and we discuss the biogenesis and function of MLOs formed through phase separation. We then present emerging proof-of-concept examples regarding the disruption of phase separation homeostasis in a diverse array of clinical conditions including neurodegenerative disorders, hearing loss, cancers, and immunological diseases. Finally, we describe the emerging discovery of chemical modulators of phase separation.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

The Proteomic Analysis of Cancer-Related Alterations in the Human Unfoldome

Many proteins lack stable 3D structures. These intrinsically disordered proteins (IDPs) or hybrid proteins containing ordered domains with intrinsically disordered protein regions (IDPRs) often carry out regulatory functions related to molecular recognition and signal transduction. IDPs/IDPRs constitute a substantial portion of the human proteome and are termed "the unfoldome". Herein, we probe the human breast cancer unfoldome and investigate relations between IDPs and key disease genes and pathways. We utilized bottom-up proteomics, MudPIT (Multidimensional Protein Identification Technology), to profile differentially expressed IDPs in human normal (MCF-10A) and breast cancer (BT-549) cell lines. Overall, we identified 2271 protein groups in the unfoldome of normal and cancer proteomes, with 148 IDPs found to be significantly differentially expressed in cancer cells. Further analysis produced annotations of 140 IDPs, which were then classified to GO (Gene Ontology) categories and pathways. In total, 65% (91 of 140) IDPs were related to various diseases, and 20% (28 of 140) mapped to cancer terms. A substantial portion of the differentially expressed IDPs contained disordered regions, confirmed by in silico characterization. Overall, our analyses suggest high levels of interactivity in the human cancer unfoldome and a prevalence of moderately and highly disordered proteins in the network.

Prediction of protein structure and AI

AlphaFold, an artificial intelligence (AI)-based tool for predicting the 3D structure of proteins, is now widely recognized for its high accuracy and versatility in the folding of human proteins. AlphaFold is useful for understanding structure-function relationships from protein 3D structure models and can serve as a template or a reference for experimental structural analysis including X-ray crystallography, NMR and cryo-EM analysis. Its use is expanding among researchers, not only in structural biology but also in other research fields. Researchers are currently exploring the full potential of AlphaFold-generated protein models. Predicting disease severity caused by missense mutations is one such application. This article provides an overview of the 3D structural modeling of AlphaFold based on deep learning techniques and highlights the challenges in predicting the pathogenicity of missense mutations.

Challenges in Brugada Syndrome Stratification: Investigating SCN5A Mutation Localization and Clinical Phenotypes

Brugada Syndrome (BrS) is a genetic heart condition linked to sudden cardiac death. Though the SCN5A gene is primarily associated with BrS, there is a lack of comprehensive studies exploring the connection between SCN5A mutation locations and the clinical presentations of the syndrome. This study aimed to address this gap and gain further understanding of the syndrome. The investigation classified 36 high-risk BrS patients based on SCN5A mutations within the transmembrane/structured (TD) and intra-domain loops (IDLs) lacking a 3D structure. We characterized the intrinsically disordered regions (IDRs) abundant in IDLs, using bioinformatics tools to predict IDRs and post-translational modifications (PTMs) in NaV1.5. Interestingly, it was found that current predictive tools often underestimate the impacts of mutations in IDLs and disordered regions. Moreover, patients with SCN5A mutations confined to IDL regions-previously deemed 'benign'-displayed clinical symptoms similar to those carrying 'damaging' variants. Our research illuminates the difficulty in stratifying patients based on SCN5A mutation locations, emphasizing the vital role of IDLs in the NaV1.5 channel's functioning and protein interactions. We advocate for caution when using predictive tools for mutation evaluation in these regions and call for the development of improved strategies in accurately assessing BrS risk.

Prediction of Quality-control Degradation Signals in Yeast Proteins

Effective proteome homeostasis is key to cellular and organismal survival, and cells therefore contain efficient quality control systems to monitor and remove potentially toxic misfolded proteins. Such general protein quality control to a large extent relies on the efficient and robust delivery of misfolded or unfolded proteins to the ubiquitin-proteasome system. This is achieved via recognition of so-called degradation motifs-degrons-that are assumed to become exposed as a result of protein misfolding. Despite their importance, the nature and sequence properties of quality-control degrons remain elusive. Here, we have used data from a yeast-based screen of 23,600 17-residue peptides to build a predictor of quality-control degrons. The resulting model, QCDPred (Quality Control Degron Prediction), achieves good accuracy using only the sequence composition of the peptides as input. Our analysis reveals that strong degrons are enriched in hydrophobic amino acids and depleted in negatively charged amino acids, in line with the expectation that they are buried in natively folded proteins. We applied QCDPred to the yeast proteome, enabling us to analyse more widely the potential effects of degrons. As an example, we show a correlation between cellular abundance and degron potential in disordered regions of proteins. Together with recent results on membrane proteins, our work suggest that the recognition of exposed hydrophobic residues is a key and generic mechanism for proteome homeostasis. QCDPred is freely available as open source code and via a web interface.

Mapping the Constrained Coding Regions in the Human Genome to Their Corresponding Proteins

Constrained Coding Regions (CCRs) in the human genome have been derived from DNA sequencing data of large cohorts of healthy control populations, available in the Genome Aggregation Database (gnomAD) [1]. They identify regions depleted of protein-changing variants and thus identify segments of the genome that have been constrained during human evolution. By mapping these DNA-defined regions from genomic coordinates onto the corresponding protein positions and combining this information with protein annotations, we have explored the distribution of CCRs and compared their co-occurrence with different protein functional features, previously annotated at the amino acid level in public databases. As expected, our results reveal that functional amino acids involved in interactions with DNA/RNA, protein-protein contacts and catalytic sites are the protein features most likely to be highly constrained for variation in the control population. More surprisingly, we also found that linear motifs, linear interacting peptides (LIPs), disorder-order transitions upon binding with other protein partners and liquid-liquid phase separating (LLPS) regions are also strongly associated with high constraint for variability. We also compared intra-species constraints in the human CCRs with inter-species conservation and functional residues to explore how such CCRs may contribute to the analysis of protein variants. As has been previously observed, CCRs are only weakly correlated with conservation, suggesting that intraspecies constraints complement interspecies conservation and can provide more information to interpret variant effects.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Identification of molecular subtypes based on liquid-liquid phase separation and cross-talk with immunological phenotype in bladder cancer

Background

Mounting evidence has demonstrated that an imbalance in liquid-liquid phase separation (LLPS) can induce alteration in the spatiotemporal coordination of biomolecular condensates, which plays a role in carcinogenesis and cachexia. However, the role of LLPS in the occurrence and progression of bladder cancer (BLCA) remains to be elucidated. Identifying the role of LLPS in carcinogenesis may aid in cancer therapeutics.

Methods

A total of 1,351 BLCA samples from six cohorts were retrieved from publicly available databases like The Cancer Genome Atlas, Gene Expression Omnibus, and ArrayExpress. The samples were divided into three distinct clusters, and their multi-dimensional heterogeneities were explored. The LLPS patterns of all patients were determined based on the LLPS-related risk score (LLPSRS), and its multifaceted landscape was depicted and experimentally validated at the multi-omics level. Finally, a cytotoxicity-related and LLPSRS-based classifier was established to predict the patient's response to immune checkpoint blockade (ICB) treatment.

Results

Three LLPS-related subtypes were identified and validated. The differences in prognosis, tumor microenvironment (TME) features, cancer hallmarks, and certain signatures of the three LLPS-related subtypes were validated. LLPSRS was calculated, which could be used as a prognostic biomarker. A close correlation was observed between clinicopathological features, genomic variations, biological mechanisms, immune infiltration in TME, chemosensitivity, and LLPSRS. Furthermore, our classifier could effectively predict immunotherapy response in patients with BLCA.

Conclusions

Our study identified a novel categorization of BLCA patients based on LLPS. The LLPSRS could predict the prognosis of patients and aid in designing personalized medicine. Further, our binary classifier could effectively predict patients' sensitivity to immunotherapy.

Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach

This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.

A large-scale polygenic risk score analysis identified candidate proteins associated with anxiety, depression and neuroticism

Psychiatric disorders and neuroticism are closely associated with central nervous system, whose proper functioning depends on efficient protein renewal. This study aims to systematically analyze the association between anxiety / depression / neuroticism and each of the 439 proteins. 47,536 pQTLs of 439 proteins in brain, plasma and cerebrospinal fluid (CSF) were collected from recent genome-wide association study. Polygenic risk scores (PRS) of the 439 proteins were then calculated using the UK Biobank cohort, including 120,729 subjects of neuroticism, 255,354 subjects of anxiety and 316,513 subjects of depression. Pearson correlation analyses were performed to evaluate the correlation between each protein and each of the mental traits by using calculated PRSs as the instrumental variables of protein. In general population, six correlations were identified in plasma and CSF such as plasma protease C1 inhibitor (C1-INH) with neuroticism score (r = - 0.011, P = 2.56 × 10^- 9) in plasma, C1-INH with neuroticism score (r = -0.010, P = 3.09 × 10^- 8) in CSF, and ERBB1 with self-reported depression (r = - 0.012, P = 4.65 × 10^- 5) in CSF. C1-INH and ERBB1 may induce neuroticism and depression by affecting brain function and synaptic development. Gender subgroup analyses found that BST1 was correlated with neuroticism score in male CSF (r = - 0.011, P = 1.80 × 10^- 5), while CNTN2 was correlated with depression score in female brain (r = - 0.013, P = 6.43 × 10^- 4). BST1 and CNTN2 may be involved in nervous system metabolism and brain health. Six common candidate proteins were associated with all three traits (P < 0.05) and were confirmed in relevant proteomic studies, such as C1-INH in plasma, CNTN2 and MSP in the brain. Our results provide novel clues for revealing the roles of proteins in the development of anxiety, depression and neuroticism.

Negative Thermal Expansion and Disorder-to-Order Collapse of an Intrinsically Disordered Protein under Marginally Denaturing Conditions

Recent work with intrinsically disordered proteins (IDPs) has projected a myriad of their survival instincts based mainly on the total charge content, the abundance of polar residues, and the paucity of hydrophobic amino acids. This work uses a plant IDP AtPP16-1 (Arabidopsis thaliana phloem protein class 16-1), whose solution NMR structure was determined by us recently, to show legitimate negative thermal expansion (NTE) of the native state. The thermal expansion continues to be negative even as the tertiary structure is perturbed by ultralow levels of urea up to 0.4 M. The NTE of these subdenatured states is called apparent NTE because they are prone to undergo conformational changes with temperature. Hydrodynamic shrinkage of the NTE IDP is also observed by dynamic light scattering (DLS) and NMR-measured global rotational correlation time (τ_c). The protein with denatured tertiary structure but otherwise preserved native-state secondary structure collapses to a dynamically rigid state. The data are mainly based on thermal coefficients of chemical shift and nuclear relaxation measured by heteronuclear NMR. The hydrodynamic shrinkage and collapse under marginally varying solvent compositions that may arise from unstable tertiary structure and dynamic disorder of chain segments across the backbone could be a generic property of IDPs.

Pan-cancer assessment of mutational landscape in intrinsically disordered hotspots reveals potential driver genes

Large-scale cancer genome sequencing has enabled the catalogs of somatic mutations; however, the mutational impact on intrinsically disordered protein regions (IDRs) has not been systematically investigated to date. Here, we comprehensively characterized the mutational landscapes of IDRs and found that IDRs have higher mutation frequencies across diverse cancers. We thus developed a computational method, ROI-Driver, to identify putative driver genes enriching IDR and domain hotspots in cancer. Numerous well-known cancer-related oncogenes or tumor suppressors that play important roles in cancer signaling regulation, development and immune response were identified at a higher resolution. In particular, the incorporation of IDR structures helps in the identification of novel potential driver genes that play central roles in human protein-protein interaction networks. Interestingly, we found that the putative driver genes with IDR hotspots were significantly enriched with predicted phase separation propensities, suggesting that IDR mutations disrupt phase separation in key cellular pathways. We also identified an appreciable number of clinically relevant genes enriching IDR mutational hotspots that exhibited differential expression patterns and are associated with cancer patient survival. In summary, combinations of mutational effects on IDRs significantly increase the sensitivity of driver detection and are likely to open new therapeutic avenues for various cancers.

Phase-Separated Subcellular Compartmentation and Related Human Diseases

In live cells, proteins and nucleic acids can associate together through multivalent interactions, and form relatively isolated phases that undertake designated biological functions and activities. In the past decade, liquid–liquid phase separation (LLPS) has gradually been recognized as a general mechanism for the intracellular organization of biomolecules. LLPS regulates the assembly and composition of dozens of membraneless organelles and condensates in cells. Due to the altered physiological conditions or genetic mutations, phase-separated condensates may undergo aberrant formation, maturation or gelation that contributes to the onset and progression of various diseases, including neurodegenerative disorders and cancers. In this review, we summarize the properties of different membraneless organelles and condensates, and discuss multiple phase separation-regulated biological processes. Based on the dysregulation and mutations of several key regulatory proteins and signaling pathways, we also exemplify how aberrantly regulated LLPS may contribute to human diseases.

Structural dynamics shape the fitness window of alanine:glyoxylate aminotransferase

The conformational landscape of a protein is constantly expanded by genetic variations that have a minimal impact on the function(s) while causing subtle effects on protein structure. The wider the conformational space sampled by these variants, the higher the probabilities to adapt to changes in environmental conditions. However, the probability that a single mutation may result in a pathogenic phenotype also increases. Here we present a paradigmatic example of how protein evolution balances structural stability and dynamics to maximize protein adaptability and preserve protein fitness. We took advantage of known genetic variations of human alanine:glyoxylate aminotransferase (AGT1), which is present as a common major allelic form (AGT-Ma) and a minor polymorphic form (AGT-Mi) expressed in 20% of Caucasian population. By integrating crystallographic studies and molecular dynamics simulations, we show that AGT-Ma is endowed with structurally unstable (frustrated) regions, which become disordered in AGT-Mi. An in-depth biochemical characterization of variants from an anticonsensus library, encompassing the frustrated regions, correlates this plasticity to a fitness window defined by AGT-Ma and AGT-Mi. Finally, co-immunoprecipitation analysis suggests that structural frustration in AGT1 could favor additional functions related to protein-protein interactions. These results expand our understanding of protein structural evolution by establishing that naturally occurring genetic variations tip the balance between stability and frustration to maximize the ensemble of conformations falling within a well-defined fitness window, thus expanding the adaptability potential of the protein.

Shell Disorder Models Detect That Omicron Has Harder Shells with Attenuation but Is Not a Descendant of the Wuhan-Hu-1 SARS-CoV-2

Before the SARS-CoV-2 Omicron variant emergence, shell disorder models (SDM) suggested that an attenuated precursor from pangolins may have entered humans in 2017 or earlier. This was based on a shell disorder analysis of SARS-CoV-1/2 and pangolin-Cov-2017. The SDM suggests that Omicron is attenuated with almost identical N (inner shell) disorder as pangolin-CoV-2017 (N-PID (percentage of intrinsic disorder): 44.8% vs. 44.9%-lower than other variants). The outer shell disorder (M-PID) of Omicron is lower than that of other variants and pangolin-CoV-2017 (5.4% vs. 5.9%). COVID-19-related CoVs have the lowest M-PIDs (hardest outer shell) among all CoVs. This is likely to be responsible for the higher contagiousness of SARS-CoV-2 and Omicron, since hard outer shell protects the virion from salivary/mucosal antimicrobial enzymes. Phylogenetic study using M reveals that Omicron branched off from an ancestor of the Wuhan-Hu-1 strain closely related to pangolin-CoVs. M, being evolutionarily conserved in COVID-19, is most ideal for COVID-19 phylogenetic study. Omicron may have been hiding among burrowing animals (e.g., pangolins) that provide optimal evolutionary environments for attenuation and increase shell hardness, which is essential for fecal-oral-respiratory transmission via buried feces. Incoming data support SDM e.g., the presence of fewer infectious particles in the lungs than in the bronchi upon infection.

Computational and experimental methods for classifying variants of unknown clinical significance

The increase in sequencing capacity, reduction in costs, and national and international coordinated efforts have led to the widespread introduction of next-generation sequencing (NGS) technologies in patient care. More generally, human genetics and genomic medicine are gaining importance for more and more patients. Some communities are already discussing the prospect of sequencing each individual's genome at time of birth. Together with digital health records, this shall enable individualized treatments and preventive measures, so-called precision medicine. A central step in this process is the identification of disease causal mutations or variant combinations that make us more susceptible for diseases. Although various technological advances have improved the identification of genetic alterations, the interpretation and ranking of the identified variants remains a major challenge. Based on our knowledge of molecular processes or previously identified disease variants, we can identify potentially functional genetic variants and, using different lines of evidence, we are sometimes able to demonstrate their pathogenicity directly. However, the vast majority of variants are classified as variants of uncertain clinical significance (VUSs) with not enough experimental evidence to determine their pathogenicity. In these cases, computational methods may be used to improve the prioritization and an increasing toolbox of experimental methods is emerging that can be used to assay the molecular effects of VUSs. Here, we discuss how computational and experimental methods can be used to create catalogs of variant effects for a variety of molecular and cellular phenotypes. We discuss the prospects of integrating large-scale functional data with machine learning and clinical knowledge for the development of accurate pathogenicity predictions for clinical applications.

Modeling the Dynamics of Protein-Protein Interfaces, How and Why?

Protein-protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein-protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein-protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein-protein interface.

An integrative multi-omics analysis based on liquid–liquid phase separation delineates distinct subtypes of lower-grade glioma and identifies a prognostic signature

Background

Emerging evidences have indicated that the aberrant liquid–liquid phase separation (LLPS) leads to the dysfunction of biomolecular condensates, thereby contributing to the tumorigenesis and progression. Nevertheless, it remains unclear whether or how the LLPS of specific molecules affects the prognosis and tumor immune microenvironment (TIME) of patients with lower-grade glioma (LGG).

Methods

We integrated the transcriptome information of 3585 LLPS-related genes to comprehensively evaluate the LLPS patterns of 423 patients with LGG in The Cancer Genome Atlas (TCGA) cohort. Then, we systematically demonstrated the differences among four LLPS subtypes based on multi-omics analyses. In addition, we constructed the LLPS-related prognostic risk score (LPRS) for individualized integrative assessment.

Results

Based on the expression profiles of 85 scaffolds, 355 regulators, and 3145 clients in LGG, we identified four LLPS subtypes, namely LS1, LS2, LS3 and LS4.

We confirmed that there were significant differences in prognosis, clinicopathological features, cancer hallmarks, genomic alterations, TIME patterns and immunotherapeutic responses among four LLPS subtypes. In addition, a prognostic signature called LPRS was constructed for individualized integrative assessment. LPRS exhibited a robust predictive capacity for prognosis of LGG patients in multiple cohorts. Moreover, LPRS was found to be correlated with clinicopathological features, cancer hallmarks, genomic alterations and TIME patterns of LGG patients. The predictive power of LPRS in response to immune checkpoint inhibitor (ICI) therapy was also prominent.

Conclusions

This study provided a novel classification of LGG patients based on LLPS. The constructed LPRS might facilitate individualized prognosis prediction and better immunotherapy options for LGG patients.

Understanding p300-transcription factor interactions using sequence variation and hybridization

The hypoxic response is central to cell function and plays a significant role in the growth and survival of solid tumours. HIF-1 regulates the hypoxic response by activating over 100...

Effect of Disease Causing Missense Mutations on Intrinsically Disordered Regions in Proteins

BACKGROUND

It is well known that disease-causing missense mutations (DCMMs) reduce the structural stability/integrity of the proteins with well-defined 3D structures, thereby impacting their molecular functions. However, it is not known in what way DCMMs affect the intrinsically disordered proteins (IDPs) that do not adopt well defined stable 3D structures.

METHODS

In order to investigate how DCMMs may impact intrinsically disordered regions (IDRs) in proteins, we undertook Molecular Dynamics (MD) based studies on three different examples of functionally important IDRs with known DCMMs. Our studies revealed that the functional impact of DCMMs is in reducing the conformational heterogeneity of IDRs, which is intrinsic and quintessential for their multi-faceted cellular roles.

RESULTS

These results are reinforced by energy landscapes of the wildtype and mutant IDRs where the former is characterized by many local minima separated by low barriers, whereas the latter are characterized by one global minimum and several local minima separated by high energy barriers. Our MD based studies also indicate that DCMMs stabilize very few structural possibilities of IDRs either by the newly formed interactions induced by the substituted side chains or by means of restricted or increased flexibilities of the backbone conformations at the mutation sites.

CONCLUSION

Furthermore, the structural possibilities stabilized by DCMMs do not support the native functional roles of the IDRs, thereby leading to disease conditions.

Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions

Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.

Cancer-Associated Mutations Perturb the Disordered Ensemble and Interactions of the Intrinsically Disordered p53 Transactivation Domain

Intrinsically disordered proteins (IDPs) are key components of regulatory networks that control crucial aspects of cell decision making. The intrinsically disordered transactivation domain (TAD) of tumor suppressor p53 mediates its interactions with multiple regulatory pathways to control the p53 homeostasis during the cellular response to genotoxic stress. Many cancer-associated mutations have been discovered in p53-TAD, but their structural and functional consequences are poorly understood. Here, by combining atomistic simulations, NMR spectroscopy, and binding assays, we demonstrate that cancer-associated mutations can significantly perturb the balance of p53 interactions with key activation and degradation regulators. Importantly, the four mutations studied in this work do not all directly disrupt the known interaction interfaces. Instead, at least three of these mutations likely modulate the disordered state of p53-TAD to perturb its interactions with regulators. Specifically, NMR and simulation analysis together suggest that these mutations can modulate the level of conformational expansion as well as rigidity of the disordered state. Our work suggests that the disordered conformational ensemble of p53-TAD can serve as a central conduit in regulating the response to various cellular stimuli at the protein-protein interaction level. Understanding how the disordered state of IDPs may be modulated by regulatory signals and/or disease associated perturbations will be essential in the studies on the role of IDPs in biology and diseases.

On the specificity of protein-protein interactions in the context of disorder

With the increased focus on intrinsically disordered proteins (IDPs) and their large interactomes, the question about their specificity — or more so on their multispecificity — arise. Here we recapitulate how specificity and multispecificity are quantified and address through examples if IDPs in this respect differ from globular proteins. The conclusion is that quantitatively, globular proteins and IDPs are similar when it comes to specificity. However, compared with globular proteins, IDPs have larger interactome sizes, a phenomenon that is further enabled by their flexibility, repetitive binding motifs and propensity to adapt to different binding partners. For IDPs, this adaptability, interactome size and a higher degree of multivalency opens for new interaction mechanisms such as facilitated exchange through trimer formation and ultra-sensitivity via threshold effects and ensemble redistribution. IDPs and their interactions, thus, do not compromise the definition of specificity. Instead, it is the sheer size of their interactomes that complicates its calculation. More importantly, it is this size that challenges how we conceptually envision, interpret and speak about their specificity.

IDRMutPred: predicting disease-associated germline nonsynonymous single nucleotide variants (nsSNVs) in intrinsically disordered regions

Motivation

Despite of the lack of folded structure, intrinsically disordered regions (IDRs) of proteins play versatile roles in various biological processes, and many nonsynonymous single nucleotide variants (nsSNVs) in IDRs are associated with human diseases. The continuous accumulation of nsSNVs resulted from the wide application of NGS has driven the development of disease-association prediction methods for decades. However, their performance on nsSNVs in IDRs remains inferior, possibly due to the domination of nsSNVs from structured regions in training data. Therefore, it is highly demanding to build a disease-association predictor specifically for nsSNVs in IDRs with better performance.

Results

We present IDRMutPred, a machine learning-based tool specifically for predicting disease-associated germline nsSNVs in IDRs. Based on 17 selected optimal features that are extracted from sequence alignments, protein annotations, hydrophobicity indices and disorder scores, IDRMutPred was trained using three ensemble learning algorithms on the training dataset containing only IDR nsSNVs. The evaluation on the two testing datasets shows that all the three prediction models outperform 17 other popular general predictors significantly, achieving the ACC between 0.856 and 0.868 and MCC between 0.713 and 0.737. IDRMutPred will prioritize disease-associated IDR germline nsSNVs more reliably than general predictors.

Availability and implementation

The software is freely available at http://www.wdspdb.com/IDRMutPred.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

A disease-causing mutation K240E disrupts ferroportin trafficking by SUMO (ferroportin SUMOylation)

Ferroportin (Fpn/IREG1/MTP1) is the only known transporter mediating iron efflux from epithelial cells and macrophages, and thus regulates how much iron is released into the circulation. Consequently, Fpn mutations are associated with haemochromatosis. Fpn itself is post-translationally regulated by hepcidin (Hepc) which induces its redistribution and degradation in a ubiquitin-dependent process. Together, the two proteins appear to be the nexus for iron homeostasis. Here we show that a rare gain-of-function mutation (K240E) that is associated with iron overload, impedes Fpn binding and subcellular trafficking by the small ubiquitin-like modifier (SUMO). Whereas wild-type Fpn is ensconced within vesicular bodies, the FpnK240E mutant appeared diffused within the cell when co-expressed with SUMO. Furthermore, compared with wild type Fpn, the sumoylation-defective mutant was constitutively-active, resulting in a lower intracellular labile iron pool than the former. These findings suggest that SUMO may regulate iron homeostasis by controlling Fpn trafficking.

•

Ferroportin (Fpn) regulates iron efflux.

•

A disease causing mutation (K240E) in a patient causes iron-overload.

•

Fpn K240 is a SUMO conjugation site important for Fpn trafficking to endosomes by SUMO.

•

The Fpn mutant K240E cannot be trafficked properly by SUMO and is a gain-of-function mutant that is constitutively active.

•

FpnK240E effluxes more iron from intracellular stores than wild type Fpn.

Design principles of gene evolution for niche adaptation through changes in protein-protein interaction networks

In contrast to fossorial and above-ground organisms, subterranean species have adapted to the extreme stresses of living underground. We analyzed the predicted protein–protein interactions (PPIs) of all gene products, including those of stress-response genes, among nine subterranean, ten fossorial, and 13 aboveground species. We considered 10,314 unique orthologous protein families and constructed 5,879,879 PPIs in all organisms using ChiPPI. We found strong association between PPI network modulation and adaptation to specific habitats, noting that mutations in genes and changes in protein sequences were not linked directly with niche adaptation in the organisms sampled. Thus, orthologous hypoxia, heat-shock, and circadian clock proteins were found to cluster according to habitat, based on PPIs rather than on sequence similarities. Curiously, "ordered" domains were preserved in aboveground species, while "disordered" domains were conserved in subterranean organisms, and confirmed for proteins in DistProt database. Furthermore, proteins with disordered regions were found to adopt significantly less optimal codon usage in subterranean species than in fossorial and above-ground species. These findings reveal design principles of protein networks by means of alterations in protein domains, thus providing insight into deep mechanisms of evolutionary adaptation, generally, and particularly of species to underground living and other confined habitats.

A variant in the RP1L1 gene in a family with occult macular dystrophy in a predicted intrinsically disordered region

SIGNIFICANCE

The responsible genetic variants for occult macular dystrophy (OMD) were found at the predicted intrinsically disordered region (IDR) of the RP1L1 gene.

PURPOSE

We examined the phenotypes and genotypes of family members from OMD. In addition, the genetic characteristics of the RP1L1 gene in OMD were investigated.

METHODS

Whole-exome sequencing was applied on two affected family members, and Sanger sequencing was performed on three members. The structural property of RP1L1 and pathogenic variants was analyzed using predictor of natural disordered regions (PONDR).

RESULTS

Two affected members showed moderate visual impairment and relative central scotoma. The spectral domain optical coherence tomography (SD-OCT) images showed an absence of the interdigitation zone (IZ) and ellipsoid zone (EZ) in one case, and an obscure EZ line in the other case. A RP1L1 variant (c.3593 C > T, p.Ser1198Phe) was identified in two affected members but not in the unaffected member. The PONDR analysis showed that the region from p.1189 to p.1248 could be predicted to be an IDR in the RP1L1 molecule. And the p. Ser1198Phe variant showed significant reduction of PONDR score.

CONCLUSIONS

Although, the major pathogenic variant of OMD is p.Arg45Trp, multiple reports indicate that the region between p.1194 and p.1201 is another hot spot of OMD. The PONDR analysis predicted that the RP1L1 molecule is one of the intrinsically disordered proteins. It is speculated that the region around p.1200 is essential for the normal function of the RP1L1 molecule, and the missense variants of that area cause the development of OMD.

Targeting Intrinsically Disordered Proteins through Dynamic Interactions

Intrinsically disordered proteins (IDPs) are over-represented in major disease pathways and have attracted significant interest in understanding if and how they may be targeted using small molecules for therapeutic purposes. While most existing studies have focused on extending the traditional structure-centric drug design strategies and emphasized exploring pre-existing structure features of IDPs for specific binding, several examples have also emerged to suggest that small molecules could achieve specificity in binding IDPs and affect their function through dynamic and transient interactions. These dynamic interactions can modulate the disordered conformational ensemble and often lead to modest compaction to shield functionally important interaction sites. Much work remains to be done on further elucidation of the molecular basis of the dynamic small molecule-IDP interaction and determining how it can be exploited for targeting IDPs in practice. These efforts will rely critically on an integrated experimental and computational framework for disordered protein ensemble characterization. In particular, exciting advances have been made in recent years in enhanced sampling techniques, Graphic Processing Unit (GPU)-computing, and protein force field optimization, which have now allowed rigorous physics-based atomistic simulations to generate reliable structure ensembles for nontrivial IDPs of modest sizes. Such de novo atomistic simulations will play crucial roles in exploring the exciting opportunity of targeting IDPs through dynamic interactions.

Modulation of p53 Transactivation Domain Conformations by Ligand Binding and Cancer-Associated Mutations

Intrinsically disordered proteins (IDPs) are important functional proteins, and their deregulation are linked to numerous human diseases including cancers. Understanding how disease-associated mutations or drug molecules can perturb the sequence-disordered ensemble-function-disease relationship of IDPs remains challenging, because it requires detailed characterization of the heterogeneous structural ensembles of IDPs. In this work, we combine the latest atomistic force field a99SB-disp, enhanced sampling technique replica exchange with solute tempering, and GPU-accelerated molecular dynamics simulations to investigate how four cancer-associated mutations, K24N, N29K/N30D, D49Y, and W53G, and binding of an anti-cancer molecule, epigallocatechin gallate (EGCG), modulate the disordered ensemble of the transactivation domain (TAD) of tumor suppressor p53. Through extensive sampling, in excess of 1.0 μs per replica, well-converged structural ensembles of wild-type and mutant p53-TAD as well as WT p53-TAD in the presence of EGCG were generated. The results reveal that mutants could induce local structural changes and affect secondary structural properties. Interestingly, both EGCG binding and N29K/N30D could also induce long-range structural reorganizations and lead to more compact structures that could shield key binding sites of p53-TAD regulators. Further analysis reveals that the effects of EGCG binding are mainly achieved through nonspecific interactions. These observations are generally consistent with on-going NMR studies and binding assays. Our studies suggest that induced conformational collapse of IDPs may be a general mechanism for shielding functional sites, thus inhibiting recognition of their targets. The current study also demonstrates that atomistic simulations provide a viable approach for studying the sequence-disordered ensemble-function-disease relationships of IDPs and developing new drug design strategies targeting regulatory IDPs.

Pathogenicity and functional impact of non-frameshifting insertion/deletion variation in the human genome

Differentiation between phenotypically neutral and disease-causing genetic variation remains an open and relevant problem. Among different types of variation, non-frameshifting insertions and deletions (indels) represent an understudied group with widespread phenotypic consequences. To address this challenge, we present a machine learning method, MutPred-Indel, that predicts pathogenicity and identifies types of functional residues impacted by non-frameshifting insertion/deletion variation. The model shows good predictive performance as well as the ability to identify impacted structural and functional residues including secondary structure, intrinsic disorder, metal and macromolecular binding, post-translational modifications, allosteric sites, and catalytic residues. We identify structural and functional mechanisms impacted preferentially by germline variation from the Human Gene Mutation Database, recurrent somatic variation from COSMIC in the context of different cancers, as well as de novo variants from families with autism spectrum disorder. Further, the distributions of pathogenicity prediction scores generated by MutPred-Indel are shown to differentiate highly recurrent from non-recurrent somatic variation. Collectively, we present a framework to facilitate the interrogation of both pathogenicity and the functional effects of non-frameshifting insertion/deletion variants. The MutPred-Indel webserver is available at http://mutpred.mutdb.org/.

Gain-of-Function Mutations: An Emerging Advantage for Cancer Biology

Advances in next-generation sequencing have identified thousands of genomic variants that perturb the normal functions of proteins, further contributing to diverse phenotypic consequences in cancer. Elucidating the functional pathways altered by loss-of-function (LOF) or gain-of-function (GOF) mutations will be crucial for prioritizing cancer-causing variants and their resultant therapeutic liabilities. In this review, we highlight the fundamental function of GOF mutations and discuss the potential mechanistic effects in the context of signaling networks. We also summarize advances in experimental and computational resources, which will dramatically help with studies on the functional and phenotypic consequences of mutations. Together, systematic investigations of the function of GOF mutations will provide an important missing piece for cancer biology and precision therapy.

Structural Perspective on Revealing and Altering Molecular Functions of Genetic Variants Linked with Diseases

Structural information of biological macromolecules is crucial and necessary to deliver predictions about the effects of mutations-whether polymorphic or deleterious (i.e., disease causing), wherein, thermodynamic parameters, namely, folding and binding free energies potentially serve as effective biomarkers. It may be emphasized that the effect of a mutation depends on various factors, including the type of protein (globular, membrane or intrinsically disordered protein) and the structural context in which it occurs. Such information may positively aid drug-design. Furthermore, due to the intrinsic plasticity of proteins, even mutations involving radical change of the structural and physico⁻chemical properties of the amino acids (native vs. mutant) can still have minimal effects on protein thermodynamics. However, if a mutation causes significant perturbation by either folding or binding free energies, it is quite likely to be deleterious. Mitigating such effects is a promising alternative to the traditional approaches of designing inhibitors. This can be done by structure-based in silico screening of small molecules for which binding to the dysfunctional protein restores its wild type thermodynamics. In this review we emphasize the effects of mutations on two important biophysical properties, stability and binding affinity, and how structures can be used for structure-based drug design to mitigate the effects of disease-causing variants on the above biophysical properties.

Predictions of Backbone Dynamics in Intrinsically Disordered Proteins Using De Novo Fragment-Based Protein Structure Predictions

Intrinsically disordaered proteins (IDPs) are a prevalent phenomenon with over 30% of human proteins estimated to have long disordered regions. Computational methods are widely used to study IDPs, however, nearly all treat disorder in a binary fashion, not accounting for the structural heterogeneity present in disordered regions. Here, we present a new de novo method, FRAGFOLD-IDP, which addresses this problem. Using 200 protein structural ensembles derived from NMR, we show that FRAGFOLD-IDP achieves superior results compared to methods which can predict related data (NMR order parameter, or crystallographic B-factor). FRAGFOLD-IDP produces very good predictions for 33.5% of cases and helps to get a better insight into the dynamics of the disordered ensembles. The results also show it is not necessary to predict the correct fold of the protein to reliably predict per-residue fluctuations. It implies that disorder is a local property and it does not depend on the fold. Our results are orthogonal to DynaMine, the only other method significantly better than the naïve prediction. We therefore combine these two using a neural network. FRAGFOLD-IDP enables better insight into backbone dynamics in IDPs and opens exciting possibilities for the design of disordered ensembles, disorder-to-order transitions, or design for protein dynamics.

The contribution of intrinsically disordered regions to protein function, cellular complexity, and human disease

In the 1960s, Christian Anfinsen postulated that the unique three-dimensional structure of a protein is determined by its amino acid sequence. This work laid the foundation for the sequence–structure–function paradigm, which states that the sequence of a protein determines its structure, and structure determines function. However, a class of polypeptide segments called intrinsically disordered regions does not conform to this postulate. In this review, I will first describe established and emerging ideas about how disordered regions contribute to protein function. I will then discuss molecular principles by which regulatory mechanisms, such as alternative splicing and asymmetric localization of transcripts that encode disordered regions, can increase the functional versatility of proteins. Finally, I will discuss how disordered regions contribute to human disease and the emergence of cellular complexity during organismal evolution.

Mapping genotype-phenotype associations of nsSNPs in coiled-coil oligomerization domains of the human proteome

We assessed the impact of disease mutations (DMs) versus polymorphisms (PYs) in coiled-coil (CC) domains in UniProt by modeling the structural and functional impact of variants in silico with the CC prediction program Multicoil. The structural impact of variants was evaluated with respect to three main metrics: the oligomerization score-to determine whether the variant is stabilizing or destabilizing-the oligomerization state, and the register-specific score. The functional impact was queried indirectly in several ways. First, we examined marginally stable CCs that were either stabilized or destabilized by the variant. Second, we looked for variants that altered the register of the wild-type CC near wild-type irregularities of likely functional importance, such as skips and stammers. Third, we searched for variants that altered the oligomerization state of the CC. DMs tended to be more destabilizing than PYs; but interestingly, PYs were more frequently associated with predicted changes in the oligomerization state. The functional impact was also queried by testing the association of CC variants with multiple phenotypes, that is, pleiotropy. Mutations in CC regions of proteins cause 155 different phenotypes and are more frequently associated with pleiotropy than proteins in general. Importantly, the CC region itself often encodes the pleiotropy.

Therapeutic Interventions of Cancers Using Intrinsically Disordered Proteins as Drug Targets: c-Myc as Model System

The concept of protein intrinsic disorder has taken the driving seat to understand regulatory proteins in general. Reports suggest that in mammals nearly 75% of signalling proteins contain long disordered regions with greater than 30 amino acid residues. Therefore, intrinsically disordered proteins (IDPs) have been implicated in several human diseases and should be considered as potential novel drug targets. Moreover, intrinsic disorder provides a huge multifunctional capability to hub proteins such as c-Myc and p53. c-Myc is the hot spot for understanding and developing therapeutics against cancers and cancer stem cells. Our past understanding is mainly based on in vitro and in vivo experiments conducted using c-Myc as whole protein. Using the reductionist approach, c-Myc oncoprotein has been divided into structured and disordered domains. A wealth of data is available dealing with the structured perspectives of c-Myc, but understanding c-Myc in terms of disordered domains has just begun. Disorderness provides enormous flexibility to proteins in general for binding to numerous partners. Here, we have reviewed the current progress on understanding c-Myc using the emerging concept of IDPs.

The role of protein intrinsic disorder in major psychiatric disorders

Although new candidate genes for Autism Spectrum Disorder (ASD), Schizophrenia (SCZ), Attention-Deficit/Hyperactivity Disorder (ADHD), and Bipolar Disorder (BD) emerged from genome-wide association studies (GWAS), their underlying molecular mechanisms remain poorly understood. Evidences of the involvement of intrinsically disordered proteins in diseases have grown in the last decade. These proteins lack tridimensional structure under physiological conditions and are involved in important cellular functions such as signaling, recognition and regulation. The aim of the present study was to identify the role and abundance of intrinsically disordered proteins in a set of psychiatric diseases and to test whether diseases are different regarding protein intrinsic disorder. Our hypothesis is that differences across psychiatric illnesses phenotypes and symptoms may arise from differences in intrinsic protein disorder content and properties of each group. A bioinformatics prediction of intrinsic disorder was performed in proteins retrieved based on top findings from GWAS, Copy Number Variation and candidate gene investigations for each disease. This approach revealed that about 80% of studied proteins presented long stretches of disorder. This amount was significantly higher than that observed in general eukaryotic proteins, and those involved in cardiovascular diseases. These results suggest that proteins with intrinsic disorder are a common feature of neurodevelopment and synaptic transmission processes which are potentially involved in the etiology of psychiatric diseases. Moreover, we identified differences between ADHD and ASD when the binary prediction of structure and putative binding sites were compared. These differences may be related to variation in symptom complexity between both diseases. © 2016 Wiley Periodicals, Inc.

The gene expression profiling of hepatocellular carcinoma by a network analysis approach shows a dominance of intrinsically disordered proteins (IDPs) between hub nodes

We have analyzed the transcriptomic data from patients with hepatocellular carcinoma (HCC) after viral HCV infection at the various stages of the disease by means of a networking analysis using the publicly available E-MTAB-950 dataset. The data was compared with those obtained in our group from HepG2 cells, a cancer cell line that lacks the viral infection. By sequential pruning of data, and also taking into account the data from cells of healthy patients as blanks, we were able to obtain a distribution of hub genes for the various stages that characterize the disease and finally, we isolated a metabolic sub-net specific to HCC alone. The general picture is that the basic organization to energetically and metabolically sustain the cells in both the normal and diseased conditions is the same, but a complex cluster of sub-networks controlled by hub genes drives the HCC progression with high metabolic flexibility and plasticity. In particular, we have extracted a sub-net of genes strictly correlated to other hub genes of the network from HepG2 cells, but specific for the HCC and mainly devoted to: (i) control at chromatin levels of cell division; (ii) control of ergastoplasmatic stress through protein degradation and misfolding; (iii) control of the immune response also through an increase of mature T-cells in the thymus. This sub-net is characterized by 26 hub genes coding for intrinsically disordered proteins with a high ability to interact with numerous molecular partners. Moreover, we have also noted that periphery molecules, that is, with one or very few interactions (e.g., cytokines or post-translational enzymes), which do not have a central role in the clusters that make up the global metabolic network, essentially have roles as information transporters. The results evidence a strong presence of intrinsically disordered proteins with key roles as hubs in the sub-networks that characterize the various stages of the disease, conferring a structural plasticity to the net nodes but an inherent functional versatility to the whole metabolic net. Thus, our present article provides a novel way of targeting the intrinsic disorder in HCC networks to dampen the cancer effects and provides new insight into the potential mechanisms of HCC. Taken together, the present findings suggest novel targets to design strategies for drug design and may support a rational intervention in the pharmacotherapy of HCC and other associated diseases.

Protein function in precision medicine: deep understanding with machine learning

Precision medicine and personalized health efforts propose leveraging complex molecular, medical and family history, along with other types of personal data toward better life. We argue that this ambitious objective will require advanced and specialized machine learning solutions. Simply skimming some low-hanging results off the data wealth might have limited potential. Instead, we need to better understand all parts of the system to define medically relevant causes and effects: how do particular sequence variants affect particular proteins and pathways? How do these effects, in turn, cause the health or disease-related phenotype? Toward this end, deeper understanding will not simply diffuse from deeper machine learning, but from more explicit focus on understanding protein function, context-specific protein interaction networks, and impact of variation on both.

Simulation of coupled folding and binding of an intrinsically disordered protein in explicit solvent with metadynamics

The C-terminal domain of measles virus nucleoprotein is an intrinsically disordered protein that could bind to the X domain (XD) of phosphoprotein P to exert its physiological function. Experiments reveal that the minimal binding unit is a 21-residue α-helical molecular recognition element (α-MoRE-MeV), which adopts a fully helical conformation upon binding to XD. Due to currently limited computing power, direct simulation of this coupled folding and binding process with atomic force field in explicit solvent cannot be achieved. In this work, two advanced sampling methods, metadynamics and parallel tempering, are combined to characterize the free energy surface of this process and investigate the underlying mechanism. Starting from an unbound and partially folded state of α-MoRE-MeV, multiple folding and binding events are observed during the simulation and the energy landscape was well estimated. The results demonstrate that the isolated α-MoRE-MeV resembles a molten globule and rapidly interconverts between random coil and multiple partially helical states in solution. The coupled folding and binding process occurs through the induced fit mechanism, with the residual helical conformations providing the initial binding sites. Upon binding, α-MoRE-MeV can easily fold into helical conformation without obvious energy barriers. Two mechanisms, namely, the system tending to adopt the structure in which the free energy of isolated α-MoRE-MeV is the minimum, and the binding energy of α-MoRE-MeV to its partner protein XD tending to the minimum, jointly dominate the coupled folding and binding process. With the advanced sampling approach, more IDP systems could be simulated and common mechanisms concerning the coupled folding and binding process could be investigated in the future.

The contribution of protein intrinsic disorder to understand the role of genetic variants uncovered by autism spectrum disorders exome studies

Several autism spectrum disorders (ASD) exome studies suggest that coding single nucleotide variants (SNVs) play an important role on ASD etiology. Usually, the pathogenic effect of missense mutations is estimated through predictors that lose accuracy for those SNVs placed in intrinsically disordered regions of protein. Here, we used bioinformatics tools to investigate the effect of mutations described in ASD published exome studies (549 mutations) in protein disorder, considering post-translational modification, PEST and Molecular Recognition Features (MoRFs) motifs. Schizophrenia and type 2 diabetes (T2D) datasets were created for comparison purposes. The frequency of mutations predicted as disordered was comparable among the three datasets (38.1% in ASD, 35.7% in schizophrenia, 46.4% in T2D). However, the frequency of SNVs predicted to lead a gain or loss of functional sites or change intrinsic disorder tendencies was higher in ASD and schizophrenia than T2D (46.9%, 36.4%, and 23.1%, respectively). The results obtained by SIFT and PolyPhen-2 indicated that 38.9% and 34.4% of the mutations predicted, respectively, as tolerated and benign showed functional alterations in disorder properties. Given the frequency of mutations placed in IDRs and their functional impact, this study suggests that alterations in intrinsic disorder properties might play a role in ASD and schizophrenia etiologies. They should be taken into consideration when researching the pathogenicity of mutations in neurodevelopmental and psychiatric diseases. Finally, mutations with functional alterations in disorder properties must be potential targets for in vitro and in vivo functional studies.

Long range recognition and selection in IDPs: the interactions of the C-terminus of p53

The C-terminal domain of p53 is an extensively studied IDP, interacting with different partners through multiple distinct conformations. To explore the interplay between preformed structural elements and intrinsic fluctuations in its folding and binding we combine extensive atomistic equilibrium and non-equilibrium simulations. We find that the free peptide segment rapidly interconverts between ordered and disordered states with significant populations of the conformations that are seen in the complexed states. The underlying global folding-binding landscape points to a synergistic mechanism in which recognition is dictated via long range electrostatic recognition which results in the formation of reactive structures as far away as 10 Å, and binding proceeds with the steering of selected conformations followed by induced folding at the target surface or within a close range.