Accurately Predicting Disordered Regions of Proteins Using Rosetta ResidueDisorder Application

Tutorial: a guide for the selection of fast and accurate computational tools for the prediction of intrinsic disorder in proteins

Intrinsic disorder is instrumental for a wide range of protein functions, and its analysis, using computational predictions from primary structures, complements secondary and tertiary structure-based approaches. In this Tutorial, we provide an overview and comparison of 23 publicly available computational tools with complementary parameters useful for intrinsic disorder prediction, partly relying on results from the Critical Assessment of protein Intrinsic Disorder prediction experiment. We consider factors such as accuracy, runtime, availability and the need for functional insights. The selected tools are available as web servers and downloadable programs, offer state-of-the-art predictions and can be used in a high-throughput manner. We provide examples and instructions for the selected tools to illustrate practical aspects related to the submission, collection and interpretation of predictions, as well as the timing and their limitations. We highlight two predictors for intrinsically disordered proteins, flDPnn as accurate and fast and IUPred as very fast and moderately accurate, while suggesting ANCHOR2 and MoRFchibi as two of the best-performing predictors for intrinsically disordered region binding. We link these tools to additional resources, including databases of predictions and web servers that integrate multiple predictive methods. Altogether, this Tutorial provides a hands-on guide to comparatively evaluating multiple predictors, submitting and collecting their own predictions, and reading and interpreting results. It is suitable for experimentalists and computational biologists interested in accurately and conveniently identifying intrinsic disorder, facilitating the functional characterization of the rapidly growing collections of protein sequences.

Predicting ion mobility collision cross sections using projection approximation with ROSIE-PARCS webserver

Ion mobility coupled to mass spectrometry informs on the shape and size of protein structures in the form of a collision cross section (CCSIM). Although there are several computational methods for predicting CCSIM based on protein structures, including our previously developed projection approximation using rough circular shapes (PARCS), the process usually requires prior experience with the command-line interface. To overcome this challenge, here we present a web application on the Rosetta Online Server that Includes Everyone (ROSIE) webserver to predict CCSIM from protein structure using projection approximation with PARCS. In this web interface, the user is only required to provide one or more PDB files as input. Results from our case studies suggest that CCSIM predictions (with ROSIE-PARCS) are highly accurate with an average error of 6.12%. Furthermore, the absolute difference between CCSIM and CCSPARCS can help in distinguishing accurate from inaccurate AlphaFold2 protein structure predictions. ROSIE-PARCS is designed with a user-friendly interface, is available publicly and is free to use. The ROSIE-PARCS web interface is supported by all major web browsers and can be accessed via this link (https://rosie.graylab.jhu.edu).

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

Ab Initio Modelling of the Structure of ToxA-like and MAX Fungal Effector Proteins

Pathogenic fungal diseases in crops are mediated by the release of effector proteins that facilitate infection. Characterising the structure of these fungal effectors is vital to understanding their virulence mechanisms and interactions with their hosts, which is crucial in the breeding of plant cultivars for disease resistance. Several effectors have been identified and validated experimentally; however, their lack of sequence conservation often impedes the identification and prediction of their structure using sequence similarity approaches. Structural similarity has, nonetheless, been observed within fungal effector protein families, creating interest in validating the use of computational methods to predict their tertiary structure from their sequence. We used Rosetta ab initio modelling to predict the structures of members of the ToxA-like and MAX effector families for which experimental structures are known to validate this method. An optimised approach was then used to predict the structures of phenotypically validated effectors lacking known structures. Rosetta was found to successfully predict the structure of fungal effectors in the ToxA-like and MAX families, as well as phenotypically validated but structurally unconfirmed effector sequences. Interestingly, potential new effector structural families were identified on the basis of comparisons with structural homologues and the identification of associated protein domains.

Investigating In Situ Expression of c-MYC and Candidate Ubiquitin-Specific Proteases in DLBCL and Assessment for Peptidyl Disruptor Molecule against c-MYC-USP37 Complex

Diffuse Large B-Cell Lymphoma (DLBCL) is the most common form of non-Hodgkin's lymphoma (NHL). Elevated expression of c-MYC in DLBCL is associated with poor prognosis of the disease. In different cancers, c-MYC has been found regulated by different ubiquitin-specific proteases (USPs), but to date, the role of USPs in c-MYC regulation has not been investigated in DLBCL. In this study, in situ co expression of c-MYC and three candidates USPs, USP28, USP36 and USP37, have been investigated in both the ABC and GCB subtypes of DLBCL. This shows that USP37 expression is positively correlated with the c-MYC expression in the ABC subtype of DLBCL. Structurally, both c-MYC and USP37 has shown large proportion of intrinsically disordered regions, minimizing their chances for full structure crystallization. Peptide array and docking simulations has shown that N-terminal region of c-MYC interacts directly with residues within and in proximity of catalytically active C19 domain of the USP37. Given the structural properties of the interaction sites in the c-MYC-USP37 complex, a peptidyl inhibitor has been designed. Molecular docking has shown that the peptide fits well in the targeted site of c-MYC, masking most of its residues involved in the binding with USP37. The findings could further be exploited to develop therapeutic interventions against the ABC subtype of DLBCL.

Prediction of Intrinsic Disorder Using Rosetta ResidueDisorder and AlphaFold2

The combination of deep learning and sequence data has transformed protein structure prediction and modeling, evidenced in the success of AlphaFold (AF). For this reason, many methods have been developed to take advantage of this success in areas where inaccurate structural modeling may limit computational predictiveness. For example, many methods have been developed to predict protein intrinsic disorder from sequence, including our Rosetta ResidueDisorder (RRD) approach. Intrinsically disordered regions in proteins are parts of the sequence that do not form ordered, folded structures under typical physiological conditions. In the original implementation of RRD, Rosetta ab initio models were generated, and disordered regions were predicted based on residue scores (disordered residues typically exist in regions of unfavorable scores). In this work, we show that by (i) replacing the ab initio modeling with AF (using the same scoring and disorder assignment approach) and (ii) updating the score function, the predictiveness improved significantly. Residues were better ranked by the order/disorder, evidenced by an improvement in receiver operating characteristic area-under-the-curve from 0.69 to 0.78 on a large (229 protein) and balanced data set (relatively even ordered versus disordered residues). Finally, the binary prediction accuracy also improved from 62% to 74% on the same data set. Our results show that the combined AF-RRD approach was as good as or better than all existing methods by these metrics (AF-RRD had the highest prediction accuracy).

Protein shape sampled by ion mobility mass spectrometry consistently improves protein structure prediction

Ion mobility (IM) mass spectrometry provides structural information about protein shape and size in the form of an orientationally-averaged collision cross-section (CCSIM). While IM data have been used with various computational methods, they have not yet been utilized to predict monomeric protein structure from sequence. Here, we show that IM data can significantly improve protein structure determination using the modelling suite Rosetta. We develop the Rosetta Projection Approximation using Rough Circular Shapes (PARCS) algorithm that allows for fast and accurate prediction of CCSIM from structure. Following successful testing of the PARCS algorithm, we use an integrative modelling approach to utilize IM data for protein structure prediction. Additionally, we propose a confidence metric that identifies near native models in the absence of a known structure. The results of this study demonstrate the ability of IM data to consistently improve protein structure prediction.

Computational Structure Prediction for Antibody-Antigen Complexes From Hydrogen-Deuterium Exchange Mass Spectrometry: Challenges and Outlook

Although computational structure prediction has had great successes in recent years, it regularly fails to predict the interactions of large protein complexes with residue-level accuracy, or even the correct orientation of the protein partners. The performance of computational docking can be notably enhanced by incorporating experimental data from structural biology techniques. A rapid method to probe protein-protein interactions is hydrogen-deuterium exchange mass spectrometry (HDX-MS). HDX-MS has been increasingly used for epitope-mapping of antibodies (Abs) to their respective antigens (Ags) in the past few years. In this paper, we review the current state of HDX-MS in studying protein interactions, specifically Ab-Ag interactions, and how it has been used to inform computational structure prediction calculations. Particularly, we address the limitations of HDX-MS in epitope mapping and techniques and protocols applied to overcome these barriers. Furthermore, we explore computational methods that leverage HDX-MS to aid structure prediction, including the computational simulation of HDX-MS data and the combination of HDX-MS and protein docking. We point out challenges in interpreting and incorporating HDX-MS data into Ab-Ag complex docking and highlight the opportunities they provide to build towards a more optimized hybrid method, allowing for more reliable, high throughput epitope identification.

Predicting Protein Conformational Disorder and Disordered Binding Sites

In the last two decades it has become increasingly evident that a large number of proteins adopt either a fully or a partially disordered conformation. Intrinsically disordered proteins are ubiquitous proteins that fulfill essential biological functions while lacking a stable 3D structure. Their conformational heterogeneity is encoded by the amino acid sequence, thereby allowing intrinsically disordered proteins or regions to be recognized based on their sequence properties. The identification of disordered regions facilitates the functional annotation of proteins and is instrumental for delineating boundaries of protein domains amenable to crystallization. This chapter focuses on the methods currently employed for predicting protein disorder and identifying intrinsically disordered binding sites.

Structural and dynamic studies of the human RNA binding protein RBM3 reveals the molecular basis of its oligomerization and RNA recognition

Human RNA-binding motif 3 protein (RBM3) is a cold-shock protein which functions in various aspects of global protein synthesis, cell proliferation and apoptosis by interacting with the components of basal translational machinery. RBM3 plays important roles in tumour progression and cancer metastasis, and also has been shown to be involved in neuroprotection and endoplasmic reticulum stress response. Here, we have solved the solution NMR structure of the N-terminal 84 residue RNA recognition motif (RRM) of RBM3. The remaining residues are rich in RGG and YGG motifs and are disordered. The RRM domain adopts a βαββαβ topology, which is found in many RNA-binding proteins. NMR-monitored titration experiments and molecular dynamic simulations show that the beta-sheet and two loops form the RNA-binding interface. Hydrogen bond, pi-pi and pi-cation are the key interactions between the RNA and the RRM domain. NMR, size exclusion chromatography and chemical cross-linking experiments show that RBM3 forms oligomers in solution, which is favoured by decrease in temperature, thus, potentially linking it to its function as a cold-shock protein. Temperature-dependent NMR studies revealed that oligomerization of the RRM domain occurs via nonspecific interactions. Overall, this study provides the detailed structural analysis of RRM domain of RBM3, its interaction with RNA and the molecular basis of its temperature-dependent oligomerization.

Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR

Hydrogen-deuterium exchange (HDX) measured by nuclear magnetic resonance (NMR) provides structural information for proteins relating to solvent accessibility and flexibility. While this structural information is beneficial, the data cannot be used exclusively to elucidate structures. However, the structural information provided by the HDX-NMR data can be supplemented by computational methods. In previous work, we developed an algorithm in Rosetta to predict structures using qualitative HDX-NMR data (categories of exchange rate). Here we expand on the effort, and utilize quantitative protection factors (PFs) from HDX-NMR for structure prediction. From observed correlations between PFs and solvent accessibility/flexibility measures, we present a scoring function to quantify the agreement with HDX data. Using a benchmark set of 10 proteins, an average improvement of 5.13 Å in root-mean-square deviation (RMSD) is observed for cases of inaccurate Rosetta predictions. Ultimately, seven out of 10 predictions are accurate without including HDX data, and nine out of 10 are accurate when using our PF-based HDX score.

Protein Structure Prediction from NMR Hydrogen-Deuterium Exchange Data

Amide hydrogen-deuterium exchange (HDX) has long been used to determine regional flexibility and binding sites in proteins; however, the data are too sparse for full structural characterization. Experiments that measure HDX rates, such as HDX-NMR, have far higher throughput compared to structure determination via X-ray crystallography, cryo-EM, or a full suite of NMR experiments. Data from HDX-NMR experiments encode information on the protein structure, making HDX a prime candidate to be supplemented by computational algorithms for protein structure prediction. We have developed a methodology to incorporate HDX-NMR data into ab initio protein structure prediction using the Rosetta software framework to predict structures based on experimental agreement. To demonstrate the efficacy of our algorithm, we examined 38 proteins with HDX-NMR data available, comparing the predicted model with and without the incorporation of HDX data into scoring. The root-mean-square deviation (rmsd, a measure of the average atomic distance between superimposed models) of the predicted model improved by 1.42 Å on average after incorporating the HDX-NMR data into scoring. The average rmsd improvement for the proteins where the selected model rmsd changed after incorporating HDX data was 3.63 Å, including one improvement of more than 11 Å and seven proteins improving by greater than 4 Å, with 12/15 proteins improving overall. Additionally, for independent verification, two proteins that were not part of the original benchmark were scored including HDX data, with a dramatic improvement of the selected model rmsd of nearly 9 Å for one of the proteins. Moreover, we have developed a confidence metric allowing us to successfully identify near-native models in the absence of a native structure. Improvement in model selection with a strong confidence measure demonstrates that protein structure prediction with HDX-NMR is a powerful tool which can be performed with minimal additional computational strain and expense.

Predicting substitutions to modulate disorder and stability in coiled-coils

Background

Coiled-coils are described as stable structural motifs, where two or more helices wind around each other. However, coiled-coils are associated with local mobility and intrinsic disorder. Intrinsically disordered regions in proteins are characterized by lack of stable secondary and tertiary structure under physiological conditions in vitro. They are increasingly recognized as important for protein function. However, characterizing their behaviour in solution and determining precisely the extent of disorder of a protein region remains challenging, both experimentally and computationally.

Results

In this work, we propose a computational framework to quantify the extent of disorder within a coiled-coil in solution and to help design substitutions modulating such disorder. Our method relies on the analysis of conformational ensembles generated by relatively short all-atom Molecular Dynamics (MD) simulations. We apply it to the phosphoprotein multimerisation domains (PMD) of Measles virus (MeV) and Nipah virus (NiV), both forming tetrameric left-handed coiled-coils. We show that our method can help quantify the extent of disorder of the C-terminus region of MeV and NiV PMDs from MD simulations of a few tens of nanoseconds, and without requiring an extensive exploration of the conformational space. Moreover, this study provided a conceptual framework for the rational design of substitutions aimed at modulating the stability of the coiled-coils. By assessing the impact of four substitutions known to destabilize coiled-coils, we derive a set of rules to control MeV PMD structural stability and cohesiveness. We therefore design two contrasting substitutions, one increasing the stability of the tetramer and the other increasing its flexibility.

Conclusions

Our method can be considered as a platform to reason about how to design substitutions aimed at regulating flexibility and stability.

A Unified De Novo Approach for Predicting the Structures of Ordered and Disordered Proteins

As recognition of the abundance and relevance of intrinsically disordered proteins (IDPs) continues to grow, demand increases for methods that can rapidly predict the conformational ensembles populated by these proteins. To date, IDP simulations have largely been dominated by molecular dynamics (MD) simulations, which require significant compute times and/or complex hardware. Recent developments in MD have afforded methods capable of simulating both ordered and disordered proteins, yet to date, accurate fold prediction from a sequence has been dominated by Monte Carlo (MC)-based methods such as Rosetta. To overcome the limitations of current approaches in IDP simulation using Rosetta while maintaining its utility for modeling folded domains, we developed PyRosetta-based algorithms that allow for the accurate de novo prediction of proteins across all degrees of foldedness along with structural ensembles of disordered proteins. Our simulations have accuracy comparable to state-of-the-art MD with vastly reduced computational demands.

Predicting Protein Complex Structure from Surface-Induced Dissociation Mass Spectrometry Data

Recently, mass spectrometry (MS) has become a viable method for elucidation of protein structure. Surface-induced dissociation (SID), colliding multiply charged protein complexes or other ions with a surface, has been paired with native MS to provide useful structural information such as connectivity and topology for many different protein complexes. We recently showed that SID gives information not only on connectivity and topology but also on relative interface strengths. However, SID has not yet been coupled with computational structure prediction methods that could use the sparse information from SID to improve the prediction of quaternary structures, i.e., how protein subunits interact with each other to form complexes. Protein-protein docking, a computational method to predict the quaternary structure of protein complexes, can be used in combination with subunit structures from X-ray crystallography and NMR in situations where it is difficult to obtain an experimental structure of an entire complex. While de novo structure prediction can be successful, many studies have shown that inclusion of experimental data can greatly increase prediction accuracy. In this study, we show that the appearance energy (AE, defined as 10% fragmentation) extracted from SID can be used in combination with Rosetta to successfully evaluate protein-protein docking poses. We developed an improved model to predict measured SID AEs and incorporated this model into a scoring function that combines the RosettaDock scoring function with a novel SID scoring term, which quantifies agreement between experiments and structures generated from RosettaDock. As a proof of principle, we tested the effectiveness of these restraints on 57 systems using ideal SID AE data (AE determined from crystal structures using the predictive model). When theoretical AEs were used, the RMSD of the selected structure improved or stayed the same in 95% of cases. When experimental SID data were incorporated on a different set of systems, the method predicted near-native structures (less than 2 Å root-mean-square deviation, RMSD, from native) for 6/9 tested cases, while unrestrained RosettaDock (without SID data) only predicted 3/9 such cases. Score versus RMSD funnel profiles were also improved when SID data were included. Additionally, we developed a confidence measure to evaluate predicted model quality in the absence of a crystal structure.

Measuring Intrinsic Disorder and Tracking Conformational Transitions Using Rosetta ResidueDisorder

Many proteins contain regions of intrinsic disorder, not folding into unique, stable conformations. Numerous experimental methods have been developed to measure the disorder of all or select residues. In the absence of experimental data, computational methods are often utilized to identify these disordered regions and thus gain a better understanding of both structure and function. Many freely available computational methods have been developed to predict regions of intrinsic disorder from the primary sequence of a protein, including our recently developed Rosetta ResidueDisorder. While these methods are very useful, they are only designed to predict intrinsic disorder from the sequence. However, it would be useful to have a method that could also measure intrinsic disorder directly from structure. Such a method might also be used to identify changes in the structure of systems that can transition from folded to unfolded or vice versa, even systems that are not intrinsically disordered. Here we extended the capabilities of Rosetta ResidueDisorder to measure the intrinsic disorder from the coordinates of a single conformation of a protein. As a proof of principle, we show that ResidueDisorder can measure the intrinsic disorder from the coordinates with a higher accuracy (69.2%) than when predicted from sequence (65.4%) using a benchmark set of 229 proteins, showing that intrinsic disorder can be measured accurately from single structures over a large range of intrinsic disorder (0-100%). Additionally, we used ResidueDisorder to analyze unfolding trajectories of 12 fast-folding, nonintrinsically disordered proteins generated using molecular dynamics (MD), specifically steered MD (SMD), high-temperature MD, and accelerated MD (aMD) as well as long-time scale folding/unfolding trajectories. Using ResidueDisorder, a clear correlation between RMSD with respect to the native structure and measured fraction of denatured residues was observed. Finally, we introduced methods to predict folding/unfolding transitions as well as a native-like structure in the absence of a crystal structure from folding/unfolding MD trajectories. Rosetta ResidueDisorder is available as an application in the Rosetta software suite with the addition of new capabilities for directly identifying denatured regions and predicting events.