The alphabet of intrinsic disorder: II. Various roles of glutamic acid in ordered and intrinsically disordered proteins

Engineering salt-tolerant Cas12f1 variants for gene-editing applications

CRISPR has revolutionized the field of genome editing in life sciences by serving as a versatile and state-of-the-art tool. Cas12f1 is a small nuclease of the bacterial immunity CRISPR system with an ideal size for cellular delivery, in contrast to CRISPR-associated (Cas) proteins like Cas9 or Cas12. However, Cas12f1 works best at low salt concentrations. In this study, we find that the plasticity of certain Cas12f1 regions (K196-Y202 and I452-L515) is negatively affected by increased salt concentrations. On this line, key protein domains (REC1, WED, Nuc, lid) that are involved in the DNA-target recognition and the activation of the catalytic RuvC domain are in turn also affected. We suggest that salt concentration should be taken in to consideration for activity assessments of Cas engineered variants, especially if the mutations are on the protospacer adjacent motif interacting domain. The results can be exploited for the engineering of Cas variants and the assessment of their activity at varying salt concentrations. We propose that the K198Q mutation can restore at great degree the compromised plasticity and could potentially lead to salt-tolerant Cas12f1 variants. The methodology can be also employed for the study of biomolecules in terms of their salinity tolerance.Communicated by Ramaswamy H. Sarma.

Digital nets conformational sampling (DNCS) - an enhanced sampling technique to explore the conformational space of intrinsically disordered peptides

We propose digital nets conformational sampling (DNCS) - an enhanced sampling technique to explore the conformational ensembles of peptides, especially intrinsically disordered peptides (IDPs). The DNCS algorithm relies on generating history-dependent samples of dihedral variables using bitwise XOR operations and binary angle measurements (BAM). The algorithm was initially studied using met-enkephalin, a highly elusive neuropeptide. The DNCS method predicted near-native structures and the energy landscape of met-enkephalin was observed to be in direct correlation with earlier studies on the neuropeptide. Clustering analysis revealed that there are only 24 low-lying conformations of the molecule. The DNCS method has then been tested for predicting optimal conformations of 42 oligopeptides of length varying from 3 to 8 residues. The closest-to-native structures of 86% of cases are near-native and 24% of them have a root mean square deviation of less than 1.00 Å with respect to their crystal structures. The results obtained reveal that the DNCS method performs well, that too in less computational time.

Understanding the Specific Implications of Amino Acids in the Antibody Development

As the demand for immunotherapy to treat and manage cancers, infectious diseases and other disorders grows, a comprehensive understanding of amino acids and their intricate role in antibody engineering has become a prime requirement. Naturally produced antibodies may not have the most suitable amino acids at the complementarity determining regions (CDR) and framework regions, for therapeutic purposes. Therefore, to enhance the binding affinity and therapeutic properties of an antibody, the specific impact of certain amino acids on the antibody's architecture must be thoroughly studied. In antibody engineering, it is crucial to identify the key amino acid residues that significantly contribute to improving antibody properties. Therapeutic antibodies with higher binding affinity and improved functionality can be achieved through modifications or substitutions with highly suitable amino acid residues. Here, we have indicated the frequency of amino acids and their association with the binding free energy in CDRs. The review also analyzes the experimental outcome of two studies that reveal the frequency of amino acids in CDRs and provides their significant correlation between the outcomes. Additionally, it discusses the various bond interactions within the antibody structure and antigen binding. A detailed understanding of these amino acid properties should assist in the analysis of antibody sequences and structures needed for designing and enhancing the overall performance of therapeutic antibodies.

Immune targeting of filarial glutaredoxin through a multi-epitope peptide-based vaccine: A reverse vaccinology approach

Despite the efforts of global programme to eliminate lymphatic filariasis (GPELF), the threat of lymphatic filariasis (LF) still looms over humanity in terms of long-term disabilities, and morbidities across the globe. In light of this situation, investigators have chosen to focus on the development of immunotherapeutics targeting the physiologically important filarial-specific proteins. Glutaredoxin (16.43 kDa) plays a pivotal role in filarial redox biology, serving as a vital contributor. In the context of the intra-host survival of filarial parasites, this antioxidant helps in mitigating the oxidative stress imposed by the host immune system. Given its significant contribution, the development of a vaccine targeting glutaredoxin holds promise as a new avenue for achieving a filaria-free world. Herein, multi-epitope-based vaccine was designed using advanced immunoinformatics approach. Initially, 4B-cell epitopes and 6 T-cell epitopes (4 MHC I and 2 MHC II) were identified from the 146 amino acid long sequence of glutaredoxin of the human filarid, Wuchereria bancrofti. Subsequent clustering of these epitopes with linker peptides finalized the vaccine structure. To boost TLR-mediated innate immunity, TLR-specific adjuvants were incorporated into the designed vaccine. After that, experimental analyses confirm the designed vaccine, Vac4 as anefficient ligand of human TLR5 to elicit protective innate immunity against filarial glutaredoxin. Immune simulation further demonstrated abundant levels of IgG and IgM as crucial contributors in triggering vaccine-induced adaptive responses in the recipients. Hence, to facilitate the validation of immunogenicity of the designed vaccine, Vac4 was cloned in silico in pET28a(+) expression vector for recombinant production. Taken together, our findings suggest that vaccine-mediated targeting of filarial glutaredoxin could be a future option for intervening LF on a global scale.

High-throughput Selection of Human de novo-emerged sORFs with High Folding Potential

De novo genes emerge from previously noncoding stretches of the genome. Their encoded de novo proteins are generally expected to be similar to random sequences and, accordingly, with no stable tertiary fold and high predicted disorder. However, structural properties of de novo proteins and whether they differ during the stages of emergence and fixation have not been studied in depth and rely heavily on predictions. Here we generated a library of short human putative de novo proteins of varying lengths and ages and sorted the candidates according to their structural compactness and disorder propensity. Using Förster resonance energy transfer combined with Fluorescence-activated cell sorting, we were able to screen the library for most compact protein structures, as well as most elongated and flexible structures. We find that compact de novo proteins are on average slightly shorter and contain lower predicted disorder than less compact ones. The predicted structures for most and least compact de novo proteins correspond to expectations in that they contain more secondary structure content or higher disorder content, respectively. Our experiments indicate that older de novo proteins have higher compactness and structural propensity compared with young ones. We discuss possible evolutionary scenarios and their implications underlying the age-dependencies of compactness and structural content of putative de novo proteins.

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

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

Thermodynamic characterization of amyloid polymorphism by microfluidic transient incomplete separation

Amyloid fibrils of proteins such as α-synuclein are a hallmark of neurodegenerative diseases and much research has focused on their kinetics and mechanisms of formation. The question as to the thermodynamic stability of such structures has received much less attention. Here, we newly utilize the principle of transient incomplete separation of species in laminar flow in combination with chemical depolymerization for the quantification of amyloid fibril stability. The relative concentrations of fibrils and monomer at equilibrium are determined through an in situ separation of these species based on their different diffusivity inside a microfluidic capillary. The method is highly sample economical, using much less than a microliter of sample per data point and its only requirement is the presence of aromatic residues (W, Y) because of its label-free nature, which makes it widely applicable. Using this method, we investigate the differences in thermodynamic stability between different fibril polymorphs of α-synuclein and quantify these differences for the first time. Importantly, we show that fibril formation can be under kinetic or thermodynamic control and that a change in solution conditions can both stabilise and destabilise amyloid fibrils. Taken together, our results establish the thermodynamic stability as a well-defined and key parameter that can contribute towards a better understanding of the physiological roles of amyloid fibril polymorphism.

A three-state mechanism for trifluoroethanol denaturation of an intrinsically disordered protein (IDP)

Relating the amino acid composition and sequence to chain folding and binding preferences of intrinsically disordered proteins (IDPs) has emerged as a huge challenge. While globular proteins have respective 3D structures that are unique to their individual functions, IDPs violate this structure-function paradigm because rather than having a well-defined structure an ensemble of rapidly interconverting disordered structures characterize an IDP. This work measures 2,2,2-trifluoroethanol (TFE)-induced equilibrium transitions of an IDP called AtPP16-1 (Arabidopsis thaliana phloem protein type 16-1) by using fluorescence, circular dichroism, infrared and nuclear magnetic resonance (NMR) methods at pH 4, 298 K. Low TFE reversibly removes the tertiary structure to produce an ensemble of obligate intermediate ($\mathrm{I}$) retaining the native-state ($\mathrm{N}$) secondary structure. The intermediate $\mathrm{I}$ is preceded by a non-obligate tryptophan-specific intermediate ${\mathrm{I}}_{\mathrm{w}}$ whose population is detectable for AtPP16-1 specifically. Accumulation of such non-obligate intermediates is discriminated according to the sequence composition of the protein. In all cases, however, a tertiary structure-unfolded general obligate intermediate $\mathrm{I}$ is indispensable. The $\mathrm{I}$ ensemble has higher helical propensity conducive to the acquisition of an exceedingly large level of α-helices by a reversible denaturation transition of $\mathrm{I}$ to the denatured state $\mathrm{D}$ as the TFE level is increased. Strikingly, it is the same $\mathrm{N}\rightleftharpoons \mathrm{I}\rightleftharpoons \mathrm{D}$ scheme typifying the TFE transitions of globular proteins. The high-energy state $\mathrm{I}$ characterized by increased helical propensity is called a universal intermediate encountered in both genera of globular and disordered proteins. Neither $\mathrm{I}$ nor $\mathrm{D}$ strictly show molten globule (MG)-like properties, dismissing the belief that TFE promotes MGs.

Expanding the repertoire of human tandem repeat RNA-binding proteins

Protein regions consisting of arrays of tandem repeats are known to bind other molecular partners, including nucleic acid molecules. Although the interactions between repeat proteins and DNA are already widely explored, studies characterising tandem repeat RNA-binding proteins are lacking. We performed a large-scale analysis of human proteins devoted to expanding the knowledge about tandem repeat proteins experimentally reported as RNA-binding molecules. This work is timely because of the release of a full set of accurate structural models for the human proteome amenable to repeat detection using structural methods. The main goal of our analysis was to build a comprehensive set of human RNA-binding proteins that contain repeats at the sequence or structure level. Our results showed that the combination of sequence and structural methods finds significantly more tandem repeat proteins than either method alone. We identified 219 tandem repeat proteins that bind RNA molecules and characterised the overlap between repeat regions and RNA-binding regions as a first step towards assessing their functional relationship. We observed differences in the characteristics of repeat regions predicted by sequence-based or structure-based methods in terms of their sequence composition, their functions and their protein domains.

Proteins released into the plant apoplast by the obligate parasitic protist Albugo selectively repress phyllosphere-associated bacteria

Biotic and abiotic interactions shape natural microbial communities. The mechanisms behind microbe-microbe interactions, particularly those protein based, are not well understood. We hypothesize that released proteins with antimicrobial activity are a powerful and highly specific toolset to shape and defend plant niches. We have studied Albugo candida, an obligate plant parasite from the protist Oomycota phylum, for its potential to modulate the growth of bacteria through release of antimicrobial proteins into the apoplast. Amplicon sequencing and network analysis of Albugo-infected and uninfected wild Arabidopsis thaliana samples revealed an abundance of negative correlations between Albugo and other phyllosphere microbes. Analysis of the apoplastic proteome of Albugo-colonized leaves combined with machine learning predictors enabled the selection of antimicrobial candidates for heterologous expression and study of their inhibitory function. We found for three candidate proteins selective antimicrobial activity against Gram-positive bacteria isolated from A. thaliana and demonstrate that these inhibited bacteria are precisely important for the stability of the community structure. We could ascribe the antibacterial activity of the candidates to intrinsically disordered regions and positively correlate it with their net charge. This is the first report of protist proteins with antimicrobial activity under apoplastic conditions that therefore are potential biocontrol tools for targeted manipulation of the microbiome.

Substrate and phosphorylation site selection by phosphoprotein phosphatases

Dynamic protein phosphorylation and dephosphorylation are essential regulatory mechanisms that ensure proper cellular signaling and biological functions. Deregulation of either reaction has been implicated in several human diseases. Here, we focus on the mechanisms that govern the specificity of the dephosphorylation reaction. Most cellular serine/threonine dephosphorylation is catalyzed by 13 highly conserved phosphoprotein phosphatase (PPP) catalytic subunits, which form hundreds of holoenzymes by binding to regulatory and scaffolding subunits. PPP holoenzymes recognize phosphorylation site consensus motifs and interact with short linear motifs (SLiMs) or structural elements distal to the phosphorylation site. We review recent advances in understanding the mechanisms of PPP site-specific dephosphorylation preference and substrate recruitment and highlight examples of their interplay in the regulation of cell division.

Overview Update: Computational Prediction of Intrinsic Disorder in Proteins

There are over 100 computational predictors of intrinsic disorder. These methods predict amino acid-level propensities for disorder directly from protein sequences. The propensities can be used to annotate putative disordered residues and regions. This unit provides a practical and holistic introduction to the sequence-based intrinsic disorder prediction. We define intrinsic disorder, explain the format of computational prediction of disorder, and identify and describe several accurate predictors. We also introduce recently released databases of intrinsic disorder predictions and use an illustrative example to provide insights into how predictions should be interpreted and combined. Lastly, we summarize key experimental methods that can be used to validate computational predictions. © 2023 Wiley Periodicals LLC.

IDP Force Fields Applied to Model PPII-Rich 33-mer Gliadin Peptides

The 33-mer gliadin peptide and its deamidated metabolite, 33-mer DGP, are the immunodominant peptides responsible for the adaptive immune response in celiac disease (CD). CD is a complex autoimmune chronic disorder triggered by gluten ingestion that affects the small intestine and affects ∼1% of the global population. The 33-mers are polyproline II-rich (PPII) and intrinsically disordered peptides (IDPs), whose structures remain elusive. We sampled the conformational ensembles of both 33-mer peptides via molecular dynamics simulations employing two force fields (FFs) (Amber ff03ws and Amber ff99SB-disp) specifically validated for other IDPs. Our results show that both FFs allow the extensive exploration of the conformational landscape, which was not possible with the standard FF GROMOS53A6 reported before. Clustering analysis of the trajectories showed that the five largest clusters (78-88% of the total structures) present elongated, semielongated, and curved conformations in both FFs. Large average radius of gyration and solvent-exposed surfaces characterized these structures. While the structures sampled are similar, the Amber ff99SB-disp trajectories explored folded conformations with a higher probability. In addition, PPII secondary structure was preserved throughout the trajectories (58-73%) together with a non-negligible content of β structures (11-23%), in agreement with previous experimental results. This work represents the initial step in studying further the interaction of these peptides with other biologically relevant molecules, which could lead to finally disclose the molecular events that lead to CD.

Probing Peptidylprolyl Bond cis/trans Status Using Distal 19 F NMR Reporters

A method for measuring peptidylprolyl bond cis-trans conformational status in peptide models is described, using 4-fluorophenylalanine (4FPhe) as a distal reporter for 19 F NMR. The %cis-Pro population was measured for peptides of the general structure Ac-X-Pro-Z-Ala-Ala-4FPhe (X and Z are proteinogenic amino acids) at pH 7.4, and provided conformational populations consistent with literature values obtained by more complex methods. This approach was applied to probe the prolyl bond status in pentapeptide models of the intrinsically disordered C-terminal region of α-synuclein, which mirrored the preferences in the Ac-X-Pro-Z-Ala-4FPhe models. Advantageously, the 19 F reporter group does not need to be adjacent to or attached to proline to provide quantifiable signals and distal 4-fluorophenylalanines can be placed so as not to influence prolyl bond conformation. Finally, we demonstrated that the prolyl bond status is not significantly affected by pH when there are ionisable amino acid residues at the carboxyl side of proline, which makes 19 F NMR an invaluable tool with which to study proline isomerism at a range of pHs and in different solvents and buffers.

Prediction of protein-protein interactions using sequences of intrinsically disordered regions

Protein-protein interactions (PPIs) play a crucial role in numerous molecular processes. Despite many efforts, mechanisms governing molecular recognition between interacting proteins remain poorly understood and it is particularly challenging to predict from sequence whether two proteins can interact. Here we present a new method to tackle this challenge using intrinsically disordered regions (IDRs). IDRs are protein segments that are functional despite lacking a single invariant three-dimensional structure. The prevalence of IDRs in eukaryotic proteins suggests that IDRs are critical for interactions. To test this hypothesis, we predicted PPIs using IDR sequences in candidate proteins in humans. Moreover, we divide the PPI prediction problem into two specific subproblems and adapt appropriate training and test strategies based on problem type. Our findings underline the importance of defining clearly the problem type and show that sequences encoding IDRs can aid in predicting specific features of the protein interaction network of intrinsically disordered proteins. Our findings further suggest that accounting for IDRs in future analyses should accelerate efforts to elucidate the eukaryotic PPI network.

Reciprocal regulatory balance within the CLEC16A-RNF41 mitophagy complex depends on an intrinsically disordered protein region

CLEC16A is an E3 ubiquitin ligase that regulates mitochondrial quality control through mitophagy and is associated with over 20 human diseases. CLEC16A forms a complex with another E3 ligase, RNF41, and a ubiquitin-specific peptidase, USP8; however, regions that regulate CLEC16A activity or the assembly of the tripartite mitophagy regulatory complex are unknown. Here, we report that CLEC16A contains an internal intrinsically disordered protein region (IDPR) that is crucial for CLEC16A function and turnover. IDPRs lack a fixed secondary structure and possess emerging, yet still equivocal roles in protein stability, interactions, and enzymatic activity. We find that the internal IDPR of CLEC16A is crucial for its degradation. CLEC16A turnover was promoted by RNF41, which binds and acts upon the internal IDPR to destabilize CLEC16A. Loss of this internal IDPR also destabilized the ubiquitin-dependent tripartite CLEC16A-RNF41-USP8 complex. Finally, the presence of an internal IDPR within CLEC16A was confirmed using NMR and circular dichroism spectroscopy. Together, our studies reveal that an IDPR is essential to control the reciprocal regulatory balance between CLEC16A and RNF41, which could be targeted to improve mitochondrial health in disease.

Disease association and comparative genomics of compositional bias in human proteins

Background: The evolutionary rate of disordered protein regions varies greatly due to the lack of structural constraints. So far, few studies have investigated the presence/absence patterns of compositional bias, indicative of disorder, across phylogenies in conjunction with human disease. In this study, we report a genome-wide analysis of compositional bias association with disease in human proteins and their taxonomic distribution. Methods: The human genome protein set provided by the Ensembl database was annotated and analysed with respect to both disease associations and the detection of compositional bias. The Uniprot Reference Proteome dataset, containing 11297 proteomes was used as target dataset for the comparative genomics of a well-defined subset of the Human Genome, including 100 characteristic, compositionally biased proteins, some linked to disease. Results: Cross-evaluation of compositional bias and disease-association in the human genome reveals a significant bias towards biased regions in disease-associated genes, with charged, hydrophilic amino acids appearing as over-represented. The phylogenetic profiling of 17 disease-associated, proteins with compositional bias across 11297 proteomes captures characteristic taxonomic distribution patterns. Conclusions: This is the first time that a combined genome-wide analysis of compositional bias, disease-association and taxonomic distribution of human proteins is reported, covering structural, functional, and evolutionary properties. The reported framework can form the basis for large-scale, follow-up projects, encompassing the entire human genome and all known gene-disease associations.

Structural Insights into ATP-Sensitive Potassium Channel Mechanics: A Role of Intrinsically Disordered Regions

Commonly used techniques, such as CryoEM or X-ray, are not able to capture the structural reorganizations of disordered regions of proteins (IDR); therefore, it is difficult to assess their functions in proteins based exclusively on experiments. To fill this gap, we used computational molecular dynamics (MD) simulation methods to capture IDR dynamics and trace biological function-related interactions in the Kir6.2/SUR1 potassium channel. This ATP-sensitive octameric complex, one of the critical elements in the insulin secretion process in human pancreatic β-cells, has four to five large, disordered fragments. Using unique MD simulations of the full Kir6.2/SUR1 channel complex, we present an in-depth analysis of the dynamics of the disordered regions and discuss the possible functions they could have in this system. Our MD results confirmed the crucial role of the N-terminus of the Kir6.2 fragment and the L0-loop of the SUR1 protein in the transfer of mechanical signals between domains that trigger insulin release. Moreover, we show that the presence of IDRs affects natural ligand binding. Our research takes us one step further toward understanding the action of this vital complex.

Chromosome Translocation t(10;19)(q26;q13) in a CIC-sarcoma

BACKGROUND/AIM

CIC-sarcomas are characterized by rearrangements of the capicua transcriptional repressor (CIC) gene on chromosome subband 19q13.2, generating chimeras in which CIC is the 5'-end partner. Most reported CIC-sarcomas have been detected using PCR amplifications together with Sanger sequencing, high throughput sequencing, and fluorescence in situ hybridization (FISH). Only a few CIC-rearranged tumors have been characterized cytogenetically. Here, we describe the cytogenetic and molecular genetic features of a CIC-sarcoma carrying a t(10;19)(q26;q13), a chromosomal rearrangement not previously detected in such neoplasms.

MATERIALS AND METHODS

A round cell sarcoma removed from the right thigh of a 57-year-old man was investigated by G-banding cytogenetics, FISH, PCR and Sanger sequencing.

RESULTS

The tumor cells had three cytogenetically related clones with the translocations t(9;18)(q22;q21) and t(10;19)(q26;q13) common to all of them. FISH with a BAC probe containing the CIC gene hybridized to the normal chromosome 19, to der(10)t(10;19), and to der(19)t(10;19). PCR using tumor cDNA as template together with Sanger sequencing detected two CIC::DUX4 fusion transcripts which both had a stop TAG codon immediately after the fusion point. Both transcripts are predicted to encode truncated CIC polypeptides lacking the carboxy terminal part of the native protein. This missing part is crucial for CIC's DNA binding capacity and interaction with other proteins.

CONCLUSION

In addition to demonstrating that CIC rearrangement in sarcomas can occur via the microscopically visible translocation t(10;19)(q26;q13), the findings in the present case provide evidence that the missing part in CIC-truncated proteins has important functions whose loss may be important in tumorigenesis.

Phase separation in viral infections

Viruses rely on the reprogramming of cellular processes to enable efficient viral replication; this often requires subcompartmentalization within the host cell. Liquid-liquid phase separation (LLPS) has emerged as a fundamental principle to organize and subdivide cellular processes, and plays an important role in viral life cycles. Despite substantial advances in the field, elucidating the exact organization and function of these organelles remains a major challenge. In this review, we summarize the biochemical basis of condensate formation, the role of LLPS during viral infection, and interplay of LLPS with innate immune responses. Finally, we discuss possible strategies and molecules to modulate LLPS during viral infections.

Evolutionary adaptation of the protein folding pathway for secretability

Secretory preproteins of the Sec pathway are targeted post-translationally and cross cellular membranes through translocases. During cytoplasmic transit, mature domains remain non-folded for translocase recognition/translocation. After translocation and signal peptide cleavage, mature domains fold to native states in the bacterial periplasm or traffic further. We sought the structural basis for delayed mature domain folding and how signal peptides regulate it. We compared how evolution diversified a periplasmic peptidyl-prolyl isomerase PpiA mature domain from its structural cytoplasmic PpiB twin. Global and local hydrogen-deuterium exchange mass spectrometry showed that PpiA is a slower folder. We defined at near-residue resolution hierarchical folding initiated by similar foldons in the twins, at different order and rates. PpiA folding is delayed by less hydrophobic native contacts, frustrated residues and a β-turn in the earliest foldon and by signal peptide-mediated disruption of foldon hierarchy. When selected PpiA residues and/or its signal peptide were grafted onto PpiB, they converted it into a slow folder with enhanced in vivo secretion. These structural adaptations in a secretory protein facilitate trafficking.

Amyloidogenesis: What Do We Know So Far?

The study of protein aggregation, and amyloidosis in particular, has gained considerable interest in recent times. Several neurodegenerative diseases, such as Alzheimer's (AD) and Parkinson's (PD) show a characteristic buildup of proteinaceous aggregates in several organs, especially the brain. Despite the enormous upsurge in research articles in this arena, it would not be incorrect to say that we still lack a crystal-clear idea surrounding these notorious aggregates. In this review, we attempt to present a holistic picture on protein aggregation and amyloids in particular. Using a chronological order of discoveries, we present the case of amyloids right from the onset of their discovery, various biophysical techniques, including analysis of the structure, the mechanisms and kinetics of the formation of amyloids. We have discussed important questions on whether aggregation and amyloidosis are restricted to a subset of specific proteins or more broadly influenced by the biophysiochemical and cellular environment. The therapeutic strategies and the significant failure rate of drugs in clinical trials pertaining to these neurodegenerative diseases have been also discussed at length. At a time when the COVID-19 pandemic has hit the globe hard, the review also discusses the plausibility of the far-reaching consequences posed by the virus, such as triggering early onset of amyloidosis. Finally, the application(s) of amyloids as useful biomaterials has also been discussed briefly in this review.

Deciphering the Alphabet of Disorder-Glu and Asp Act Differently on Local but Not Global Properties

Compared to folded proteins, the sequences of intrinsically disordered proteins (IDPs) are enriched in polar and charged amino acids. Glutamate is one of the most enriched amino acids in IDPs, while the chemically similar amino acid aspartate is less enriched. So far, the underlying functional differences between glutamates and aspartates in IDPs remain poorly understood. In this study, we examine the differential effects of aspartate and glutamates in IDPs by comparing the function and conformational ensemble of glutamate and aspartate variants of the disordered protein Dss1, using a range of assays, including interaction studies, nuclear magnetic resonance spectroscopy, small-angle X-ray scattering and molecular dynamics simulation. First, we analyze the sequences of the rapidly growing database of experimentally verified IDPs (DisProt) and show that glutamate enrichment is not caused by a taxonomy bias in IDPs. From analyses of local and global structural properties as well as cell growth and protein-protein interactions using a model acidic IDP from yeast and three Glu/Asp variants, we find that while the Glu/Asp variants support similar function and global dimensions, the variants differ in their binding affinities and population of local transient structural elements. We speculate that these local structural differences may play roles in functional diversity, where glutamates can support increased helicity, important for folding and binding, while aspartates support extended structures and form helical caps, as well as playing more relevant roles in, e.g., transactivation domains and ion-binding.

Fluorinated Protein and Peptide Materials for Biomedical Applications

Fluorination represents one of the most powerful modern design strategies to impart biomacromolecules with unique functionality, empowering them for widespread application in the biomedical realm. However, the properties of fluorinated protein materials remain unpredictable due to the heavy context-dependency of the surrounding atoms influenced by fluorine’s strong electron-withdrawing tendencies. This review aims to discern patterns and elucidate design principles governing the biochemical synthesis and rational installation of fluorine into protein and peptide sequences for diverse biomedical applications. Several case studies are presented to deconvolute the overgeneralized fluorous stabilization effect and critically examine the duplicitous nature of the resultant enhanced chemical and thermostability as it applies to use as biomimetic therapeutics, drug delivery vehicles, and bioimaging modalities.

Structural and mutational analysis of MazE6-operator DNA complex provide insights into autoregulation of toxin-antitoxin systems

Of the 10 paralogs of MazEF Toxin-Antitoxin system in Mycobacterium tuberculosis, MazEF6 plays an important role in multidrug tolerance, virulence, stress adaptation and Non Replicative Persistant (NRP) state establishment. The solution structures of the DNA binding domain of MazE6 and of its complex with the cognate operator DNA show that transcriptional regulation occurs by binding of MazE6 to an 18 bp operator sequence bearing the TANNNT motif (-10 region). Kinetics and thermodynamics of association, as determined by NMR and ITC, indicate that the nMazE6-DNA complex is of high affinity. Residues in N-terminal region of MazE6 that are key for its homodimerization, DNA binding specificity, and the base pairs in the operator DNA essential for the protein-DNA interaction, have been identified. It provides a basis for design of chemotherapeutic agents that will act via disruption of TA autoregulation, leading to cell death.

The dimeric MazE6 antitoxin binds to a specific sequence in its cognate operator DNA for autoregulation, and the key residues for dimerization and DNA binding are identified.

Low Complexity Induces Structure in Protein Regions Predicted as Intrinsically Disordered

There is increasing evidence that many intrinsically disordered regions (IDRs) in proteins play key functional roles through interactions with other proteins or nucleic acids. These interactions often exhibit a context-dependent structural behavior. We hypothesize that low complexity regions (LCRs), often found within IDRs, could have a role in inducing local structure in IDRs. To test this, we predicted IDRs in the human proteome and analyzed their structures or those of homologous sequences in the Protein Data Bank (PDB). We then identified two types of simple LCRs within IDRs: regions with only one (polyX or homorepeats) or with only two types of amino acids (polyXY). We were able to assign structural information from the PDB more often to these LCRs than to the surrounding IDRs (polyX 61.8% > polyXY 50.5% > IDRs 39.7%). The most frequently observed polyX and polyXY within IDRs contained E (Glu) or G (Gly). Structural analyses of these sequences and of homologs indicate that polyEK regions induce helical conformations, while the other most frequent LCRs induce coil structures. Our work proposes bioinformatics methods to help in the study of the structural behavior of IDRs and provides a solid basis suggesting a structuring role of LCRs within them.

Compositional Bias of Intrinsically Disordered Proteins and Regions and Their Predictions

Intrinsically disordered regions (IDRs) carry out many cellular functions and vary in length and placement in protein sequences. This diversity leads to variations in the underlying compositional biases, which were demonstrated for the short vs. long IDRs. We analyze compositional biases across four classes of disorder: fully disordered proteins; short IDRs; long IDRs; and binding IDRs. We identify three distinct biases: for the fully disordered proteins, the short IDRs and the long and binding IDRs combined. We also investigate compositional bias for putative disorder produced by leading disorder predictors and find that it is similar to the bias of the native disorder. Interestingly, the accuracy of disorder predictions across different methods is correlated with the correctness of the compositional bias of their predictions highlighting the importance of the compositional bias. The predictive quality is relatively low for the disorder classes with compositional bias that is the most different from the “generic” disorder bias, while being much higher for the classes with the most similar bias. We discover that different predictors perform best across different classes of disorder. This suggests that no single predictor is universally best and motivates the development of new architectures that combine models that target specific disorder classes.

The BCC7 Protein Contributes to the Toxoplasma Basal Pole by Interfacing between the MyoC Motor and the IMC Membrane Network

T. gondii is a eukaryotic parasite that has evolved a stage called tachyzoite which multiplies in host cells by producing two daughter cells internally. These nascent tachyzoites bud off their mother and repeat the division process until the expanding progenies escape to settle and multiply in other host cells. Over these intra- and extra-cellular phases, the tachyzoite maintains an essential apicobasal polarity that emerges through a unique bidirectional budding process of the elongating cells. This process requires the assembly of several molecular complexes that, at the nascent pole, encompass structural and myosin motor elements. To characterize a recently identified basal pole marker named BCC7 with respect to the posterior myosin J and myosin C motors, we used conventional biochemistry as well as advanced proteomic and in silico analysis in conjunction with live and super resolution microscopy of transgenic fluorescent tachyzoites. We document that BCC7 forms a ribbed ring below which myosin C motor entities distribute regularly. In addition, we identified-among 13 BCC7 putative partners-two novel and five known members of the inner membrane complex (IMC) family which ends at the apical side of the ring. Therefore, BCC7 could assist the stabilization of the IMC plaques and contribute to the parasite biomechanical properties.

Notch ankyrin domain: evolutionary rise of a thermodynamic sensor

Notch signalling pathway plays a key role in metazoan biology by contributing to resolution of binary decisions in the life cycle of cells during development. Outcomes such as proliferation/differentiation dichotomy are resolved by transcriptional remodelling that follows a switch from Notch^on to Notch^off state, characterised by dissociation of Notch intracellular domain (NICD) from DNA-bound RBPJ. Here we provide evidence that transitioning to the Notch^off state is regulated by heat flux, a phenomenon that aligns resolution of fate dichotomies to mitochondrial activity. A combination of phylogenetic analysis and computational biochemistry was utilised to disclose structural adaptations of Notch1 ankyrin domain that enabled function as a sensor of heat flux. We then employed DNA-based micro-thermography to measure heat flux during brain development, followed by analysis in vitro of the temperature-dependent behaviour of Notch1 in mouse neural progenitor cells. The structural capacity of NICD to operate as a thermodynamic sensor in metazoans stems from characteristic enrichment of charged acidic amino acids in β-hairpins of the ankyrin domain that amplify destabilising inter-residue electrostatic interactions and render the domain thermolabile. The instability emerges upon mitochondrial activity which raises the perinuclear and nuclear temperatures to 50 °C and 39 °C, respectively, leading to destabilization of Notch1 transcriptional complex and transitioning to the Notch^off state. Notch1 functions a metazoan thermodynamic sensor that is switched on by intercellular contacts, inputs heat flux as a proxy for mitochondrial activity in the Notch^on state via the ankyrin domain and is eventually switched off in a temperature-dependent manner.

LAG3 associates with TCR-CD3 complexes and suppresses signaling by driving co-receptor-Lck dissociation

LAG3 is an inhibitory receptor that is highly expressed on exhausted T cells. Although LAG3-targeting immunotherapeutics are currently in clinical trials, how LAG3 inhibits T cell function remains unclear. Here, we show that LAG3 moved to the immunological synapse and associated with the T cell receptor (TCR)-CD3 complex in CD4+ and CD8+ T cells, in the absence of binding to major histocompatibility complex class II-its canonical ligand. Mechanistically, a phylogenetically conserved, acidic, tandem glutamic acid-proline repeat in the LAG3 cytoplasmic tail lowered the pH at the immune synapse and caused dissociation of the tyrosine kinase Lck from the CD4 or CD8 co-receptor, which resulted in a loss of co-receptor-TCR signaling and limited T cell activation. These observations indicated that LAG3 functioned as a signal disruptor in a major histocompatibility complex class II-independent manner, and provide insight into the mechanism of action of LAG3-targeting immunotherapies.

AlphaFold2: A Role for Disordered Protein/Region Prediction?

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

Intra-protein interactions of SARS-CoV-2 and SARS: a bioinformatic analysis for plausible explanation regarding stability, divergency, and severity

The current nightmare for the whole world is COVID-19. The occurrence of concentrated pneumonia cases in Wuhan city, Hubei province of China, was first reported on December 30, 2019. SARS-CoV first disclosed in 2002 but had not outspread worldwide. After 18 years, in 2020, it reemerged and outspread worldwide as SARS-CoV-2 (COVID-19), as the most dangerous virus-creating disease in the world. Is it possible to create a favorable evolution within the short time (18 years)? If possible, then what are those properties or factors that are changed in SARS-CoV-2 to make it undefeated? What are the fundamental differences between SARS-CoV-2 and SARS? The study is one of the initiatives to find out all those queries. Here, four types of protein sequences from SARS-CoV-2 and SARS were retrieved from the database to study their physicochemical and structural properties. Results showed that charged residues are playing a pivotal role in SARS-CoV-2 evolution and contribute to the helix stabilization. The formation of the cyclic salt bridge and other intra-protein interactions specially network aromatic-aromatic interaction also play the crucial role in SAS-CoV-2. This comparative study will help to understand the evolution from SARS to SARS-CoV-2 and helpful in protein engineering.

Intrinsically Disordered Proteins: Perspective on COVID-19 Infection and Drug Discovery

Since the beginning of the COVID-19 pandemic caused by SARS-CoV-2, millions of patients have been diagnosed and many of them have died from the disease worldwide. The identification of novel therapeutic targets are of utmost significance for prevention and treatment of COVID-19. SARS-CoV-2 is a single-stranded RNA virus with a 30 kb genome packaged into a membrane-enveloped virion, transcribing several tens of proteins. The belief that the amino acid sequence of proteins determines their 3D structure which, in turn, determines their function has been a central principle of molecular biology for a long time. Recently, it has been increasingly realized, however, that there is a large group of proteins that lack a fixed or ordered 3D structure, yet they exhibit important biological activities─so-called intrinsically disordered proteins and protein regions (IDPs/IDRs). Disordered regions in viral proteins are generally associated with viral infectivity and pathogenicity because they endow the viral proteins the ability to easily and promiscuously bind to host proteins; therefore, the proteome of SARS-CoV-2 has been thoroughly examined for intrinsic disorder. It has been recognized that, in fact, the SARS-CoV-2 proteome exhibits significant levels of structural order, with only the nucleocapsid (N) structural protein and two of the nonstructural proteins being highly disordered. The spike (S) protein of SARS-CoV-2 exhibits significant levels of structural order, yet its predicted percentage of intrinsic disorder is still higher than that of the spike protein of SARS-CoV. Noteworthy, however, even though IDPs/IDRs are not common in the SARS-CoV-2 proteome, the existing ones play major roles in the functioning and virulence of the virus and are thus promising drug targets for rational antiviral drug design. Presented here is a COVID-19 perspective on the intrinsically disordered proteins, summarizing recent results on the SARS-CoV-2 proteome disorder features, their physiological and pathological relevance, and their prominence as prospective drug target sites.

Unifying coarse-grained force fields for folded and disordered proteins

Liquid-liquid phase separation drives the formation of biological condensates that play essential roles in transcriptional regulation and signal sensing. Computational modeling could provide high-resolution structural characterizations of these condensates and help uncover physicochemical interactions that dictate their stability. However, many protein molecules involved in phase separation often contain multiple ordered domains connected with flexible, structureless linkers. Simulating such proteins necessitates force fields with consistent accuracy for both folded and disordered proteins. We provide a critical review of existing coarse-grained force fields for disordered proteins and highlight the challenges in their application to folded proteins. After discussing existing algorithms for force field parameterization, we propose an optimization strategy that should lead to computer models with improved transferability across protein types.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

Matrin3: Disorder and ALS Pathogenesis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by the degeneration of both upper and lower motor neurons in the brain and spinal cord. ALS is associated with protein misfolding and inclusion formation involving RNA-binding proteins, including TAR DNA-binding protein (TDP-43) and fused in sarcoma (FUS). The 125-kDa Matrin3 is a highly conserved nuclear DNA/RNA-binding protein that is implicated in many cellular processes, including binding and stabilizing mRNA, regulating mRNA nuclear export, modulating alternative splicing, and managing chromosomal distribution. Mutations in MATR3, the gene encoding Matrin3, have been identified as causal in familial ALS (fALS). Matrin3 lacks a prion-like domain that characterizes many other ALS-associated RNA-binding proteins, including TDP-43 and FUS, however, our bioinformatics analyses and preliminary studies document that Matrin3 contains long intrinsically disordered regions that may facilitate promiscuous interactions with many proteins and may contribute to its misfolding. In addition, these disordered regions in Matrin3 undergo numerous post-translational modifications, including phosphorylation, ubiquitination and acetylation that modulate the function and misfolding of the protein. Here we discuss the disordered nature of Matrin3 and review the factors that may promote its misfolding and aggregation, two elements that might explain its role in ALS pathogenesis.

Chromatin-contact atlas reveals disorder-mediated protein interactions and moonlighting chromatin-associated RBPs

RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1459 chromatin-associated proteins, ∼48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 403 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about 40% of the caRBPs bind chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin.

Partial structure, dampened mobility, and modest impact of a His tag in the SARS-CoV-2 Nsp2 C-terminal region

Intrinsically disordered proteins (IDPs) play essential roles in regulating physiological processes in eukaryotic cells. Many viruses use their own IDPs to "hack" these processes to deactivate host defenses and promote viral growth. Thus, viral IDPs are attractive drug targets. While IDPs are hard to study by X-ray crystallography or cryo-EM, atomic level information on their conformational preferences and dynamics can be obtained using NMR spectroscopy. SARS-CoV-2 Nsp2, whose C-terminal region (CtR) is predicted to be disordered, interacts with human proteins that regulate translation initiation and endosome vesicle sorting. Molecules that block these interactions could be valuable leads for drug development. The 13Cβ and backbone 13CO, 1HN, 13Cα, and 15N nuclei of Nsp2's 45-residue CtR were assigned and used to characterize its structure and dynamics in three contexts; namely: (1) retaining an N-terminal His tag, (2) without the His tag and with an adventitious internal cleavage, and (3) lacking both the His tag and the internal cleavage. Two five-residue segments adopting a minor extended population were identified. Overall, the dynamic behavior is midway between a completely rigid and a fully flexible chain. Whereas the presence of an N-terminal His tag and internal cleavage stiffen and loosen, respectively, neighboring residues, they do not affect the tendency of two regions to populate extended conformations.

Identification of a Chemotherapeutic Lead Molecule for the Potential Disruption of the FAM72A-UNG2 Interaction to Interfere with Genome Stability, Centromere Formation, and Genome Editing

Family with sequence similarity 72 A (FAM72A) is a pivotal mitosis-promoting factor that is highly expressed in various types of cancer. FAM72A interacts with the uracil-DNA glycosylase UNG2 to prevent mutagenesis by eliminating uracil from DNA molecules through cleaving the N-glycosylic bond and initiating the base excision repair pathway, thus maintaining genome integrity. In the present study, we determined a specific FAM72A-UNG2 heterodimer protein interaction using molecular docking and dynamics. In addition, through in silico screening, we identified withaferin B as a molecule that can specifically prevent the FAM72A-UNG2 interaction by blocking its cell signaling pathways. Our results provide an excellent basis for possible therapeutic approaches in the clinical treatment of cancer.

Expanding the Disorder-Function Paradigm in the C-Terminal Tails of Erbbs

ErbBs are receptor tyrosine kinases involved not only in development, but also in a wide variety of diseases, particularly cancer. Their extracellular, transmembrane, juxtamembrane, and kinase folded domains were described extensively over the past 20 years, structurally and functionally. However, their whole C-terminal tails (CTs) following the kinase domain were only described at atomic resolution in the last 4 years. They were shown to be intrinsically disordered. The CTs are known to be tyrosine-phosphorylated when the activated homo- or hetero-dimers of ErbBs are formed. Their phosphorylation triggers interaction with phosphotyrosine binding (PTB) or Src Homology 2 (SH2) domains and activates several signaling pathways controling cellular motility, proliferation, adhesion, and apoptosis. Beyond this passive role of phosphorylated domain and site display for partners, recent structural and function studies unveiled active roles in regulation of phosphorylation and interaction: the CT regulates activity of the kinase domain; different phosphorylation states have different compaction levels, potentially modulating the succession of phosphorylation events; and prolines have an important role in structure, dynamics, and possibly regulatory interactions. Here, we review both the canonical role of the disordered CT domains of ErbBs as phosphotyrosine display domains and the recent findings that expand the known range of their regulation functions linked to specific structural and dynamic features.

Differential regulation of mRNA fate by the human Ccr4-Not complex is driven by coding sequence composition and mRNA localization

BACKGROUND

Regulation of protein output at the level of translation allows for a rapid adaptation to dynamic changes to the cell's requirements. This precise control of gene expression is achieved by complex and interlinked biochemical processes that modulate both the protein synthesis rate and stability of each individual mRNA. A major factor coordinating this regulation is the Ccr4-Not complex. Despite playing a role in most stages of the mRNA life cycle, no attempt has been made to take a global integrated view of how the Ccr4-Not complex affects gene expression.

RESULTS

This study has taken a comprehensive approach to investigate post-transcriptional regulation mediated by the Ccr4-Not complex assessing steady-state mRNA levels, ribosome position, mRNA stability, and protein production transcriptome-wide. Depletion of the scaffold protein CNOT1 results in a global upregulation of mRNA stability and the preferential stabilization of mRNAs enriched for G/C-ending codons. We also uncover that mRNAs targeted to the ER for their translation have reduced translational efficiency when CNOT1 is depleted, specifically downstream of the signal sequence cleavage site. In contrast, translationally upregulated mRNAs are normally localized in p-bodies, contain disorder-promoting amino acids, and encode nuclear localized proteins. Finally, we identify ribosome pause sites that are resolved or induced by the depletion of CNOT1.

CONCLUSIONS

We define the key mRNA features that determine how the human Ccr4-Not complex differentially regulates mRNA fate and protein synthesis through a mechanism linked to codon composition, amino acid usage, and mRNA localization.

Bi-allelic variants in HOPS complex subunit VPS41 cause cerebellar ataxia and abnormal membrane trafficking

Membrane trafficking is a complex, essential process in eukaryotic cells responsible for protein transport and processing. Deficiencies in vacuolar protein sorting (VPS) proteins, key regulators of trafficking, cause abnormal intracellular segregation of macromolecules and organelles and are linked to human disease. VPS proteins function as part of complexes such as the homotypic fusion and vacuole protein sorting (HOPS) tethering complex, composed of VPS11, VPS16, VPS18, VPS33A, VPS39 and VPS41. The HOPS-specific subunit VPS41 has been reported to promote viability of dopaminergic neurons in Parkinson’s disease but to date has not been linked to human disease. Here, we describe five unrelated families with nine affected individuals, all carrying homozygous variants in VPS41 that we show impact protein function. All affected individuals presented with a progressive neurodevelopmental disorder consisting of cognitive impairment, cerebellar atrophy/hypoplasia, motor dysfunction with ataxia and dystonia, and nystagmus. Zebrafish disease modelling supports the involvement of VPS41 dysfunction in the disorder, indicating lysosomal dysregulation throughout the brain and providing support for cerebellar and microglial abnormalities when vps41 was mutated. This provides the first example of human disease linked to the HOPS-specific subunit VPS41 and suggests the importance of HOPS complex activity for cerebellar function.

Energy landscapes of protein aggregation and conformation switching in intrinsically disordered proteins

The protein folding problem was apparently solved recently by the advent of a deep learning method for protein structure prediction called AlphaFold. However, this program is not able to make predictions about the protein folding pathways. Moreover, it only treats about half of the human proteome, as the remaining proteins are intrinsically disordered or contain disordered regions. By definition these proteins differ from natively folded proteins and do not adopt a properly folded structure in solution. However these intrinsically disordered proteins (IDPs) also systematically differ in amino acid composition and uniquely often become folded upon binding to an interaction partner. These factors preclude solving IDP structures by current machine-learning methods like AlphaFold, which also cannot solve the protein aggregation problem, since this meta-folding process can give rise to different aggregate sizes and structures. An alternative computational method is provided by molecular dynamics simulations that already successfully explored the energy landscapes of IDP conformational switching and protein aggregation in multiple cases. These energy landscapes are very different from those of 'simple' protein folding, where one energy funnel leads to a unique protein structure. Instead, the energy landscapes of IDP conformational switching and protein aggregation feature a number of minima for different competing low-energy structures. In this review, I discuss the characteristics of these multifunneled energy landscapes in detail, illustrated by molecular dynamics simulations that elucidated the underlying conformational transitions and aggregation processes.

Priming mycobacterial ESX-secreted protein B to form a channel-like structure

ESX-1 is a major virulence factor of Mycobacterium tuberculosis, a secretion machinery directly involved in the survival of the microorganism from the immune system defence. It disrupts the phagosome membrane of the host cell through a contact-dependent mechanism. Recently, the structure of the inner-membrane core complex of the homologous ESX-3 and ESX-5 was resolved; however, the elements involved in the secretion through the outer membrane or those acting on the host cell membrane are unknown. Protein substrates might form this missing element. Here, we describe the oligomerisation process of the ESX-1 substrate EspB, which occurs upon cleavage of its C-terminal region and is favoured by an acidic environment. Cryo-electron microscopy data shows that quaternary structure of EspB is conserved across slow growing species, but not in the fast growing M. smegmatis. EspB assembles into a channel with dimensions and characteristics suitable for the transit of ESX-1 substrates, as shown by the presence of another EspB trapped within. Our results provide insight into the structure and assembly of EspB, and suggests a possible function as a structural element of ESX-1.

Deep conservation of ribosome stall sites across RNA processing genes

The rate of translation can vary depending on the mRNA template. During the elongation phase the ribosome can transiently pause or permanently stall. A pause can provide the nascent protein with the time to fold or be transported, while stalling can serve as quality control and trigger degradation of aberrant mRNA and peptide. Ribosome profiling has allowed for the genome-wide detection of such pauses and stalls, but due to library-specific biases, these predictions are often unreliable. Here, we take advantage of the deep conservation of protein synthesis machinery, hypothesizing that similar conservation could exist for functionally important locations of ribosome slowdown, here collectively called stall sites. We analyze multiple ribosome profiling datasets from phylogenetically diverse eukaryotes: yeast, fruit fly, zebrafish, mouse and human to identify conserved stall sites. We find thousands of stall sites across multiple species, with the enrichment of proline, glycine and negatively charged amino acids around conserved stalling. Many of the sites are found in RNA processing genes, suggesting that stalling might have a conserved role in RNA metabolism. In summary, our results provide a rich resource for the study of conserved stalling and indicate possible roles of stalling in gene regulation.

Consistent Force Field Captures Homologue-Resolved HP1 Phase Separation

Many proteins have been shown to function via liquid-liquid phase separation. Computational modeling could offer much needed structural details of protein condensates and reveal the set of molecular interactions that dictate their stability. However, the presence of both ordered and disordered domains in these proteins places a high demand on the model accuracy. Here, we present an algorithm to derive a coarse-grained force field, MOFF, which can model both ordered and disordered proteins with consistent accuracy. It combines maximum entropy biasing, least-squares fitting, and basic principles of energy landscape theory to ensure that MOFF recreates experimental radii of gyration while predicting the folded structures for globular proteins with lower energy. The theta temperature determined from MOFF separates ordered and disordered proteins at 300 K and exhibits a strikingly linear relationship with amino acid sequence composition. We further applied MOFF to study the phase behavior of HP1, an essential protein for post-translational modification and spatial organization of chromatin. The force field successfully resolved the structural difference of two HP1 homologues despite their high sequence similarity. We carried out large-scale simulations with hundreds of proteins to determine the critical temperature of phase separation and uncover multivalent interactions that stabilize higher-order assemblies. In all, our work makes significant methodological strides to connect theories of ordered and disordered proteins and provides a powerful tool for studying liquid-liquid phase separation with near-atomistic details.

Discovery of Arginine Methylation, Phosphorylation, and Their Co-occurrence in Condensate-Associated Proteins in Saccharomyces cerevisiae

The formation of condensates in membraneless organelles is thought to be driven by protein phase separation. Arginine methylation and serine/threonine phosphorylation are important in the phase separation process; however, these post-translational modifications are often present in intrinsically disordered regions that are difficult to analyze with standard proteomic techniques. To understand their presence and co-occurrence in condensate-associated proteins, here, we use a multiprotease and multi-tandem mass spectrometry (MS/MS) fragmentation approach, coupled with heavy methyl stable isotope labeling of amino acids in cell culture (SILAC) and phospho- or methyl-peptide enrichment. For Saccharomyces cerevisiae, we report a 50% increase in the known arginine methylproteome, involving 15 proteins that are all condensate-associated. Importantly, some of these proteins have arginine methylation on all predicted sites-providing evidence that this modification can be pervasive. We explored whether arginine-methylated, condensate-associated proteins are also phosphorylated and found 12 such proteins to carry phosphorylated serine or threonine. In Npl3, Ded1, and Sbp1, single peptides were found to carry both modifications, indicating a co-occurrence in close proximity and on the same protein molecule. These co-modifications occur in regions of disorder, whereas arginine methylation is typically on regions of disorder that are also basic. For phosphorylation, its association with charged regions of condensate-associated proteins was less consistent, although some regions with multisite phosphorylation sites were strongly acidic. We conclude that arginine-methylated proteins associated with condensates are typically also modified with protein phosphorylation.

New genotypes of Helicobacter Pylori VacA d-region identified from global strains

BACKGROUND

Pathogenesis of Helicobacter Pylori (HP) vacuolating toxin A (vacA) depends on polymorphic diversity within the signal (s), middle (m), intermediate (i), deletion (d) and c-regions. These regions show distinct allelic diversity. The s-region, m-region and the c-region (a 15 bp deletion at the 3'-end region of the p55 domain of the vacA gene) exist as 2 types (s1, s2, m1, m2, c1 and c2), while the i-region has 3 allelic types (i1, i2 and i3). The locus of d-region of the vacA gene has also been classified into 2 genotypes, namely d1 and d2. We investigated the "d-region"/"loop region" through bioinformatics, to predict its properties and relation to disease. One thousand two hundred fifty-nine strains from the NCBI nucleotide database and the dryad database with complete vacA sequences were included in the study. The sequences were aligned using BioEdit and analyzed using Lasergene and BLAST. The secondary structure and physicochemical properties of the region were predicted using PredictProtein.

RESULTS

We identified 31 highly polymorphic genotypes in the "d-region", with a mean length of 34 amino acids (9 ~ 55 amino acids). We further classified the 31 genotypes into 3 main types, namely K-type (strains starting with the KDKP motif in the "d-region"), Q-type (strains starting with the KNQT motif), and E-type (strains starting with the ESKT motif) respectively. The most common type, K-type, is more prevalent in cancer patients (80.87%) and is associated with the s1i1m1c1 genotypes (P < .01). Incidentally, a new region expressing sequence diversity (2 aa deletion) at the C-terminus of the p55 domain of vacA was identified during bioinformatics analysis.

CONCLUSIONS

Prediction of secondary structures shows that the "d-region" adopts a loop conformation and is a disordered region.

Molecular simulations of IDPs: From ensemble generation to IDP interactions leading to disorder-to-order transitions

Intrinsically disordered proteins (IDPs) lack a well-defined three-dimensional structure but do exhibit some dynamical and structural ordering. The structural plasticity of IDPs indicates that entropy-driven motions are crucial for their function. Many IDPs undergo function-related disorder-to-order transitions upon by their interaction with specific binding partners. Approaches that are based on both experimental and theoretical tools enable the biophysical characterization of IDPs. Molecular simulations provide insights into IDP structural ensembles and disorder-to-order transition mechanisms. However, such studies depend strongly on the chosen force field parameters and simulation techniques. In this chapter, we provide an overview of IDP characteristics, review all-atom force fields recently developed for IDPs, and present molecular dynamics-based simulation methods that allow IDP ensemble generation as well as the characterization of disorder-to-order transitions. In particular, we introduce metadynamics, replica exchange molecular dynamics simulations, and also kinetic models resulting from Markov State modeling, and provide various examples for the successful application of these simulation methods to IDPs.