Comprehensive analysis of NMR data using advanced line shape fitting

Structural basis of substrate recognition and allosteric activation of the proapoptotic mitochondrial HtrA2 protease

The mitochondrial serine protease HtrA2 is a human homolog of the Escherichia coli Deg-proteins exhibiting chaperone and proteolytic roles. HtrA2 is involved in both apoptotic regulation via its ability to degrade inhibitor-of-apoptosis proteins (IAPs), as well as in cellular maintenance as part of the cellular protein quality control machinery, by preventing the possible toxic accumulation of aggregated proteins. In this study, we use advanced solution NMR spectroscopy methods combined with biophysical characterization and biochemical assays to elucidate the crucial role of the substrate recognizing PDZ domain. This domain regulates the protease activity of HtrA2 by triggering an intricate allosteric network involving the regulatory loops of the protease domain. We further show that divalent metal ions can both positively and negatively modulate the activity of HtrA2, leading to a refined model of HtrA2 regulation within the apoptotic pathway.

The Listeria monocytogenes persistence factor ClpL is a potent stand-alone disaggregase

Heat stress can cause cell death by triggering the aggregation of essential proteins. In bacteria, aggregated proteins are rescued by the canonical Hsp70/AAA+ (ClpB) bi-chaperone disaggregase. Man-made, severe stress conditions applied during, e.g., food processing represent a novel threat for bacteria by exceeding the capacity of the Hsp70/ClpB system. Here, we report on the potent autonomous AAA+ disaggregase ClpL from Listeria monocytogenes that provides enhanced heat resistance to the food-borne pathogen enabling persistence in adverse environments. ClpL shows increased thermal stability and enhanced disaggregation power compared to Hsp70/ClpB, enabling it to withstand severe heat stress and to solubilize tight aggregates. ClpL binds to protein aggregates via aromatic residues present in its N-terminal domain (NTD) that adopts a partially folded and dynamic conformation. Target specificity is achieved by simultaneous interactions of multiple NTDs with the aggregate surface. ClpL shows remarkable structural plasticity by forming diverse higher assembly states through interacting ClpL rings. NTDs become largely sequestered upon ClpL ring interactions. Stabilizing ring assemblies by engineered disulfide bonds strongly reduces disaggregation activity, suggesting that they represent storage states.

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

Molecular switching in transcription through splicing and proline-isomerization regulates stress responses in plants

The Arabidopsis thaliana DREB2A transcription factor interacts with the negative regulator RCD1 and the ACID domain of subunit 25 of the transcriptional co-regulator mediator (Med25) to integrate stress signals for gene expression, with elusive molecular interplay. Using biophysical and structural analyses together with high-throughput screening, we reveal a bivalent binding switch in DREB2A containing an ACID-binding motif (ABS) and the known RCD1-binding motif (RIM). The RIM is lacking in a stress-induced DREB2A splice variant with retained transcriptional activity. ABS and RIM bind to separate sites on Med25-ACID, and NMR analyses show a structurally heterogeneous complex deriving from a DREB2A-ABS proline residue populating cis- and trans-isomers with remote impact on the RIM. The cis-isomer stabilizes an α-helix, while the trans-isomer may introduce energetic frustration facilitating rapid exchange between activators and repressors. Thus, DREB2A uses a post-transcriptionally and post-translationally modulated switch for transcriptional regulation.

Optimal control pulses for the 1.2-GHz (28.2-T) NMR spectrometers

The ability to measure nuclear magnetic resonance (NMR) spectra with a large sample volume is crucial for concentration-limited biological samples to attain adequate signal-to-noise (S/N) ratio. The possibility to measure with a 5-mm cryoprobe is currently absent at the 1.2-GHz NMR instruments due to the exceedingly high radio frequency power demands, which is four times compared to 600-MHz instruments. Here, we overcome the high-power demands by designing optimal control (OC) pulses with up to 20 times lower power requirements than currently necessary at a 1.2-GHz spectrometer. We show that multidimensional biomolecular NMR experiments constructed using these OC pulses can bestow improvement in the S/N ratio of up to 26%. With the expected power limitations of a 5-mm cryoprobe, we observe an enhancement in the S/N ratio of more than 240% using our OC sequences. This motivates the development of a cryoprobe with a larger volume than the current 3-mm cryoprobes.

The Prenucleation Equilibrium of the Parathyroid Hormone Determines the Critical Aggregation Concentration and Amyloid Fibril Nucleation

Nucleation and growth of amyloid fibrils were found to only occur in supersaturated solutions above a critical concentration (ccrit ). The biophysical meaning of ccrit remained mostly obscure, since typical low values of ccrit in the sub-μM range hamper investigations of potential oligomeric states and their structure. Here, we investigate the parathyroid hormone PTH84 as an example of a functional amyloid fibril forming peptide with a comparably high ccrit of 67±21 μM. We describe a complex concentration dependent prenucleation ensemble of oligomers of different sizes and secondary structure compositions and highlight the occurrence of a trimer and tetramer at ccrit as possible precursors for primary fibril nucleation. Furthermore, the soluble state found in equilibrium with fibrils adopts to the prenucleation state present at ccrit . Our study sheds light onto early events of amyloid formation directly related to the critical concentration and underlines oligomer formation as a key feature of fibril nucleation. Our results contribute to a deeper understanding of the determinants of supersaturated peptide solutions. In the current study we present a biophysical approach to investigate ccrit of amyloid fibril formation of PTH84 in terms of secondary structure, cluster size and residue resolved intermolecular interactions during oligomer formation. Throughout the investigated range of concentrations (1 μM to 500 μM) we found different states of oligomerization with varying ability to contribute to primary fibril nucleation and with a concentration dependent equilibrium. In this context, we identified the previously described ccrit of PTH84 to mark a minimum concentration for the formation of homo-trimers/tetramers. These investigations allowed us to characterize molecular interactions of various oligomeric states that are further converted into elongation competent fibril nuclei during the lag phase of a functional amyloid forming peptide.

Barriers to resolution in 1H NMR of rotating solids

The role of 1H solid-state NMR in structure elucidation of solids is becoming more preponderant, particularly as faster magic-angle spinning rates (MAS) become available which improve 1H detected assignment strategies. However, current 1H spectral resolution is still relatively poor, with linewidths of typically a few hundred Hz, even at the fastest rates available today. Here we detail and assess the factors limiting proton linewidths and line shapes in MAS experiments with five different samples, exemplifying the different sources of broadening that affect the residual linewidth. We disentangle the different contributions through one- and two-dimensional experiments: by using dilution to identify the contribution of ABMS; by using extensive deuteration to identify the dipolar contributions; and by using variable MAS rates to determine the ratio between homogeneous and inhomogeneous components. We find that the overall widths and the nature of the different contributions to the linewidths can vary very considerably. While we find that faster spinning always yields narrower lines and longer coherence lifetimes, we also find that for some resonances the dipolar contribution is no longer dominant at 100 kHz MAS. When the inhomogeneous sources of broadening, such as ABMS and chemical shift disorder, are dominant, two-dimensional 1H-1H correlation experiments yield better resolution for assignment. Particularly the extraction of the antidiagonal of a 2D peak will remove any correlated inhomogeneous broadening, giving substantially narrower 1H linewidths.

Structure and dynamics of the mitochondrial DNA-compaction factor Abf2 from S. cerevisiae

Mitochondria are essential organelles that produce most of the energy via the oxidative phosphorylation (OXPHOS) system in all eukaryotic cells. Several essential subunits of the OXPHOS system are encoded by the mitochondrial genome (mtDNA) despite its small size. Defects in mtDNA maintenance and expression can lead to severe OXPHOS deficiencies and are amongst the most common causes of mitochondrial disease. The mtDNA is packaged as nucleoprotein structures, referred to as nucleoids. The conserved mitochondrial proteins, ARS-binding factor 2 (Abf2) in the budding yeast Saccharomyces cerevisiae (S. cerevisiae) and mitochondrial transcription factor A (TFAM) in mammals, are nucleoid associated proteins (NAPs) acting as condensing factors needed for packaging and maintenance of the mtDNA. Interestingly, gene knockout of Abf2 leads, in yeast, to the loss of mtDNA and respiratory growth, providing evidence for its crucial role. On a structural level, the condensing factors usually contain two DNA binding domains called high-mobility group boxes (HMG boxes). The co-operating mechanical activities of these domains are crucial in establishing the nucleoid architecture by bending the DNA strands. Here we used advanced solution NMR spectroscopy methods to characterize the dynamical properties of Abf2 together with its interaction with DNA. We find that the two HMG-domains react notably different to external cues like temperature and salt, indicating distinct functional properties. Biophysical characterizations show the pronounced difference of these domains upon DNA-binding, suggesting a refined picture of the Abf2 functional cycle.

Applications of Solution NMR Spectroscopy in Quality Assessment and Authentication of Bovine Milk

Nuclear magnetic resonance (NMR) spectroscopy is emerging as a promising technique for the analysis of bovine milk, primarily due to its non-destructive nature, minimal sample preparation requirements, and comprehensive approach to untargeted milk analysis. These inherent strengths of NMR make it a formidable complementary tool to mass spectrometry-based techniques in milk metabolomic studies. This review aims to provide a comprehensive overview of the applications of NMR techniques in the quality assessment and authentication of bovine milk. It will focus on the experimental setup and data processing techniques that contribute to achieving accurate and highly reproducible results. The review will also highlight key studies that have utilized commonly used NMR methodologies in milk analysis, covering a wide range of application fields. These applications include determining milk animal species and feeding regimes, as well as assessing milk nutritional quality and authenticity. By providing an overview of the diverse applications of NMR in milk analysis, this review aims to demonstrate the versatility and significance of NMR spectroscopy as an invaluable tool for milk and dairy metabolomics research and hence, for assessing the quality and authenticity of bovine milk.

Arginine anchor points govern H3 tail dynamics

Chromatin is dynamically reorganized spatially and temporally, and the post-translational modification of histones is a key component of this regulation. The basic subunit of chromatin is the nucleosome core particle, consisting of two copies each of the histones H2A, H2B, H3, and H4 around which ∼147 base pairs of DNA wrap. The intrinsically disordered histone termini, or tails, protrude from the core and are heavily post-translationally modified. Previous studies have shown that the histone tails exist in dynamic ensembles of DNA-bound states within the nucleosome. Histone tail interactions with DNA are involved in nucleosome conformation and chromatin organization. Charge-modulating histone post-translational modifications (PTMs) are poised to perturb the dynamic interactions between histone tails and DNA. Arginine side chains form favorable interactions with DNA and are sites of charge-modulating PTMs such as citrullination. Our current focus is on the H3 tail, the longest histone tail. Four arginine residues are relatively evenly spaced along the H3 tail sequence, suggesting multivalent interactions with DNA poised for regulation by PTMs. In this study, we use NMR nuclear spin relaxation experiments to investigate the contribution of arginine residues to H3 tail dynamics within the nucleosome core particle. By neutralizing arginine via mutation to glutamine, we begin to work towards a comprehensive understanding of the contribution of individual residues to H3 tail dynamics. We find that neutralization of arginine residues results in increased regional mobility of the H3 tails, with implications for understanding the direct effects of arginine citrullination. Altogether, these studies support a role for dynamics within the histone language and emphasize the importance of charge-modulating histone PTMs in regulating chromatin dynamics, starting at the level of the basic subunit of chromatin.

Morphology-Dependent Interactions between α-Synuclein Monomers and Fibrils

Amyloid fibrils may adopt different morphologies depending on the solution conditions and the protein sequence. Here, we show that two chemically identical but morphologically distinct α-synuclein fibrils can form under identical conditions. This was observed by nuclear magnetic resonance (NMR), circular dichroism (CD), and fluorescence spectroscopy, as well as by cryo-transmission electron microscopy (cryo-TEM). The results show different surface properties of the two morphologies, A and B. NMR measurements show that monomers interact differently with the different fibril surfaces. Only a small part of the N-terminus of the monomer interacts with the fibril surface of morphology A, compared to a larger part of the monomer for morphology B. Differences in ThT binding seen by fluorescence titrations, and mesoscopic structures seen by cryo-TEM, support the conclusion of the two morphologies having different surface properties. Fibrils of morphology B were found to have lower solubility than A. This indicates that fibrils of morphology B are thermodynamically more stable, implying a chemical potential of fibrils of morphology B that is lower than that of morphology A. Consequently, at prolonged incubation time, fibrils of morphology B remained B, while an initially monomorphic sample of morphology A gradually transformed to B.

Upstream of N-Ras C-terminal cold shock domains mediate poly(A) specificity in a novel RNA recognition mode and bind poly(A) binding protein

RNA binding proteins (RBPs) often engage multiple RNA binding domains (RBDs) to increase target specificity and affinity. However, the complexity of target recognition of multiple RBDs remains largely unexplored. Here we use Upstream of N-Ras (Unr), a multidomain RBP, to demonstrate how multiple RBDs orchestrate target specificity. A crystal structure of the three C-terminal RNA binding cold-shock domains (CSD) of Unr bound to a poly(A) sequence exemplifies how recognition goes beyond the classical ππ-stacking in CSDs. Further structural studies reveal several interaction surfaces between the N-terminal and C-terminal part of Unr with the poly(A)-binding protein (pAbp). All interactions are validated by mutational analyses and the high-resolution structures presented here will guide further studies to understand how both proteins act together in cellular processes.

Cosolvent Dimethyl Sulfoxide Influences Protein-Ligand Binding Kinetics via Solvent Viscosity Effects: Revealing the Success Rate of Complex Formation Following Diffusive Protein-Ligand Encounter

Protein-ligand-exchange kinetics determines the duration of biochemical signals and consequently plays an important role in drug design. Binding studies commonly require solubilization of designed ligands in solvents such as dimethyl sulfoxide (DMSO), resulting in residual amounts of DMSO following titration of solubilized ligands into aqueous protein samples. Therefore, it is critical to establish whether DMSO influences protein-ligand binding. Here, we address the general and indirect effect of DMSO on protein-ligand binding caused by solvent viscosity, which is strongly dependent on the relative concentrations of DMSO and water. As a model system, we studied the binding of a drug-like ligand to the carbohydrate recognition domain of galectin-3 in the presence of variable amounts of DMSO. We used isothermal titration calorimetry to characterize binding thermodynamics and ¹⁵N NMR relaxation to monitor kinetics. The binding enthalpy is not affected, but we observe a subtle trend of increasingly unfavorable entropy of binding, and consequently decreased affinity, with increasing DMSO concentration. The increasing concentration of DMSO results in a reduced association rate of binding, while the dissociation rate is less affected. The observed association rate is inversely proportional to the viscosity of the DMSO-water mixture, as expected from theory, but significantly reduced from the diffusion-controlled limit. By comparing the viscosity dependence of the observed association rate with that of the theoretical diffusion-controlled association rate, we estimate the success rate of productive complex formation following an initial encounter of proteins and ligands, showing that only one out of several hundred binding "attempts" are successful.

Structure and dynamics of the von Willebrand Factor C6 domain

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.

Chaperoning of the histone octamer by the acidic domain of DNA repair factor APLF

Nucleosome assembly requires the coordinated deposition of histone complexes H3-H4 and H2A-H2B to form a histone octamer on DNA. In the current paradigm, specific histone chaperones guide the deposition of first H3-H4 and then H2A-H2B. Here, we show that the acidic domain of DNA repair factor APLF (APLF^AD) can assemble the histone octamer in a single step and deposit it on DNA to form nucleosomes. The crystal structure of the APLF^AD-histone octamer complex shows that APLF^AD tethers the histones in their nucleosomal conformation. Mutations of key aromatic anchor residues in APLF^AD affect chaperone activity in vitro and in cells. Together, we propose that chaperoning of the histone octamer is a mechanism for histone chaperone function at sites where chromatin is temporarily disrupted.

COLMARq: A Web Server for 2D NMR Peak Picking and Quantitative Comparative Analysis of Cohorts of Metabolomics Samples

Highly quantitative metabolomics studies of complex biological mixtures are facilitated by the resolution enhancement afforded by 2D NMR spectra such as 2D ¹³C-¹H HSQC spectra. Here, we describe a new public web server, COLMARq, for the semi-automated analysis of sets of 2D HSQC spectra of cohorts of samples. The workflow of COLMARq includes automated peak picking using the deep neural network DEEP Picker, quantitative cross-peak volume extraction by numerical fitting using Voigt Fitter, the matching of corresponding cross-peaks across cohorts of spectra, peak volume normalization between different spectra, database query for metabolite identification, and basic univariate and multivariate statistical analyses of the results. COLMARq allows the analysis of cross-peaks that belong to both known and unknown metabolites. After a user has uploaded cohorts of 2D ¹³C-¹H HSQC and optionally 2D ¹H-¹H TOCSY spectra in their preferred format, all subsequent steps on the web server can be performed fully automatically, allowing manual editing if needed and the sessions can be saved for later use. The accuracy, versatility, and interactive nature of COLMARq enables quantitative metabolomics analysis, including biomarker identification, of a broad range of complex biological mixtures as is illustrated for cohorts of samples from bacterial cultures of Pseudomonas aeruginosa in both its biofilm and planktonic states.

Observation of Arginine Side-Chain Motions Coupled to the Global Conformational Exchange Process in Deubiquitinase A

Coupled motions have been demonstrated to be functionally important in a number of enzymes. Noncovalent side-chain interactions play essential roles in coordinating the motions across different structural elements in a protein. However, most of the dynamic studies of proteins are focused on backbone amides or methyl groups in the side chains and little is known about the polar and charged side chains. We have previously characterized the conformational dynamics of deubiquitinase A (DUBA), an isopeptidase, on the microsecond-to-millisecond (μs-ms) time scales with the amide ¹H Carr-Purcell-Meiboom-Gill (CPMG) experiment. We detected a global conformational exchange process on a time scale of approximately 200 μs, which involves most of the structural elements in DUBA, including the active site and the substrate binding interface. Here, we extend our previous study on backbone amides to the arginine side-chain N^ε-H^ε groups using a modified ¹H CPMG pulse sequence that can efficiently detect both backbone amide and arginine side-chain N^ε-H^ε signals in a single experiment. We found that the side chains of three arginines display motions on the same time scale as the backbone amides. Mutations of two of the three arginines to alanines result in a decrease in enzyme activity. One of these two arginines is located in a loop involved in substrate binding. This loop is not visible in the backbone amide-detected experiments due to excess line broadening induced by motions on the μs-ms time scales. These results clearly demonstrate that the motions of some arginine side chains are coupled to the global conformational exchange process and provide an additional probe for motions in a functionally important loop that did not yield visible backbone amide signals, suggesting the value of side-chain experiments on DUBA. The modified ¹H CPMG pulse sequence allows the simultaneous characterization of backbone and arginine side-chain dynamics without any increase in data acquisition time and can be applied to the dynamic studies of any protein that displays measurable amide ¹H relaxation dispersion.

Structure determination of high-energy states in a dynamic protein ensemble

Macromolecular function frequently requires that proteins change conformation into high-energy states^1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts⁵ (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion⁶ (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

Structural Analysis of the Menangle Virus P Protein Reveals a Soft Boundary between Ordered and Disordered Regions

The paramyxoviral phosphoprotein (P protein) is the non-catalytic subunit of the viral RNA polymerase, and coordinates many of the molecular interactions required for RNA synthesis. All paramyxoviral P proteins oligomerize via a centrally located coiled-coil that is connected to a downstream binding domain by a dynamic linker. The C-terminal region of the P protein coordinates interactions between the catalytic subunit of the polymerase, and the viral nucleocapsid housing the genomic RNA. The inherent flexibility of the linker is believed to facilitate polymerase translocation. Here we report biophysical and structural characterization of the C-terminal region of the P protein from Menangle virus (MenV), a bat-borne paramyxovirus with zoonotic potential. The MenV P protein is tetrameric but can dissociate into dimers at sub-micromolar protein concentrations. The linker is globally disordered and can be modeled effectively as a worm-like chain. However, NMR analysis suggests very weak local preferences for alpha-helical and extended beta conformation exist within the linker. At the interface between the disordered linker and the structured C-terminal binding domain, a gradual disorder-to-order transition occurs, with X-ray crystallographic analysis revealing a dynamic interfacial structure that wraps the surface of the binding domain.

Structure and dynamics of the quaternary hunchback mRNA translation repression complex

A key regulatory process during Drosophila development is the localized suppression of the hunchback mRNA translation at the posterior, which gives rise to a hunchback gradient governing the formation of the anterior-posterior body axis. This suppression is achieved by a concerted action of Brain Tumour (Brat), Pumilio (Pum) and Nanos. Each protein is necessary for proper Drosophila development. The RNA contacts have been elucidated for the proteins individually in several atomic-resolution structures. However, the interplay of all three proteins during RNA suppression remains a long-standing open question. Here, we characterize the quaternary complex of the RNA-binding domains of Brat, Pum and Nanos with hunchback mRNA by combining NMR spectroscopy, SANS/SAXS, XL/MS with MD simulations and ITC assays. The quaternary hunchback mRNA suppression complex comprising the RNA binding domains is flexible with unoccupied nucleotides functioning as a flexible linker between the Brat and Pum-Nanos moieties of the complex. Moreover, the presence of the Pum-HD/Nanos-ZnF complex has no effect on the equilibrium RNA binding affinity of the Brat RNA binding domain. This is in accordance with previous studies, which showed that Brat can suppress mRNA independently and is distributed uniformly throughout the embryo.

Rapid measurement of heteronuclear transverse relaxation rates using non-uniformly sampled R1ρ accordion experiments

Multidimensional, heteronuclear NMR relaxation methods are used extensively to characterize the dynamics of biological macromolecules. Acquisition of relaxation datasets on proteins typically requires significant measurement time, often several days. Accordion spectroscopy offers a powerful means to shorten relaxation rate measurements by encoding the "relaxation dimension" into the indirect evolution period in multidimensional experiments. Time savings can also be achieved by non-uniform sampling (NUS) of multidimensional NMR data, which is used increasingly to improve spectral resolution or increase sensitivity per unit time. However, NUS is not commonly implemented in relaxation experiments, because most reconstruction algorithms are inherently nonlinear, leading to problems when estimating signal intensities, relaxation rate constants and their error bounds. We have previously shown how to avoid these shortcomings by combining accordion spectroscopy with NUS, followed by data reconstruction using sparse exponential mode analysis, thereby achieving a dramatic decrease in the total length of longitudinal relaxation experiments. Here, we present the corresponding transverse relaxation experiment, taking into account the special considerations required for its successful implementation in the framework of the accordion-NUS approach. We attain the highest possible precision in the relaxation rate constants by optimizing the NUS scheme with respect to the Cramér-Rao lower bound of the variance of the estimated parameter, given the total number of sampling points and the spectrum-specific signal characteristics. The resulting accordion-NUS R 1 ρ relaxation experiment achieves comparable precision in the parameter estimates compared to conventional CPMG (Carr-Purcell-Meiboom-Gill) R 2 or spin-lock R 1 ρ experiments while saving an order of magnitude in experiment time.

15N CPMG Relaxation Dispersion for the Investigation of Protein Conformational Dynamics on the µs-ms Timescale

Protein conformational dynamics play fundamental roles in regulation of enzymatic catalysis, ligand binding, allostery, and signaling, which are important biological processes. Understanding how the balance between structure and dynamics governs biological function is a new frontier in modern structural biology and has ignited several technical and methodological developments. Among these, CPMG relaxation dispersion solution NMR methods provide unique, atomic-resolution information on the structure, kinetics, and thermodynamics of protein conformational equilibria occurring on the µs-ms timescale. Here, the study presents detailed protocols for acquisition and analysis of a ¹⁵N relaxation dispersion experiment. As an example, the pipeline for the analysis of the µs-ms dynamics in the C-terminal domain of bacteria Enzyme I is shown.

Cholesterol-containing lipid nanodiscs promote an α-synuclein binding mode that accelerates oligomerization

Dysregulation of the biosynthesis of cholesterol and other lipids has been implicated in many neurological diseases, including Parkinson's disease. Misfolding of α-synuclein (α-Syn), the main actor in Parkinson's disease, is associated with changes in a lipid environment. However, the exact molecular mechanisms underlying cholesterol effect on α-Syn binding to lipids as well as α-Syn oligomerization and fibrillation remain elusive, as does the relative importance of cholesterol compared to other factors. We probed the interactions and fibrillation behaviour of α-Syn using styrene-maleic acid nanodiscs, containing zwitterionic and anionic lipid model systems with and without cholesterol. Surface plasmon resonance and thioflavin T fluorescence assays were employed to monitor α-Syn binding, as well as fibrillation in the absence and presence of membrane models. 1 H-15 N-correlated NMR was used to monitor the fold of α-Syn in response to nanodisc binding, determining individual residue apparent affinities for the nanodisc-contained bilayers. The addition of cholesterol inhibited α-Syn interaction with lipid bilayers and, however, significantly promoted α-Syn fibrillation, with a more than a 20-fold reduction of lag times before fibrillation onset. When α-Syn bilayer interactions were analysed at an individual residue level by solution-state NMR, we observed two different effects of cholesterol. In nanodiscs made of DOPC, the addition of cholesterol modulated the NAC part of α-Syn, leading to stronger interaction of this region with the lipid bilayer. In contrast, in the nanodiscs comprising DOPC, DOPE and DOPG, the NAC part was mostly unaffected by the presence of cholesterol, while the binding of the N and the C termini was both inhibited.

Conformational Dynamics of Deubiquitinase A and Functional Implications

Deubiquitinase A (DUBA) belongs to the ovarian tumor family of deubiquitinating enzymes and was initially identified as a negative regulator of type I interferons, whose overproduction has been linked to autoimmune diseases. The deubiquitinating activity of DUBA is positively regulated by phosphorylation at a single serine residue, S177, which results in minimal structural changes. We have previously shown that phosphorylation induces a two-state conformational equilibrium observed only in the active form of DUBA, highlighting the functional importance of DUBA dynamics. Here, we report the conformational dynamics of DUBA on the microsecond-to-millisecond time scales characterized by nuclear magnetic resonance relaxation dispersion experiments. We found that motions on these time scales are highly synchronized in the phosphorylated and nonphosphorylated DUBA. Despite the overall similarity of these two forms, different dynamic properties were observed in helix α1 and the neighboring regions, including residue S177, which likely contribute to the activation of DUBA by phosphorylation. Moreover, our data suggest that transient unfolding of helix α6 drives the global conformational process and that mutations can be introduced to modulate this process, which provides a basis for future studies to define the exact functional roles of motions in DUBA activation and substrate specificity.

Gain-of-Function Variant p.Pro2555Arg of von Willebrand Factor Increases Aggregate Size through Altering Stem Dynamics

The multimeric plasma glycoprotein (GP) von Willebrand factor (VWF) is best known for recruiting platelets to sites of injury during primary hemostasis. Generally, mutations in the VWF gene lead to loss of hemostatic activity and thus the bleeding disorder von Willebrand disease. By employing cone and platelet aggregometry and microfluidic assays, we uncovered a platelet GPIIb/IIIa-dependent prothrombotic gain of function (GOF) for variant p.Pro2555Arg, located in the C4 domain, leading to an increase in platelet aggregate size. We performed complementary biophysical and structural investigations using circular dichroism spectra, small-angle X-ray scattering, nuclear magnetic resonance spectroscopy, molecular dynamics simulations on the single C4 domain, and dimeric wild-type and p.Pro2555Arg constructs. C4-p.Pro2555Arg retained the overall structural conformation with minor populations of alternative conformations exhibiting increased hinge flexibility and slow conformational exchange. The dimeric protein becomes disordered and more flexible. Our data suggest that the GOF does not affect the binding affinity of the C4 domain for GPIIb/IIIa. Instead, the increased VWF dimer flexibility enhances temporal accessibility of platelet-binding sites. Using an interdisciplinary approach, we revealed that p.Pro2555Arg is the first VWF variant, which increases platelet aggregate size and shows a shear-dependent function of the VWF stem region, which can become hyperactive through mutations. Prothrombotic GOF variants of VWF are a novel concept of a VWF-associated pathomechanism of thromboembolic events, which is of general interest to vascular health but not yet considered in diagnostics. Thus, awareness should be raised for the risk they pose. Furthermore, our data implicate the C4 domain as a novel antithrombotic drug target.

How directed evolution reshapes the energy landscape in an enzyme to boost catalysis

The advent of biocatalysts designed computationally and optimized by laboratory evolution provides an opportunity to explore molecular strategies for augmenting catalytic function. Applying a suite of nuclear magnetic resonance, crystallography, and stopped-flow techniques to an enzyme designed for an elementary proton transfer reaction, we show how directed evolution gradually altered the conformational ensemble of the protein scaffold to populate a narrow, highly active conformational ensemble and accelerate this transformation by nearly nine orders of magnitude. Mutations acquired during optimization enabled global conformational changes, including high-energy backbone rearrangements, that cooperatively organized the catalytic base and oxyanion stabilizer, thus perfecting transition-state stabilization. The development of protein catalysts for many chemical transformations could be facilitated by explicitly sampling conformational substates during design and specifically stabilizing productive substates over all unproductive conformations.

NMR fragment screening reveals a novel small molecule binding site near the catalytic surface of the disulfide-dithiol oxidoreductase enzyme DsbA from Burkholderia pseudomallei

The presence of suitable cavities or pockets on protein structures is a general criterion for a therapeutic target protein to be classified as 'druggable'. Many disease-related proteins that function solely through protein-protein interactions lack such pockets, making development of inhibitors by traditional small-molecule structure-based design methods much more challenging. The 22 kDa bacterial thiol oxidoreductase enzyme, DsbA, from the gram-negative bacterium Burkholderia pseudomallei (BpsDsbA) is an example of one such target. The crystal structure of oxidized BpsDsbA lacks well-defined surface pockets. BpsDsbA is required for the correct folding of numerous virulence factors in B. pseudomallei, and genetic deletion of dsbA significantly attenuates B. pseudomallei virulence in murine infection models. Therefore, BpsDsbA is potentially an attractive drug target. Herein we report the identification of a small molecule binding site adjacent to the catalytic site of oxidized BpsDsbA. ¹H^N CPMG relaxation dispersion NMR measurements suggest that the binding site is formed transiently through protein dynamics. Using fragment-based screening, we identified a small molecule that binds at this site with an estimated affinity of K_D ~ 500 µM. This fragment inhibits BpsDsbA enzymatic activity in vitro. The binding mode of this molecule has been characterized by NMR data-driven docking using HADDOCK. These data provide a starting point towards the design of more potent small molecule inhibitors of BpsDsbA.

Resolving overlapped signals with automated FitNMR analytical peak modeling

Nuclear magnetic resonance (NMR) is a valuable tool for determining the structures of molecules and probing their dynamics. A longstanding problem facing both small-molecule and macromolecular NMR is overlapped signals in crowded spectra. To address this, we have developed a method that extracts peak features by fitting analytically derived models of NMR lineshapes. The approach takes into account the effects of truncation, apodization, and the resulting artifacts, while avoiding systematic errors that have affected other models. Even severely overlapped peaks, beyond the point of coalescence, can be distinguished in both simulated and experimental data. We show that the method can measure unresolved backbone scalar couplings directly from a 2D proton-nitrogen spectrum of a de novo designed mini protein. The algorithm is implemented in the FitNMR open-source R package and can be used to analyze nearly any type of single or multidimensional data from small molecules or biomolecules.

Pseudo-RNA-Binding Domains Mediate RNA Structure Specificity in Upstream of N-Ras

RNA-binding proteins (RBPs) commonly feature multiple RNA-binding domains (RBDs), which provide these proteins with a modular architecture. Accumulating evidence supports that RBP architectural modularity and adaptability define the specificity of their interactions with RNA. However, how multiple RBDs recognize their cognate single-stranded RNA (ssRNA) sequences in concert remains poorly understood. Here, we use Upstream of N-Ras (Unr) as a model system to address this question. Although reported to contain five ssRNA-binding cold-shock domains (CSDs), we demonstrate that Unr includes an additional four CSDs that do not bind RNA (pseudo-RBDs) but are involved in mediating RNA tertiary structure specificity by reducing the conformational heterogeneity of Unr. Disrupting the interactions between canonical and non-canonical CSDs impacts RNA binding, Unr-mediated translation regulation, and the Unr-dependent RNA interactome. Taken together, our studies reveal a new paradigm in protein-RNA recognition, where interactions between RBDs and pseudo-RBDs select RNA tertiary structures, influence RNP assembly, and define target specificity.

NMR insights into the pre-amyloid ensemble and secretion targeting of the curli subunit CsgA

The biofilms of Enterobacteriaceae are fortified by assembly of curli amyloid fibres on the cell surface. Curli not only provides structural reinforcement, but also facilitates surface adhesion. To prevent toxic intracellular accumulation of amyloid precipitate, secretion of the major curli subunit, CsgA, is tightly regulated. In this work, we have employed solution state NMR spectroscopy to characterise the structural ensemble of the pre-fibrillar state of CsgA within the bacterial periplasm, and upon recruitment to the curli pore, CsgG, and the secretion chaperone, CsgE. We show that the N-terminal targeting sequence (N) of CsgA binds specifically to CsgG and that its subsequent sequestration induces a marked transition in the conformational ensemble, which is coupled to a preference for CsgE binding. These observations lead us to suggest a sequential model for binding and structural rearrangement of CsgA at the periplasmic face of the secretion machinery.

Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional ¹⁹F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional ¹⁹F spectra to those from two-dimensional ¹H-¹⁵N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional ¹⁹F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of ¹⁹F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.

Methyl TROSY spectroscopy: A versatile NMR approach to study challenging biological systems

A major goal in structural biology is to unravel how molecular machines function in detail. To that end, solution-state NMR spectroscopy is ideally suited as it is able to study biological assemblies in a near natural environment. Based on methyl TROSY methods, it is now possible to record high-quality data on complexes that are far over 100 kDa in molecular weight. In this review, we discuss the theoretical background of methyl TROSY spectroscopy, the information that can be extracted from methyl TROSY spectra and approaches that can be used to assign methyl resonances in large complexes. In addition, we touch upon insights that have been obtained for a number of challenging biological systems, including the 20S proteasome, the RNA exosome, molecular chaperones and G-protein-coupled receptors. We anticipate that methyl TROSY methods will be increasingly important in modern structural biology approaches, where information regarding static structures is complemented with insights into conformational changes and dynamic intermolecular interactions.

Structure, dynamics and roX2-lncRNA binding of tandem double-stranded RNA binding domains dsRBD1,2 of Drosophila helicase Maleless

Maleless (MLE) is an evolutionary conserved member of the DExH family of helicases in Drosophila. Besides its function in RNA editing and presumably siRNA processing, MLE is best known for its role in remodelling non-coding roX RNA in the context of X chromosome dosage compensation in male flies. MLE and its human orthologue, DHX9 contain two tandem double-stranded RNA binding domains (dsRBDs) located at the N-terminal region. The two dsRBDs are essential for localization of MLE at the X-territory and it is presumed that this involves binding roX secondary structures. However, for dsRBD1 roX RNA binding has so far not been described. Here, we determined the solution NMR structure of dsRBD1 and dsRBD2 of MLE in tandem and investigated its role in double-stranded RNA (dsRNA) binding. Our NMR and SAXS data show that both dsRBDs act as independent structural modules in solution and are canonical, non-sequence-specific dsRBDs featuring non-canonical KKxAXK RNA binding motifs. NMR titrations combined with filter binding experiments and isothermal titration calorimetry (ITC) document the contribution of dsRBD1 to dsRNA binding in vitro. Curiously, dsRBD1 mutants in which dsRNA binding in vitro is strongly compromised do not affect roX2 RNA binding and MLE localization in cells. These data suggest alternative functions for dsRBD1 in vivo.

Intrinsically disordered protein as carbon nanotube dispersant: How dynamic interactions lead to excellent colloidal stability

The rich pool of protein conformations combined with the dimensions and properties of carbon nanotubes create new possibilities in functional materials and nanomedicine. Here, the intrinsically disordered protein α-synuclein is explored as a dispersant of single-walled carbon nanotubes (SWNTs) in water. We use a range of spectroscopic methods to quantify the amount of dispersed SWNT and to elucidate the binding mode of α-synuclein to SWNT. The dispersion ability of α-synuclein is good even with mild sonication and the obtained dispersion is very stable over time. The whole polypeptide chain is involved in the interaction accompanied by a fraction of the chain changing into a helical structure upon binding. Similar to other dispersants, we observe that only a small fraction (15-20%) of α-synuclein is adsorbed on the SWNT surface with an average residence time below 10 ms.

Extreme Nonuniform Sampling for Protein NMR Dynamics Studies in Minimal Time

NMR spectroscopy is an extraordinarily rich source of quantitative dynamics of proteins in solution using spin relaxation or chemical exchange saturation transfer (CEST) experiments. However, ¹⁵N-CEST measurements require prolonged multidimensional, so-called pseudo-3D HSQC experiments where the pseudo dimension is a radio frequency offset Δω of a weak ¹⁵N saturation field. Nonuniform sampling (NUS) approaches have the potential to significantly speed up these measurements, but they also carry the risk of introducing serious artifacts and the systematic optimization of nonuniform sampling schedules has remained elusive. It is demonstrated here how this challenge can be addressed by using fitted cross-peaks of a reference 2D HSQC experiment as footprints, which are subsequently used to reconstruct cross-peak amplitudes of a pseudo-3D data set as a function of Δω by a linear least-squares fit. It is shown for protein Im7 how the approach can yield highly accurate CEST profiles based on an absolutely minimally sampled (AMS) data set allowing a speed-up of a factor 20-30. Spectrum-specific optimized nonuniform sampling (SONUS) schemes based on the Cramer-Rao lower bound metric were critical to achieve such a performance, revealing also more general properties of optimal sampling schedules. This is the first systematic exploration and optimization of NUS schedules for the dramatic speed-up of quantitative multidimensional NMR measurements that minimize unwanted errors.

Switching the Post-translational Modification of Translation Elongation Factor EF-P

Tripeptides with two consecutive prolines are the shortest and most frequent sequences causing ribosome stalling. The bacterial translation elongation factor P (EF-P) relieves this arrest, allowing protein biosynthesis to continue. A seven amino acids long loop between beta-strands β3/β4 is crucial for EF-P function and modified at its tip by lysylation of lysine or rhamnosylation of arginine. Phylogenetic analyses unveiled an invariant proline in the -2 position of the modification site in EF-Ps that utilize lysine modifications such as Escherichia coli. Bacteria with the arginine modification like Pseudomonas putida on the contrary have selected against it. Focusing on the EF-Ps from these two model organisms we demonstrate the importance of the β3/β4 loop composition for functionalization by chemically distinct modifications. Ultimately, we show that only two amino acid changes in E. coli EF-P are needed for switching the modification strategy from lysylation to rhamnosylation.

NMR chemical shift assignment of the C-terminal region of the Menangle virus phosphoprotein

Menangle virus is a bat-borne paramyxovirus with zoonotic potential. The single-stranded RNA genome of the virus is encapsidated in a helical nucleocapsid which is the template for both transcription and genome replication. Each of these operations is performed by the viral RNA polymerase. The phosphoprotein is the non-catalytic subunit of the polymerase, and its C-terminal region enables the polymerase to engage with the nucleocapsid. Here, we report the ¹H, ¹⁵N, and ¹³C chemical shift assignments of the C-terminal region (amino acids 267-388) of the Menangle virus phosphoprotein. This region has a bipartite character, with a highly flexible and structurally disordered sequence preceding a structured nucleocapsid-binding domain. NMR chemical shift assignment will enable the detailed characterization of the dynamic behavior of the phosphoprotein, and its functional linkage with polymerase translocation.

Time-domain signal modelling in multidimensional NMR experiments for estimation of relaxation parameters

We present a model-based method for estimation of relaxation parameters from time-domain NMR data specifically suitable for processing data in popular 2D phase-sensitive experiments. Our model is formulated in terms of commutative bicomplex algebra, which allows us to use the complete information available in an NMR signal acquired with principles of quadrature detection without disregarding any of its dimensions. Compared to the traditional intensity-analysis method, our model-based approach offers an important advantage for the analysis of overlapping peaks and is robust over a wide range of signal-to-noise ratios. We assess its performance with simulated experiments and then apply it for determination of [Formula: see text], [Formula: see text], and [Formula: see text] relaxation rates in datasets of a protein with more than 100 cross peaks.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

Rescue of conformational dynamics in enzyme catalysis by directed evolution

Rational design and directed evolution have proved to be successful approaches to increase catalytic efficiencies of both natural and artificial enzymes. Protein dynamics is recognized as important, but due to the inherent flexibility of biological macromolecules it is often difficult to distinguish which conformational changes are directly related to function. Here, we use directed evolution on an impaired mutant of the proline isomerase CypA and identify two second-shell mutations that partially restore its catalytic activity. We show both kinetically, using NMR spectroscopy, and structurally, by room-temperature X-ray crystallography, how local perturbations propagate through a large allosteric network to facilitate conformational dynamics. The increased catalysis selected for in the evolutionary screen is correlated with an accelerated interconversion between the two catalytically essential conformational sub-states, which are both captured in the high-resolution X-ray ensembles. Our data provide a glimpse of an evolutionary trajectory and show how subtle changes can fine-tune enzyme function.