The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Quantitative analysis of the slow exchange process by 19F NMR in the presence of scalar and dipolar couplings: applications to the ribose 2'-19F probe in nucleic acids

Solution NMR spectroscopy is a particularly powerful technique for characterizing the functional dynamics of biomolecules, which is typically achieved through the quantitative characterization of chemical exchange processes via the measurement of spin relaxation rates. In addition to the conventional nuclei such as 15N and 13C, which are abundant in biomolecules, fluorine-19 (19F) has recently garnered attention and is being widely used as a site-specific spin probe. While 19F offers the advantages of high sensitivity and low background, it can be susceptible to artifacts in quantitative relaxation analyses due to a multitude of dipolar and scalar coupling interactions with nearby 1H spins. In this study, we focused on the ribose 2'-19F spin probe in nucleic acids and investigated the effects of 1H-19F spin interactions on the quantitative characterization of slow exchange processes on the millisecond time scale. We demonstrated that the 1H-19F dipolar coupling can significantly affect the interpretation of 19F chemical exchange saturation transfer (CEST) experiments when 1H decoupling is applied, while the 1H-19F interactions have a lesser impact on Carr-Purcell-Meiboom-Gill relaxation dispersion applications. We also proposed a modified CEST scheme to alleviate these artifacts along with experimental verifications on self-complementary RNA systems. The theoretical framework presented in this study can be widely applied to various 19F spin systems where 1H-19F interactions are operative, further expanding the utility of 19F relaxation-based NMR experiments.

Increasing the accuracy of exchange parameters reporting on slow dynamics by performing CEST experiments with 'high' B1 fields

Over the last decade chemical exchange saturation transfer (CEST) NMR methods have emerged as powerful tools to characterize biomolecular conformational dynamics occurring between a visible major state and 'invisible' minor states. The ability of the CEST experiment to detect these minor states, and provide precise exchange parameters, hinges on using appropriate B1 field strengths during the saturation period. Typically, a pair of B1 fields with ω1 (=2πB1) values around the exchange rate kex are chosen. Here we show that the transverse relaxation rate of the minor state resonance (R2,B) also plays a crucial role in determining the B1 fields that lead to the most informative datasets. Using [Formula: see text]  ≥ kex, to guide the choice of B1, instead of kex, leads to data wherefrom substantially more accurate exchange parameters can be derived. The need for higher B1 fields, guided by K, is demonstrated by studying the conformational exchange in two mutants of the 71 residue FF domain with kex ∼ 11 s-1 and ∼ 72 s-1, respectively. In both cases analysis of CEST datasets recorded using B1 field values guided by kex lead to imprecise exchange parameters, whereas using B1 values guided by K resulted in precise site-specific exchange parameters. The conclusions presented here will be valuable while using CEST to study slow processes at sites with large intrinsic relaxation rates, including carbonyl sites in small to medium sized proteins, amide 15N sites in large proteins and when the minor state dips are broadened due to exchange among the minor states.

Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( k e x , F I 1 ~ 550 s-1, k e x , F I 2 ~ 1200 s-1, k e x , I 1 I 2 ~ 5000 s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15 °C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Picosecond Dynamics of a Small Molecule in Its Bound State with an Intrinsically Disordered Protein

Intrinsically disordered proteins (IDPs) are highly dynamic biomolecules that rapidly interconvert among many structural conformations. These dynamic biomolecules are involved in cancers, neurodegeneration, cardiovascular illnesses, and viral infections. Despite their enormous therapeutic potential, IDPs have generally been considered undruggable because of their lack of classical long-lived binding pockets for small molecules. Currently, only a few instances are known where small molecules have been observed to interact with IDPs, and this situation is further exacerbated by the limited sensitivity of experimental techniques to detect such binding events. Here, using experimental nuclear magnetic resonance (NMR) spectroscopy 19F transverse spin-relaxation measurements, we discovered that a small molecule, 5-fluoroindole, interacts with the disordered domains of non-structural protein 5A from hepatitis C virus with a Kd of 260 ± 110 μM. Our analysis also allowed us to determine the rotational correlation times (τc) for the free and bound states of 5-fluoroindole. In the free state, we observed a rotational correlation time of 27.0 ± 1.3 ps, whereas in the bound state, τc only increased to 46 ± 10 ps. Our findings imply that it is possible for small molecules to engage with IDPs in exceptionally dynamic ways, in sharp contrast to the rigid binding modes typically exhibited when small molecules bind to well-defined binding pockets within structured proteins.

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

19F NMR Untersuchung des Konformationsaustauschs mehrerer Zustände im synthetischen Neomycin-bindenden Riboschalter

Der synthetische Neomycin‐bindende Riboschalter interagiert mit seinem Liganden Neomycin sowie mit den verwandten Antibiotika Ribostamycin und Paromomycin. Die Bindung dieser Aminoglykoside induziert sehr ähnliche Grundzustandsstrukturen in der RNA, allerdings kann nur Neomycin die Initiierung der Translation effizient unterdrücken. Der molekulare Ursprung dieser Unterschiede wurde auf Unterschiede in der Dynamik der Ligand‐Riboschalter‐Komplexe zurückgeführt. In diesem Artikel kombinieren wir fünf komplementäre fluorbasierte NMR‐Methoden, um die Dynamik der drei Riboschalter‐Komplexe im Sekunden‐ bis Mikrosekundenbereich genau zu quantifizieren. Unsere Daten offenbaren komplexe Austauschprozesse mit bis zu vier strukturell unterschiedlichen Zuständen. Wir interpretieren unsere Ergebnisse in einem Modell, das ein Zusammenspiel zwischen verschiedenen chemischen Gruppen in den Antibiotika und spezifischen Basen im Riboschalter zeigt. Allgemeiner unterstreichen unsere Daten das Potenzial von 19F NMR‐Methoden, komplexe Austauschprozesse mit mehreren angeregten Zuständen zu charakterisieren.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

Optimizing frequency sampling in CEST experiments

For the past decade chemical exchange saturation transfer (CEST) experiments have been successfully applied to study exchange processes in biomolecules involving sparsely populated, transiently formed conformers. Initial implementations focused on extensive sampling of the CEST frequency domain, requiring significant measurement times. Here we show that the lengthy sampling schemes often used are not optimal and that reduced frequency sampling schedules can be developed without a priori knowledge of the exchange parameters, that only depend on the chosen B₁ field, and, to a lesser extent, on the intrinsic transverse relaxation rates of ground state spins. The reduced sampling approach described here can be used synergistically with other methods for reducing measurement times such as those that excite multiple frequencies in the CEST dimension simultaneously, or make use of non-uniform sampling of indirectly detected time domains, to further decrease measurement times. The proposed approach is validated by analysis of simulated and experimental datasets.

Clustered tree regression to learn protein energy change with mutated amino acid

Accurate and effective prediction of mutation-induced protein energy change remains a great challenge and of great interest in computational biology. However, high resource consumption and insufficient structural information of proteins severely limit the experimental techniques and structure-based prediction methods. Here, we design a structure-independent protocol to accurately and effectively predict the mutation-induced protein folding free energy change with only sequence, physicochemical and evolutionary features. The proposed clustered tree regression protocol is capable of effectively exploiting the inherent data patterns by integrating unsupervised feature clustering by K-means and supervised tree regression using XGBoost, and thus enabling fast and accurate protein predictions with different mutations, with an average Pearson correlation coefficient of 0.83 and an average root-mean-square error of 0.94kcal/mol. The proposed sequence-based method not only eliminates the dependence on protein structures, but also has potential applications in protein predictions with rare structural information.

NMR spectroscopy, excited states and relevance to problems in cell biology - transient pre-nucleation tetramerization of huntingtin and insights into Huntington's disease

Solution nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for analyzing three-dimensional structure and dynamics of macromolecules at atomic resolution. Recent advances have exploited the unique properties of NMR in exchanging systems to detect, characterize and visualize excited sparsely populated states of biological macromolecules and their complexes, which are only transient. These states are invisible to conventional biophysical techniques, and play a key role in many processes, including molecular recognition, protein folding, enzyme catalysis, assembly and fibril formation. All the NMR techniques make use of exchange between sparsely populated NMR-invisible and highly populated NMR-visible states to transfer a magnetization property from the invisible state to the visible one where it can be easily detected and quantified. There are three classes of NMR experiments that rely on differences in distance, chemical shift or transverse relaxation (molecular mass) between the NMR-visible and -invisible species. Here, I illustrate the application of these methods to unravel the complex mechanism of sub-millisecond pre-nucleation oligomerization of the N-terminal region of huntingtin, encoded by exon-1 of the huntingtin gene, where CAG expansion leads to Huntington's disease, a fatal autosomal-dominant neurodegenerative condition. I also discuss how inhibition of tetramerization blocks the much slower (by many orders of magnitude) process of fibril formation.