Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy

Large dynamics of a phase separating arginine-glycine-rich domain revealed via nuclear and electron spins

Liquid-liquid phase separation is the key process underlying formation of membrane-less compartments in cells. A highly dynamic cellular body with rapid component exchange is Cajal body (CB), which supports the extensive compositional dynamics of the RNA splicing machinery, spliceosome. Here, we select an arginine-glycine (RG)-rich segment of coilin, the major component of CB, establish its RNA-induced phase separation, and through combined use of nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) probes, interrogate its dynamics within the crowded interior of formed droplets. Taking advantage of glycine-based singlet-states, we show that glycines retain a large level of sub-nanoseconds dynamics inside the coilin droplets. Furthermore, the continuous-wave (CW) and electron-electron dipolar (PELDOR) and electron-nucleus hyperfine coupling EPR data (HYSCORE) support the RNA-induced formation of dynamic coilin droplets with high coilin peptide concentrations. The combined NMR and EPR data reveal the high dynamics of the RG-rich coilin within droplets and suggest its potential role in the large dynamics of CBs.

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.

Structural and dynamical determinants of a β-sheet-enriched intermediate involved in amyloid fibrillar assembly of human prion protein

The conformational conversion of the cellular prion protein (PrP^C) into a misfolded, aggregated and infectious scrapie isoform is associated with prion disease pathology and neurodegeneration. Despite the significant number of experimental and theoretical studies the molecular mechanism regulating this structural transition is still poorly understood. Here, via Nuclear Magnetic Resonance (NMR) methodologies we investigate at the atomic level the mechanism of the human HuPrP(90–231) thermal unfolding and characterize the conformational equilibrium between its native structure and a β-enriched intermediate state, named β-PrPI. By comparing the folding mechanisms of metal-free and Cu²⁺-bound HuPrP(23–231) and HuPrP(90–231) we show that the coupling between the N- and C-terminal domains, through transient electrostatic interactions, is the key molecular process in tuning long-range correlated μs–ms dynamics that in turn modulate the folding process. Moreover, via thioflavin T (ThT)-fluorescence fibrillization assays we show that β-PrPI is involved in the initial stages of PrP fibrillation, overall providing a clear molecular description of the initial phases of prion misfolding. Finally, we show by using Real-Time Quaking-Induced Conversion (RT-QuIC) that the β-PrPI acts as a seed for the formation of amyloid aggregates with a seeding activity comparable to that of human infectious prions.

The N-ter domain in HuPrP regulates the folding mechanism by tuning the long-range μs–ms dynamics. Removal of the N-ter domain triggers the formation of a stable β-enriched intermediate state inducing amyloid aggregates with HuPrP^Sc seeding activity.

NMR Relaxation Dispersion Methods for the Structural and Dynamic Analysis of Quickly Interconverting, Low-Populated Conformational Substates

Most biomolecular processes involve proteins shuttling among different conformational states, particularly from highly populated ground states to the lowly populated excited states that determine the interconversion rates and biological function, and which are invisible to most structural biology techniques. These structural transitions are rare and relatively fast: happen in the millisecond-microsecond timescale (ms-μs). NMR spectroscopy can access these timescales via relaxation dispersion techniques (RD-NMR). The exchange parameters extracted from RD-NMR experiments provide pivotal information on these otherwise invisible states that reports on key properties of the high free energy, reactive regions of the protein's energy landscape, including the mechanisms of folding/unfolding and of the interconversion between active and inactive states. Here, we describe a simple, step-by-step protocol to carry out RD-NMR experiments on proteins to detect the existence of such conformational substates and characterize their structural properties (chemical shifts).

Model free analysis of experimental residual dipolar couplings in small organic compounds

Residual dipolar couplings (RDCs) contain information on the relative arrangement and dynamics of internuclear spin vectors in chemical compounds. Classically, RDC data is analyzed by fitting to structure models, while model-free approaches (MFA) directly relate RDCs to the corresponding internuclear vectors. The recently introduced software TITANIA implements the MFA and extracts structure and dynamics parameters directly from RDCs to facilitate de novo structure refinement for small organic compounds. Encouraged by our previous results on simulated data, we herein focus on the prerequisites and challenges faced when using purely experimental data for this approach. These concern mainly the fact that not all couplings are accessible in all media, leading to voids in the RDC matrix and the concomitant effects on the structure refinement. It is shown that RDC data sets obtained experimentally from currently available alignment media and measurement methods are of sufficient quality to allow relative configuration determination even when the relative configuration of the analyte is completely unknown.

TITANIA: Model-Free Interpretation of Residual Dipolar Couplings in the Context of Organic Compounds

Residual dipolar couplings (RDCs) become increasingly important as additional NMR parameters in the structure elucidation of organic compounds but are usually used in fitting procedures to discriminate between (computed) structures that are in accordance with RDCs and others that can be ruled out. Thus, the determination of configurations requires prior structural information. The direct use of RDCs as restraints to construct structures based on RDCs has only recently begun also in organic compounds. No protocol has been published though that uses the vector and dynamics information available in multialignment data sets directly for the joint determination of conformation and configuration of organic compounds. This is proposed in the current study. We show that by employing these data, even a flat or random start structure converges into the correctly configured structure when employing multiple alignment data sets in our iterative procedure. The requirements in terms of the number of RDCs and alignment media necessary are discussed in detail.

Arginine Side-Chain Hydrogen Exchange: Quantifying Arginine Side-Chain Interactions in Solution

The rate with which labile backbone hydrogen atoms in proteins exchange with the solvent has long been used to probe protein interactions in aqueous solutions. Arginine, an essential amino acid found in many interaction interfaces, is capable of an impressive range of interactions via its guanidinium group. The hydrogen exchange rate of the guanidinium hydrogens therefore becomes an important measure to quantify side-chain interactions. Herein we present an NMR method to quantify the hydrogen exchange rates of arginine side-chain 1 Hϵ protons and thus present a method to gauge the strength of arginine side-chain interactions. The method employs 13 C-detection and the one-bond deuterium isotope shift observed for 15 Nϵ to generate two exchanging species in 1 H2 O/2 H2 O mixtures. An application to the protein T4 Lysozyme is shown, where protection factors calculated from the obtained exchange rates correlate well with the interactions observed in the crystal structure. The methodology presented provides an important step towards characterising interactions of arginine side-chains in enzymes, in phase separation, and in protein interaction interfaces in general.

Utilizing dipole-dipole cross-correlated relaxation for the measurement of angles between pairs of opposing CαHα-CαHα bonds in anti-parallel β-sheets

Dipole-dipole cross-correlated relaxation (CCR) between two spin pairs is rich with macromolecular structural and dynamic information on inter-nuclear bond vectors. Measurement of short range dipolar CCR rates has been demonstrated for a variety of inter-nuclear vector spin pairs in proteins and nucleic acids, where the multiple quantum coherence necessary for observing the CCR rate is created by through-bond scalar coupling. In principle, CCR rates can be measured for any pair of inter-nuclear vectors where coherence can be generated between one spin of each spin pair, regardless of both the distance between the two spin pairs and the distance of the two spins forming the multiple quantum coherence. In practice, however, long range CCR (lrCCR) rates are challenging to measure due to difficulties in linking spatially distant spin pairs. By utilizing through-space relaxation allowed coherence transfer (RACT), we have developed a new method for the measurement of lrCCR rates involving CαHα bonds on opposing anti-parallel β-strands. The resulting lrCCR rates are straightforward to interpret since only the angle between the two vectors modulates the strength of the interference effect. We applied our lrCCR measurement to the third immunoglobulin-binding domain of the streptococcal protein G (GB3) and utilize published NMR ensembles and static NMR/X-ray structures to highlight the relationship between the lrCCR rates and the CαHα-CαHα inter-bond angle and bond mobility. Furthermore, we employ the lrCCR rates to guide the selection of sub-ensembles from the published NMR ensembles for enhancing the structural and dynamic interpretation of the data. We foresee this methodology for measuring lrCCR rates as improving the generation of structural ensembles by providing highly accurate details concerning the orientation of CαHα bonds on opposing anti-parallel β-strands.

Label-free NMR-based dissociation kinetics determination

Understanding the dissociation of molecules is the basis to modulate interactions of biomedical interest. Optimizing drugs for dissociation rates is found to be important for their efficacy, selectivity, and safety. Here, we show an application of the high-power relaxation dispersion (RD) method to the determination of the dissociation rates of weak binding ligands from receptors. The experiment probes proton RD on the ligand and, therefore, avoids the need for any isotopic labeling. The large ligand excess eases the detection significantly. Importantly, the use of large spin-lock fields allows the detection of faster dissociation rates than other relaxation approaches. Moreover, this experimental approach allows to access directly the off-rate of the binding process without the need for analyzing a series of samples with increasing ligand saturation. The validity of the method is shown with small molecule interactions using two macromolecules, bovine serum albumin and tubulin heterodimers.

Forbidden Coherence Transfer of 19F Nuclei to Quantitatively Measure the Dynamics of a CF₃-Containing Ligand in Receptor-Bound States

The dynamic property of a ligand in the receptor-bound state is an important metric to characterize the interactions in the ligand–receptor interface, and the development of an experimental strategy to quantify the amplitude of motions in the bound state is of importance to introduce the dynamic aspect into structure-guided drug development (SGDD). Fluorine modifications are frequently introduced at the hit-to-lead optimization stage to enhance the binding potency and other characteristics of a ligand. However, the effects of fluorine modifications are generally difficult to predict, owing to the pleiotropic nature of the interactions. In this study, we report an NMR-based approach to experimentally evaluate the local dynamics of trifluoromethyl (CF₃)-containing ligands in the receptor-bound states. For this purpose, the forbidden coherence transfer (FCT) analysis, which has been used to study the dynamics of methyl moieties in proteins, was extended to the ¹⁹F nuclei of CF₃-containing ligands. By applying this CF₃–FCT analysis to a model interaction system consisting of a ligand, AST-487, and a receptor, p38α, we successfully quantified the amplitude of the CF₃ dynamics in the p38α-bound state. The strategy would bring the CF₃-containing ligands within the scope of dynamic SGDD to improve the affinity and specificity for the drug-target receptors.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

High-power (1)H composite pulse decoupling provides artifact free exchange-mediated saturation transfer (EST) experiments

Exchange-mediated saturation transfer (EST) provides critical information regarding dynamics of molecules. In typical applications EST is studied by either scanning a wide range of (15)N chemical shift offsets where the applied (15)N irradiation field strength is on the order of hundreds of Hertz or, scanning a narrow range of (15)N chemical shift offsets where the applied (15)N irradiation field-strength is on the order of tens of Hertz during the EST period. The (1)H decoupling during the EST delay is critical as incomplete decoupling causes broadening of the EST profile, which could possibly result in inaccuracies of the extracted kinetic parameters and transverse relaxation rates. Currently two different (1)H decoupling schemes have been employed, intermittently applied 180° pulses and composite-pulse-decoupling (CPD), for situations where a wide range, or narrow range of (15)N chemical shift offsets are scanned, respectively. We show that high-power CPD provides artifact free EST experiments, which can be universally implemented regardless of the offset range or irradiation field-strengths.

Kinetics of the Antibody Recognition Site in the Third IgG-Binding Domain of Protein G

Protein dynamics occurring on a wide range of timescales play a crucial role in governing protein function. Particularly, motions between the globular rotational correlation time (τc ) and 40 μs (supra-τc window), strongly influence molecular recognition. This supra-τc window was previously hidden, owing to a lack of experimental methods. Recently, we have developed a high-power relaxation dispersion (RD) experiment for measuring kinetics as fast as 4 μs. For the first time, this method, performed under super-cooled conditions, enabled us to detect a global motion in the first β-turn of the third IgG-binding domain of protein G (GB3), which was extrapolated to 371±115 ns at 310 K. Furthermore, the same residues show the plasticity in the model-free residual dipolar coupling (RDC) order parameters and in an ensemble encoding the supra-τc dynamics. This β-turn is involved in antibody binding, exhibiting the potential link of the observed supra-τc motion with molecular recognition.

Speeding-up exchange-mediated saturation transfer experiments by Fourier transform

Protein motions over various time scales are crucial for protein function. NMR relaxation dispersion experiments play a key role in explaining these motions. However, the study of slow conformational changes with lowly populated states remained elusive. The recently developed exchange-mediated saturation transfer experiments allow the detection and characterization of such motions, but require extensive measurement time. Here we show that, by making use of Fourier transform, the total acquisition time required to measure an exchange-mediated saturation transfer profile can be reduced by twofold in case that one applies linear prediction. In addition, we demonstrate that the analytical solution for R1ρ experiments can be used for fitting the exchange-mediated saturation transfer profile. Furthermore, we show that simultaneous analysis of exchange-mediated saturation transfer profiles with two different radio-frequency field strengths is required for accurate and precise characterization of the exchange process and the exchanging states.

Population shuffling between ground and high energy excited states

Stochastic processes powered by thermal energy lead to protein motions traversing time-scales from picoseconds to seconds. Fundamental to protein functionality is the utilization of these dynamics for tasks such as catalysis, folding, and allostery. A hierarchy of motion is hypothesized to connect and synergize fast and slow dynamics toward performing these essential activities. Population shuffling predicts a "top-down" temporal hierarchy, where slow time-scale conformational interconversion leads to a shuffling of the free energy landscape for fast time-scale events. Until now, population shuffling was only applied to interconverting ground states. Here, we extend the framework of population shuffling to be applicable for a system interconverting between low energy ground and high energy excited states, such as the SH3 domain mutants G48M and A39V/N53P/V55L from the Fyn tyrosine kinase, providing another tool for accessing the structural dynamics of high energy excited states. Our results indicate that the higher energy gauche - rotameric state for the leucine χ2 dihedral angle contributes significantly to the distribution of rotameric states in both the major and minor forms of the SH3 domain. These findings are corroborated with unrestrained molecular dynamics (MD) simulations on both the major and minor states of the SH3 domain demonstrating high correlations between experimental and back-calculated leucine χ2 rotameric populations. Taken together, we demonstrate how fast time-scale rotameric side-chain population distributions can be extracted from slow time-scale conformational exchange data further extending the scope and the applicability of the population shuffling model.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, "conservation of experimental data in analysis of Fourier" (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b 6 f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.

Organic synthesis: march of the machines

Organic synthesis is changing; in a world where budgets are constrained and the environmental impacts of practice are scrutinized, it is increasingly recognized that the efficient use of human resource is just as important as material use. New technologies and machines have found use as methods for transforming the way we work, addressing these issues encountered in research laboratories by enabling chemists to adopt a more holistic systems approach in their work. Modern developments in this area promote a multi-disciplinary approach and work is more efficient as a result. This Review focuses on the concepts, procedures and methods that have far-reaching implications in the chemistry world. Technologies have been grouped as topics of opportunity and their recent applications in innovative research laboratories are described.

Kinetic modulation of a disordered protein domain by phosphorylation

Phosphorylation is a major post-translational mechanism of regulation that frequently targets disordered protein domains, but it remains unclear how phosphorylation modulates disordered states of proteins. Here we determine the kinetics and energetics of a disordered protein domain the kinase-inducible domain (KID) of the transcription factor CREB and that of its phosphorylated form pKID, using high-throughput molecular dynamic simulations. We identify the presence of a metastable, partially ordered state with a 60-fold slowdown in conformational kinetics that arises due to phosphorylation, kinetically stabilizing residues known to participate in an early binding intermediate. We show that this effect is only partially reconstituted by mutation to glutamate, indicating that the phosphate is uniquely required for the long-lived state to arise. This mechanism of kinetic modulation could be important for regulation beyond conformational equilibrium shifts.

Local isotropic diffusion approximation for coupled internal and overall molecular motions in NMR spin relaxation

The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].