ShereKhan--calculating exchange parameters in relaxation dispersion data from CPMG experiments

Modulation of Primary and Secondary Processes in Tau Fibril Formation by Salt-Induced Dynamics

The initial stages of amyloid fibrilization begin with the monomers populating aggregation-prone conformers. Characterization of such aggregation-prone conformers is crucial in the study of neurodegenerative diseases. The current study characterizes the aggregation pathway of two tau protein constructs that have been recently demonstrated to form Alzheimer's (AD) fibril structures with divalent ions and chronic traumatic encephalopathy (CTE) fibril structures with monovalent ions. The results highlight the involvement of identical residues in both the primary and secondary processes of both AD and CTE fibril propagation. Nuclear magnetic resonance relaxation experiments reveal increased flexibility of the motifs 321KCGS within R3 and 364PGGGN within R4 in the presence of MgCl2/NaCl, correlating with faster aggregation kinetics and indicating efficient primary nucleation. Notably, the seeded aggregation kinetics of the tau monomers in the presence and absence of metal ions are strikingly different. This correlates with the overall sign of the 15N-ΔR2 profile specifying the dominant mechanism involved in the process of aggregation.

Epigenetic CpG duplex marks probed by an evolved DNA reader via a well-tempered conformational plasticity

5-methylcytosine (mC) and its TET-oxidized derivatives exist in CpG dyads of mammalian DNA and regulate cell fate, but how their individual combinations in the two strands of a CpG act as distinct regulatory signals is poorly understood. Readers that selectively recognize such novel 'CpG duplex marks' could be versatile tools for studying their biological functions, but their design represents an unprecedented selectivity challenge. By mutational studies, NMR relaxation, and MD simulations, we here show that the selectivity of the first designer reader for an oxidized CpG duplex mark hinges on precisely tempered conformational plasticity of the scaffold adopted during directed evolution. Our observations reveal the critical aspect of defined motional features in this novel reader for affinity and specificity in the DNA/protein interaction, providing unexpected prospects for further design progress in this novel area of DNA recognition.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Conformational Dynamics of the Hepatitis B Virus Pre-genomic RNA on Multiple Time Scales: Implications for Viral Replication

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication.

Identification of Structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

Structural Insights into the Unique Modes of Relaxin-Binding and Tethered-Agonist Mediated Activation of RXFP1 and RXFP2

Our poor understanding of the mechanism by which the peptide-hormone H2 relaxin activates its G protein coupled receptor, RXFP1 and the related receptor RXFP2, has hindered progress in its therapeutic development. Both receptors possess large ectodomains, which bind H2 relaxin, and contain an N-terminal LDLa module that is essential for receptor signaling and postulated to be a tethered agonist. Here, we show that a conserved motif (GDxxGWxxxF), C-terminal to the LDLa module, is critical for receptor activity. Importantly, this motif adopts different structures in RXFP1 and RXFP2, suggesting distinct activation mechanisms. For RXFP1, the motif is flexible, weakly associates with the LDLa module, and requires H2 relaxin binding to stabilize an active conformation. Conversely, the GDxxGWxxxF motif in RXFP2 is more closely associated with the LDLa module, forming an essential binding interface for H2 relaxin. These differences in the activation mechanism will aid drug development targeting these receptors.

RING NMR dynamics: software for analysis of multiple NMR relaxation experiments

Molecular motions are fundamental to the existence of life, and NMR spectroscopy remains one of the most useful and powerful methods to measure their rates and molecular characteristics. Multiple experimental methods are available for measuring the NMR relaxation properties and these can require different methods for extracting model parameters. We present here a new software application, RING NMR Dynamics, that is designed to support analysis of multiple relaxation types. The initial release of RING NMR Dynamics supports the analysis of exponential decay experiments such as T₁ and T₂, as well as CEST and R₂ and R_1ρ relaxation dispersion. The software runs on multiple operating systems in both a command line mode and a user-friendly GUI that allows visualizing and simulating relaxation data. Interaction with another program, NMRFx Analyst, allows drilling down from the derived relaxation parameters to the raw spectral data.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe-S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol-1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

19F NMR reveals the conformational properties of free thrombin and its zymogen precursor prethrombin-2

The conformational properties of trypsin-like proteases and their zymogen forms remain controversial because of a lack of sufficient information on their free forms. Specifically, it is unclear whether the free protease is zymogen-like and shifts to its mature form upon a ligand-induced fit or exists in multiple conformations in equilibrium from which the ligand selects the optimal fit via conformational selection. Here we report the results of ¹⁹F NMR measurements that reveal the conformational properties of a protease and its zymogen precursor in the free form. Using the trypsin-like, clotting protease thrombin as a relevant model system, we show that its conformation is quite different from that of its direct zymogen precursor prethrombin-2 and more similar to that of its fully active Na⁺-bound form. The results cast doubts on recent hypotheses that free thrombin is zymogen-like and transitions to protease-like forms upon ligand binding. Rather, they validate the scenario emerged from previous findings of X-ray crystallography and rapid kinetics supporting a pre-existing equilibrium between open (E) and closed (E*) forms of the active site. In this scenario, prethrombin-2 is more dynamic and exists predominantly in the E* form, whereas thrombin is more rigid and exists predominantly in the E form. Ligand binding to thrombin takes place exclusively in the E form without significant changes in the overall conformation. In summary, these results disclose the structural architecture of the free forms of thrombin and prethrombin-2, consistent with an E*-E equilibrium and providing no evidence that free thrombin is zymogen-like.

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.

Determinants of oligosaccharide specificity of the carbohydrate-binding modules of AMP-activated protein kinase

We solved the structures of β1- and β2-carbohydrate-binding modules (CBMs) of AMP-activated protein kinase (AMPK) bound to a branched carbohydrate. The additional threonine within the β2-module allows it to bind single α1,6-branched carbohydrates, such as partially degraded glycogen, with greater affinity.

An exact solution for R2,eff in CPMG experiments in the case of two site chemical exchange

The Carr-Purcell-Meiboom-Gill (CPMG) experiment is widely used to quantitatively analyse the effects of chemical exchange on NMR spectra. In a CPMG experiment, the effective transverse relaxation rate, R2,eff, is typically measured as a function of the pulse frequency, νCPMG. Here, an exact expression for how R2,eff varies with νCPMG is derived for the commonly encountered scenario of two-site chemical exchange of in-phase magnetisation. This result, summarised in Appendix A, generalises a frequently used equation derived by Carver and Richards, published in 1972. The expression enables more rapid analysis of CPMG data by both speeding up calculation of R2,eff over numerical methods by a factor of ca. 130, and yields exact derivatives for use in data analysis. Moreover, the derivation provides insight into the physical principles behind the experiment.

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

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

Enhanced accuracy of kinetic information from CT-CPMG experiments by transverse rotating-frame spectroscopy

Micro-to-millisecond motions of proteins transmit pivotal signals for protein function. A powerful technique for the measurement of these motions is nuclear magnetic resonance spectroscopy. One of the most widely used methodologies for this purpose is the constant-time Carr-Purcell-Meiboom-Gill (CT-CPMG) relaxation dispersion experiment where kinetic and structural information can be obtained at atomic resolution. Extraction of accurate kinetics determined from CT-CPMG data requires refocusing frequencies that are much larger than the nuclei's exchange rate between states. We investigated the effect when fast processes are probed by CT-CPMG experiments via simulation and show that if the intrinsic relaxation rate (R(CT-CPMG)(2,0)) is not known a priori the extraction of accurate kinetics is hindered. Errors on the order of 50 % in the exchange rate are attained when processes become fast, but are minimized to 5 % with a priori (CT-CPMG)(2,0)) information. To alleviate this shortcoming, we developed an experimental scheme probing (CT-CPMG)(2,0)) with large amplitude spin-lock fields, which specifically contains the intrinsic proton longitudinal Eigenrelaxation rate. Our approach was validated with ubiquitin and the Oscillatoria agardhii agglutinin (OAA). For OAA, an underestimation of 66 % in the kinetic rates was observed if (CT-CPMG)(2,0)) is not included during the analysis of CT-CPMG data and result in incorrect kinetics and imprecise amplitude information. This was overcome by combining CT-CPMG with (CT-CPMG)(2,0)) measured with a high power R1ρ experiment. In addition, the measurement of (CT-CPMG)(2,0)) removes the ambiguities in choosing between different models that describe CT-CPMG data.