Quantification of Compactness and Local Order in the Ensemble of the Intrinsically Disordered Protein FCP1

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

13C Direct Detected NMR for Challenging Systems

Thanks to recent improvements in NMR spectrometer hardware and pulse sequence design, modern 13C NMR has become a useful tool for biomolecular applications. The complete assignment of a protein can be accomplished by using 13C detected multinuclear experiments and it can provide unique information relevant for the study of a variety of different biomolecules including paramagnetic proteins and intrinsically disordered proteins. A wide range of NMR observables can be measured, concurring to the structural and dynamic characterization of a protein in isolation, as part of a larger complex, or even inside a living cell. We present the different properties of 13C with respect to 1H, which provide the rationale for the experiments developed and their application, the technical aspects that need to be faced, and the many experimental variants designed to address different cases. Application areas where these experiments successfully complement proton NMR are also described.

Influence of electrostatic forces on the association kinetics and conformational ensemble of an intrinsically disordered protein

Recent work has revealed that the association of a disordered region of a protein with a folded binding partner can occur as rapidly as association between two folded proteins. This is the case for the phosphatase calcineurin (CaN) and its association with its activator calmodulin. Calmodulin binds to the intrinsically disordered regulatory domain of CaN. Previous studies have shown that electrostatic steering can accelerate the binding of folded proteins with disordered ligands. Given that electrostatic forces are strong determinants of disordered protein ensembles, the relationship between electrostatics, conformational ensembles, and quaternary interactions is unclear. Here, we employ experimental approaches to explore the impact of electrostatic interactions on the association of calmodulin with the disordered regulatory region of CaN. We find that estimated association rate constants of calmodulin with our chosen calmodulin-substrates are within the diffusion-limited regime. The association rates are dependent on the ionic strength, indicating that favorable electrostatic forces increase the rate of association. Further, we show that charged amino acids outside the calmodulin-binding site modulate the binding rate. Conformational ensembles obtained from computer simulations suggest that electrostatic interactions within the regulatory domain might bias the conformational ensemble such that the calmodulin binding region is readily accessible. Given the prevalence of charged residues in disordered protein chains, our findings are likely relevant to many protein-protein interactions.

Connecting Coil-to-Globule Transitions to Full Phase Diagrams for Intrinsically Disordered Proteins

Phase separation is thought to underlie spatial and temporal organization that is required for controlling biochemical reactions in cells. Multivalence of interaction motifs, also known as stickers, is a defining feature of proteins that drive phase separation. Intrinsically disordered proteins with stickers uniformly distributed along the linear sequence can serve as scaffold molecules that drive phase separation. The sequence-intrinsic contributions of disordered proteins to phase separation can be discerned by computing or measuring sequence-specific phase diagrams. These help to delineate the combinations of protein concentration and a suitable control parameter, such as temperature, that support phase separation. Here, we present an approach that combines detailed simulations with a numerical adaptation of an analytical Gaussian cluster theory to enable the calculation of sequence-specific phase diagrams. Our approach leverages the known equivalence between the driving forces for single-chain collapse in dilute solutions and the driving forces for phase separation in concentrated solutions. We demonstrate the application of the theory-aided computations through calculation of phase diagrams for a set of archetypal intrinsically disordered low-complexity domains. We also leverage theories to compute sequence-specific percolation lines and thereby provide a thermodynamic framework for hardening transitions that have been observed for many biomolecular condensates.

Information theoretic measures for quantifying sequence-ensemble relationships of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) contribute to a multitude of functions. De novo design of IDPs should open the door to modulating functions and phenotypes controlled by these systems. Recent design efforts have focused on compositional biases and specific sequence patterns as the design features. Analysis of the impact of these designs on sequence-function relationships indicates that individual sequence/compositional parameters are insufficient for describing sequence-function relationships in IDPs. To remedy this problem, we have developed information theoretic measures for sequence–ensemble relationships (SERs) of IDPs. These measures rely on prior availability of statistically robust conformational ensembles derived from all atom simulations. We show that the measures we have developed are useful for comparing sequence-ensemble relationships even when sequence is poorly conserved. Based on our results, we propose that de novo designs of IDPs, guided by knowledge of their SERs, should provide improved insights into their sequence–ensemble–function relationships.

Investigating the ferric ion binding site of magnetite biomineralisation protein Mms6

The biomineralization protein Mms6 has been shown to be a major player in the formation of magnetic nanoparticles both within the magnetosomes of magnetotactic bacteria and as an additive in synthetic magnetite precipitation assays. Previous studies have highlighted the ferric iron binding capability of the protein and this activity is thought to be crucial to its mineralizing properties. To understand how this protein binds ferric ions we have prepared a series of single amino acid substitutions within the C-terminal binding region of Mms6 and have used a ferric binding assay to probe the binding site at the level of individual residues which has pinpointed the key residues of E44, E50 and R55 involved in Mms6 ferric binding. No aspartic residues bound ferric ions. A nanoplasmonic sensing experiment was used to investigate the unstable EER44, 50,55AAA triple mutant in comparison to native Mms6. This suggests a difference in interaction with iron ions between the two and potential changes to the surface precipitation of iron oxide when the pH is increased. All-atom simulations suggest that disruptive mutations do not fundamentally alter the conformational preferences of the ferric binding region. Instead, disruption of these residues appears to impede a sequence-specific motif in the C-terminus critical to ferric ion binding.

Improvements to the ABSINTH Force Field for Proteins Based on Experimentally Derived Amino Acid Specific Backbone Conformational Statistics

We present an improved version of the ABSINTH implicit solvation model and force field paradigm (termed ABSINTH-C) by incorporating a grid-based term that bootstraps against experimentally derived and computationally optimized conformational statistics for blocked amino acids. These statistics provide high-resolution descriptions of the intrinsic backbone dihedral angle preferences for all 20 amino acids. The original ABSINTH model generates Ramachandran plots that are too shallow in terms of the basin structures and too permissive in terms of dihedral angle preferences. We bootstrap against the reference optimized landscapes and incorporate CMAP-like residue-specific terms that help us reproduce the intrinsic dihedral angle preferences of individual amino acids. These corrections that lead to ABSINTH-C are achieved by balancing the incorporation of the new residue-specific terms with the accuracies inherent to the original ABSINTH model. We demonstrate the robustness of ABSINTH-C through a series of examples to highlight the preservation of accuracies as well as examples that demonstrate the improvements. Our efforts show how the recent experimentally derived and computationally optimized coil-library landscapes can be used as a touchstone for quantifying errors and making improvements to molecular mechanics force fields.

Solution Ensemble of the C-Terminal Domain from the Transcription Factor Pdx1 Resembles an Excluded Volume Polymer

The pancreatic and duodenal homeobox 1 (Pdx1) is an essential pancreatic transcription factor. The C-terminal intrinsically disordered domain of Pdx1 (Pdx1-C) has a heavily biased amino acid composition; most notably, 18 of 83 residues are proline, including a hexaproline cluster near the middle of the chain. For these reasons, Pdx1-C is an attractive target for structure characterization, given the availability of suitable methods. To determine the solution ensembles of disordered proteins, we have developed a suite of ¹³C direct-detect NMR experiments that provide high spectral quality, even in the presence of strong proline enrichment. Here, we have extended our suite of NMR experiments to include four new pulse programs designed to record backbone residual dipolar couplings in a ¹³C,¹⁵N-CON detection format. Using our NMR strategy, in combination with small-angle X-ray scattering measurements and Monte Carlo simulations, we have determined that Pdx1-C is extended in solution, with a radius of gyration and internal scaling similar to that of an excluded volume polymer, and a subtle tendency toward a collapsed structure to the N-terminal side of the hexaproline sequence. This structure leaves Pdx1-C exposed for interactions with trans-regulatory co-factors that contribute with Pdx1 to transcription control in the cell.

Control of transcriptional activity by design of charge patterning in the intrinsically disordered RAM region of the Notch receptor

Intrinsically disordered regions (IDRs) play important roles in proteins that regulate gene expression. A prominent example is the intracellular domain of the Notch receptor (NICD), which regulates the transcription of Notch-responsive genes. The NICD sequence includes an intrinsically disordered RAM region and a conserved ankyrin (ANK) domain. The 111-residue RAM region mediates bivalent interactions of NICD with the transcription factor CSL. Although the sequence of RAM is poorly conserved, the linear patterning of oppositely charged residues shows minimal variation. The conformational properties of polyampholytic IDRs are governed as much by linear charge patterning as by overall charge content. Here, we used sequence design to assess how changing the charge patterning within RAM affects its conformational properties, the affinity of NICD to CSL, and Notch transcriptional activity. Increased segregation of oppositely charged residues leads to linear decreases in the global dimensions of RAM and decreases the affinity of a construct including a C-terminal ANK domain (RAMANK) for CSL. Increasing charge segregation from WT RAM sharply decreases transcriptional activation for all permutants. Activation also decreases for some, but not all, permutants with low charge segregation, although there is considerable variation. Our results suggest that the RAM linker is more than a passive tether, contributing local and/or long-range sequence features that modulate interactions within NICD and with downstream components of the Notch pathway. We propose that sequence features within IDRs have evolved to ensure an optimal balance of sequence-encoded conformational properties, interaction strengths, and cellular activities.

Visualizing single-stranded nucleic acids in solution

Single-stranded nucleic acids (ssNAs) are ubiquitous in many key cellular functions. Their flexibility limits both the number of high-resolution structures available, leaving only a small number of protein-ssNA crystal structures, while forcing solution investigations to report ensemble averages. A description of the conformational distributions of ssNAs is essential to more fully characterize biologically relevant interactions. We combine small angle X-ray scattering (SAXS) with ensemble-optimization methods (EOM) to dynamically build and refine sets of ssNA structures. By constructing candidate chains in representative dinucleotide steps and refining the models against SAXS data, a broad array of structures can be obtained to match varying solution conditions and strand sequences. In addition to the distribution of large scale structural parameters, this approach reveals, for the first time, intricate details of the phosphate backbone and underlying strand conformations. Such information on unperturbed strands will critically inform a detailed understanding of an array of problems including protein-ssNA binding, RNA folding and the polymer nature of NAs. In addition, this scheme, which couples EOM selection with an iteratively refining pool to give confidence in the underlying structures, is likely extendable to the study of other flexible systems.

Application of NMR to studies of intrinsically disordered proteins

The prevalence of intrinsically disordered protein regions, particularly in eukaryotic proteins, and their clear functional advantages for signaling and gene regulation have created an imperative for high-resolution structural and mechanistic studies. NMR spectroscopy has played a central role in enhancing not only our understanding of the intrinsically disordered native state, but also how that state contributes to biological function. While pathological functions associated with protein aggregation are well established, it has recently become clear that disordered regions also mediate functionally advantageous assembly into high-order structures that promote the formation of membrane-less sub-cellular compartments and even hydrogels. Across the range of functional assembly states accessed by disordered regions, post-translational modifications and regulatory macromolecular interactions, which can also be investigated by NMR spectroscopy, feature prominently. Here we will explore the many ways in which NMR has advanced our understanding of the physical-chemical phase space occupied by disordered protein regions and provide prospectus for the future role of NMR in this emerging and exciting field.

Computational and theoretical advances in studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are increasingly realized to play diverse biological roles, ranging from molecular signaling to the formation of membraneless organelles. Their high degree of disorder makes them more challenging to study using the techniques of conventional structural biology, because any observable will be averaged over a heterogeneous ensemble of structures. Molecular simulations and theory are therefore a natural complement to experiment for studying the structure, dynamics and function of IDPs. The diverse time and length scales relevant to the roles played by IDPs require flexibility in the techniques applied. Here, I summarize some of the developments in simulation and theory in recent years, which have been driven by the desire to better capture IDP properties at different time- and length-scales. I also provide an outlook for how methods can be improved in the future and emerging problems which may be addressed by theory and simulation.