Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.     

High-resolution X-ray and NMR structures of the SMN Tudor domain: conformational variation in the binding site for symmetrically dimethylated arginine residues

The SMN protein, which is linked to spinal muscular atrophy (SMA), plays an important role in the assembly of the spliceosomal small nuclear ribonucleoprotein complexes. This function requires binding of SMN to the arginine-glycine (RG) rich C-terminal tails of the Sm proteins, which contain symmetrically dimethylated arginine residues (sDMA) in vivo. Using NMR titrations, we show that the SMN Tudor domain recognizes these sDMAs in the methylated RG repeats. Upon complex formation a cluster of conserved aromatic residues in the SMN Tudor domain interacts with the sDMA methyl groups. We present two high resolution structures of the uncomplexed SMN Tudor domain, a 1.8A crystal structure and an NMR structure that has been refined against a large number of backbone and side-chain residual dipolar couplings. The backbone conformation of both structures is very similar, however, differences are observed for the cluster of conserved aromatic side-chains in the sDMA binding pocket. In order to validate these variations we introduce a novel application of residual dipolar couplings for aromatic rings. We show that structural information can be derived from aromatic ring residual dipolar couplings, even in the presence of internal motions such as ring flipping. These residual dipolar couplings and ring current shifts independently confirm that the SMN Tudor domain adopts two different conformations in the sDMA binding pocket. The observed structural variations may play a role for the recognition of sDMAs.     

Amplitudes and time scales of picosecond-to-microsecond motion in proteins studied by solid-state NMR: a critical evaluation of experimental approaches and application to crystalline ubiquitin

Solid-state NMR provides insight into protein motion over time scales ranging from picoseconds to seconds. While in solution state the methodology to measure protein dynamics is well established, there is currently no such consensus protocol for measuring dynamics in solids. In this article, we perform a detailed investigation of measurement protocols for fast motions, i.e. motions ranging from picoseconds to a few microseconds, which is the range covered by dipolar coupling and relaxation experiments. We perform a detailed theoretical investigation how dipolar couplings and relaxation data can provide information about amplitudes and time scales of local motion. We show that the measurement of dipolar couplings is crucial for obtaining accurate motional parameters, while systematic errors are found when only relaxation data are used. Based on this realization, we investigate how the REDOR experiment can provide such data in a very accurate manner. We identify that with accurate rf calibration, and explicit consideration of rf field inhomogeneities, one can obtain highly accurate absolute order parameters. We then perform joint model-free analyses of 6 relaxation data sets and dipolar couplings, based on previously existing, as well as new data sets on microcrystalline ubiquitin. We show that nanosecond motion can be detected primarily in loop regions, and compare solid-state data to solution-state relaxation and RDC analyses. The protocols investigated here will serve as a useful basis towards the establishment of a routine protocol for the characterization of ps–μs motions in proteins by solid-state NMR.     

Arginine kinase: joint crystallographic and NMR RDC analyses link substrate-associated motions to intrinsic flexibility

The phosphagen kinase family, including creatine and arginine kinases (AKs), catalyzes the reversible transfer of a “high-energy” phosphate between ATP and a phosphoguanidino substrate. They have become a model for the study of both substrate-induced conformational change and intrinsic protein dynamics. Prior crystallographic studies indicated large substrate-induced domain rotations, but differences among a recent set of AK structures were interpreted as a plastic deformation. Here, the structure of Limulus substrate-free AK is refined against high-resolution crystallographic data and compared quantitatively with NMR chemical shifts and residual dipolar couplings (RDCs). This demonstrates the feasibility of this type of RDC analysis of proteins that are large by NMR standards (42 kDa) and illuminates the solution structure, free from crystal-packing constraints. Detailed comparison of the 1.7 Å resolution substrate-free crystal structure against the 1.7 Å transition-state analog complex shows large substrate-induced domain motions that can be broken down into movements of smaller quasi-rigid bodies. The solution-state structure of substrate-free AK is most consistent with an equilibrium of substrate-free and substrate-bound structures, with the substrate-free form dominating, but with varying displacements of the quasi-rigid groups. Rigid-group rotations evident from the crystal structures are about axes previously associated with intrinsic millisecond dynamics using NMR relaxation dispersion. Thus, “substrate-induced” motions are along modes that are intrinsically flexible in the substrate-free enzyme and likely involve some degree of conformational selection.     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.     

Novel NMR tools to study structure and dynamics of biomembranes

Nuclear magnetic resonance (NMR) studies on biomembranes have benefited greatly from introduction of magic angle spinning (MAS) NMR techniques. Improvements in MAS probe technology, combined with the higher magnetic field strength of modern instruments, enables almost liquid-like resolution of lipid resonances. The cross-relaxation rates measured by nuclear Overhauser enhancement spectroscopy (NOESY) provide new insights into conformation and dynamics of lipids with atomic-scale resolution. The data reflect the tremendous motional disorder in the lipid matrix. Transfer of magnetization by spin diffusion along the proton network of lipids is of secondary relevance, even at a long NOESY mixing time of 300 ms. MAS experiments with re-coupling of anisotropic interactions, like the 13C-(1)H dipolar couplings, benefit from the excellent resolution of 13C shifts that enables assignment of the couplings to specific carbon atoms. The traditional 2H NMR experiments on deuterated lipids have higher sensitivity when conducted on oriented samples at higher magnetic field strength. A very large number of NMR parameters from lipid bilayers is now accessible, providing information about conformation and dynamics for every lipid segment. The NMR methods have the sensitivity and resolution to study lipid-protein interaction, lateral lipid organization, and the location of solvents and drugs in the lipid matrix.     

Weak alignment offers new NMR opportunities to study protein structure and dynamics

Protein solution nuclear magnetic resonance (NMR) can be conducted in a slightly anisotropic environment, where the orientational distribution of the proteins is no longer random. In such an environment, the large one-bond internuclear dipolar interactions no longer average to zero and report on the average orientation of the corresponding vectors relative to the magnetic field. The desired very weak ordering, on the order of 10(-3), can be induced conveniently by the use of aqueous nematic liquid crystalline suspensions or by anisotropically compressed hydrogels. The resulting residual dipolar interactions are scaled down by three orders of magnitude relative to their static values, but nevertheless can be measured at high accuracy. They are very precise reporters on the average orientation of bonds relative to the molecular alignment frame, and they can be used in a variety of ways to enrich our understanding of protein structure and function. Applications to date have focused primarily on validation of structures, determined by NMR, X-ray crystallography, or homology modeling, and on refinement of structures determined by conventional NMR approaches. Although de novo structure determination on the basis of dipolar couplings suffers from a severe multiple minimum problem, related to the degeneracy of dipolar coupling relative to inversion of the internuclear vector, a number of approaches can address this problem and potentially can accelerate the NMR structure determination process considerably. In favorable cases, where large numbers of dipolar couplings can be measured, inconsistency between measured values can report on internal motions.     

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.     

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R₁ and R₂), and relaxation dispersion (R_1ρ) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (∼30 μs) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.     

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.     

Solution NMR conformation of glycosaminoglycans

Nuclear magnetic resonance (NMR) spectroscopy has been giving a pivotal contribution to the progress of glycomics, mostly by elucidating the structural, dynamical, conformational and intermolecular binding aspects of carbohydrates. Particularly in the field of conformation, NOE resonances, scalar couplings, residual dipolar couplings, and chemical shift anisotropy offsets have been the principal NMR parameters utilized. Molecular dynamics calculations restrained by NMR-data input are usually employed in conjunction to generate glycosidic bond dihedral angles. Glycosaminoglycans (GAGs) are a special class of sulfated polysaccharides extensively studied worldwide. Besides regulating innumerous physiological processes, these glycans are also widely explored in the global market as either clinical or nutraceutical agents. The conformational aspects of GAGs are key regulators to the quality of interactions with the functional proteins involved in biological events. This report discusses the solution conformation of each GAG type analyzed by one or more of the above-mentioned methods.     

Rapid determination of protein folds using residual dipolar couplings

Over the next few years, various genome projects will sequence many new genes and yield many new gene products. Many of these products will have no known function and little, if any, sequence homology to existing proteins. There is reason to believe that a rapid determination of a protein fold, even at low resolution, can aid in the identification of function and expedite the determination of structure at higher resolution. Recently devised NMR methods of measuring residual dipolar couplings provide one route to the determination of a fold. They do this by allowing the alignment of previously identified secondary structural elements with respect to each other. When combined with constraints involving loops connecting elements or other short-range experimental distance information, a fold is produced. We illustrate this approach to protein fold determination on (15)N-labeled Eschericia coli acyl carrier protein using a limited set of (15)N-(1)H and (1)H-(1)H dipolar couplings. We also illustrate an approach using a more extended set of heteronuclear couplings on a related protein, (13)C, (15)N-labeled NodF protein from Rhizobium leguminosarum.     

Backbone dynamics of the human CC chemokine eotaxin: fast motions, slow motions, and implications for receptor binding.

Eotaxin is a member of the chemokine family of about 40 proteins that induce cell migration. Eotaxin binds the CC chemokine receptor CCR3 that is highly expressed by eosinophils, and it is considered important in the pathology of chronic respiratory disorders such as asthma. The high resolution structure of eotaxin is known. The 74 amino acid protein has two disulfide bridges and shows a typical chemokine fold comprised of a core of three antiparallel beta-strands and an overlying alpha-helix. In this paper, we report the backbone dynamics of eotaxin determined through 15N-T1, T2, and [1H]-15N nuclear Overhauser effect heteronuclear multidimensional NMR experiments. This is the first extensive study of the dynamics of a chemokine derived from 600, 500, and 300 MHz NMR field strengths. From the T1, T2, and NOE relaxation data, parameters that describe the internal motions of eotaxin were derived using the Lipari-Szabo model free analysis. The most ordered regions of the protein correspond to the known secondary structure elements. However, surrounding the core, the regions known to be functionally important in chemokines show a range of motions on varying timescales. These include extensive subnanosecond to picosecond motions in the N-terminus, C-terminus, and the N-loop succeeding the disulfides. Analysis of rotational diffusion anisotropy of eotaxin and chemical exchange terms at multiple fields also allowed the confident identification of slow conformational exchange through the "30s" loop, disulfides, and adjacent residues. In addition, we show that these motions may be attenuated in the dimeric form of a synthetic eotaxin. The structure and dynamical basis for eotaxin receptor binding is discussed in light of the dynamics data.     

Real-time pure shift <Superscript>15</Superscript>N HSQC of proteins: a real improvement in resolution and sensitivity

Spectral resolution in proton NMR spectroscopy is reduced by the splitting of resonances into multiplets due to the effect of homonuclear scalar couplings. Although these effects are often hidden in protein NMR spectroscopy by low digital resolution and routine apodization, behind the scenes homonuclear scalar couplings increase spectral overcrowding. The possibilities for biomolecular NMR offered by new pure shift NMR methods are illustrated here. Both resolution and sensitivity are improved, without any increase in experiment time. In these experiments, free induction decays are collected in short bursts of data acquisition, with durations short on the timescale of J-evolution, interspersed with suitable refocusing elements. The net effect is real-time (t ₂) broadband homodecoupling, suppressing the multiplet structure caused by proton–proton interactions. The key feature of the refocusing elements is that they discriminate between the resonances of active (observed) and passive (coupling partner) spins. This can be achieved either by using band-selective refocusing or by the BIRD element, in both cases accompanied by a nonselective 180° proton pulse. The latter method selects the active spins based on their one-bond heteronuclear J-coupling to ¹⁵N, while the former selects a region of the ¹H spectrum. Several novel pure shift experiments are presented, and the improvements in resolution and sensitivity they provide are evaluated for representative samples: the N-terminal domain of PGK; ubiquitin; and two mutants of the small antifungal protein PAF. These new experiments, delivering improved sensitivity and resolution, have the potential to replace the current standard HSQC experiments.     

Correlated Inter-Domain Motions in Adenylate Kinase

Correlated inter-domain motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. Here we characterize at structural level the inter-domain coupling in a multidomain enzyme, Adenylate Kinase (AK), using computational methods that exploit the shape information encoded in residual dipolar couplings (RDCs) measured under steric alignment by nuclear magnetic resonance (NMR). We find experimental evidence for a multi-state equilibrium distribution along the opening/closing pathway of Adenylate Kinase, previously proposed from computational work, in which inter-domain interactions disfavour states where only the AMP binding domain is closed. In summary, we provide a robust experimental technique for study of allosteric regulation in AK and other enzymes.     

Weak alignment NMR: a hawk-eyed view of biomolecular structure

Imposing a very slight deviation from the isotropic random distribution of macromolecules in solution in an NMR sample tube permits the measurement of residual internuclear dipolar couplings (RDCs). Such interactions are very sensitive functions of the time-averaged orientation of the corresponding internuclear vectors and thereby offer highly precise structural information. In recent years, advances have been made both in the technology to measure RDCs and in the computational procedures that integrate this information in the structure determination process. The exceptional precision with which RDCs can be measured under weakly aligned conditions is also starting to reveal the mostly, but not universally, subtle effects of internal protein dynamics. Importantly, RDCs potentially can reveal motions taking place on a timescale slower than rotational diffusion and analysis is uniquely sensitive to the direction of motion, not just its amplitude.     

Breaking symmetry in the structure determination of (large) symmetric protein dimers

We demonstrate a novel methodology to disrupt the symmetry in the NMR spectra of homodimers. A paramagnetic probe is introduced sub-stoichiometrically to create an asymmetric system with the paramagnetic probe residing on only one monomer within the dimer. This creates sufficient magnetic anisotropy for resolution of symmetry-related overlapped resonances and, consequently, detection of pseudocontact shifts and residual dipolar couplings specific to each monomeric component. These pseudocontact shifts can be readily incorporated into existing structure refinement calculations and enable determination of monomer orientation within the dimeric protein. This methodology can be widely used for solution structure determination of symmetric dimers.     

Mapping backbone dynamics in solution with site-directed spin labeling: GCN4-58 bZip free and bound to DNA

In site-directed spin labeling, a nitroxide-containing side chain is introduced at selected sites in a protein. The EPR spectrum of the labeled protein encodes information about the motion of the nitroxide on the nanosecond time scale, which has contributions from the rotary diffusion of the protein, from internal motions in the side chain, and from backbone fluctuations. In the simplest model for the motion of noninteracting (surface) side chains, the contribution from the internal motion is sequence independent, as is that from protein rotary diffusion. Hence, differences in backbone motions should be revealed by comparing the sequence-dependent motions of nitroxides at structurally homologous sites. To examine this model, nitroxide side chains were introduced, one at a time, along the GCN4-58 bZip sequence, for which NMR (15)N relaxation experiments have identified a striking gradient of backbone mobility along the DNA-binding region [Bracken et al. (1999) J. Mol. Biol. 285, 2133]. Spectral simulation techniques and a simple line width measure were used to extract dynamical parameters from the EPR spectra, and the results reveal a mobility gradient similar to that observed in NMR relaxation, indicating that side chain motions mirror backbone motions. In addition, the sequence-dependent side chain dynamics were analyzed in the DNA/protein complex, which has not been previously investigated by NMR relaxation methods. As anticipated, the backbone motions are damped in the DNA-bound state, although a gradient of motion persists with residues at the DNA-binding site being the most highly ordered, similar to those of helices on globular proteins.