Mass Spectrometry of RNA-Binding Proteins during Liquid-Liquid Phase Separation Reveals Distinct Assembly Mechanisms and Droplet Architectures.

Liquid-liquid phase separation (LLPS) of heterogeneous ribonucleoproteins (hnRNPs) drives the formation of membraneless organelles, but structural information about their assembled states is still lacking. Here, we address this challenge through a combination of protein engineering, native ion mobility mass spectrometry, and molecular dynamics simulations. We used an LLPS-compatible spider silk domain and pH changes to control the self-assembly of the hnRNPs FUS, TDP-43, and hCPEB3, which are implicated in neurodegeneration, cancer, and memory storage. By releasing the proteins inside the mass spectrometer from their native assemblies, we could monitor conformational changes associated with liquid-liquid phase separation. We find that FUS monomers undergo an unfolded-to-globular transition, whereas TDP-43 oligomerizes into partially disordered dimers and trimers. hCPEB3, on the other hand, remains fully disordered with a preference for fibrillar aggregation over LLPS. The divergent assembly mechanisms revealed by ion mobility mass spectrometry of soluble protein species that exist under LLPS conditions suggest structurally distinct complexes inside liquid droplets that may impact RNA processing and translation depending on biological context.

Double Mutant of Chymotrypsin Inhibitor 2 Stabilized through Increased Conformational Entropy.

The conformational heterogeneity of a folded protein can affect not only its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here, we have examined whether changed native-state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side chains with molecular dynamics simulations to quantify the native-state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: a reduction in conformational dynamics (and entropy estimated from this) of the native state of the L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native-state entropy might be needed to improve the predictions of the effect of mutations on protein stability.

On the specificity of protein-protein interactions in the context of disorder.

With the increased focus on intrinsically disordered proteins (IDPs) and their large interactomes, the question about their specificity - or more so on their multispecificity - arise. Here we recapitulate how specificity and multispecificity are quantified and address through examples if IDPs in this respect differ from globular proteins. The conclusion is that quantitatively, globular proteins and IDPs are similar when it comes to specificity. However, compared with globular proteins, IDPs have larger interactome sizes, a phenomenon that is further enabled by their flexibility, repetitive binding motifs and propensity to adapt to different binding partners. For IDPs, this adaptability, interactome size and a higher degree of multivalency opens for new interaction mechanisms such as facilitated exchange through trimer formation and ultra-sensitivity via threshold effects and ensemble redistribution. IDPs and their interactions, thus, do not compromise the definition of specificity. Instead, it is the sheer size of their interactomes that complicates its calculation. More importantly, it is this size that challenges how we conceptually envision, interpret and speak about their specificity.

Charge Interactions in a Highly Charge-Depleted Protein.

Electrostatic forces are important for protein folding and are favored targets of protein engineering. However, interactions between charged residues are difficult to study because of the complex network of interactions found in most proteins. We have designed a purposely simple system to investigate this problem by systematically introducing individual and pairs of charged and titratable residues in a protein otherwise free of such residues. We used constant pH molecular dynamics simulations, NMR spectroscopy, and thermodynamic double mutant cycles to probe the structure and energetics of the interaction between the charged residues. We found that the partial burial of surface charges contributes to a shift in pKa value, causing an aspartate to titrate in the neutral pH range. Additionally, the interaction between pairs of residues was found to be highly context dependent, with some pairs having no apparent preferential interaction, while other pairs would engage in coupled titration forming a highly stabilized salt bridge. We find good agreement between experiments and simulations and use the simulations to rationalize our observations and to provide a detailed mechanistic understanding of the electrostatic interactions.

Software for reconstruction of nonuniformly sampled NMR data.

Nonuniform sampling (NUS) of multidimensional NMR experiments is a powerful tool to obtain high-resolution spectra with less instrument time. With NUS, only a subset of the points needed for conventional Fourier transformation is recorded, and sophisticated algorithms are needed to reconstruct the missing data points. During the last decade, several software packages implementing the reconstruction algorithms have emerged and been refined and now result in spectra of almost similar quality as spectra from conventionally recorded and processed data. However, from the number of literature references to the reconstruction methods, many more multidimensional NMR spectra could presumably be recorded with NUS. To help researchers considering to start using NUS for their NMR experiments, we here review 13 different reconstruction methods found in five software packages (CambridgeCS, hmsIST, MddNMR, NESTA-NMR, and SMILE). We have compared how the methods run with the provided example scripts for reconstructing a nonuniform sampled three-dimensional 15 N-NOESY-HSQC at sampling densities from 5% to 50%. Overall, the spectra are all of similar quality above 20% sampling density. Thus, without any particular knowledge about the details of the reconstruction algorithms, significant reduction in the experiment time can be achieved. Below 20% sampling density, the intensities of particular weak peaks start being affected. MddNMR's IST with virtual echo and the SMILE algorithms still reproduced the spectra with the highest accuracy of peak intensities.

The three-dimensional structure of an H-superfamily conotoxin reveals a granulin fold arising from a common ICK cysteine framework.

Venomous marine cone snails produce peptide toxins (conotoxins) that bind ion channels and receptors with high specificity and therefore are important pharmacological tools. Conotoxins contain conserved cysteine residues that form disulfide bonds that stabilize their structures. To gain structural insight into the large, yet poorly characterized conotoxin H-superfamily, we used NMR and CD spectroscopy along with MS-based analyses to investigate H-Vc7.2 from Conus victoriae, a peptide with a VI/VII cysteine framework. This framework has CysI-CysIV/CysII-CysV/CysIII-CysVI connectivities, which have invariably been associated with the inhibitor cystine knot (ICK) fold. However, the solution structure of recombinantly expressed and purified H-Vc7.2 revealed that although it displays the expected cysteine connectivities, H-Vc7.2 adopts a different fold consisting of two stacked ß-hairpins with opposing ß-strands connected by two parallel disulfide bonds, a structure homologous to the N-terminal region of the human granulin protein. Using structural comparisons, we subsequently identified several toxins and nontoxin proteins with this "mini-granulin" fold. These findings raise fundamental questions concerning sequence-structure relationships within peptides and proteins and the key determinants that specify a given fold.

Global analysis of protein stability by temperature and chemical denaturation.

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 µl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 µg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

Structure of the competence pilus major pilin ComGC in Streptococcus pneumoniae.

Type IV pili are important virulence factors on the surface of many pathogenic bacteria and have been implicated in a wide range of diverse functions, including attachment, twitching motility, biofilm formation, and horizontal gene transfer. The respiratory pathogen Streptococcus pneumoniae deploys type IV pili to take up DNA during transformation. These "competence pili" are composed of the major pilin protein ComGC and exclusively assembled during bacterial competence, but their biogenesis remains unclear. Here, we report the high resolution NMR structure of N-terminal truncated ComGC revealing a highly flexible and structurally divergent type IV pilin. It consists of only three α-helical segments forming a well-defined electronegative cavity and confined electronegative and hydrophobic patches. The structure is particularly flexible between the first and second α-helix with the first helical part exhibiting slightly slower dynamics than the rest of the pilin, suggesting that the first helix is involved in forming the pilus structure core and that parts of helices two and three are primarily surface-exposed. Taken together, our results provide the first structure of a type IV pilin protein involved in the formation of competence-induced pili in Gram-positive bacteria and corroborate the remarkable structural diversity among type IV pilin proteins.

Behaviour of intrinsically disordered proteins in protein-protein complexes with an emphasis on fuzziness.

Intrinsically disordered proteins (IDPs) do not, by themselves, fold into a compact globular structure. They are extremely dynamic and flexible, and are typically involved in signalling and transduction of information through binding to other macromolecules. The reason for their existence may lie in their malleability, which enables them to bind several different partners with high specificity. In addition, their interactions with other macromolecules can be regulated by a variable amount of chemically diverse post-translational modifications. Four kinetically and energetically different types of complexes between an IDP and another macromolecule are reviewed: (1) simple two-state binding involving a single binding site, (2) avidity, (3) allovalency and (4) fuzzy binding; the last three involving more than one site. Finally, a qualitative definition of fuzzy binding is suggested, examples are provided, and its distinction to allovalency and avidity is highlighted and discussed.

The Pathogenic A2V Mutant Exhibits Distinct Aggregation Kinetics, Metal Site Structure, and Metal Exchange of the Cu2+ -Aß Complex.

A prominent current hypothesis is that impaired metal ion homeostasis may contribute to Alzheimer's disease (AD). We elucidate the interaction of Cu2+ with wild-type (WT) Aß1-40 and the genetic variants A2T and A2V which display increasing pathogenicity as A2T

A Soluble, Folded Protein without Charged Amino Acid Residues.

Charges are considered an integral part of protein structure and function, enhancing solubility and providing specificity in molecular interactions. We wished to investigate whether charged amino acids are indeed required for protein biogenesis and whether a protein completely free of titratable side chains can maintain solubility, stability, and function. As a model, we used a cellulose-binding domain from Cellulomonas fimi, which, among proteins of more than 100 amino acids, presently is the least charged in the Protein Data Bank, with a total of only four titratable residues. We find that the protein shows a surprising resilience toward extremes of pH, demonstrating stability and function (cellulose binding) in the pH range from 2 to 11. To ask whether the four charged residues present were required for these properties of this protein, we altered them to nontitratable ones. Remarkably, this chargeless protein is produced reasonably well in Escherichia coli, retains its stable three-dimensional structure, and is still capable of strong cellulose binding. To further deprive this protein of charges, we removed the N-terminal charge by acetylation and studied the protein at pH 2, where the C-terminus is effectively protonated. Under these conditions, the protein retains its function and proved to be both soluble and have a reversible folding-unfolding transition. To the best of our knowledge, this is the first time a soluble, functional protein with no titratable side chains has been produced.

Amyloid-ß and α-Synuclein Decrease the Level of Metal-Catalyzed Reactive Oxygen Species by Radical Scavenging and Redox Silencing.

The formation of reactive oxygen species (ROS) is linked to the pathogenesis of neurodegenerative diseases. Here we have investigated the effect of soluble and aggregated amyloid-ß (Aß) and α-synuclein (αS), associated with Alzheimer's and Parkinson's diseases, respectively, on the Cu(2+)-catalyzed formation of ROS in vitro in the presence of a biological reductant. We find that the levels of ROS, and the rate by which ROS is generated, are significantly reduced when Cu(2+) is bound to Aß or αS, particularly when they are in their oligomeric or fibrillar forms. This effect is attributed to a combination of radical scavenging and redox silencing mechanisms. Our findings suggest that the increase in ROS associated with the accumulation of aggregated Aß or αS does not result from a particularly ROS-active form of these peptides, but rather from either a local increase of Cu(2+) and other ROS-active metal ions in the aggregates or as a downstream consequence of the formation of the pathological amyloid structures.

Non-uniform sampling of NMR relaxation data.

The use of non-uniform sampling of NMR spectra may give significant reductions in the data acquisition time. For quantitative experiments such as the measurement of spin relaxation rates, non-uniform sampling is however not widely used as inaccuracies in peak intensities may lead to errors in the extracted dynamic parameters. By systematic reducing the coverage of the Nyquist grid of (15)N Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion datasets for four different proteins and performing a full data analysis of the resulting non-uniform sampled datasets, we have compared the performance of the multi-dimensional decomposition and iterative re-weighted least-squares algorithms in reconstructing spectra with accurate peak intensities. As long as a single fully sampled spectrum is included in a series of otherwise non-uniform sampled two-dimensional spectra, multi-dimensional decomposition reconstructs the non-uniform sampled spectra with high accuracy. For two of the four analyzed datasets, a coverage of only 20% results in essentially the same results as the fully sampled data. As exemplified by other data, such a low coverage is in general not enough to produce reliable results. We find that a coverage level not compromising the final results can be estimated by recording a single full two-dimensional spectrum and reducing the spectrum quality in silico.

Protein interacting with C-kinase 1 (PICK1) binding promiscuity relies on unconventional PSD-95/discs-large/ZO-1 homology (PDZ) binding modes for nonclass II PDZ ligands.

PDZ domain proteins control multiple cellular functions by governing assembly of protein complexes. It remains unknown why individual PDZ domains can bind the extreme C terminus of very diverse binding partners and maintain selectivity. By employing NMR spectroscopy, together with molecular modeling, mutational analysis, and fluorescent polarization binding experiments, we identify here three structural mechanisms explaining why the PDZ domain of PICK1 selectively binds >30 receptors, transporters, and kinases. Class II ligands, including the dopamine transporter, adopt a canonical binding mode with promiscuity obtained via differential packing in the binding groove. Class I ligands, such as protein kinase Cα, depend on residues upstream from the canonical binding sequence that are likely to interact with flexible loop residues of the PDZ domain. Finally, we obtain evidence that the unconventional ligand ASIC1a has a dual binding mode involving a canonical insertion and a noncanonical internal insertion with the two C-terminal residues forming interactions outside the groove. Together with an evolutionary analysis, the data show how unconventional binding modes might evolve for a protein recognition domain to expand the repertoire of functionally important interactions.

Aggregation-Prone Amyloid-ß⋠Cu(II) Species Formed on the Millisecond Timescale under Mildly Acidic Conditions.

Metal ions and their interaction with the amyloid beta (Aß) peptide might be key elements in the development of Alzheimer's disease. In this work the effect of Cu(II) on the aggregation of Aß is explored on a timescale from milliseconds to days, both at physiological pH and under mildly acidic conditions, by using stopped-flow kinetic measurements (fluorescence and light-scattering), (1) H NMR relaxation and ThT fluorescence. A minimal reaction model that relates the initial Cu(II) binding and Aß folding with downstream aggregation is presented. We demonstrate that a highly aggregation prone Aß⋅Cu(II) species is formed on the sub-second timescale at mildly acidic pH. This observation might be central to the molecular origin of the known detrimental effect of acidosis in Alzheimer's disease.

relax: the analysis of biomolecular kinetics and thermodynamics using NMR relaxation dispersion data.

UNLABELLED: Nuclear magnetic resonance (NMR) is a powerful tool for observing the motion of biomolecules at the atomic level. One technique, the analysis of relaxation dispersion phenomenon, is highly suited for studying the kinetics and thermodynamics of biological processes. Built on top of the relax computational environment for NMR dynamics is a new dispersion analysis designed to be comprehensive, accurate and easy-to-use. The software supports more models, both numeric and analytic, than current solutions. An automated protocol, available for scripting and driving the graphical user interface (GUI), is designed to simplify the analysis of dispersion data for NMR spectroscopists. Decreases in optimization time are granted by parallelization for running on computer clusters and by skipping an initial grid search by using parameters from one solution as the starting point for another -using analytic model results for the numeric models, taking advantage of model nesting, and using averaged non-clustered results for the clustered analysis. AVAILABILITY AND IMPLEMENTATION: The software relax is written in Python with C modules and is released under the GPLv3+ license. Source code and precompiled binaries for all major operating systems are available from http://www.nmr-relax.com. CONTACT: edward@nmr-relax.com.

The pKa value and accessibility of cysteine residues are key determinants for protein substrate discrimination by glutaredoxin.

The enzyme glutaredoxin catalyzes glutathione exchange, but little is known about its interaction with protein substrates. Very different proteins are substrates in vitro, and the enzyme seems to have low requirements for specific protein interactions. Here we present a systematic investigation of the interaction between human glutaredoxin 1 and glutathionylated variants of a single model protein. Thus, single cysteine variants of acyl-coenzyme A binding protein were produced creating a set of substrates in the same protein background. The rate constants for deglutathionylation differ by more than 2 orders of magnitude between the best (k1 = 1.75 × 10(5) M(-1) s(-1)) and the worst substrate (k1 = 4 × 10(2) M(-1) s(-1)). The pKa values of the substrate cysteine residues were determined by NMR spectroscopy and found to vary from 8.2 to 9.9. Rates of glutaredoxin 1-catalyzed deglutathionylation were assessed with respect to substrate cysteine pKa values, cysteine residue accessibility, local stability, and backbone dynamics. Good substrates are characterized by a combination of high accessibility of the glutathionylated site and low pKa of the cysteine residue.

Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection.

Protein dynamics on the microsecond-millisecond time scales often play a critical role in biological function. NMR relaxation dispersion experiments are powerful approaches for investigating biologically relevant dynamics with site-specific resolution, as shown by a growing number of publications on enzyme catalysis, protein folding, ligand binding, and allostery. To date, the majority of studies has probed the backbone amides or side-chain methyl groups, while experiments targeting other sites have been used more sparingly. Aromatic side chains are useful probes of protein dynamics, because they are over-represented in protein binding interfaces, have important catalytic roles in enzymes, and form a sizable part of the protein interior. Here we present an off-resonance R 1ρ experiment for measuring microsecond to millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimization. Using selective excitation and inversion of the narrow component of the (13)C doublet, the experiment achieves significant sensitivity enhancement in terms of both signal intensity and the fractional contribution from exchange to transverse relaxation; additional signal enhancement is achieved by optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-selected R 1ρ experiment by measuring exchange parameters for Y23 in bovine pancreatic trypsin inhibitor at a temperature of 328 K, where the ring flip is in the fast exchange regime with a mean waiting time between flips of 320 µs. The determined chemical shift difference matches perfectly with that measured from the NMR spectrum at lower temperatures, where separate peaks are observed for the two sites. We further show that potentially complicating effects of strong scalar coupling between protons (Weininger et al. in J Phys Chem B 117: 9241-9247, 2013b) can be accounted for using a simple expression, and provide recommendations for data acquisition when the studied system exhibits this behavior. The present method extends the repertoire of relaxation methods tailored for aromatic side chains by enabling studies of faster processes and improved control over artifacts due to strong coupling.

Helical propensity in an intrinsically disordered protein accelerates ligand binding.

Many intrinsically disordered proteins fold upon binding to other macromolecules. The secondary structure present in the well-ordered complex is often formed transiently in the unbound state. The consequence of such transient structure for the binding process is, however, not clear. The activation domain of the activator for thyroid hormone and retinoid receptors (ACTR) is intrinsically disordered and folds upon binding to the nuclear coactivator binding domain (NCBD) of the CREB binding protein. A number of mutants was designed that selectively perturbs the amount of secondary structure in unbound ACTR without interfering with the intermolecular interactions between ACTR and NCBD. Using NMR spectroscopy and fluorescence-monitored stopped-flow kinetic measurements we show that the secondary structure content in helix 1 of ACTR indeed influences the binding kinetics. The results thus support the notion of preformed secondary structure as an important determinant for molecular recognition in intrinsically disordered proteins.