CheSPI: chemical shift secondary structure population inference

Dynamic networks connect the USP14 active site region with the proteasome interaction surface

Ubiquitin-specific protease 14 (USP14) is a member of the USP family responsible for the catalytic removal of ubiquitin (Ub) from proteins directed to the proteasome, implicated in the pathogenesis of neurodegeneration and cancer. Crystallography and cryo-EM analysis have identified loop regions crucial for the deubiquitinase activity of USP14, specifically those involved in Ub and proteasome binding. However, the structural changes in USP14 upon ligand binding to these regions are minimal, indicating significant yet uncharacterized dynamic contributions to its function. In this study, through structural and dynamical NMR experiments and functional evaluation, we demonstrate that small mutations designed to impact Ub binding and catalytic activity without disturbing the USP structure display both local and long-range effects. The affected residues connect the catalytic site and the Ub binding region with the proteasome interaction surface through a network of loops, which show varied dynamics on the ps-ms time scale. Collectively, our findings experimentally reveal different aspects of dynamic connections within USP14, suggesting the presence of allosteric networks that link enzyme activity with regulatory function. The identification of coupled clusters of possible allostery participants in the free USP domain provides new insights into the dynamic regulation of USP14, with potential implications for understanding its role in cellular processes.

Assignment of the N-terminal domain of mouse cGAS

Cyclic GMP-AMP synthase (cGAS) is a DNA-sensing enzyme that is a member of the nucleotidyltransferase (NTase) family and functions as a DNA sensor. The protein is comprised of a catalytic NTase core domain and an unstructured hypervariable N-terminal domain (NTD) that was reported to increase protein activity by providing an additional DNA-binding surface. We report nearly complete 1H, 15N, and 13C backbone chemical-shift assignments of mouse cGAS NTD (residues 5-146), obtained with a set of 3D and 4D solution NMR experiments. Analysis of the chemical-shift values confirms that the NTD is intrinsically disordered. These resonance assignments can provide the basis for further studies such as activation by DNA and protein-protein interactions.

Structural Plasticity as a Driver of the Maturation of Pro-Interleukin-18

Dynamics are often critical for biomolecular function. Herein we explore the role of motion in driving the maturation process of pro-IL-18, a potent pro-inflammatory cytokine that is cleaved by caspases-1 and -4 to generate the mature form of the protein. An NMR dynamics study of pro-IL-18, probing time scales over 12 orders of magnitude and focusing on 1H, 13C, and 15N spin probes along the protein backbone and amino-acid side chains, reveals a plastic structure, with millisecond time scale dynamics occurring in a pair of β-strands, β1 and β*, that show large structural variations in a comparison of caspase-free and bound pro-IL-18 states. Fits of the relaxation data to a three-site model of exchange showed that the ground state secondary structure is maintained in the excited conformers, with the side chain of I48 that undergoes a buried-to-exposed conformational change in the caspase-free to -bound transition of pro-IL-18, sampling a more extensive range of torsion angles in one of the excited states characterized, suggesting partial unpacking in this region. Hydrogen exchange measurements establish the occurrence of an additional process, whereby strands β1 and β* locally unfold. Our data are consistent with a hierarchy of dynamic events that likely prime pro-IL-18 for facile caspase binding.

Transient Structural Properties of the Rho GDP-Dissociation Inhibitor

Rho GTPases, master spatial regulators of a wide range of cellular processes, are orchestrated by complex formation with guanine nucleotide dissociation inhibitors (RhoGDIs). These have been thought to possess an unstructured N-terminus that inhibits nucleotide exchange of their client upon binding/folding. Via NMR analyses, molecular dynamics simulations, and biochemical assays, we reveal instead pertinent structural properties transiently maintained both, in the presence and absence of the client, imposed onto the terminus context-specifically by modulating interactions with the surface of the folded C-terminal domain. These observations revise the long-standing textbook picture of the GTPases' mechanism of membrane extraction. Rather than by a disorder-to-order transition upon binding of an inhibitory peptide, the intricate and highly selective extraction process of RhoGTPases is orchestrated via a dynamic ensemble bearing preformed transient structural properties, suitably modulated by the specific surrounding along the multi-step process.

Control of G protein-coupled receptor function via membrane-interacting intrinsically disordered C-terminal domains

G protein-coupled receptors (GPCRs) control intracellular signaling cascades via agonist-dependent coupling to intracellular transducers including heterotrimeric G proteins, GPCR kinases (GRKs), and arrestins. In addition to their critical interactions with the transmembrane core of active GPCRs, all three classes of transducers have also been reported to interact with receptor C-terminal domains (CTDs). An underexplored aspect of GPCR CTDs is their possible role as lipid sensors given their proximity to the membrane. CTD-membrane interactions have the potential to control the accessibility of key regulatory CTD residues to downstream effectors and transducers. Here, we report that the CTDs of two closely related family C GPCRs, metabotropic glutamate receptor 2 (mGluR2) and mGluR3, bind to membranes and that this interaction can regulate receptor function. We first characterize CTD structure with NMR spectroscopy, revealing lipid composition-dependent modes of membrane binding. Using molecular dynamics simulations and structure-guided mutagenesis, we then identify key conserved residues and cancer-associated mutations that modulate CTD-membrane binding. Finally, we provide evidence that mGluR3 transducer coupling is controlled by CTD-membrane interactions in live cells, which may be subject to regulation by CTD phosphorylation and changes in membrane composition. This work reveals an additional mechanism of GPCR modulation, suggesting that CTD-membrane binding may be a general regulatory mode throughout the broad GPCR superfamily.

Transient interdomain interactions in free USP14 shape its conformational ensemble

The deubiquitinase (DUB) ubiquitin-specific protease 14 (USP14) is a dual domain protein that plays a regulatory role in proteasomal degradation and has been identified as a promising therapeutic target. USP14 comprises a conserved USP domain and a ubiquitin-like (Ubl) domain separated by a 25-residue linker. The enzyme activity of USP14 is autoinhibited in solution, but is enhanced when bound to the proteasome, where the Ubl and USP domains of USP14 bind to the Rpn1 and Rpt1/Rpt2 units, respectively. No structure of full-length USP14 in the absence of proteasome has yet been presented, however, earlier work has described how transient interactions between Ubl and USP domains in USP4 and USP7 regulate DUB activity. To better understand the roles of the Ubl and USP domains in USP14, we studied the Ubl domain alone and in full-length USP14 by nuclear magnetic resonance spectroscopy and used small angle x-ray scattering and molecular modeling to visualize the entire USP14 protein ensemble. Jointly, our results show how transient interdomain interactions between the Ubl and USP domains of USP14 predispose its conformational ensemble for proteasome binding, which may have functional implications for proteasome regulation and may be exploited in the design of future USP14 inhibitors.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

The pathogen-encoded signalling receptor Tir exploits host-like intrinsic disorder for infection

The translocated intimin receptor (Tir) is an essential type III secretion system (T3SS) effector of attaching and effacing pathogens contributing to the global foodborne disease burden. Tir acts as a cell-surface receptor in host cells, rewiring intracellular processes by targeting multiple host proteins. We investigated the molecular basis for Tir's binding diversity in signalling, finding that Tir is a disordered protein with host-like binding motifs. Unexpectedly, also are several other T3SS effectors. By an integrative approach, we reveal that Tir dimerises via an antiparallel OB-fold within a highly disordered N-terminal cytosolic domain. Also, it has a long disordered C-terminal cytosolic domain partially structured at host-like motifs that bind lipids. Membrane affinity depends on lipid composition and phosphorylation, highlighting a previously unrecognised host interaction impacting Tir-induced actin polymerisation and cell death. Furthermore, multi-site tyrosine phosphorylation enables Tir to engage host SH2 domains in a multivalent fuzzy complex, consistent with Tir's scaffolding role and binding promiscuity. Our findings provide insights into the intracellular Tir domains, highlighting the ability of T3SS effectors to exploit host-like protein disorder as a strategy for host evasion.

Disordered-to-ordered transitions in assembly factors allow the complex II catalytic subunit to switch binding partners

Complex II (CII) activity controls phenomena that require crosstalk between metabolism and signaling, including neurodegeneration, cancer metabolism, immune activation, and ischemia-reperfusion injury. CII activity can be regulated at the level of assembly, a process that leverages metastable assembly intermediates. The nature of these intermediates and how CII subunits transfer between metastable complexes remains unclear. In this work, we identify metastable species containing the SDHA subunit and its assembly factors, and we assign a preferred temporal sequence of appearance of these species during CII assembly. Structures of two species show that the assembly factors undergo disordered-to-ordered transitions without the appearance of significant secondary structure. The findings identify that intrinsically disordered regions are critical in regulating CII assembly, an observation that has implications for the control of assembly in other biomolecular complexes.

Backbone 1H, 15N and 13C resonance assignments of the 27kDa fluorescent protein mCherry

mCherry is one of the most successfully applied monomeric red fluorescent proteins (RFPs) for in vivo and in vitro imaging. However, questions pertaining to the photostability of the RFPs remain and rational further engineering of their photostability requires information about the fluorescence quenching mechanism in solution. To this end, NMR spectroscopic investigations might be helpful, and we present the near-complete backbone NMR chemical shift assignment to aid in this pursuit.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

SAM1 domain of SASH1 harbors distinctive structural heterogeneity

The sterile alpha motif (SAM) domains are among the most versatile protein domains in biology, and the variety of the oligomerization states contribute to their diverse roles in many diseases. A better understanding of the structure and dynamics of various SAM domains will provide a scientific basis for drug development targeting them. Here, we used SEC-MALS, HPLC, NMR, and other biophysical techniques to characterize the structural features and dynamics of the SAM1 domain in SASH1. SASH1 is a scaffold protein belonging to the same family as SASH3. Unlike the dimerization seen in SASH3's SAM domain, our SEC-MALS and SE-HPLC showed that SAM1 exists primarily as a less compact monomer with a minor oligomer. NMR assignment, relaxation, and exchange experiments revealed the presence of both a disordered monomer and a more structured oligomer with multiple timescale exchange regimes in solution. Mutagenesis and SE-HPLC showed that D663A/T664K substitutions in SAM1 increased its oligomerization. In sum, this study is the first to characterize a disordered structure for a SAM domain, provides additional evidence and framework for the diversity of SAM domains, and identifies a region in SAM1 as a potential starting point to further characterize the structural mechanism of oligomerization of the domain.

A glutamine-based single α-helix scaffold to target globular proteins

The binding of intrinsically disordered proteins to globular ones can require the folding of motifs into α-helices. These interactions offer opportunities for therapeutic intervention but their modulation with small molecules is challenging because they bury large surfaces. Linear peptides that display the residues that are key for binding can be targeted to globular proteins when they form stable helices, which in most cases requires their chemical modification. Here we present rules to design peptides that fold into single α-helices by instead concatenating glutamine side chain to main chain hydrogen bonds recently discovered in polyglutamine helices. The resulting peptides are uncharged, contain only natural amino acids, and their sequences can be optimized to interact with specific targets. Our results provide design rules to obtain single α-helices for a wide range of applications in protein engineering and drug design.

NMR resonance assignment of the N-terminal GTPase domain of human Miro2 Bound to GTP

Miro2 and Miro1 are mitochondrial-associated proteins critical for regulating mitochondrial movement within the cell. Both Miro1 and Miro2 have roles in promoting neuron function, but recently Miro2 has been shown to have additional roles in response to nutrient starvation in tumor cells. Miro1 and 2 consist of two small GTPase domains flanking a pair of EF-hands. The N-terminal GTPase (nGTPase) domain is responsible for initiating mitochondrial trafficking and interactions with GCN1 in prostate cancer. The crystal structure of Miro1 nGTPase bound to GTP has been solved. However, no structural data is available for the nGTPase domain of Miro2. To better understand the similarities and differences in the functions of Miro1 and Miro2, we have initiated structural studies of Miro2. Here we report the backbone NMR chemical shift assignments of a 22 KDa construct of the nGTPase domain of Miro2 bound to GTP that includes residues 1-180 of the full-length protein. We affirm that the overall secondary structure of this complex closely resembles that of Miro1 nGTPase bound to GTP. Minor variations in the overall structures can be attributed to crystal packing interactions in the structure of Miro1. These NMR studies will form the foundation for future work identifying the specific interaction sites between Miro2 and its cellular binding partners.

The Ni(II)-Binding Activity of the Intrinsically Disordered Region of Human NDRG1, a Protein Involved in Cancer Development

Nickel exposure is associated with tumors of the respiratory tract such as lung and nasal cancers, acting through still-uncharacterized mechanisms. Understanding the molecular basis of nickel-induced carcinogenesis requires unraveling the mode and the effects of Ni(II) binding to its intracellular targets. A possible Ni(II)-binding protein and a potential focus for cancer treatment is hNDRG1, a protein induced by Ni(II) through the hypoxia response pathway, whose expression correlates with higher cancer aggressiveness and resistance to chemotherapy in lung tissue. The protein sequence contains a unique C-terminal sequence of 83 residues (hNDRG1*C), featuring a three-times-repeated decapeptide, involved in metal binding, lipid interaction and post-translational phosphorylation. In the present work, the biochemical and biophysical characterization of unmodified hNDRG1*C was performed. Bioinformatic analysis assigned it to the family of the intrinsically disordered regions and the absence of secondary and tertiary structure was experimentally proven by circular dichroism and NMR. Isothermal titration calorimetry revealed the occurrence of a Ni(II)-binding event with micromolar affinity. Detailed information on the Ni(II)-binding site and on the residues involved was obtained in an extensive NMR study, revealing an octahedral paramagnetic metal coordination that does not cause any major change of the protein backbone, which is coherent with CD analysis. hNDRG1*C was found in a monomeric form by light-scattering experiments, while the full-length hNDRG1 monomer was found in equilibrium between the dimer and tetramer, both in solution and in human cell lines. The results are the first essential step for understanding the cellular function of hNDRG1*C at the molecular level, with potential future applications to clarify its role and the role of Ni(II) in cancer development.

IDPConformerGenerator: A Flexible Software Suite for Sampling the Conformational Space of Disordered Protein States

The power of structural information for informing biological mechanisms is clear for stable folded macromolecules, but similar structure-function insight is more difficult to obtain for highly dynamic systems such as intrinsically disordered proteins (IDPs) which must be described as structural ensembles. Here, we present IDPConformerGenerator, a flexible, modular open-source software platform for generating large and diverse ensembles of disordered protein states that builds conformers that obey geometric, steric, and other physical restraints on the input sequence. IDPConformerGenerator samples backbone phi (φ), psi (ψ), and omega (ω) torsion angles of relevant sequence fragments from loops and secondary structure elements extracted from folded protein structures in the RCSB Protein Data Bank and builds side chains from robust Monte Carlo algorithms using expanded rotamer libraries. IDPConformerGenerator has many user-defined options enabling variable fractional sampling of secondary structures, supports Bayesian models for assessing the agreement of IDP ensembles for consistency with experimental data, and introduces a machine learning approach to transform between internal and Cartesian coordinates with reduced error. IDPConformerGenerator will facilitate the characterization of disordered proteins to ultimately provide structural insights into these states that have key biological functions.

Solution NMR backbone assignments of disordered Olduvai protein domain CON1 employing Hα-detected experiments

Olduvai protein domains, encoded by the NBPF gene family, are responsible for the largest increase in copy number of any protein-coding region in the human genome. This has spawned various genetics studies which have linked these domains to human brain development and divergence from our primate ancestors, as well as currently relevant cognitive diseases such as schizophrenia and autism spectrum disorder (ASD). There are six separate Olduvai domains which together form the majority of the various protein products of the NBPF genes. The six domains include three conserved domains (CON1-3), and three human-lineage-specific domains (HLS1-3) which occur in triplet. Here, we present the solution nuclear magnetic resonance backbone assignments for the CON1 domain, which has been linked to the severity of ASD. The data confirm that CON1 is an intrinsically disordered protein (IDP). Additionally, we use innovative Hα-detected experiments which allow us to not only assign the Hα atoms and N atoms of proline residues, but also to assign residues where HN-experiments suffered from peak overlap or broadening.