A computational study on mitogenome-encoded proteins of Pavo cristatus and Pavo muticus identifies key genetic variations with functional implications.

BACKGROUND: The Pavo cristatus population, native to the Indian subcontinent, is thriving well in India. However, the Pavo muticus population, native to the tropical forests of Southeast Asia, has reduced drastically and has been categorised as an endangered group. To understand the probable genetic factors associated with the decline of P. muticus, we compared the mitogenome-encoded proteins (13 proteins) between these two species. RESULTS: Our data revealed that the most frequent variant between these two species was mtND1, which had an alteration in 9.57% residues, followed by mtND5 and mtATP6. We extended our study on the rest of the proteins and observed that cytochrome c oxidase subunits 1, 2, and 3 do not have any change. The 3-dimensional structure of all 13 proteins was modeled using the Phyre2 programme. Our data show that most of the proteins are alpha helical, and the variations observed in P. muticus reside on the surface of the respective proteins. The effect of variation on protein function was also predicted, and our results show that amino acid substitution in mtND1 at 14 sites could be deleterious. Similarly, destabilising changes were observed in mtND1, 2, 3, 4, 5, and 6 and mtATP6-8 due to amino acid substitution in P. muticus. Furthermore, protein disorder scores were considerably altered in mtND1, 2, and 5 of P. muticus. CONCLUSIONS: The results presented here strongly suggest that variations in mitogenome-encoded proteins of P. cristatus and P. muticus may alter their structure and functions. Subsequently, these variations could alter energy production and may correlate with the decline in the population of P. muticus.

Critical assessment of protein intrinsic disorder prediction (CAID) - Results of round 2.

Protein intrinsic disorder (ID) is a complex and context-dependent phenomenon that covers a continuum between fully disordered states and folded states with long dynamic regions. The lack of a ground truth that fits all ID flavors and the potential for order-to-disorder transitions depending on specific conditions makes ID prediction challenging. The CAID2 challenge aimed to evaluate the performance of different prediction methods across different benchmarks, leveraging the annotation provided by the DisProt database, which stores the coordinates of ID regions when there is experimental evidence in the literature. The CAID2 challenge demonstrated varying performance of different prediction methods across different benchmarks, highlighting the need for continued development of more versatile and efficient prediction software. Depending on the application, researchers may need to balance performance with execution time when selecting a predictor. Methods based on AlphaFold2 seem to be good ID predictors but they are better at detecting absence of order rather than ID regions as defined in DisProt. The CAID2 predictors can be freely used through the CAID Prediction Portal, and CAID has been integrated into OpenEBench, which will become the official platform for running future CAID challenges.

Possible functional proximity of various organisms based on the bioinformatics analysis of their taste receptors.

Taste is one of the essential senses in providing the organism a faithful representation of the external world. Taste perception is responsible for basic food and drink appraisal and bestows the organism with valuable discriminatory power. Umami and sweet are "good" tastes that promote consumption of nutritive food, whereas bitter and sour are "bad" tastes that alert the organism to toxins and low pH, promoting rejection of foods containing harmful substances. Not every animal has the same sense of taste as humans. Variation in the taste receptor genes contributes to inter and intra organism differences of taste (sweet/bitter) sensation and preferences. Therefore a deeper understanding was needed to comprehend taste perception by various vertebrates and accordingly elucidate a possible proximity among them. In this study, a total 20 Type-1 (sweet) and 189 Type-2 (bitter) taste receptor complete-amino acid sequences were taken from the 20 vertebrate organisms (18 mammalian, 1 Aves, and 1 amphibian). Among 10 primates, 8 including humans were very close based on genomics of taste receptors and rodent organisms viz. the rat and mouse were away from them. This investigation throws light on the similitude and dissimilitude of perception of sweet and bitter taste among 20 different organisms, steered by quantitative analysis of their genomic data. Furthermore, it enlightened that ligand binding affinity of sweet/bitter taste molecules in the taste receptors of any proximal pair of organisms would be similar.

Deep learning in prediction of intrinsic disorder in proteins.

Intrinsic disorder prediction is an active area that has developed over 100 predictors. We identify and investigate a recent trend towards the development of deep neural network (DNN)-based methods. The first DNN-based method was released in 2013 and since 2019 deep learners account for majority of the new disorder predictors. We find that the 13 currently available DNN-based predictors are diverse in their topologies, sizes of their networks and the inputs that they utilize. We empirically show that the deep learners are statistically more accurate than other types of disorder predictors using the blind test dataset from the recent community assessment of intrinsic disorder predictions (CAID). We also identify several well-rounded DNN-based predictors that are accurate, fast and/or conveniently available. The popularity, favorable predictive performance and architectural flexibility suggest that deep networks are likely to fuel the development of future disordered predictors. Novel hybrid designs of deep networks could be used to adequately accommodate for diversity of types and flavors of intrinsic disorder. We also discuss scarcity of the DNN-based methods for the prediction of disordered binding regions and the need to develop more accurate methods for this prediction.

Emergence of unique SARS-CoV-2 ORF10 variants and their impact on protein structure and function.

The devastating impact of the ongoing coronavirus disease 2019 (COVID-19) on public health, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has made targeting the COVID-19 pandemic a top priority in medical research and pharmaceutical development. Surveillance of SARS-CoV-2 mutations is essential for the comprehension of SARS-CoV-2 variant diversity and their impact on virulence and pathogenicity. The SARS-CoV-2 open reading frame 10 (ORF10) protein interacts with multiple human proteins CUL2, ELOB, ELOC, MAP7D1, PPT1, RBX1, THTPA, TIMM8B, and ZYG11B expressed in lung tissue. Mutations and co-occurring mutations in the emerging SARS-CoV-2 ORF10 variants are expected to impact the severity of the virus and its associated consequences. In this article, we highlight 128 single mutations and 35 co-occurring mutations in the unique SARS-CoV-2 ORF10 variants. The possible predicted effects of these mutations and co-occurring mutations on the secondary structure of ORF10 variants and host protein interactomes are presented. The findings highlight the possible effects of mutations and co-occurring mutations on the emerging 140 ORF10 unique variants from secondary structure and intrinsic protein disorder perspectives.

Protein feature analysis of heat shock induced ubiquitination sites reveals preferential modification site localization.

Protein aggregation is indicative of failing protein quality control systems. These systems are responsible for the refolding or degradation of aberrant and misfolded proteins. Heat stress can cause proteins to misfold, triggering cellular responses including a marked increase in the ubiquitination of proteins. This response has been characterized in yeast, however more studies are needed within mammalian cells. Herein, we examine proteins that become ubiquitinated during heat shock in human tissue culture cells using diGly enrichment coupled with mass spectrometry. A majority of these proteins are localized in the nucleus or cytosol. Proteins which are conjugated under stress display longer sequence lengths, more interaction partners, and more hydrophobic patches than controls but do not show lower melting temperatures. Furthermore, heat-induced conjugation sites occur less frequently in disordered regions and are closer to hydrophobic patches than other ubiquitination sites; perhaps providing novel insight into the molecular mechanism mediating this response. Nuclear and cytosolic pools of modified proteins appear to have different protein features. Using a pulse-SILAC approach, we found that both long-lived and newly-synthesized proteins are conjugated under stress. Modified long-lived proteins are predominately nuclear and were distinct from newly-synthesized proteins, indicating that different pathways may mediate the heat-induced increase of polyubiquitination. SIGNIFICANCE: The maintenance of protein homeostasis requires a balance of protein synthesis, folding, and degradation. Under stress conditions, the cell must rapidly adapt by increasing its folding capacity to eliminate aberrant proteins. A major pathway for proteolysis is mediated by the ubiquitin proteasome system. While increased ubiquitination after heat stress was observed over 30 years ago, it remains unclear which proteins are conjugated during heat shock in mammalian cells and by what means this conjugation occurs. In this study, we combined SILAC-based mass spectrometry with computational analyses to reveal features associated to proteins ubiquitinated while under heat shock. Interestingly, we found that conjugation sites induced by the stress are less often located within disordered regions and more often located near hydrophobic patches. Our study showcases how proteomics can reveal distinct feature associated to a cohort of proteins that are modified post translationally and how the ubiquitin conjugation sites are preferably selected in these conditions. Our work opens a new path for delineating the molecular mechanisms leading to the heat stress response and the regulation of protein homeostasis.

Mapping the nucleolar proteome reveals a spatiotemporal organization related to intrinsic protein disorder.

The nucleolus is essential for ribosome biogenesis and is involved in many other cellular functions. We performed a systematic spatiotemporal dissection of the human nucleolar proteome using confocal microscopy. In total, 1,318 nucleolar proteins were identified; 287 were localized to fibrillar components, and 157 were enriched along the nucleoplasmic border, indicating a potential fourth nucleolar subcompartment: the nucleoli rim. We found 65 nucleolar proteins (36 uncharacterized) to relocate to the chromosomal periphery during mitosis. Interestingly, we observed temporal partitioning into two recruitment phenotypes: early (prometaphase) and late (after metaphase), suggesting phase-specific functions. We further show that the expression of MKI67 is critical for this temporal partitioning. We provide the first proteome-wide analysis of intrinsic protein disorder for the human nucleolus and show that nucleolar proteins in general, and mitotic chromosome proteins in particular, have significantly higher intrinsic disorder level compared to cytosolic proteins. In summary, this study provides a comprehensive and essential resource of spatiotemporal expression data for the nucleolar proteome as part of the Human Protein Atlas.

Distinct Features of Stress Granule Proteins Predict Localization in Membraneless Organelles.

Recently generated proteomic data provides unprecedented insight into stress granule composition and stands as fruitful ground for further analysis. Stress granules are stress-induced biological assemblies that are of keen interest due to being linked to both long-term cell viability and a variety of protein aggregation-based diseases. Herein, we compile recently published stress granule composition data, formed specifically through heat and oxidative stress, for both mammalian (Homo sapiens) and yeast (Saccharomyces cerevisiae) cells. Interrogation of the data reveals that stress granule proteins are enriched in features that favor protein liquid-liquid phase separation, being highly disordered, soluble, and abundant while maintaining a high level of protein-protein interactions under basal conditions. Furthermore, these "stress granuleomes" are shown to be enriched for multidomained, RNA-binding proteins with increased potential for post-translational modifications. Findings are consistent with the notion that stress granule formation is driven by protein liquid-liquid phase separation. Furthermore, stress granule proteins appear poised near solubility limits while possessing the ability to dynamically alter their phase behavior in response to external threat. Interestingly, several features, such as protein disorder, are more prominent among stress granule proteins that share homologs between yeast and mammalian systems also found within stress-induced foci. We culminate results from our stress granule analysis into novel predictors for granule incorporation and validate the mammalian predictor's performance against multiple types of membraneless condensates and by colocalization microscopy.

Challenges in the Structural-Functional Characterization of Multidomain, Partially Disordered Proteins CBP and p300: Preparing Native Proteins and Developing Nanobody Tools.

The structural and functional characterization of large multidomain signaling proteins containing long disordered linker regions represents special methodological and conceptual challenges. These proteins show extreme structural heterogeneity and have complex posttranslational modification patterns, due to which traditional structural biology techniques provide results that are often difficult to interpret. As demonstrated through the example of two such multidomain proteins, CREB-binding protein (CBP) and its paralogue, p300, even the expression and purification of such proteins are compromised by their extreme proteolytic sensitivity and structural heterogeneity. In this chapter, we describe the effective expression of CBP and p300 in a eukaryotic host, Sf9 insect cells, followed by their tandem affinity purification based on two terminal tags to ensure their structural integrity. The major focus of this chapter is on the development of novel accessory tools, single-domain camelid antibodies (nanobodies), for structural-functional characterization. Specific nanobodies against full-length CBP and p300 can specifically target their different regions and can be used for their marking, labeling, and structural stabilization in a broad range of in vitro and in vivo studies. Here, we describe four high-affinity nanobodies binding to the KIX and the HAT domains, either mimicking known interacting partners or revealing new functionally relevant conformations. As immunization of llamas results in nanobody libraries with a great sequence variation, deep sequencing and interaction analysis with different regions of the proteins provide a novel approach toward developing a panel of specific nanobodies.

A subset of functional adaptation mutations alter propensity for α-helical conformation in the intrinsically disordered glucocorticoid receptor tau1core activation domain.

BACKGROUND: Adaptive mutations that alter protein functionality are enriched within intrinsically disordered protein regions (IDRs), thus conformational flexibility correlates with evolvability. Pre-structured motifs (PreSMos) with transient propensity for secondary structure conformation are believed to be important for IDR function. The glucocorticoid receptor tau1core transcriptional activation domain (GR tau1core) domain contains three α-helical PreSMos in physiological buffer conditions. METHODS: Sixty change-of-function mutants affecting the intrinsically disordered 58-residue GR tau1core were studied using disorder prediction and molecular dynamics simulations. RESULTS: Change-of-function mutations were partitioned into seven clusters based on their effect on IDR predictions and gene activation activity. Some mutations selected from clusters characterized by mutations altering the IDR prediction score, altered the apparent stability of the α-helical form of one of the PreSMos in molecular dynamics simulations, suggesting PreSMo stabilization or destabilization as strategies for functional adaptation. Indeed all tested gain-of-function mutations affecting this PreSMo were associated with increased stability of the α-helical PreSMo conformation, suggesting that PreSMo stabilization may be the main mechanism by which adaptive mutations can increase the activity of this IDR type. Some mutations did not appear to affect PreSMo stability. CONCLUSIONS: Changes in PreSMo stability account for the effects of a subset of change-of-function mutants affecting the GR tau1core IDR. GENERAL SIGNIFICANCE: Long IDRs occur in about 50% of human proteins. They are poorly characterized despite much recent attention. Our results suggest the importance of a subtle balance between PreSMo stability and IDR activity, which may provide a novel target for future pharmaceutical intervention.

Nuclear and nucleolar activity of linker histone variant H1.0.

Histone H1.0 belongs to the class of linker histones (H1), although it is substantially distinct from other histone H1 family members. The differences can be observed in the chromosomal location and organization of the histone H1.0 encoding gene, as well as in the length and composition of its amino acid chain. Whereas somatic (H1.1-H1.5) histone H1 variants are synthesized in the cell cycle S-phase, histone H1.0 is synthesized throughout the cell cycle. By replacing somatic H1 variants during cell maturation, histone H1.0 is gradually deposited in low dividing cells and achieves the highest level of expression in the terminally differentiated cells. Compared to other differentiation-specific H1 histone (H5) characteristic for unique tissue and organisms, the distribution of histone H1.0 remains non-specific. Classic investigations emphasize that histone H1.0 is engaged in the organization of nuclear chromatin accounting for formation and maintenance of its nucleosomal and higher-order structure, and thus influences gene expression. However, the recent data confirmed histone H1.0 peculiar localization in the nucleolus and unexpectedly revealed its potential for regulation of nucleolar, RNA-dependent, activity via interaction with other proteins. According to such findings, histone H1.0 participates in the formation of gene-coded information through its control at both transcriptional and translational levels. In order to reappraise the biological significance of histone H1.0, both aspects of its activity are presented in this review.

Evolutionary volatile Cysteines and protein disorder in the fast evolving tunicate Oikopleura dioica.

Cysteine (Cys) is regarded as the most conservative amino acid in nature, something that does not occur in the tunicate Oikopleura dioica, where this amino acid is one of the fastest evolving. In this work we analyze some of the causes of this intriguing absence of conservation. Considering the well-known stabilizing role of Cys, it was first investigated whether the lack of conservation was accompanied by an increase in intrinsic protein disorder. In contrast to expectations, it was found that O. dioica is the chordate that has the lowest levels of intrinsic disorder, while vertebrates (represented by Bos taurus) contain the most disordered proteins. Oikopleura proteins are shorter than their homologs in other Chordates (Ciona and B. taurus proteins are respectively 11% and 18% longer). This process of protein shortening was more intense in intrinsic disordered regions. As a result proteins became not only shorter but also more compact. It is also reported here that the conservation/divergence behavior of Cys depends on whether they are located in ordered or disordered regions. In the four species analyzed, disordered Cys are majorly (> 75%) not conserved at all. Ordered Cys instead, are much more free to diverge in Oikopleura than in the other chordates. We hypothesize that the preferential deletion of disordered regions resulted in a decreased protein disorder and a direct elimination (by deletion) of many ancestral Cys. Besides, the alterations (shortening or complete elimination) of some disordered regions (loops/random coils) probably promoted further Cys evolutionary volatility, because some ancestral Cys (and other amino acids which play a role in stability like Trp) located outside deleted regions became redundant due to the loss of their stabilizing partners.

The high-risk HPV16 E7 oncoprotein mediates interaction between the transcriptional coactivator CBP and the retinoblastoma protein pRb.

The oncoprotein E7 from human papillomavirus (HPV) strains that confer high cancer risk mediates cell transformation by deregulating host cellular processes and activating viral gene expression through recruitment of cellular proteins such as the retinoblastoma protein (pRb) and the cyclic-AMP response element binding binding protein (CBP) and its paralog p300. Here we show that the intrinsically disordered N-terminal region of E7 from high-risk HPV16 binds the TAZ2 domain of CBP with greater affinity than E7 from low-risk HPV6b. HPV E7 and the tumor suppressor p53 compete for binding to TAZ2. The TAZ2 binding site in E7 overlaps the LxCxE motif that is crucial for interaction with pRb. While TAZ2 and pRb compete for binding to a monomeric E7 polypeptide, the full-length E7 dimer mediates an interaction between TAZ2 and pRb by promoting formation of a ternary complex. Cell-based assays show that expression of full-length HPV16 E7 promotes increased pRb acetylation and that this response depends both on the presence of CBP/p300 and on the ability of E7 to form a dimer. These observations suggest a model for the oncogenic effect of high-risk HPV16 E7. The disordered region of one E7 molecule in the homodimer interacts with the pocket domain of pRb, while the same region of the other E7 molecule binds the TAZ2 domain of CBP/p300. Through its ability to dimerize, E7 recruits CBP/p300 and pRb into a ternary complex, bringing the histone acetyltransferase domain of CBP/p300 into proximity to pRb and promoting acetylation, leading to disruption of cell cycle control.

An overview of the sequence features of N- and C-terminal segments of the human chemokine receptors.

Chemokine receptors play a crucial role in the cellular signaling enrolling extracellular ligands chemotactic proteins which recruit immune cells. They possess seven trans-membrane helices, an extracellular N-terminal region with three extracellular hydrophilic loops being important for search and recognition of specific ligand(s), and an intracellular C-terminal region with three intracellular loops that couple G-proteins. Although the functional aspects of the terminal segments of the extra-and intra-cellular G proteins are universally identified, the molecular basis on which they rest are still unclear because they are not definable by means of X-rays due to their high mobility and are not easy to study in the membrane. The purpose of this work is to define which physical-chemical properties of the terminal segments of the human chemokine receptors are at the basis of their functional mechanisms. Therefore, we have evaluated their physical-chemical properties in terms of amino acid composition, local flexibility, disorder propensity, net charge distribution and putative sites of post-translational modifications. Our results support the conclusion that all 19 C-terminal and N-terminal segments of human chemokine receptors are very flexible due to the systematic presence of intrinsic disorder. Although, the purpose of this plasticity clearly appears that of controlling and modulating the binding of ligands, we provide evidence that the overlap of linearly charged stretches, intrinsic disorder and post-translational modification sites, consistently found in these motives, is a necessary feature to exert the function. The role of the intrinsic disorder has been discussed considering the structural information coming from intrinsically disordered model compounds which support the view that the chemokine terminals have to be considered as strong polyampholytes or polyelectrolytes where conformational ensembles and structural transitions between them are modulated by charge fraction variations. Also the role of post-translational modifications has been found coherent with this view because, changing the charge fraction, they guide structural transitions between ensembles. Moreover, we have also considered our results from an evolutionary point of view in order to understand if the features found in humans were also present in other species. Our data evidenced that the structural features of the human terminals of the chemokine receptors were shared and evolutionarily conserved particularly among mammals. This means that the various organisms not only tolerate but select intrinsic disorder for the terminal regions of their receptors, reflecting constraints that point to molecular recognition. In conclusion the terminal segments of chemokine receptors must be considered as strong polyampholytes where the charge fraction variations induced by post-translational modifications are the driving physico-chemical feature able to adapt the conformations of the terminal segments to their functions.

Protein domain definition should allow for conditional disorder.

Proteins are often classified in a binary fashion as either structured or disordered. However this approach has several deficits. Firstly, protein folding is always conditional on the physiochemical environment. A protein which is structured in some circumstances will be disordered in others. Secondly, it hides a fundamental asymmetry in behavior. While all structured proteins can be unfolded through a change in environment, not all disordered proteins have the capacity for folding. Failure to accommodate these complexities confuses the definition of both protein structural domains and intrinsically disordered regions. We illustrate these points with an experimental study of a family of small binding domains, drawn from the RNA polymerase of mumps virus and its closest relatives. Assessed at face value the domains fall on a structural continuum, with folded, partially folded, and near unstructured members. Yet the disorder present in the family is conditional, and these closely related polypeptides can access the same folded state under appropriate conditions. Any heuristic definition of the protein domain emphasizing conformational stability divides this domain family in two, in a way that makes no biological sense. Structural domains would be better defined by their ability to adopt a specific tertiary structure: a structure that may or may not be realized, dependent on the circumstances. This explicitly allows for the conditional nature of protein folding, and more clearly demarcates structural domains from intrinsically disordered regions that may function without folding.

Hydrogel formation by multivalent IDPs: A reincarnation of the microtrabecular lattice?

Based on high-voltage electron microscopic (HVEM) data of fixed cultured cells, an elaborate three-dimensional network of filaments, including and interconnecting other elements of the cytoskeleton, was observed in cells some half a century ago. Despite many attempts and comparative studies, this "microtrabecular lattice" (MTL) of the cytoplasmic ground substance could not be established as a genuine component of the eukaryotic cell, and is mostly considered today as a sample-preparation artifact of protein adherence and cross-linking to the cytoskeleton. Here we elaborate on the provocative idea that recent observations of hydrogel-forming phase transitions of repetitive regions of intrinsically disordered proteins (IDPs) bear resemblance in creation, organization and physical appearance to the MTL. We review this phenomenon in detail, and suggest that phase transitions of actin regulatory proteins, neurofilament side-arms and other proteins could generate non-uniform spatial distribution of cytoplasmic material in the vicinity of the cytoskeleton that might even give rise to fixation phenomena resembling the MTL. Whether such hydrogel formation by IDPs is a general physical phenomenon, will remain to be seen, nevertheless, the underlying organizational principle provokes novel experimental studies to uncover the ensuing higher-level regulation of cell physiology, in which the despised and long-forgotten concept of MTL might give some interesting leads.