Measurement of eight scalar and dipolar couplings for methine-methylene pairs in proteins and nucleic acids

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide ¹⁵N-¹H RDCs, to ¹H-¹H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

The Exact NOE as an Alternative in Ensemble Structure Determination

The structure-function paradigm is increasingly replaced by the structure-dynamics-function paradigm. All protein activity is steered by the interplay between enthalpy and entropy. Conformational dynamics serves as a proxy of conformational entropy. Therefore, it is essential to study not only the average conformation but also the spatial sampling of a protein on all timescales. To this purpose, we have established a protocol for determining multiple-state ensembles of proteins based on exact nuclear Overhauser effects (eNOEs). We have recently extended our previously reported eNOE data set for the protein GB3 by a very large set of backbone and side-chain residual dipolar couplings and three-bond J couplings. Here, we demonstrate that at least four structural states are required to represent the complete data set by dissecting the contributions to the CYANA target function, which quantifies restraint violations in structure calculation. We present a four-state ensemble of GB3, which largely preserves the characteristics obtained from eNOEs only. Due to the abundance of the input data, the ensemble and χ(1) angles in particular are well suited for cross-validation of the input data and comparison to x-ray structures. Principal component analysis is used to automatically identify and validate relevant states of the ensembles. Overall, our findings suggest that eNOEs are a valuable alternative to traditional NMR probes in spatial elucidation of proteins.

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp(3) atom, and (iii) reduced 1-4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric-isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14-C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14-C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric-isothermal MD simulations 2-4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490-1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Compiled data set of exact NOE distance limits, residual dipolar couplings and scalar couplings for the protein GB3

We compiled an NMR data set consisting of exact nuclear Overhauser enhancement (eNOE) distance limits, residual dipolar couplings (RDCs) and scalar (J) couplings for GB3, which forms one of the largest and most diverse data set for structural characterization of a protein to date. All data have small experimental errors, which are carefully estimated. We use the data in the research article Vogeli et al., 2015, Complementarity and congruence between exact NOEs and traditional NMR probes for spatial decoding of protein dynamics, J. Struct. Biol., 191, 3, 306–317, doi:10.1016/j.jsb.2015.07.008 [1] for cross-validation in multiple-state structural ensemble calculation. We advocate this set to be an ideal test case for molecular dynamics simulations and structure calculations.

ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB

Molecular mechanics is powerful for its speed in atomistic simulations, but an accurate force field is required. The Amber ff99SB force field improved protein secondary structure balance and dynamics from earlier force fields like ff99, but weaknesses in side chain rotamer and backbone secondary structure preferences have been identified. Here, we performed a complete refit of all amino acid side chain dihedral parameters, which had been carried over from ff94. The training set of conformations included multidimensional dihedral scans designed to improve transferability of the parameters. Improvement in all amino acids was obtained as compared to ff99SB. Parameters were also generated for alternate protonation states of ionizable side chains. Average errors in relative energies of pairs of conformations were under 1.0 kcal/mol as compared to QM, reduced 35% from ff99SB. We also took the opportunity to make empirical adjustments to the protein backbone dihedral parameters as compared to ff99SB. Multiple small adjustments of φ and ψ parameters were tested against NMR scalar coupling data and secondary structure content for short peptides. The best results were obtained from a physically motivated adjustment to the φ rotational profile that compensates for lack of ff99SB QM training data in the β-ppII transition region. Together, these backbone and side chain modifications (hereafter called ff14SB) not only better reproduced their benchmarks, but also improved secondary structure content in small peptides and reproduction of NMR χ1 scalar coupling measurements for proteins in solution. We also discuss the Amber ff12SB parameter set, a preliminary version of ff14SB that includes most of its improvements.

Stereospecific assignments in proteins using exact NOEs

Recently developed methods to measure distances in proteins with high accuracy by "exact" nuclear Overhauser effects (eNOEs) make it possible to determine stereospecific assignments, which are particularly important to fully exploit the accuracy of the eNOE distance measurements. Stereospecific assignments are determined by comparing the eNOE-derived distances to protein structure bundles calculated without stereospecific assignments, or an independently determined crystal structure. The absolute and relative CYANA target function difference upon swapping the stereospecific assignment of a diastereotopic group yields the respective stereospecific assignment. We applied the method to the eNOE data set that has recently been obtained for the third immunoglobulin-binding domain of protein G (GB3). The 884 eNOEs provide relevant data for 47 of the total of 75 diastereotopic groups. Stereospecific assignments could be established for 45 diastereotopic groups (96 %) using the X-ray structure, or for 27 diastereotopic groups (57 %) using structures calculated with the eNOE data set without stereospecific assignments, all of which are in agreement with those determined previously. The latter case is relevant for structure determinations based on eNOEs. The accuracy of the eNOE distance measurements is crucial for making stereospecific assignments because applying the same method to the traditional NOE data set for GB3 with imprecise upper distance bounds yields only 13 correct stereospecific assignments using the X-ray structure or 2 correct stereospecific assignments using NMR structures calculated without stereospecific assignments.

The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins

Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

Improved side-chain torsion potentials for the Amber ff99SB protein force field

Recent advances in hardware and software have enabled increasingly long molecular dynamics (MD) simulations of biomolecules, exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, we further advance simulation accuracy by improving the amino acid side-chain torsion potentials of the Amber ff99SB force field. First, we used simulations of model alpha-helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, we optimized the side-chain torsion potentials of these residues to match new, high-level quantum-mechanical calculations. Finally, we used microsecond-timescale MD simulations in explicit solvent to validate the resulting force field against a large set of experimental NMR measurements that directly probe side-chain conformations. The new force field, which we have termed Amber ff99SB-ILDN, exhibits considerably better agreement with the NMR data. Proteins 2010. © 2010 Wiley-Liss, Inc.

Side chain: backbone projections in aromatic and ASX residues from NMR cross-correlated relaxation

The measurements of cross-correlated relaxation rates between H(N)-N and C(beta)-C(gamma) intraresidual and sequential dipolar interactions is demonstrated in ASN, ASP and aromatic residues. The experiment can be used for deuterated samples and no additional knowledge such as Karplus parametrizations is required for the analysis. The data constitutes a new type of information since no other method relates the C(beta)-C(gamma) bond to H(N)-N. Using this method the dominant populations of rotamer states of chi 1 can be readily cross checked provided that phi or psi are known. In addition, dynamics on all timescales can be probed. As opposed to standard dynamics analysis of isolated bonds, the presented observables depend on relative dynamics with an interesting prospect to analyze correlated fluctuations of the two torsion angles phi or psi with chi 1. Experimental rates are compared to single conformer and ensemble representations of GB3 and ubiquitin. In particular, it is found that the recently published ubiquitin ensemble 2k39 improves the agreement obtained for 1UBQ. In general, however, input data restricting ASX and aromatic side chains in structure calculation is sparse highlighting the need for new NMR observables.

An efficient protocol for the complete incorporation of methyl-protonated alanine in perdeuterated protein

A strategy for the introduction of ((1)H,(13)C-methyl)-alanine into perdeuterated proteins is described. Specific protonation of alanine methyl groups to a level of 95% can be achieved by overexpressing proteins in M9/D(2)O based bacterial growth medium supplemented with 800 mg/l of 2-[(2)H], 3-[(13)C] L: -alanine. However, though simple, this approach results in undesired, non-specific background labeling due to isotope scrambling via different amino acid metabolic pathways. Following a careful analysis of known metabolic pathways we found that co-addition of perdeuterated forms of alpha-ketoisovalerate-d(7), succinate-d(4) and L: -isoleucine-d(10) with labeled L: -alanine, reduces undesired background labeling to <1%. When combined with recently developed methyl TROSY experiments, this methyl-specific labeling protocol permits the acquisition of excellent quality correlation spectra of alanine methyl groups in high molecular weight proteins. Our cost effective strategy offers a significant enhancement in the level of incorporation of methyl-labeled alanine in overexpressed proteins over previously reported methods.

Solution structure determination of oligoureas using methylene spin state selective NMR at 13C natural abundance

Ability of N,N'-linked oligoureas containing proteinogenic side chains to adopt a stable helix conformation in solution has been described recently. NMR as well as circular dichroism (CD) spectroscopies were employed to gain insight into their specific fold. It is herein proposed to extend the structural information available on these peptidomimetics by an advantageous use of a methylene spin state selective NMR experiment. Homodecoupling provided by the pulse scheme made it possible to readily measure conformation-dependent (3)J(HH) constants that are difficult if not impossible to obtain with standard NMR experiments. Adding those couplings to the NMR restraints improved the quality of the structure calculations significantly, as judged by a ca 30% decrease of the root mean square deviation (RMSD) obtained over an ensemble of 20 structures. Moreover, accurate determination of individual (1)J(CH) couplings within each methylene group revealed uniform values throughout the oligourea sequence, with (1)J(CH) systematically slightly larger for the pro-S hydrogen than for the pro-R. As shown in this study, the methylene spin state selective NMR experiment displays a good intrinsic sensitivity and could therefore provide valuable structural information at (13)C natural abundance for peptidomimetic molecules and foldamers bearing diastereotopic methylene protons.

Liquid-crystal NMR structure of HIV TAR RNA bound to its SELEX RNA aptamer reveals the origins of the high stability of the complex

Transactivation-response element (TAR) is a stable stem-loop structure of HIV RNA, which plays a crucial role during the life cycle of the virus. The apical loop of TAR acts as a binding site for essential cellular cofactors required for the replication of HIV. High-affinity aptamers directed against the apical loop of TAR have been identified by the SELEX approach. The RNA aptamers with the highest affinity for TAR fold as hairpins and form kissing complexes with the targeted RNA through loop-loop interactions. The aptamers with the strongest binding properties all possess a GA base pair combination at the loop-closing position. Using liquid-crystal NMR methodology, we have obtained a structural model in solution of a TAR-aptamer kissing complex with an unprecedented accuracy. This high-resolution structure reveals that the GA base pair is unilaterally shifted toward the 5' strand and is stabilized by a network of intersugar hydrogen bonds. This specific conformation of the GA base pair allows for the formation of two supplementary stable base-pair interactions. By systematic permutations of the loop-closing base pair, we establish that the identified atomic interactions, which form the basis for the high stability of the complex, are maintained in several other kissing complexes. This study rationalizes the stabilizing role of the loop-closing GA base pairs in kissing complexes and may help the development or improvement of drugs against RNA loops of viruses or pathogens as well as the conception of biochemical tools targeting RNA hairpins involved in important biological functions.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

P.E.HSQC: a simple experiment for simultaneous and sign-sensitive measurement of (1JCH+DCH) and (2JHH+DHH) couplings

The angular information content of residual dipolar couplings between nuclei of fixed distance makes the accurate and sign-sensitive measurement of (1JCH+DCH) and (2JHH+DHH) couplings highly desirable. Experiments published so far are typically highly specialized for the effective measurement of a subset of couplings. The P.E.HSQC presented here, is an E.COSY based experiment which allows the simultaneous measurement of all heteronuclear and homonuclear couplings within CH, CH2, and CH3 groups in a single spectrum with the necessary precision and sign information. The simplicity of the approach and the absence of artefacts like phase distortions due to antiphase evolution make it ideally suited for coupling determination of organic molecules at natural abundance.

13C direct detected NMR increases the detectability of residual dipolar couplings

13C direct detection is becoming an increasingly efficient approach to identify signals of residues that escape detection in 1H detected experiments. Pulse sequences have been developed to obtain 1H partially recoupled experiments for the measurement of the 1JHalphaCalpha and 1JHN couplings with the same resolution available in conventional 1H detected experiments. A consistent set of backbone rdc obtained without any 1H-based experiment has been obtained and shown to be effective for protein solution structure determination.

Simultaneous determination of one- and two-bond scalar and residual dipolar couplings between 13C', 13Calpha and 15N spins in proteins

Four simple and sensitive HNCO-based methods for measurement of 1J(C'Calpha), 1J(NCalpha) and 2J(NCalpha) coupling constants in protein main chains are presented. Three of these experiments enable the simultaneous measurement of 1J(C'Calpha), 1J(NCalpha) and 2J(NCalpha) couplings. Exploitation of the E.COSY principle provides excellent dispersion of cross peaks in the resulting 3D spectra. The couplings can be retrieved with good accuracy from peak-to-peak separations. Karplus parameterizations are provided for 1J(NCalpha) and 2J(NCalpha), obtained from a nearly complete set of couplings of human ubiquitin. In addition, feasibility of the proposed methodology for measuring several residual dipolar couplings (RDCs) simultaneously is assessed.

Optimum spin-state selection for all multiplicities in the acquisition dimension of the HSQC experiment

Most conventional heteronuclear spin-state-selective (S(3)) NMR experiments only work for a specific multiplicity, typically IS spin systems. Here, we introduce a general and efficient IPAP strategy to achieve S(3) editing simultaneously for all multiplicities in the acquisition dimension of the HSQC experiment. Complementary in-phase (HSQC-IP) and anti-phase (HSQC-AP) data are separately recorded with a simple phase exchange of two 90 degrees proton pulses involved in the mixing process of the F2-coupled sensitivity-improved HSQC pulse sequence. Additive and subtractive linear combination of these IP/AP data generates simplified F2-alpha/beta-spin-edited HSQC subspectra for all IS, I(2)S, and I(3)S spin systems and combines enhanced and optimized sensitivity with excellent tolerance to unwanted cross-talk contributions over a considerable range of coupling constants. Practical aspects such as pulse phase settings, transfer efficiency dependence, inter-pulse delay optimization, and percentage of cross-talk are theoretically analyzed and discussed as a function of each I(n)S multiplicity. Particular emphasis on the features associated to spin-editing in diastereotopic I(2)S spin systems and application to the measurement of long-range proton-carbon coupling constants are also provided.

Measurement of one-bond 13Calpha-1Halpha residual dipolar coupling constants in proteins by selective manipulation of CalphaHalpha spins

We have developed new 2D and 3D experiments for the measurement of C(alpha)-H(alpha) residual dipolar coupling constants in (13)C and (15)N labelled proteins. Two experiments, 2D (HNCO)-(J-CA)NH and 3D (HN)CO-(J-CA)NH, sample the C(alpha)-H(alpha) splitting by means of C(alpha) magnetization, while 2D (J-HACACO)NH and 3D J-HA(CACO)NH use H(alpha) magnetization to achieve a similar result. In the 2D experiments the coupling evolution is superimposed on the evolution of the (15)N chemical shifts and the IPAP principle is used to obtain (1)H-(15)N HSQC-like spectra from which the splitting is determined. The use of a third dimension in 3D experiments reduces spectral overlap to the point where use of an IPAP scheme may not be necessary. The length of the sampling interval in the J-dimension of these experiments is dictated solely by the relaxation properties of C(alpha) or H(alpha) nuclei. This was made possible by the use of C(alpha) selective pulses in combination with either a DPFGSE or modified BIRD pulses. Inclusion of these pulse sequence elements in the J-evolution periods removes unwanted spin-spin interactions. This allows prolonged sampling periods ( approximately 25 ms) yielding higher precision C(alpha)-H(alpha) splitting determination than is achievable with existing frequency based methods.

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.