Divided-Evolution-Based Pulse Scheme for Quantifying Exchange Processes in Proteins: Powerful Complement to Relaxation Dispersion Experiments

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.

Investigation of the Low-Populated Excited States of the HIV-1 Nucleocapsid Domain

The nucleocapsid domain (NCd), located at the C-terminus of the HIV-1 Gag protein, is involved in numerous stages of the replication cycle, such as the packaging of the viral genome and reverse transcription. It exists under different forms through the viral life cycle, depending on the processing of Gag by the HIV-1 protease. NCd is constituted of two adjacent zinc knuckles (ZK1 and ZK2), separated by a flexible linker and flanked by disordered regions. Here, conformational equilibria between a major and two minor states were highlighted exclusively in ZK2, by using CPMG and CEST NMR experiments. These minor states appear to be temperature dependent, and their populations are highest at physiological temperature. These minor states are present both in NCp7, the mature form of NCd, and in NCp9 and NCp15, the precursor forms of NCd, with increased populations. The role of these minor states in the targeting of NCd by drugs and its binding properties is discussed.

NMR Chemical Exchange Measurements Reveal That N6-Methyladenosine Slows RNA Annealing

N⁶-Methyladenosine (m⁶A) is an abundant epitranscriptomic modification that plays important roles in many aspects of RNA metabolism. While m⁶A is thought to mainly function by recruiting reader proteins to specific RNA sites, the modification can also reshape RNA-protein and RNA-RNA interactions by altering RNA structure mainly by destabilizing base pairing. Little is known about how m⁶A and other epitranscriptomic modifications might affect the kinetic rates of RNA folding and other conformational transitions that are also important for cellular activity. Here, we used NMR R_1ρ relaxation dispersion and chemical exchange saturation transfer to noninvasively and site-specifically measure nucleic acid hybridization kinetics. The methodology was validated on two DNA duplexes and then applied to examine how a single m⁶A alters the hybridization kinetics in two RNA duplexes. The results show that m⁶A minimally impacts the rate constant for duplex dissociation, changing k_off by ∼1-fold but significantly slows the rate of duplex annealing, decreasing k_on by ∼7-fold. A reduction in the annealing rate was observed robustly for two different sequence contexts at different temperatures, both in the presence and absence of Mg²⁺. We propose that rotation of the N⁶-methyl group from the preferred syn conformation in the unpaired nucleotide to the energetically disfavored anti conformation required for Watson-Crick pairing is responsible for the reduced annealing rate. The results help explain why in mRNA m⁶A slows down tRNA selection and more generally suggest that m⁶A may exert cellular functions by reshaping the kinetics of RNA conformational transitions.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

Measuring the signs of the methyl 1H chemical shift differences between major and 'invisible' minor protein conformational states using methyl 1H multi-quantum spectroscopy

Carr-Purcell-Meiboom-Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and 'invisible' minor states. The exchange rate(s) ([Formula: see text]), population(s) of the minor state(s) and the absolute value of the chemical shift difference [Formula: see text] (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of [Formula: see text] that is required to reconstruct the spectrum of the 'invisible' minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl [Formula: see text] CPMG experiment (Yuwen in Angew Chem 55:11490-11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of [Formula: see text] is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme ([Formula: see text]) and (2) the [Formula: see text] kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus ([Formula: see text]).

Determining rotational dynamics of the guanidino group of arginine side chains in proteins by carbon-detected NMR

A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine to characterise the interactions mediated by arginine side chains.

Arginine residues are imperative for many active sites and protein-interaction interfaces. A new NMR-based method is presented to determine the rotational dynamics around the N_ε–C_ζ bond of arginine side chains. An application to a 19 kDa protein shows that the strengths of interactions involving arginine side chains can be characterised.

Mechanisms of amyloid formation revealed by solution NMR

Amyloid fibrils are proteinaceous elongated aggregates involved in more than fifty human diseases. Recent advances in electron microscopy and solid state NMR have allowed the characterization of fibril structures to different extents of refinement. However, structural details about the mechanism of fibril formation remain relatively poorly defined. This is mainly due to the complex, heterogeneous and transient nature of the species responsible for assembly; properties that make them difficult to detect and characterize in structural detail using biophysical techniques. The ability of solution NMR spectroscopy to investigate exchange between multiple protein states, to characterize transient and low-population species, and to study high molecular weight assemblies, render NMR an invaluable technique for studies of amyloid assembly. In this article we review state-of-the-art solution NMR methods for investigations of: (a) protein dynamics that lead to the formation of aggregation-prone species; (b) amyloidogenic intrinsically disordered proteins; and (c) protein-protein interactions on pathway to fibril formation. Together, these topics highlight the power and potential of NMR to provide atomic level information about the molecular mechanisms of one of the most fascinating problems in structural biology.

NMR studies of the dynamics of high-spin nitrophorins: comparative studies of NP4 and NP2 at close to physiological pH

The β-barrel nitrophorin (NP) heme proteins are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands. NO is bound to iron of the NPs and is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1–NP4, which have sequence similarities in two pairs, NP1 and NP4 (90% identical) and NP2 and NP3 (80% identical). The available crystal structures of NP4 have been used to propose that pH-dependent changes in the conformation of two loops between adjacent β-strands at the front opening of the protein, the A–B and G–H loops, determine the rate of NO release. At pH 7.3, NP4 releases NO 17 times faster than NP2 does. In this work, the aqua complexes of NP4 and NP2 have been investigated by nuclear magnetic resonance (NMR) relaxation measurements to probe the pico- to nanosecond and micro- to millisecond time scale motions at two pH values, 6.5 and 7.3. It is found that NP4-OH₂ is fairly rigid and only residues in the loop regions show dynamics at pH 6.5; at pH 7.3, much more dynamics of the loops and most of the β-strands are observed while the α-helices remain fairly rigid. In comparison, NP2-OH₂ shows much less dynamics, albeit somewhat more than that of the previously reported NP2-NO complex [Muthu, D., Berry, R. E., Zhang, H., and Walker, F. A. (2013) Biochemistry 52, 7910–7925]. The reasons for this major difference between NP4 and NP2 are discussed.

NMR studies of the dynamics of nitrophorin 2 bound to nitric oxide

The Rhodnius nitrophorins are β-barrel proteins of the lipocalin fold with a heme protruding from the open end of the barrel. They are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands, where NO is bound to iron. NO is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1-NP4. At pH 7.3, NP4 releases NO 17 times faster than NP2 does, as measured by stopped-flow kinetics. A number of crystal structures of the least abundant protein, NP4, are available. These structures have been used to propose that two loops between adjacent β-strands at the front opening of the protein, the A-B and G-H loops, determine the rate of NO release. To learn how the protein loops contribute to the release of NO for each of the nitrophorins, the dynamics of these proteins are being studied in our laboratory. In this work, the NP2-NO complex has been investigated by nuclear magnetic resonance relaxation measurements to probe the picosecond-to-nanosecond and microsecond-to-millisecond time scale motions at three pH values, 5.0, 6.5, and 7.3. It is found that at pH 5.0 and 6.5, the NP2-NO complex is rigid and only a few residues in the loop regions show dynamics, while at pH 7.3, somewhat more dynamics, particularly of the A-B loop, are observed. Comparison to other lipocalins shows that all are relatively rigid, and that the dynamics of lipocalins in general are much more subtle than those of mainly α-helical proteins.

Heteronuclear Adiabatic Relaxation Dispersion (HARD) for quantitative analysis of conformational dynamics in proteins

NMR relaxation methods probe biomolecular motions over a wide range of timescales. In particular, the rotating frame spin-lock R(1ρ) and Carr-Purcell-Meiboom-Gill (CPMG) R(2) experiments are commonly used to characterize μs to ms dynamics, which play a critical role in enzyme folding and catalysis. In an effort to complement these approaches, we introduced the Heteronuclear Adiabatic Relaxation Dispersion (HARD) method, where dispersion in rotating frame relaxation rate constants (longitudinal R(1ρ) and transverse R(2ρ)) is created by modulating the shape and duration of adiabatic full passage (AFP) pulses. Previously, we showed the ability of the HARD method to detect chemical exchange dynamics in the fast exchange regime (k(ex)∼10(4)-10(5) s(-1)). In this article, we show the sensitivity of the HARD method to slower exchange processes by measuring R(1ρ) and R(2ρ) relaxation rates for two soluble proteins (ubiquitin and 10C RNA ligase). One advantage of the HARD method is its nominal dependence on the applied radio frequency field, which can be leveraged to modulate the dispersion in the relaxation rate constants. In addition, we also include product operator simulations to define the dynamic range of adiabatic R(1ρ) and R(2ρ) that is valid under all exchange regimes. We conclude from both experimental observations and simulations that this method is complementary to CPMG-based and rotating frame spin-lock R(1ρ) experiments to probe conformational exchange dynamics for biomolecules. Finally, this approach is germane to several NMR-active nuclei, where relaxation rates are frequency-offset independent.

Increasing the exchange time-scale that can be probed by CPMG relaxation dispersion NMR

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy has emerged as a valuable tool to characterize conformational exchange between major and minor states in a large variety of biomolecules. The window of exchange that is amenable for study, corresponding to rates on the order of 2000 s(-1) or less, is limiting, however. Here we show that a combined analysis of both amide (15)N and (1)H(N) CPMG profiles and major state exchange induced (15)N chemical shift changes leads to significant increases in the exchange time scale for which accurate exchange parameters and chemical shift differences between the interconverting states can be obtained. The utility of the approach is illustrated with examples involving a pair of protein systems that are in the moderately fast exchange regime. In these cases the analysis of dispersion profiles alone is not sufficient to obtain robust measures of exchange parameters and chemical shift differences. Inclusion of major state exchange induced (15)N chemical shift changes measured in ((15)N-(1)H(N)) HMQC and HSQC data sets in addition to the (15)N and (1)H(N) dispersion profiles in the analysis "breaks" the correlation in parameters, allowing accurate values to be obtained. The approach is straightforward to implement and makes use of HMQC/HSQC data sets that are recorded as a matter of routine to obtain chemical shifts of the excited state. It promises to increase the range of exchanging systems involving low populated, transiently formed excited states that can be studied by relaxation dispersion NMR.