Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution

High-performance virtual screening by targeting a high-resolution RNA dynamic ensemble

Dynamic ensembles hold great promise in advancing RNA-targeted drug discovery. Here we subjected the transactivation response element (TAR) RNA from human immunodeficiency virus type-1 to experimental high-throughput screening against ~100,000 drug-like small molecules. Results were augmented with 170 known TAR-binding molecules and used to generate sublibraries optimized for evaluating enrichment when virtually screening a dynamic ensemble of TAR determined by combining NMR spectroscopy data and molecular dynamics simulations. Ensemble-based virtual screening scores molecules with an area under the receiver operator characteristic curve of ~0.85-0.94 and with ~40-75% of all hits falling within the top 2% of scored molecules. The enrichment decreased significantly for ensembles generated from the same molecular dynamics simulations without input NMR data and for other control ensembles. The results demonstrate that experimentally determined RNA ensembles can significantly enrich libraries with true hits and that the degree of enrichment is dependent on the accuracy of the ensemble.

Overview of refinement procedures within REFMAC5: utilizing data from different sources

Refinement is a process that involves bringing into agreement the structural model, available prior knowledge and experimental data. To achieve this, the refinement procedure optimizes a posterior conditional probability distribution of model parameters, including atomic coordinates, atomic displacement parameters (B factors), scale factors, parameters of the solvent model and twin fractions in the case of twinned crystals, given observed data such as observed amplitudes or intensities of structure factors. A library of chemical restraints is typically used to ensure consistency between the model and the prior knowledge of stereochemistry. If the observation-to-parameter ratio is small, for example when diffraction data only extend to low resolution, the Bayesian framework implemented in REFMAC5 uses external restraints to inject additional information extracted from structures of homologous proteins, prior knowledge about secondary-structure formation and even data obtained using different experimental methods, for example NMR. The refinement procedure also generates the `best' weighted electron-density maps, which are useful for further model (re)building. Here, the refinement of macromolecular structures using REFMAC5 and related tools distributed as part of the CCP4 suite is discussed.

NMR Spectroscopy to Study MAP Kinase Binding to MAP Kinase Phosphatases

NMR spectroscopy and other solution methods are increasingly being used to obtain novel insights into the mechanisms by which MAPK regulatory proteins bind and direct the activity of MAPKs. Here, we describe how interactions between the MAPK p38α and its regulatory proteins are studied using NMR spectroscopy, isothermal titration calorimetry, and small angle X-ray scattering (SAXS).

Unequivocal determination of caulamidines A and B: application and validation of new tools in the structure elucidation tool box

Newly described NMR experimental approaches can provide valuable structural details and a complementary means of structure verification.

Ambiguities and errors in the structural assignment of organic molecules hinder both drug discovery and total synthesis efforts. Newly described NMR experimental approaches can provide valuable structural details and a complementary means of structure verification. The caulamidines are trihalogenated alkaloids from a marine bryozoan with an unprecedented structural scaffold. Their unique carbon and nitrogen framework was deduced by conventional NMR methods supplemented by new experiments that define 2-bond heteronuclear connectivities, reveal very long-range connectivity data, or visualize the ^35,37Cl isotopic effect on chlorinated carbons. Computer-assisted structural elucidation (CASE) analysis of the spectroscopic data for caulamidine A provided only one viable structural alternative. Anisotropic NMR parameters, specifically residual dipolar coupling and residual chemical shift anisotropy data, were measured for caulamidine A and compared to DFT-calculated values for the proposed structure, the CASE-derived alternative structure, and two energetically feasible stereoisomers. Anisotropy-based NMR experiments provide a global, orthogonal means to verify complex structures free from investigator bias. The anisotropic NMR data were fully consistent with the assigned structure and configuration of caulamidine A. Caulamidine B has the same heterocyclic scaffold as A but a different composition and pattern of halogen substitution. Caulamidines A and B inhibited both wild-type and drug-resistant strains of the malaria parasite Plasmodium falciparum at low micromolar concentrations, yet were nontoxic to human cells.

Perspectives on paramagnetic NMR from a life sciences infrastructure

The effects arising in NMR spectroscopy because of the presence of unpaired electrons, collectively referred to as "paramagnetic NMR" have attracted increasing attention over the last decades. From the standpoint of the structural and mechanistic biology, paramagnetic NMR provides long range restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements through NMR and X-ray data. These restraints also provide information on structure rearrangements and conformational variability in biomolecular systems. Theoretical improvements in quantum chemistry calculations can nowadays allow for accurate calculations of the paramagnetic data from a molecular structural model, thus providing a tool to refine the metal coordination environment by matching the paramagnetic effects observed far away from the metal. Furthermore, the availability of an improved technology (higher fields and faster magic angle spinning) has promoted paramagnetic NMR applications in the fast-growing area of biomolecular solid-state NMR. Major improvements in dynamic nuclear polarization have been recently achieved, especially through the exploitation of the Overhauser effect occurring through the contact-driven relaxation mechanism: the very large enhancement of the ¹³C signal observed in a variety of liquid organic compounds at high fields is expected to open up new perspectives for applications of solution NMR.

Progressive Stereo Locking (PSL): A Residual Dipolar Coupling Based Force Field Method for Determining the Relative Configuration of Natural Products and Other Small Molecules

Establishing the relative configuration of a bioactive natural product represents the most challenging part in determining its structure. Residual dipolar couplings (RDCs) are sensitive probes of the relative spatial orientation of internuclear vectors. We adapted a force field structure calculation methodology to allow free sampling of both R and S configurations of the stereocenters of interest. The algorithm uses a floating alignment tensor in a simulated annealing protocol to identify the conformations and configurations that best fit experimental RDC and distance restraints (from NOE and J-coupling data). A unique configuration (for rigid molecules) or a very small number of configurations (for less rigid molecules) of the structural models having the lowest chiral angle energies and reasonable magnitudes of the alignment tensor are provided as the best predictions of the unknown configuration. For highly flexible molecules, the progressive locking of their stereocenters into their statistically dominant R or S state dramatically reduces the number of possible relative configurations. The result is verified by checking that the same configuration is obtained by initiating the locking from different regions of the molecule. For all molecules tested having known configurations (with conformations ranging from mostly rigid to highly flexible), the method accurately determined the correct configuration.

Integrating NOE and RDC using sum-of-squares relaxation for protein structure determination

We revisit the problem of protein structure determination from geometrical restraints from NMR, using convex optimization. It is well-known that the NP-hard distance geometry problem of determining atomic positions from pairwise distance restraints can be relaxed into a convex semidefinite program (SDP). However, often the NOE distance restraints are too imprecise and sparse for accurate structure determination. Residual dipolar coupling (RDC) measurements provide additional geometric information on the angles between atom-pair directions and axes of the principal-axis-frame. The optimization problem involving RDC is highly non-convex and requires a good initialization even within the simulated annealing framework. In this paper, we model the protein backbone as an articulated structure composed of rigid units. Determining the rotation of each rigid unit gives the full protein structure. We propose solving the non-convex optimization problems using the sum-of-squares (SOS) hierarchy, a hierarchy of convex relaxations with increasing complexity and approximation power. Unlike classical global optimization approaches, SOS optimization returns a certificate of optimality if the global optimum is found. Based on the SOS method, we proposed two algorithms-RDC-SOS and RDC-NOE-SOS, that have polynomial time complexity in the number of amino-acid residues and run efficiently on a standard desktop. In many instances, the proposed methods exactly recover the solution to the original non-convex optimization problem. To the best of our knowledge this is the first time SOS relaxation is introduced to solve non-convex optimization problems in structural biology. We further introduce a statistical tool, the Cramér-Rao bound (CRB), to provide an information theoretic bound on the highest resolution one can hope to achieve when determining protein structure from noisy measurements using any unbiased estimator. Our simulation results show that when the RDC measurements are corrupted by Gaussian noise of realistic variance, both SOS based algorithms attain the CRB. We successfully apply our method in a divide-and-conquer fashion to determine the structure of ubiquitin from experimental NOE and RDC measurements obtained in two alignment media, achieving more accurate and faster reconstructions compared to the current state of the art.

Structural basis for ligand binding to an enzyme by a conformational selection pathway

Proteins can bind target molecules through either induced fit or conformational selection pathways. In the conformational selection model, a protein samples a scarcely populated high-energy state that resembles a target-bound conformation. In enzymatic catalysis, such high-energy states have been identified as crucial entities for activity and the dynamic interconversion between ground states and high-energy states can constitute the rate-limiting step for catalytic turnover. The transient nature of these states has precluded direct observation of their properties. Here, we present a molecular description of a high-energy enzyme state in a conformational selection pathway by an experimental strategy centered on NMR spectroscopy, protein engineering, and X-ray crystallography. Through the introduction of a disulfide bond, we succeeded in arresting the enzyme adenylate kinase in a closed high-energy conformation that is on-pathway for catalysis. A 1.9-Å X-ray structure of the arrested enzyme in complex with a transition state analog shows that catalytic sidechains are properly aligned for catalysis. We discovered that the structural sampling of the substrate free enzyme corresponds to the complete amplitude that is associated with formation of the closed and catalytically active state. In addition, we found that the trapped high-energy state displayed improved ligand binding affinity, compared with the wild-type enzyme, demonstrating that substrate binding to the high-energy state is not occluded by steric hindrance. Finally, we show that quenching of fast time scale motions observed upon ligand binding to adenylate kinase is dominated by enzyme-substrate interactions and not by intramolecular interactions resulting from the conformational change.

Recent excitements in protein NMR: Large proteins and biologically relevant dynamics

The advent of Transverse Relaxation Optimized SpectroscopY (TROSY) and perdeuteration allowed biomolecular NMR spectroscopists to overcome the size limitation barrier (approx. 20 kDa) in de novo structure determination of proteins. The utility of these techniques was immediately demonstrated on large proteins and protein complexes (e.g. GroELGroES, ClpP protease, Hsp90-p53, 20S proteasome, etc.). Further, recent methodological developments such as Residual Dipolar Couplings and Paramagnetic Relaxation Enhancement allowed accurate measurement of long-range structural restraints. Additionally, Carr-Purcell-Meiboom-Gill (CPMG), rotating frame relaxation experiments (R1(rho)) and saturation transfer experiments (CEST and DEST) created never-before accessibility to the (mu)s-ms timescale dynamic parameters that led to the deeper understanding of biological processes. Meanwhile, the excitement in the field continued with a series of developments in the fast data acquisition methods allowing rapid structural studies on less stable proteins. This review aims to discuss important developments in the field of biomolecular NMR spectroscopy in the recent past, i.e., in the post TROSY era. These developments not only gave access to the structural studies of large protein assemblies, but also revolutionized tools in the arsenal of today's biomolecular NMR and point to a bright future of biomolecular NMR spectroscopy.

Designing DNA Nanotube Liquid Crystals as a Weak-Alignment Medium for NMR Structure Determination of Membrane Proteins

Thirty percent of the human proteome is composed of membrane proteins that can perform a wide range of cellular functions and communications. They represent the core of modern medicine as the targets of about 50 % of all prescription pharmaceuticals. However, elucidating the structure of membrane proteins has represented a constant challenge, even in the modern era. To date, only a few hundred high-resolution structural models of membrane proteins are available. This chapter describes the emergence of DNA nanotechnology as a powerful tool for the structural characterization of membrane protein using solution-state nuclear magnetic resonance (NMR) spectroscopy. Here, we detail the large-scale synthesis of detergent-resistant DNA nanotubes that can be assembled into a dilute liquid crystal to be used as a weak-alignment media in solution NMR structure determination of membrane proteins.

Perspectives in the application of residual dipolar couplings in the structure elucidation of weakly aligned small molecules

This perspective article aims to review the general methodology in the application of residual dipolar couplings (RDCs) in the structure elucidation of small molecules and give the author's view on challenges for future applications. Recent improvements in the availability of alignment media, new pulse sequences for the measurement of couplings and improvements in the analysis software have garnered widespread interest in the technique. However, further generalization is needed in order to make RDC analysis into a truly "routine" method. Copyright © 2016 John Wiley & Sons, Ltd.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Extraction of distance restraints from pure shift NOE experiments

NMR techniques incorporating pure shift methods to improve signal resolution have recently attracted much attention, owing to their potential use in studies of increasingly complex molecular systems. Extraction of frequencies from these simplified spectra enables easier structure determination, but only a few of the methods presented provide structural parameters derived from signal integral measurements. In particular, for quantification of the nuclear Overhauser effect (NOE) it is highly desirable to utilize pure shift techniques where signal overlap normally prevents accurate signal integration, to enable measurement of a larger number of interatomic distances. However, robust methods for the measurement of interatomic distances using the recently developed pure shift techniques have not been reported to date. In this work we discuss some of the factors determining the accuracy of measurements of signal integrals in interferogram-based Zangger-Sterk (ZS) pure shift NMR experiments. The ZS broadband homodecoupling technique is used in different experiments designed for quantitative NOE determination from pure shift spectra. It is shown that the techniques studied can be used for quantitative extraction of NOE-derived distance restraints, as exemplified for the test case of strychnine.

Structure of a low-population binding intermediate in protein-RNA recognition

The interaction of the HIV-1 protein transactivator of transcription (Tat) and its cognate transactivation response element (TAR) RNA transactivates viral transcription and represents a paradigm for the widespread occurrence of conformational rearrangements in protein-RNA recognition. Although the structures of free and bound forms of TAR are well characterized, the conformations of the intermediates in the binding process are still unknown. By determining the free energy landscape of the complex using NMR residual dipolar couplings in replica-averaged metadynamics simulations, we observe two low-population intermediates. We then rationally design two mutants, one in the protein and another in the RNA, that weaken specific nonnative interactions that stabilize one of the intermediates. By using surface plasmon resonance, we show that these mutations lower the release rate of Tat, as predicted. These results identify the structure of an intermediate for RNA-protein binding and illustrate a general strategy to achieve this goal with high resolution.

Deciphering solution scattering data with experimentally guided molecular dynamics simulations

Time-resolved X-ray solution scattering is an increasingly popular method to measure conformational changes in proteins. Extracting structural information from the resulting difference X-ray scattering data is a daunting task. We present a method in which the limited but precious information encoded in such scattering curves is combined with the chemical knowledge of molecular force fields. The molecule of interest is then refined toward experimental data using molecular dynamics simulation. Therefore, the energy landscape is biased toward conformations that agree with experimental data. We describe and verify the method, and we provide an implementation in GROMACS.

Solution NMR Spectroscopy for the Study of Enzyme Allostery

Allostery is a ubiquitous biological regulatory process in which distant binding sites within a protein or enzyme are functionally and thermodynamically coupled. Allosteric interactions play essential roles in many enzymological mechanisms, often facilitating formation of enzyme-substrate complexes and/or product release. Thus, elucidating the forces that drive allostery is critical to understanding the complex transformations of biomolecules. Currently, a number of models exist to describe allosteric behavior, taking into account energetics as well as conformational rearrangements and fluctuations. In the following Review, we discuss the use of solution NMR techniques designed to probe allosteric mechanisms in enzymes. NMR spectroscopy is unequaled in its ability to detect structural and dynamical changes in biomolecules, and the case studies presented herein demonstrate the range of insights to be gained from this valuable method. We also provide a detailed technical discussion of several specialized NMR experiments that are ideally suited for the study of enzymatic allostery.

Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Solving structures or structural ensembles of large macromolecular systems in solution poses a challenging problem. While NMR provides structural information at atomic resolution, increased spectral complexity, chemical shift overlap, and short transverse relaxation times (associated with slow tumbling) render application of the usual techniques that have been so successful for medium sized systems (<50 kDa) difficult. Solution X-ray scattering, on the other hand, is not limited by molecular weight but only provides low resolution structural information related to the overall shape and size of the system under investigation. Here we review how combining atomic resolution structures of smaller domains with sparse experimental data afforded by NMR residual dipolar couplings (which yield both orientational and shape information) and solution X-ray scattering data in rigid-body simulated annealing calculations provides a powerful approach for investigating the structural aspects of conformational dynamics in large multidomain proteins. The application of this hybrid methodology is illustrated for the 128 kDa dimer of bacterial Enzyme I which exists in a variety of open and closed states that are sampled at various points in the catalytic cycles, and for the capsid protein of the human immunodeficiency virus.

Pf1 filamentous phage as an alignment tool for generating local and global structural information in nucleic acids

Abstract Pf1 filamentous phage represent a simple versatile method for generating partially ordered macromolecules in solution. The phage allow tunable degrees of alignment of macromolecules under a wide range of temperature and solvent conditions. The negatively charged phage are ideal for aligning negatively charged nucleic acids and these phage-nucleic acid solutions are stable indefinitely. We have used Pf1 phage to align various DNA and RNA molecules in solution for measurement of dipolar coupling interactions. These dipolar couplings can be used to improve the local structure of nucleic acids. More importantly they also contain information on the global structure, such as DNA bending, which presently cannot be obtained by standard NMR methods. The principles involved in using Pf1 phage to generate solutions of partially order macromolecules will be discussed. The use of (1)H-(1)H, (1)H-(13)C and (1)H-(15)N dipolar couplings for generating angle constraints for structure refinement of nucleic acids will also be discussed.

Site-Specific Internal Motions in GB1 Protein Microcrystals Revealed by 3D ²H-¹³C-¹³C Solid-State NMR Spectroscopy

(2)H quadrupolar line shapes deliver rich information about protein dynamics. A newly designed 3D (2)H-(13)C-(13)C solid-state NMR magic angle spinning (MAS) experiment is presented and demonstrated on the microcrystalline β1 immunoglobulin binding domain of protein G (GB1). The implementation of (2)H-(13)C adiabatic rotor-echo-short-pulse-irradiation cross-polarization (RESPIRATION CP) ensures the accuracy of the extracted line shapes and provides enhanced sensitivity relative to conventional CP methods. The 3D (2)H-(13)C-(13)C spectrum reveals (2)H line shapes for 140 resolved aliphatic deuterium sites. Motional-averaged (2)H quadrupolar parameters obtained from the line-shape fitting identify side-chain motions. Restricted side-chain dynamics are observed for a number of polar residues including K13, D22, E27, K31, D36, N37, D46, D47, K50, and E56, which we attribute to the effects of salt bridges and hydrogen bonds. In contrast, we observe significantly enhanced side-chain flexibility for Q2, K4, K10, E15, E19, N35, N40, and E42, due to solvent exposure and low packing density. T11, T16, and T17 side chains exhibit motions with larger amplitudes than other Thr residues due to solvent interactions. The side chains of L5, V54, and V29 are highly rigid because they are packed in the core of the protein. High correlations were demonstrated between GB1 side-chain dynamics and its biological function. Large-amplitude side-chain motions are observed for regions contacting and interacting with immunoglobulin G (IgG). In contrast, rigid side chains are primarily found for residues in the structural core of the protein that are absent from protein binding and interactions.

How to tackle protein structural data from solution and solid state: An integrated approach

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

An encodable lanthanide binding tag with reduced size and flexibility for measuring residual dipolar couplings and pseudocontact shifts in large proteins

Metal ions serve important roles in structural biology applications from long-range perturbations seen in magnetic resonance experiments to electron-dense signatures in X-ray crystallography data; however, the metal ion must be secured in a molecular framework to achieve the maximum benefit. Polypeptide-based lanthanide-binding tags (LBTs) represent one option that can be directly encoded within a recombinant protein expression construct. However, LBTs often exhibit significant mobility relative to the target molecule. Here we report the characterization of improved LBTs sequences for insertion into a protein loop. These LBTs were inserted to connect two parallel alpha helices of an immunoglobulin G (IgG)-binding Z domain platform. Variants A and B bound Tb(3+) with high affinity (0.70 and 0.13 μM, respectively) and displayed restricted LBT motion. Compared to the parent construct, the metal-bound A experienced a 2.5-fold reduction in tag motion as measured by magnetic field-induced residual dipolar couplings and was further studied in a 72.2 kDa complex with the human IgG1 fragment crystallizable (IgG1 Fc) glycoprotein. The appearance of both pseudo-contact shifts (-0.221 to 0.081 ppm) and residual dipolar couplings (-7.6 to 14.3 Hz) of IgG1 Fc resonances in the IgG1 Fc:(variant A:Tb(3+))2 complex indicated structural restriction of the LBT with respect to the Fc. These studies highlight the applicability of improved LBT sequences with reduced mobility to probe the structure of macromolecular systems.

Analysis of the solution conformations of T4 lysozyme by paramagnetic NMR spectroscopy

Paramagnetic data show that the average structure of T4-lysozyme in solution is more open than its crystal structure.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Nitrogen detected TROSY at high field yields high resolution and sensitivity for protein NMR

Detection of (15)N in multidimensional NMR experiments of proteins has sparsely been utilized because of the low gyromagnetic ratio (γ) of nitrogen and the presumed low sensitivity of such experiments. Here we show that selecting the TROSY components of proton-attached (15)N nuclei (TROSY (15)NH) yields high quality spectra in high field magnets (>600 MHz) by taking advantage of the slow (15)N transverse relaxation and compensating for the inherently low (15)N sensitivity. The (15)N TROSY transverse relaxation rates increase modestly with molecular weight but the TROSY gain in peak heights depends strongly on the magnetic field strength. Theoretical simulations predict that the narrowest line width for the TROSY (15)NH component can be obtained at 900 MHz, but sensitivity reaches its maximum around 1.2 GHz. Based on these considerations, a (15)N-detected 2D (1)H-(15)N TROSY-HSQC ((15)N-detected TROSY-HSQC) experiment was developed and high-quality 2D spectra were recorded at 800 MHz in 2 h for 1 mM maltose-binding protein at 278 K (τc ~ 40 ns). Unlike for (1)H detected TROSY, deuteration is not mandatory to benefit (15)N detected TROSY due to reduced dipolar broadening, which facilitates studies of proteins that cannot be deuterated, especially in cases where production requires eukaryotic expression systems. The option of recording (15)N TROSY of proteins expressed in H2O media also alleviates the problem of incomplete amide proton back exchange, which often hampers the detection of amide groups in the core of large molecular weight proteins that are expressed in D2O culture media and cannot be refolded for amide back exchange. These results illustrate the potential of (15)NH-detected TROSY experiments as a means to exploit the high resolution offered by high field magnets near and above 1 GHz.

Ensembles of a small number of conformations with relative populations

In our previous work, we proposed a new way to represent protein native states, using ensembles of a small number of conformations with relative Populations, or ESP in short. Using Ubiquitin as an example, we showed that using a small number of conformations could greatly reduce the potential of overfitting and assigning relative populations to protein ensembles could significantly improve their quality. To demonstrate that ESP indeed is an excellent alternative to represent protein native states, in this work we compare the quality of two ESP ensembles of Ubiquitin with several well-known regular ensembles or average structure representations. Extensive amount of significant experimental data are employed to achieve a thorough assessment. Our results demonstrate that ESP ensembles, though much smaller in size comparing to regular ensembles, perform equally or even better sometimes in all four different types of experimental data used in the assessment, namely, the residual dipolar couplings, residual chemical shift anisotropy, hydrogen exchange rates, and solution scattering profiles. This work further underlines the significance of having relative populations in describing the native states.

Structural and Dynamic Basis for Low-Affinity, High-Selectivity Binding of L-Glutamine by the Glutamine Riboswitch

Naturally occurring L-glutamine riboswitches occur in cyanobacteria and marine metagenomes, where they reside upstream of genes involved in nitrogen metabolism. By combining X-ray, NMR, and MD, we characterized an L-glutamine-dependent conformational transition in the Synechococcus elongatus glutamine riboswitch from tuning fork to L-shaped alignment of stem segments. This transition generates an open ligand-binding pocket with L-glutamine selectivity enforced by Mg(2+)-mediated intermolecular interactions. The transition also stabilizes the P1 helix through a long-range "linchpin" Watson-Crick G-C pair-capping interaction, while melting a short helix below P1 potentially capable of modulating downstream readout. NMR data establish that the ligand-free glutamine riboswitch in Mg(2+) solution exists in a slow equilibrium between flexible tuning fork and a minor conformation, similar, but not identical, to the L-shaped bound conformation. We propose that an open ligand-binding pocket combined with a high conformational penalty for forming the ligand-bound state provide mechanisms for reducing binding affinity while retaining high selectivity.

Applications of high dimensionality experiments to biomolecular NMR

High dimensionality NMR experiments facilitate resonance assignment and precise determination of spectral parameters such as coupling constants. Sparse non-uniform sampling enables acquisition of experiments of high dimensionality with high resolution in acceptable time. In this review we present and compare some significant applications of NMR experiments of dimensionality higher than three in the field of biomolecular studies in solution.

Measuring Residual Dipolar Couplings in Excited Conformational States of Nucleic Acids by CEST NMR Spectroscopy

Nucleic acids undergo structural transitions to access sparsely populated and transiently lived conformational states--or excited conformational states--that play important roles in diverse biological processes. Despite ever-increasing detection of these functionally essential states, 3D structure determination of excited states (ESs) of RNA remains elusive. This is largely due to challenges in obtaining high-resolution structural constraints in these ESs by conventional structural biology approaches. Here, we present nucleic-acid-optimized chemical exchange saturation transfer (CEST) NMR spectroscopy for measuring residual dipolar couplings (RDCs), which provide unique long-range angular constraints in ESs of nucleic acids. We demonstrate these approaches on a fluoride riboswitch, where one-bond (13)C-(1)H RDCs from both base and sugar moieties provide direct structural probes into an ES of the ligand-free riboswitch.

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.

Unveiling Inherent Degeneracies in Determining Population-Weighted Ensembles of Interdomain Orientational Distributions Using NMR Residual Dipolar Couplings: Application to RNA Helix Junction Helix Motifs

A growing number of studies employ time-averaged experimental data to determine dynamic ensembles of biomolecules. While it is well-known that different ensembles can satisfy experimental data to within error, the extent and nature of these degeneracies, and their impact on the accuracy of the ensemble determination remains poorly understood. Here, we use simulations and a recently introduced metric for assessing ensemble similarity to explore degeneracies in determining ensembles using NMR residual dipolar couplings (RDCs) with specific application to A-form helices in RNA. Various target ensembles were constructed representing different domain-domain orientational distributions that are confined to a topologically restricted (<10%) conformational space. Five independent sets of ensemble averaged RDCs were then computed for each target ensemble and a "sample and select" scheme used to identify degenerate ensembles that satisfy RDCs to within experimental uncertainty. We find that ensembles with different ensemble sizes and that can differ significantly from the target ensemble (by as much as ∑Ω ∼ 0.4 where ∑Ω varies between 0 and 1 for maximum and minimum ensemble similarity, respectively) can satisfy the ensemble averaged RDCs. These deviations increase with the number of unique conformers and breadth of the target distribution, and result in significant uncertainty in determining conformational entropy (as large as 5 kcal/mol at T = 298 K). Nevertheless, the RDC-degenerate ensembles are biased toward populated regions of the target ensemble, and capture other essential features of the distribution, including the shape. Our results identify ensemble size as a major source of uncertainty in determining ensembles and suggest that NMR interactions such as RDCs and spin relaxation, on their own, do not carry the necessary information needed to determine conformational entropy at a useful level of precision. The framework introduced here provides a general approach for exploring degeneracies in ensemble determination for different types of experimental data.

Information content of long-range NMR data for the characterization of conformational heterogeneity

Long-range NMR data, namely residual dipolar couplings (RDCs) from external alignment and paramagnetic data, are becoming increasingly popular for the characterization of conformational heterogeneity of multidomain biomacromolecules and protein complexes. The question addressed here is how much information is contained in these averaged data. We have analyzed and compared the information content of conformationally averaged RDCs caused by steric alignment and of both RDCs and pseudocontact shifts caused by paramagnetic alignment, and found that, despite the substantial differences, they contain a similar amount of information. Furthermore, using several synthetic tests we find that both sets of data are equally good towards recovering the major state(s) in conformational distributions.

Molecular Dynamics of Biomolecules through Direct Analysis of Dipolar Couplings

Residual dipolar couplings (RDCs) are important probes in structural biology, but their analysis is often complicated by the determination of an alignment tensor or its associated assumptions. We here apply the maximum entropy principle to derive a tensor-free formalism which allows for direct, dynamic analysis of RDCs and holds the classic tensor formalism as a special case. Specifically, the framework enables us to robustly analyze data regardless of whether a clear separation of internal and overall dynamics is possible. Such a separation is often difficult in the core subjects of current structural biology, which include multidomain and intrinsically disordered proteins as well as nucleic acids. We demonstrate the method is tractable and self-consistent and generalizes to data sets comprised of observations from multiple different alignment conditions.

Paramagnetic relaxation enhancement of membrane proteins by incorporation of the metal-chelating unnatural amino acid 2-amino-3-(8-hydroxyquinolin-3-yl)propanoic acid (HQA)

The use of paramagnetic constraints in protein NMR is an active area of research because of the benefits of long-range distance measurements (>10 Å). One of the main issues in successful execution is the incorporation of a paramagnetic metal ion into diamagnetic proteins. The most common metal ion tags are relatively long aliphatic chains attached to the side chain of a selected cysteine residue with a chelating group at the end where it can undergo substantial internal motions, decreasing the accuracy of the method. An attractive alternative approach is to incorporate an unnatural amino acid that binds metal ions at a specific site on the protein using the methods of molecular biology. Here we describe the successful incorporation of the unnatural amino acid 2-amino-3-(8-hydroxyquinolin-3-yl)propanoic acid (HQA) into two different membrane proteins by heterologous expression in E. coli. Fluorescence and NMR experiments demonstrate complete replacement of the natural amino acid with HQA and stable metal chelation by the mutated proteins. Evidence of site-specific intra- and inter-molecular PREs by NMR in micelle solutions sets the stage for the use of HQA incorporation in solid-state NMR structure determinations of membrane proteins in phospholipid bilayers.

Integral membrane protein structure determination using pseudocontact shifts

Obtaining enough experimental restraints can be a limiting factor in the NMR structure determination of larger proteins. This is particularly the case for large assemblies such as membrane proteins that have been solubilized in a membrane-mimicking environment. Whilst in such cases extensive deuteration strategies are regularly utilised with the aim to improve the spectral quality, these schemes often limit the number of NOEs obtainable, making complementary strategies highly beneficial for successful structure elucidation. Recently, lanthanide-induced pseudocontact shifts (PCSs) have been established as a structural tool for globular proteins. Here, we demonstrate that a PCS-based approach can be successfully applied for the structure determination of integral membrane proteins. Using the 7TM α-helical microbial receptor pSRII, we show that PCS-derived restraints from lanthanide binding tags attached to four different positions of the protein facilitate the backbone structure determination when combined with a limited set of NOEs. In contrast, the same set of NOEs fails to determine the correct 3D fold. The latter situation is frequently encountered in polytopical α-helical membrane proteins and a PCS approach is thus suitable even for this particularly challenging class of membrane proteins. The ease of measuring PCSs makes this an attractive route for structure determination of large membrane proteins in general.

Using modern tools to probe the structure-function relationship of fatty acid synthases

Fatty acid biosynthesis is essential to life and represents one of the most conserved pathways in nature, preserving the same handful of chemical reactions across all species. Recent interest in the molecular details of the de novo fatty acid synthase (FAS) has been heightened by demand for renewable fuels and the emergence of multidrug-resistant bacterial strains. Central to FAS is the acyl carrier protein (ACP), a protein chaperone that shuttles the growing acyl chain between catalytic enzymes within the FAS. Human efforts to alter fatty acid biosynthesis for oil production, chemical feedstock, or antimicrobial purposes has been met with limited success, due in part to a lack of detailed molecular information behind the ACP-partner protein interactions inherent to the pathway. This review will focus on recently developed tools for the modification of ACP and analysis of protein-protein interactions, such as mechanism-based crosslinking, and the studies exploiting them. Discussion specific to each enzymatic domain will focus first on mechanism and known inhibitors, followed by available structures and known interactions with ACP. Although significant unknowns remain, new understandings of the intricacies of FAS point to future advances in manipulating this complex molecular factory.

NMR Methods for the Study of Instrinsically Disordered Proteins Structure, Dynamics, and Interactions: General Overview and Practical Guidelines

Thanks to recent improvements in NMR instrumentation, pulse sequence design, and sample preparation, a panoply of new NMR tools has become available for atomic resolution characterization of intrinsically disordered proteins (IDPs) that are optimized for the particular chemical and spectroscopic properties of these molecules. A wide range of NMR observables can now be measured on increasingly complex IDPs that report on their structural and dynamic properties in isolation, as part of a larger complex, or even inside an entire living cell. Herein we present basic NMR concepts, as well as optimised tools available for the study of IDPs in solution. In particular, the following sections are discussed hereafter: a short introduction to NMR spectroscopy and instrumentation (Sect. 3.1), the effect of order and disorder on NMR observables (Sect. 3.2), particular challenges and bottlenecks for NMR studies of IDPs (Sect. 3.3), 2D HN and CON NMR experiments: the fingerprint of an IDP (Sect. 3.4), tools for overcoming major bottlenecks of IDP NMR studies (Sect. 3.5), 13C detected experiments (Sect. 3.6), from 2D to 3D: from simple snapshots to site-resolved characterization of IDPs (Sect. 3.7), sequential NMR assignment: 3D experiments (Sect. 3.8), high-dimensional NMR experiments (nD, with n>3) (Sect. 3.9) and conclusions and perspectives (Sect. 3.10).

Ensemble Calculation for Intrinsically Disordered Proteins Using NMR Parameters

Intrinsically disordered proteins (IDPs) perform their function despite their lack of well-defined tertiary structure. Residual structure has been observed in IDPs, commonly described as transient/dynamic or expressed in terms of fractional populations. In order to understand how the protein primary sequence dictates the dynamic and structural properties of IDPs and in general to understand how IDPs function, atomic-level descriptions are needed. Nuclear magnetic resonance spectroscopy provides information about local and long-range structure in IDPs at amino acid specific resolution and can be used in combination with ensemble descriptions to represent the dynamic nature of IDPs. In this chapter we describe sample-and-select approaches for ensemble modelling of local structural propensities in IDPs with specific emphasis on validation of these ensembles.

A self-assembled, π-stacked complex as a finely-tunable magnetic aligner for biomolecular NMR applications

The 1D aggregates of π-stacked self-assembled complexes worked as magnetic aligners upon an analyte protein to observe RDC in NMR.

FANTEN: a new web-based interface for the analysis of magnetic anisotropy-induced NMR data

Pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) arising from the presence of paramagnetic metal ions in proteins as well as RDCs due to partial orientation induced by external orienting media are nowadays routinely measured as a part of the NMR characterization of biologically relevant systems. PCSs and RDCs are becoming more and more popular as restraints (1) to determine and/or refine protein structures in solution, (2) to monitor the extent of conformational heterogeneity in systems composed of rigid domains which can reorient with respect to one another, and (3) to obtain structural information in protein-protein complexes. The use of both PCSs and RDCs proceeds through the determination of the anisotropy tensors which are at the origin of these NMR observables. A new user-friendly web tool, called FANTEN (Finding ANisotropy TENsors), has been developed for the determination of the anisotropy tensors related to PCSs and RDCs and has been made freely available through the WeNMR ( http://fanten-enmr.cerm.unifi.it:8080 ) gateway. The program has many new features not available in other existing programs, among which the possibility of a joint analysis of several sets of PCS and RDC data and the possibility to perform rigid body minimizations.

Improvements to REDCRAFT: a software tool for simultaneous characterization of protein backbone structure and dynamics from residual dipolar couplings

Within the past two decades, there has been an increase in the acquisition of residual dipolar couplings (RDC) for investigations of biomolecular structures. Their use however is still not as widely adopted as the traditional methods of structure determination by NMR, despite their potential for extending the limits in studies that examine both the structure and dynamics of biomolecules. This is in part due to the difficulties associated with the analysis of this information-rich data type. The software analysis tool REDCRAFT was previously introduced to address some of these challenges. Here we describe and evaluate a number of additional features that have been incorporated in order to extend its computational and analytical capabilities. REDCRAFT's more traditional enhancements integrate a modified steric collision term, as well as structural refinement in the rotamer space. Other, non-traditional improvements include: the filtering of viable structures based on relative order tensor estimates, decimation of the conformational space based on structural similarity, and forward/reverse folding of proteins. Utilizing REDCRAFT's newest features we demonstrate de-novo folding of proteins 1D3Z and 1P7E to within less than 1.6 Å of the corresponding X-ray structures, using as many as four RDCs per residue and as little as two RDCs per residue, in two alignment media. We also show the successful folding of a structure to less than 1.6 Å of the X-ray structure using {C(i-1)-N(i), N(i)-H(i), and C(i-1)-H(i)} RDCs in one alignment medium, and only {N(i)-H(i)} in the second alignment medium (a set of data which can be collected on deuterated samples). The program is available for download from our website at http://ifestos.cse.sc.edu .

Insights into domain-domain motions in proteins and RNA from solution NMR

Many multidomain proteins and ribonucleic acids consist of domains that autonomously fold and that are linked together by flexible junctions. This architectural design allows domains to sample a wide range of positions with respect to one another, yet do so in a way that retains structural specificity, since the number of sampled conformations remains extremely small compared to the total conformations that would be sampled if the domains were connected by an infinitely long linker. This “tuned” flexibility in interdomain conformation is in turn used in many biochemical processes.

There is great interest in characterizing the dynamic properties of multidomain systems, and moving beyond conventional descriptions in terms of static structures, toward the characterization of population-weighted ensembles describing a distribution of many conformations sampled in solution. There is also great interest in understanding the design principles and underlying physical and chemical interactions that specify the nature of interdomain flexibility. NMR spectroscopy is one of the most powerful techniques for characterizing motions in complex biomolecules and has contributed greatly toward our basic understanding of dynamics in proteins and nucleic acids and its role in folding, recognition, and signaling.

Here, we review methods that have been developed in our laboratories to address these challenges. Our approaches are based on the ability of one domain of the molecule to self-align in a magnetic field, or to dominate the overall orientation of the molecule, so that the conformational freedom of other domains can be assessed by their degree of alignment induced by the aligned part. In turn, this self-alignment ability can be intrinsic or can be caused by tagging appropriate constructs to the molecule of interest. In general, self-alignment is due to magnetic susceptibility anisotropy. Nucleic acids with elongated helices have this feature, as well as several paramagnetic metal centers that can be found in, or attached to, a protein domain.

Multiple paramagnetic effects through a tagged reporter protein

Paramagnetic effects provide unique information about the structure and dynamics of biomolecules. We developed a method in which the lanthanoid tag is not directly attached to the protein of interest, but instead to a "reporter" protein, which binds and then transmits paramagnetic information to the target. The designed method allows access to a large number of paramagnetic restraints and residual dipolar couplings produced from independent molecular alignments in high-molecular-weight proteins with unknown 3D structure.