The effects of residual dipolar coupling on carnosine in proton muscle spectra

Carnosine, an MR-visible dipeptide in human muscle, is well characterized by two peaks at ~8 and ~7 ppm from C2 and C4 imidazole protons. Like creatine and other metabolites, carnosine is subject to residual dipolar coupling in the anisotropic environment of muscle fibers, but the effects have not been studied extensively. Single-voxel TE 30-32 PRESS spectra from three different 3T studies were acquired from gastrocnemius medialis and soleus muscles in the human lower leg. In these studies, carnosine T2 values were measured, and spectra were obtained at three different foot angles. LCModel was used to fit the carnosine peaks with a basis set that was generated using shaped RF pulses and included a range of dipolar couplings affecting the C4 peak. A seven-parameter analytic expression was used to fit the CH2 doublets of creatine. It incorporated an optimized "effective TE" value to model the effect of shaped RF pulses. The fits confirm that the triplet C4 peak of carnosine is dipolar coupled to a pair of CH2 protons, with no need to include a contribution from a separate pool of freely rotating uncoupled carnosine. Moreover, the couplings experienced by carnosine C4 protons and creatine CH2 protons are strongly correlated (R2  =  0.88, P<0.001), exhibiting a similar 3cos2 θ - 1 dependence on the angle θ between fiber orientation and B0. T2 values for the singlet C2 peak of gastrocnemius carnosine are inversely proportional to the C4 dipolar coupling strength (R2  = 0.97, P < 0.001), which in turn is a function of foot orientation. This dependence indicates that careful positioning of the foot while acquiring lower leg muscle spectra is important to obtain reproducible carnosine concentrations. As proton magnetic resonance spectroscopy of carnosine is currently used to non-invasively estimate the muscle fiber typology, these results have important implications in sport science.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

Residual Dipolar Coupling Based Conformational Analysis Allows the Configurational Assessment of Steroids with up to Eight Stereocenters

Residual dipolar couplings (RDCs) induced by anisotropic media have been proved as a powerful tool for the structure elucidation of organic molecules in solution in nuclear magnetic resonance (NMR) based analysis. The value of dipolar couplings to solve complex conformational and configurational problems represents indeed an appealing analytical tool for the pharmaceutical industry particularly focusing on the stereochemistry characterization of NCEs since the early phase of the drug development process. In our work, RDCs were used for the conformational and configurational study of synthetic steroids with multiple stereocenters - prednisone and beclomethasone dipropionate (BDP) -. For both molecules the correct relative configuration was identified among all the possible diastereoisomers (32 and 128 respectively) arising from the compounds stereogenic carbons. Only for prednisone the use of additional experimental data (i. e. rOes) was necessary to resolve the right stereochemical structure.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

Cross-Linked Poly-4-Acrylomorpholine: A Flexible and Reversibly Compressible Aligning Gel for Anisotropic NMR Analysis of Peptides and Small Molecules in Water

Determination of the solution conformation of both small organic molecules and peptides in water remains a substantial hurdle in using NMR solution conformations to guide drug design due to the lack of easy to use alignment media. Herein we report the design of a flexible compressible chemically cross-linked poly-4-acrylomorpholine gel that can be used for the alignment of both small molecules and cyclic peptides in water. To test the new gel, residual dipolar couplings (RDCs) and J-coupling constants were used in the configurational analysis of strychnine hydrochloride, a molecule that has been studied extensively in organic solvents as well as a small cyclic peptide that is known to form an α-helix in water. The conformational ensembles for each molecule with the best fit to the data are reported. Identification of minor conformers in water that cannot easily be determined by conventional NOE measurements will facilitate the use of RDC experiments in structure-based drug design.

Unveiling the "invisible" druggable conformations of GDP-bound inactive Ras

The prevalent view on whether Ras is druggable has gradually changed in the recent decade with the discovery of effective inhibitors binding to cryptic sites unseen in the native structures. Despite the promising advances, therapeutics development toward higher potency and specificity is challenged by the elusive nature of these binding pockets. Here we derive a conformational ensemble of guanosine diphosphate (GDP)-bound inactive Ras by integrating spin relaxation-validated atomistic simulation with NMR chemical shifts and residual dipolar couplings, which provides a quantitative delineation of the intrinsic dynamics up to the microsecond timescale. The experimentally informed ensemble unequivocally demonstrates the preformation of both surface-exposed and buried cryptic sites in Ras•GDP, advocating design of inhibition by targeting the transient druggable conformers that are invisible to conventional experimental methods. The viability of the ensemble-based rational design has been established by retrospective testing of the ability of the Ras•GDP ensemble to identify known ligands from decoys in virtual screening.

Conformational Space Sampled by Domain Reorientation of Linear Diubiquitin Reflected in Its Binding Mode for Target Proteins

Linear polyubiquitin chains regulate diverse signaling proteins, in which the chains adopt various conformations to recognize different target proteins. Thus, the structural plasticity of the chains plays an important role in controlling the binding events. Herein, paramagnetic NMR spectroscopy is employed to explore the conformational space sampled by linear diubiquitin, a minimal unit of linear polyubiquitin, in its free state. Rigorous analysis of the data suggests that, regarding the relative positions of the ubiquitin units, particular regions of conformational space are preferentially sampled by the molecule. By combining these results with further data collected for charge-reversal derivatives of linear diubiquitin, structural insights into the factors underlying the binding events of linear diubiquitin are obtained.

Practical aspects of oligopeptide AAKLVFF as an alignment medium for the measurements of residual dipolar coupling of organic molecules

Practical aspects of the oligopeptide AAKLVFF as an alignment medium are discussed, including large-scale synthesis of the oligopeptide, detailed description of preparation of the alignment medium, and acquisition of the RDCs. The resulting orienting medium is stable and highly homogeneous with tunable alignment strength in methanol.

Probing Long-Range Anisotropic Interactions: a General and Sign-Sensitive Strategy to Measure 1 H-1 H Residual Dipolar Couplings as a Key Advance for Organic Structure Determination

Residual dipolar couplings (RDCs) are amongst the most powerful NMR parameters for organic structure elucidation. In order to maximize their effectiveness in increasingly complex cases such as flexible compounds, a maximum of RDCs between nuclei sampling a large distribution of orientations is needed, including sign information. For this, the easily accessible one-bond ¹ H-¹³ C RDCs alone often fall short. Long-range ¹ H-¹ H RDCs are both abundant and typically sample highly complementary orientations, but accessing them in a sign-sensitive way has been severely obstructed due to the overflow of ¹ H-¹ H couplings. Here, we present a generally applicable strategy that allows the measurement of a large number of ¹ H-¹ H RDCs, including their signs, which is based on a combination of an improved PSYCHEDELIC method and a new selective constant-time β-COSY experiment. The potential of ¹ H-¹ H RDCs to better determine molecular alignment and to discriminate between enantiomers and diastereomers is demonstrated.

Protein structure prediction assisted with sparse NMR data in CASP13

CASP13 has investigated the impact of sparse NMR data on the accuracy of protein structure prediction. NOESY and ¹⁵ N-¹ H residual dipolar coupling data, typical of that obtained for ¹⁵ N,¹³ C-enriched, perdeuterated proteins up to about 40 kDa, were simulated for 11 CASP13 targets ranging in size from 80 to 326 residues. For several targets, two prediction groups generated models that are more accurate than those produced using baseline methods. Real NMR data collected for a de novo designed protein were also provided to predictors, including one data set in which only backbone resonance assignments were available. Some NMR-assisted prediction groups also did very well with these data. CASP13 also assessed whether incorporation of sparse NMR data improves the accuracy of protein structure prediction relative to nonassisted regular methods. In most cases, incorporation of sparse, noisy NMR data results in models with higher accuracy. The best NMR-assisted models were also compared with the best regular predictions of any CASP13 group for the same target. For six of 13 targets, the most accurate model provided by any NMR-assisted prediction group was more accurate than the most accurate model provided by any regular prediction group; however, for the remaining seven targets, one or more regular prediction method provided a more accurate model than even the best NMR-assisted model. These results suggest a novel approach for protein structure determination, in which advanced prediction methods are first used to generate structural models, and sparse NMR data is then used to validate and/or refine these models.

NMR Analysis of Apo Glutamine-Binding Protein Exposes Challenges in the Study of Interdomain Dynamics

Glutamine-binding protein (GlnBP) displays an apo, "open" and a holo, "closed" crystal form, mutually related by a rigid-body reorientation of its domains. A fundamental question about such large-scale conformational transitions, whether the closed state exists in the absence of ligand, is controversial in the case of GlnBP. NMR observations have indicated no evidence of the closed form, whereas experimentally validated computations have suggested a remarkable ca. 40 % population. Herein, a paramagnetic NMR strategy designed to detect the putative apo-closed species shows that a major population of the latter is highly improbable. Further, NMR residual dipolar couplings collected under three anisotropic conditions do not reveal differential domain alignment and establish that the average solution conformation is satisfied by the apo-open crystal structure. Our results indicate that the computational prediction of large-scale interdomain motions is not trivial and may lead to erroneous conclusions without proper experimental validation.

Continuous Interdomain Orientation Distributions Reveal Components of Binding Thermodynamics

The flexibility of biological macromolecules is an important structural determinant of function. Unfortunately, the correlations between different motional modes are poorly captured by discrete ensemble representations. Here, we present new ways to both represent and visualize correlated interdomain motions. Interdomain motions are determined directly from residual dipolar couplings, represented as a continuous conformational distribution, and visualized using the disk-on-sphere representation. Using the disk-on-sphere representation, features of interdomain motions, including correlations, are intuitively visualized. The representation works especially well for multidomain systems with broad conformational distributions.This analysis also can be extended to multiple probability density modes, using a Bingham mixture model. We use this new paradigm to study the interdomain motions of staphylococcal protein A, which is a key virulence factor contributing to the pathogenicity of Staphylococcus aureus. We capture the smooth transitions between important states and demonstrate the utility of continuous distribution functions for computing the reorientational components of binding thermodynamics. Such insights allow for the dissection of the dynamic structural components of functionally important intermolecular interactions.

Solution structure of a reverse transcriptase recognition site of a LINE RNA from zebrafish

Long interspersed nuclear element (LINE) is known to be transposed by reverse transcription using its RNA transcript. Recognition of the 3' stem-loop of LINE RNA by its reverse transcriptase (RT) is an important step of the retrotransposition. Our previous study revealed that the second G residue (G8) in the GGAUA loop of a 17mer LINE RNA from eel, UnaL2-17, is recognized by its RT and the U residue (U10) in the same loop is required to maintain the loop structure (Baba S, Kajikawa M, Okada N, Kawai G. Solution structure of an RNA stem-loop derived from the 3' conserved region of eel LINE UnaL2. RNA 2004;10:1380-1387). ZfL2-2, a LINE from zebrafish, has the same 3' stem-loop with UnaL2 and ZfL2-1 has similar but distinct 3' stem-loop with an insertion which can form an additional stem-loop. Here, we determined the solution structure of the 34mer RT recognition site of the LINE RNA (ZfL2-1-34). It was found that ZfL2-1-34 forms a hairpin with an internal loop, the tertiary structure of which is superimposed with that of ZfL2-2. It is noted that A10 and the inserted stem-loop, starting with A12, in ZfL2-1-34 located at the positions corresponding to those of G8 and U10, respectively, in UnaL2-17. These results strongly suggest that the two LINEs share the similar recognition mechanism and the A10 in ZfL2-1-34 is the determinant recognized by its RT.

A Self-Assembled Oligopeptide as a Versatile NMR Alignment Medium for the Measurement of Residual Dipolar Couplings in Methanol

Residual dipolar coupling (RDC) is a powerful structural parameter for the determination of the constitution, conformation, and configuration of organic molecules. Herein, we report the first liquid crystal-based orienting medium that is compatible with MeOH, thus enabling RDC acquisitions of a wide range of intermediate to polar organic molecules. The liquid crystals were produced from self-assembled oligopeptide nanotubes (AAKLVFF), which are stable at very low concentrations. The presented alignment medium is highly homogeneous, and the size of RDCs can be scaled with the concentration of the peptide. To assess the accuracy of the RDC measurement by employing this new medium, seven bioactive natural products from different classes were chosen and analyzed. The straightforward preparation of the anisotropic alignment sample will offer a versatile and robust protocol for the routine RDC measurement of natural products.

Conformational Heterogeneity in the Activation Mechanism of Bax

Bax is known for its pro-apoptotic role within the mitochondrial pathway of apoptosis. However, the mechanism for transitioning Bax from cytosolic to membrane-bound oligomer remains elusive. Previous nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) studies defined monomeric Bax as conformationally homogeneous. Yet it has recently been proposed that monomeric Bax exists in equilibrium with a minor state that is distinctly different from its NMR structure. Here, we revisited the structural analysis of Bax using methods uniquely suited for unveiling "invisible" states of proteins, namely, NMR paramagnetic relaxation enhancements and EPR double electron-electron resonance (DEER). Additionally we examined the effect of glycerol, the co-solvent of choice in DEER studies, on the structure of Bax using NMR chemical-shift perturbations and residual dipolar couplings. Based on our combined NMR and EPR results, Bax is a conformationally homogeneous protein prior to its activation.

Direct Experimental Evidence for Halogen-Aryl π Interactions in Solution from Molecular Torsion Balances

We dissected halogen-aryl π interactions experimentally using a bicyclic N-arylimide based molecular torsion balances system, which is based on the influence of the non-bonded interaction on the equilibria between folded and unfolded states. Through comparison of balances modulated by higher halogens with fluorine balances, we determined the magnitude of the halogen-aryl π interactions in our unimolecular systems to be larger than -5.0 kJ mol^-1 , which is comparable with the magnitude estimated in the biomolecular systems. Our study provides direct experimental evidence of halogen-aryl π interactions in solution, which until now have only been revealed in the solid state and evaluated theoretically by quantum-mechanical calculations.

Elucidation of the Enantiodiscrimination Properties of a Nonracemic Chiral Alignment Medium through Gel-based Capillary Electrochromatography: Separation of the Mefloquine Stereoisomers

Enantiodiscrimination and enantioseparation are two highly important processes in chemistry, often performed by using NMR spectroscopy and chromatography. For a better understanding of the mechanistic details, the same system should be studied by both methods. In addition, isotropic and anisotropic NMR parameters should be obtained, the latter using alignment media so that residual dipolar couplings and chemical-shift anisotropies can be measured. Consequently, a chiral alignment medium was used for the first time in chiral gel-based capillary electrochromatography with the four stereoisomers of the antimalaria drug mefloquine as test compounds. Chromatographic data verify that enantiodiscrimination obtained with this alignment gel is caused by differences in the equilibrium constants related to associate formation. Hence, the chromatographic separation provides physicochemical data that form a basis for the understanding and optimization of alignment processes, and vice versa.

A hybrid NMR/SAXS-based approach for discriminating oligomeric protein interfaces using Rosetta

Oligomeric proteins are important targets for structure determination in solution. While in most cases the fold of individual subunits can be determined experimentally, or predicted by homology-based methods, protein-protein interfaces are challenging to determine de novo using conventional NMR structure determination protocols. Here we focus on a member of the bet-V1 superfamily, Aha1 from Colwellia psychrerythraea. This family displays a broad range of crystallographic interfaces none of which can be reconciled with the NMR and SAXS data collected for Aha1. Unlike conventional methods relying on a dense network of experimental restraints, the sparse data are used to limit conformational search during optimization of a physically realistic energy function. This work highlights a new approach for studying minor conformational changes due to structural plasticity within a single dimeric interface in solution.

Conformation of inhibitor-free HIV-1 protease derived from NMR spectroscopy in a weakly oriented solution

Flexibility of the glycine-rich flaps is known to be essential for catalytic activity of the HIV-1 protease, but their exact conformations at the different stages of the enzymatic pathway remain subject to much debate. Although hundreds of crystal structures of protease-inhibitor complexes have been solved, only about a dozen inhibitor-free protease structures have been reported. These latter structures reveal a large diversity of flap conformations, ranging from closed to semi-open to wide open. To evaluate the average structure in solution, we measured residual dipolar couplings (RDCs) and compared these to values calculated for crystal structures representative of the closed, semi-open, and wide-open states. The RDC data clearly indicate that the inhibitor-free protease, on average, adopts a closed conformation in solution that is very similar to the inhibitor-bound state. By contrast, a highly drug-resistant protease mutant, PR20, adopts the wide-open flap conformation.

Measuring membrane protein bond orientations in nanodiscs via residual dipolar couplings

Membrane proteins are involved in numerous vital biological processes. To understand membrane protein functionality, accurate structural information is required. Usually, structure determination and dynamics of membrane proteins are studied in micelles using either solution state NMR or X-ray crystallography. Even though invaluable information has been obtained by this approach, micelles are known to be far from ideal mimics of biological membranes often causing the loss or decrease of membrane protein activity. Recently, nanodiscs, which are composed of a lipid bilayer surrounded by apolipoproteins, have been introduced as a more physiological alternative than micelles for NMR investigations on membrane proteins. Here, we show that membrane protein bond orientations in nanodiscs can be obtained by measuring residual dipolar couplings (RDCs) with the outer membrane protein OmpX embedded in nanodiscs using Pf1 phage as an alignment medium. The presented collection of membrane protein RDCs in nanodiscs represents an important step toward more comprehensive structural and dynamical NMR-based investigations of membrane proteins in a natural bilayer environment.

Indirect use of deuterium in solution NMR studies of protein structure and hydrogen bonding

A description of the utility of deuteration in protein NMR is provided with an emphasis on quantitative evaluation of the effects of deuteration on a number of NMR parameters of proteins: (1) chemical shifts, (2) scalar coupling constants, (3) relaxation properties (R1 and R2 rates) of nuclei directly attached to one or more deuterons as well as protons of methyl groups in a highly deuterated environment, (4) scalar relaxation of 15N and 13C nuclei in 15N-D and 13C-D spin systems as a measure of hydrogen bonding strength, and (5) NOE-based applications of deuteration in NMR studies of protein structure. The discussion is restricted to the 'indirect' use of deuterium in the sense that the description of NMR parameters and properties of the nuclei affected by nearby deuterons (15N, 13C, 1H) is provided rather than those of deuterium itself.

Partial alignment and measurement of residual dipolar couplings of proteins under high hydrostatic pressure

High-pressure NMR spectroscopy has emerged as a complementary approach for investigating various structural and thermodynamic properties of macromolecules. Noticeably absent from the array of experimental restraints that have been employed to characterize protein structures at high hydrostatic pressure is the residual dipolar coupling, which requires the partial alignment of the macromolecule of interest. Here we examine five alignment media that are commonly used at ambient pressure for this purpose. We find that the spontaneous alignment of Pf1 phage, d(GpG) and a C12E5/n-hexnanol mixture in a magnetic field is preserved under high hydrostatic pressure. However, DMPC/DHPC bicelles and collagen gel are found to be unsuitable. Evidence is presented to demonstrate that pressure-induced structural changes can be identified using the residual dipolar coupling.

CHARMM36 all-atom additive protein force field: validation based on comparison to NMR data

Protein structure and dynamics can be characterized on the atomistic level with both nuclear magnetic resonance (NMR) experiments and molecular dynamics (MD) simulations. Here, we quantify the ability of the recently presented CHARMM36 (C36) force field (FF) to reproduce various NMR observables using MD simulations. The studied NMR properties include backbone scalar couplings across hydrogen bonds, residual dipolar couplings (RDCs) and relaxation order parameter, as well as scalar couplings, RDCs, and order parameters for side-chain amino- and methyl-containing groups. It is shown that the C36 FF leads to better correlation with experimental data compared to the CHARMM22/CMAP FF and suggest using C36 in protein simulations. Although both CHARMM FFs contains the same nonbond parameters, our results show how the changes in the internal parameters associated with the peptide backbone via CMAP and the χ1 and χ2 dihedral parameters leads to improved treatment of the analyzed nonbond interactions. This highlights the importance of proper treatment of the internal covalent components in modeling nonbond interactions with molecular mechanics FFs.

Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach

Amyloid β (Aβ) peptides are a primary component of fibrils and oligomers implicated in the etiology of Alzheimer's disease (AD). However, the intrinsic flexibility of these peptides has frustrated efforts to investigate the secondary and tertiary structure of Aβ monomers, whose conformational landscapes directly contribute to the kinetics and thermodynamics of Aβ aggregation. In this work, de novo replica exchange molecular dynamics (REMD) simulations on the microseconds-per-replica timescale are used to characterize the structural ensembles of Aβ42, Aβ40, and M35-oxidized Aβ42, three physiologically relevant isoforms with substantially different aggregation properties. J-coupling data calculated from the REMD trajectories were compared to corresponding NMR-derived values acquired through two different pulse sequences, revealing that all simulations converge on the order of hundreds of nanoseconds-per-replica toward ensembles that yield good agreement with experiment. Though all three Aβ species adopt highly heterogeneous ensembles, these are considerably more structured compared to simulations on shorter timescales. Prominent in the C-terminus are antiparallel β-hairpins between L17-A21, A30-L36, and V39-I41, similar to oligomer and fibril intrapeptide models that expose these hydrophobic side chains to solvent and may serve as hotspots for self-association. Compared to reduced Aβ42, the absence of a second β-hairpin in Aβ40 and the sampling of alternate β topologies by M35-oxidized Aβ42 may explain the reduced aggregation rates of these forms. A persistent V24-K28 bend motif, observed in all three species, is stabilized by buried backbone to side-chain hydrogen bonds with D23 and a cross-region salt bridge between E22 and K28, highlighting the role of the familial AD-linked E22 and D23 residues in Aβ monomer folding. These characterizations help illustrate the conformational landscapes of Aβ monomers at atomic resolution and provide insight into the early stages of Aβ aggregation pathways.

Solid-state ¹³C NMR reveals annealing of raft-like membranes containing cholesterol by the intrinsically disordered protein α-Synuclein

Misfolding and aggregation of the intrinsically disordered protein α-Synuclein (αS) in Lewy body plaques are characteristic markers of late-stage Parkinson's disease. It is well established that membrane binding is initiated at the N-terminus of the protein and affects biasing of conformational ensembles of αS. However, little is understood about the effect of αS on the membrane lipid bilayer. One hypothesis is that intrinsically disordered αS alters the structural properties of the membrane, thereby stabilizing the bilayer against fusion. Here, we used two-dimensional (13)C separated local-field NMR to study interaction of the wild-type α-Synuclein (wt-αS) or its N-terminal (1-25) amino acid sequence (N-αS) with a cholesterol-enriched ternary membrane system. This lipid bilayer mimics cellular raft-like domains in the brain that are proposed to be involved in neuronal membrane fusion. The two-dimensional dipolar-recoupling pulse sequence DROSS (dipolar recoupling on-axis with scaling and shape preservation) was implemented to measure isotropic (13)C chemical shifts and (13)C-(1)H residual dipolar couplings under magic-angle spinning. Site-specific changes in NMR chemical shifts and segmental order parameters indicate that both wt-αS and N-αS bind to the membrane interface and change lipid packing within raft-like membranes. Mean-torque modeling of (13)C-(1)H NMR order parameters shows that αS induces a remarkable thinning of the bilayer (≈6Å), accompanied by an increase in phospholipid cross-sectional area (≈10Å(2)). This perturbation is characterized as membrane annealing and entails structural remodeling of the raft-like liquid-ordered phase. We propose this process is implicated in regulation of synaptic membrane fusion that may be altered by aggregation of αS in Parkinson's disease.

Investigation of the Polymeric Properties of α-Synuclein and Comparison with NMR Experiments: A Replica Exchange Molecular Dynamics Study

Intrinsically disordered proteins (IDPs) have been shown to be involved in a number of cellular functions, in addition to their predominance in diseased states. α-synuclein may be described as one such IDP implicated in the pathology of Parkinson's disease. Understanding the conformational characteristics of the monomeric state of α-synuclein is necessary for understanding the role of the monomer conformation in aggregation. Polymer theories have been applied to investigate the statistical properties of homopolymeric IDPs. Here we use Replica Exchange Molecular Dynamics (REMD) simulations using temperature as a proxy for solvent quality to examine how well these theories developed for homopolymeric chains describe heteropolymeric α-synuclein. Our results indicate that α-synuclein behaves like a homopolymer at the extremes of solvent quality, while in the intermediate solvent regime, the uneven distribution of charged residues along the sequence strongly influences the conformations adopted by the chain. We refine the ensemble extracted from the REMD simulations of α-synuclein, which shows the best qualitative agreement with experiment, by fitting to the experimental NMR Residual Dipolar Couplings (RDCs) and Paramagnetic Relaxation Enhancements (PREs). Our results demonstrate that the detailed shape of the RDC patterns are sensitive to the angular correlations that are local in sequence while longer range anti-correlations which arise from packing constraints affect the RDC magnitudes.

Backbone resonance assignment and order tensor estimation using residual dipolar couplings

An NMR investigation of proteins with known X-ray structures is of interest in a number of endeavors. Performing these studies through nuclear magnetic resonance (NMR) requires the costly step of resonance assignment. The prevalent assignment strategy does not make use of existing structural information and requires uniform isotope labeling. Here we present a rapid and cost-effective method of assigning NMR data to an existing structure-either an X-ray or computationally modeled structure. The presented method, Exhaustively Permuted Assignment of RDCs (EPAR), utilizes unassigned residual dipolar coupling (RDC) data that can easily be obtained by NMR spectroscopy. The algorithm uses only the backbone N-H RDCs from multiple alignment media along with the amino acid type of the RDCs. It is inspired by previous work from Zweckstetter and provides several extensions. We present results on 13 synthetic and experimental datasets from 8 different structures, including two homodimers. Using just two alignment media, EPAR achieves an average assignment accuracy greater than 80%. With three media, the average accuracy is higher than 94%. The algorithm also outputs a prediction of the assignment accuracy, which has a correlation of 0.77 to the true accuracy. This prediction score can be used to establish the needed confidence in assignment accuracy.

Structural NMR of protein oligomers using hybrid methods

Solving structures of native oligomeric protein complexes using traditional high-resolution NMR techniques remains challenging. However, increased utilization of computational platforms, and integration of information from less traditional NMR techniques with data from other complementary biophysical methods, promises to extend the boundary of NMR-applicable targets. This article reviews several of the techniques capable of providing less traditional and complementary structural information. In particular, the use of orientational constraints coming from residual dipolar couplings and residual chemical shift anisotropy offsets are shown to simplify the construction of models for oligomeric complexes, especially in cases of weak homo-dimers. Combining this orientational information with interaction site information supplied by computation, chemical shift perturbation, paramagnetic surface perturbation, cross-saturation and mass spectrometry allows high resolution models of the complexes to be constructed with relative ease. Non-NMR techniques, such as mass spectrometry, EPR and small angle X-ray scattering, are also expected to play increasingly important roles by offering alternative methods of probing the overall shape of the complex. Computational platforms capable of integrating information from multiple sources in the modeling process are also discussed in the article. And finally a new, detailed example on the determination of a chemokine tetramer structure will be used to illustrate how a non-traditional approach to oligomeric structure determination works in practice.

Engineering encodable lanthanide-binding tags into loop regions of proteins

Lanthanide-binding tags (LBTs) are valuable tools for investigation of protein structure, function, and dynamics by NMR spectroscopy, X-ray crystallography, and luminescence studies. We have inserted LBTs into three different loop positions (denoted L, R, and S) of the model protein interleukin-1β (IL1β) and varied the length of the spacer between the LBT and the protein (denoted 1−3). Luminescence studies demonstrate that all nine constructs bind Tb3+ tightly in the low nanomolar range. No significant change in the fusion protein occurs from insertion of the LBT, as shown by two X-ray crystallographic structures of the IL1β-S1 and IL1β-L3 constructs and for the remaining constructs by comparing the 1H−15N heteronuclear single-quantum coherence NMR spectra with that of the wild-type IL1β. Additionally, binding of LBT-loop IL1β proteins to their native binding partner in vitro remains unaltered. X-ray crystallographic phasing was successful using only the signal from the bound lanthanide. Large residual dipolar couplings (RDCs) could be determined by NMR spectroscopy for all LBT-loop constructs and revealed that the LBT-2 series were rigidly incorporated into the interleukin-1β structure. The paramagnetic NMR spectra of loop-LBT mutant IL1β-R2 were assigned and the Δχ tensor components were calculated on the basis of RDCs and pseudocontact shifts. A structural model of the IL1β-R2 construct was calculated using the paramagnetic restraints. The current data provide support that encodable LBTs serve as versatile biophysical tags when inserted into loop regions of proteins of known structure or predicted via homology modeling.

Simultaneous structure and dynamics of a membrane protein using REDCRAFT: membrane-bound form of Pf1 coat protein

A strategy for simultaneous study of the structure and internal dynamics of a membrane protein is described using the REDCRAFT algorithm. The membrane-bound form of the Pf1 major coat protein (mbPf1) was used as an example. First, synthetic data is utilized to validate the simultaneous study of structure and dynamics with REDCRAFT using dihedral restraints and backbone N-H RDCs from two different alignments. Subsequently, the validated analysis is applied to experimental data and confirms that REDCRAFT produces meaningful structures from sparse RDC data. Furthermore, simulated data from a two-state jump motion is used to illustrate the necessity for simultaneous consideration of structure and dynamics. Disregarding internal dynamics during the course of structure determination is shown to produce an average-state that is not related to the two intermediate states. During the analysis of RDC data from the dynamic model, REDCRAFT appropriately identifies the region separating the static and dynamic domains of the protein. Finally, analysis of experimental data strongly suggests the existence of internal motion between the amphipathic and the transmembrane helices of the membrane-bound form of the protein. The ability to perform fragmented structure determination of each domain without a priori assumption of the order tensors allows an independent determination of the order tensors, which yields a more comprehensive description of protein structure and dynamics and is particularly relevant to the study of membrane proteins.

A device for the measurement of residual chemical shift anisotropy and residual dipolar coupling in soluble and membrane-associated proteins

Residual dipolar coupling (RDC) and residual chemical shift anisotropy (RCSA) report on orientational properties of a dipolar bond vector and a chemical shift anisotropy principal axis system, respectively. They can be highly complementary in the analysis of backbone structure and dynamics in proteins as RCSAs generally include a report on vectors out of a peptide plane while RDCs usually report on in-plane vectors. Both RDC and RCSA average to zero in isotropic solutions and require partial orientation in a magnetic field to become observable. While the alignment and measurement of RDC has become routine, that of RCSA is less common. This is partly due to difficulties in providing a suitable isotopic reference spectrum for the measurement of the small chemical shift offsets coming from RCSA. Here we introduce a device (modified NMR tube) specifically designed for accurate measurement of reference and aligned spectra for RCSA measurements, but with a capacity for RDC measurements as well. Applications to both soluble and membrane anchored proteins are illustrated.

Major groove width variations in RNA structures determined by NMR and impact of 13C residual chemical shift anisotropy and 1H-13C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by (1)H- (1)H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic (1)H- (13)C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.

The solution structure of the Mg2+ form of soybean calmodulin isoform 4 reveals unique features of plant calmodulins in resting cells

Soybean calmodulin isoform 4 (sCaM4) is a plant calcium-binding protein, regulating cellular responses to the second messenger Ca(2+). We have found that the metal ion free (apo-) form of sCaM4 possesses a half unfolded structure, with the N-terminal domain unfolded and the C-terminal domain folded. This result was unexpected as the apo-forms of both soybean calmodulin isoform 1 (sCaM1) and mammalian CaM (mCaM) are fully folded. Because of the fact that free Mg(2+) ions are always present at high concentrations in cells (0.5-2 mM), we suggest that Mg(2+) should be bound to sCaM4 in nonactivated cells. CD studies revealed that in the presence of Mg(2+) the initially unfolded N-terminal domain of sCaM4 folds into an alpha-helix-rich structure, similar to the Ca(2+) form. We have used the NMR backbone residual dipolar coupling restraints (1)D(NH), (1)D(C alpha H alpha), and (1)D(C'C alpha) to determine the solution structure of the N-terminal domain of Mg(2+)-sCaM4 (Mg(2+)-sCaM4-NT). Compared with the known structure of Ca(2+)-sCaM4, the structure of the Mg(2+)-sCaM4-NT does not fully open the hydrophobic pocket, which was further confirmed by the use of the fluorescent probe ANS. Tryptophan fluorescence experiments were used to study the interactions between Mg(2+)-sCaM4 and CaM-binding peptides derived from smooth muscle myosin light chain kinase and plant glutamate decarboxylase. These results suggest that Mg(2+)-sCaM4 does not bind to Ca(2+)-CaM target peptides and therefore is functionally similar to apo-mCaM. The Mg(2+)- and apo-structures of the sCaM4-NT provide unique insights into the structure and function of some plant calmodulins in resting cells.

Improvement and analysis of computational methods for prediction of residual dipolar couplings

We describe a new, computationally efficient method for computing the molecular alignment tensor based on the molecular shape. The increase in speed is achieved by re-expressing the problem as one of numerical integration, rather than a simple uniform sampling (as in the PALES method), and by using a convex hull rather than a detailed representation of the surface of a molecule. This method is applicable to bicelles, PEG/hexanol, and other alignment media that can be modeled by steric restrictions introduced by a planar barrier. This method is used to further explore and compare various representations of protein shape by an equivalent ellipsoid. We also examine the accuracy of the alignment tensor and residual dipolar couplings (RDC) prediction using various ab initio methods. We separately quantify the inaccuracy in RDC prediction caused by the inaccuracy in the orientation and in the magnitude of the alignment tensor, concluding that orientation accuracy is much more important in accurate prediction of RDCs.

Efficient and accurate estimation of relative order tensors from lambda-maps

The rapid increase in the availability of RDC data from multiple alignment media in recent years has necessitated the development of more sophisticated analyses that extract the RDC data's full information content. This article presents an analysis of the distribution of RDCs from two media (2D-RDC data), using the information obtained from a lambda-map. This article also introduces an efficient algorithm, which leverages these findings to extract the order tensors for each alignment medium using unassigned RDC data in the absence of any structural information. The results of applying this 2D-RDC analysis method to synthetic and experimental data are reported in this article. The relative order tensor estimates obtained from the 2D-RDC analysis are compared to order tensors obtained from the program REDCAT after using assignment and structural information. The final comparisons indicate that the relative order tensors estimated from the unassigned 2D-RDC method very closely match the results from methods that require assignment and structural information. The presented method is successful even in cases with small datasets. The results of analyzing experimental RDC data for the protein 1P7E are presented to demonstrate the potential of the presented work in accurately estimating the principal order parameters from RDC data that incompletely sample the RDC space. In addition to the new algorithm, a discussion of the uniqueness of the solutions is presented; no more than two clusters of distinct solutions have been shown to satisfy each lambda-map.

Improved accuracy of 15N-1H scalar and residual dipolar couplings from gradient-enhanced IPAP-HSQC experiments on protonated proteins

The presence of dipole-dipole cross-correlated relaxation as well as unresolved E.COSY effects adversely impacts the accuracy of (1)J(NH) splittings measured from gradient-enhanced IPAP-HSQC spectra. For isotropic samples, the size of the systematic errors caused by these effects depends on the values of (2)J(NHalpha), (3)J(NHbeta) and (3)J(HNHalpha). Insertion of band-selective (1)H decoupling pulses in the IPAP-HSQC experiment eliminates these systematic errors and for the protein GB3 yields (1)J(NH) splittings that agree to within a root-mean-square difference of 0.04 Hz with values measured for perdeuterated GB3. Accuracy of the method is also highlighted by a good fit to the GB3 structure of the (1)H-(15)N RDCs extracted from the minute differences in (1)J(NH) splitting measured at 500 and 750 MHz (1)H frequencies, resulting from magnetic susceptibility anisotropy. A nearly complete set of (2)J(NHalpha) couplings was measured in GB3 in order to evaluate whether the impact of cross-correlated relaxation is dominated by the (15)N-(1)H(alpha) or (15)N-(1)H(beta) dipolar interaction. As expected, we find that (2)J(NHalpha) < or = 2 Hz, with values in the alpha-helix (0.86 +/- 0.52 Hz) slightly larger than in beta-sheet (0.66 +/- 0.26 Hz). Results indicate that under isotropic conditions, N-H(N)/N-H(beta) cross-correlated relaxation often dominates. Unresolved E.COSY effects under isotropic conditions involve (3)J(HNHalpha) and J(NHalpha), but when weakly aligned any aliphatic proton proximate to both N and H(N) can contribute.

Estimation of relative order tensors, and reconstruction of vectors in space using unassigned RDC data and its application

Advances in NMR instrumentation and pulse sequence design have resulted in easier acquisition of Residual Dipolar Coupling (RDC) data. However, computational and theoretical analysis of this type of data has continued to challenge the international community of investigators because of their complexity and rich information content. Contemporary use of RDC data has required a-priori assignment, which significantly increases the overall cost of structural analysis. This article introduces a novel algorithm that utilizes unassigned RDC data acquired from multiple alignment media (nD-RDC, n3) for simultaneous extraction of the relative order tensor matrices and reconstruction of the interacting vectors in space. Estimation of the relative order tensors and reconstruction of the interacting vectors can be invaluable in a number of endeavors. An example application has been presented where the reconstructed vectors have been used to quantify the fitness of a template protein structure to the unknown protein structure. This work has other important direct applications such as verification of the novelty of an unknown protein and validation of the accuracy of an available protein structure model in drug design. More importantly, the presented work has the potential to bridge the gap between experimental and computational methods of structure determination.

Rapid classification of protein structure models using unassigned backbone RDCs and probability density profile analysis (PDPA)

A method of identifying the best structural model for a protein of unknown structure from a list of structural candidates using unassigned 15N1H residual dipolar coupling (RDC) data and probability density profile analysis (PDPA) is described. Ten candidate structures have been obtained for the structural genomics target protein PF2048.1 using ROBETTA. 15N1H residual dipolar couplings have been measured from NMR spectra of the protein in two alignment media and these data have been analyzed using PDPA to rank the models in terms of their ability to represent the actual structure. A number of advantages in using this method to characterize a protein structure become apparent. RDCs can easily and rapidly be acquired, and without the need for assignment, the cost and duration of data acquisition is greatly reduced. The approach is quite robust with respect to imprecise and missing data. In the case of PF2048.1, a 79 residue protein, only 58 and 55 of the total RDC data were observed. The method can accelerate structure determination at higher resolution using traditional NMR spectroscopy by providing a starting point for the addition of NOEs and other NMR structural data.

Protein backbone motions viewed by intraresidue and sequential HN-Halpha residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.

RDC-assisted modeling of symmetric protein homo-oligomers

Protein oligomerization serves an important function in biological processes, yet solving structures of protein oligomers has always been a challenge. For solution NMR, the challenge arises both from the increased size of these systems and, in the case of homo-oligomers, from ambiguities in assignment of intra- as opposed to intersubunit NOEs. In this study, we present a residual dipolar coupling (RDC)-assisted method for constructing models of homo-oligomers with purely rotational symmetry. Utilizing the fact that one of the principal axes of the tensor describing the alignment needed for RDC measurement is always parallel to the oligomer symmetry axis, it is possible to greatly restrict possible models for the oligomer. Here, it is shown that, if the monomer structure is known, all allowed dimer models can be constructed using a grid search algorithm and evaluated based on RDC simulations and the quality of the interface between the subunits. Using the Bacillus subtilis protein YkuJ as an example, it is shown that the evaluation criteria based on just two sets of NH RDCs are very selective and can unambiguously produce a model in good agreement with an existing X-ray structure of YkuJ.

Structure of tumor suppressor p53 and its intrinsically disordered N-terminal transactivation domain

Proteins with intrinsically disordered domains are implicated in a vast range of biological processes, especially in cell signaling and regulation. Having solved the quaternary structure of the folded domains in the tumor suppressor p53 by a multidisciplinary approach, we have now determined the average ensemble structure of the intrinsically disordered N-terminal transactivation domain (TAD) by using residual dipolar couplings (RDCs) from NMR spectroscopy and small-angle x-ray scattering (SAXS). Remarkably, not only were we able to measure RDCs of the isolated TAD, but we were also able to do so for the TAD in both the full-length tetrameric p53 protein and in its complex with a specific DNA response element. We determined the orientation of the TAD ensemble relative to the core domain, found that the TAD was stiffer in the proline-rich region (residues 64-92), which has a tendency to adopt a polyproline II (PPII) structure, and projected the TAD away from the core. We located the TAD in SAXS experiments on a complex between tetrameric p53 and four Taz2 domains that bind tightly to the TAD (residues 1-57) and acted as "reporters." The p53-Taz2 complex was an extended cross-shaped structure. The quality of the SAXS data enabled us to model the disordered termini and the folded domains in the complex with DNA. The core domains enveloped the response element in the center of the molecule, with the Taz2-bound TADs projecting outward from the core.

Conformational analysis of a dermatan sulfate-derived tetrasaccharide by NMR, molecular modeling, and residual dipolar couplings

The solution conformation behavior of a dermatan-derived tetrasaccharide--Delta HexA-(1-->3)-GalNAc4S-beta-(1-->4)-IdoA-alpha-(1-->3)-red-GalNAc4S (S is a sulfate group)--has been explored by means of NMR spectroscopy, especially by NOE-based conformational analysis. The tetrasaccharide was present as four species, two of which are chemically different in the anomeric orientation of the reducing 2-deoxy-2-acetamido-galactose (red-GalNAc) residue, while the other two are the result of different conformations of the iduronic acid (IdoA) unit. The two alpha-beta-interconverting anomers were present in a 0.6:1 ratio. Ring conformations have been defined by analysis of (3)J(H,H) coupling constants and interresidual NOE contacts. Both 2-deoxy-2-acetamido-galactose (GalNAc) residues were found in the (4)C(1) chair conformation, the unsaturated uronic acid (Delta-Hex A) adopts a strongly predominant half-chair (1)H(2) conformation, while the IdoA residue exists either in the (1)C(4) chair or in the (2)S(0) skewed boat geometries, in a 4:1 ratio. There is a moderate flexibility of Phi and Psi torsions as suggested by nuclear Overhauser effects (NOEs), molecular modeling (MM), and molecular dynamics (MD) studies. This was further investigated by residual dipolar couplings (RDCs). One-bond C--H RDCs ((1)D(C,H)) and long-range H-H ((3)D(H,H)) RDCs were measured for the tetrasaccharide in a phage solution and interpreted in combination with restrained MD simulation. The RDC-derived data substantially confirmed the validity of the conformer distribution resulting from the NOE-derived simulations, but allowed an improved definition of the conformational behavior of the oligosaccharides in solution. In summary, the data show a moderate flexibility of the four tetrasaccharide species at the central glycosidic linkage. Differences in the shapes of species with the IdoA in skew and in chair conformations and in the distribution of the sulfate groups have also been highlighted.

Structure of a self-splicing group II intron catalytic effector domain 5: parallels with spliceosomal U6 RNA

Domain 5 (D5) is absolutely required for all catalytic functions of group II introns. Here we describe the solution NMR structure, electrostatic calculations, and detailed magnesium ion-binding surface of D5 RNA from the Pylaiella littoralis large ribosomal RNA intron (D5-PL). The overall structure consists of a hairpin capped by a GNRA tetraloop. The stem is divided into lower and upper helices of 8 and 5 bp, respectively, separated by an internal bulge. The D5-PL internal bulge nucleotides stack into the helical junction, resulting in a coupling between the bulge A25 and the closing base pair (G8-C27) of the lower helix. Comparison of the D5-PL structure to previously reported related structures indicates that our structure is most similar, in the helical regions, to the crystal structure of D5 from yeast Ai5gamma (D5-Ai5gamma) and the NMR structure of the U6 snRNA stem-loop region. Our structure differs in many respects from both the NMR and X-ray structures of D5-Ai5gamma in the bulge region. Electrostatic calculations and NMR chemical shift perturbation analyses reveal magnesium ion-binding sites in the tetraloop, internal bulge, and the AGC triad in the lower stem. Our results suggest that the structure, electrostatic environment, and the magnesium ion-binding sites within the tetraloop, bulge, and triad regions are conserved features of the splicing machinery of both the group II introns and the spliceosome that are likely key for catalytic function.

Solution structure of (gamma)S-crystallin by molecular fragment replacement NMR

The solution structure of murine gammaS-crystallin (gammaS) has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of dipolar couplings, recorded in different alignment media, and supplemented by a small number of NOE distance restraints. gammaS consists of two topologically similar domains, arranged with an approximate twofold symmetry, and each domain shows close structural homology to closely related (approximately 50% sequence identity) domains found in other members of the gamma-crystallin family. Each domain consists of two four-strand "Greek key" beta-sheets. Although the domains are tightly anchored to one another by the hydrophobic surfaces of the two inner Greek key motifs, the N-arm, the interdomain linker and several turn regions show unexpected flexibility and disorder in solution. This may contribute entropic stabilization to the protein in solution, but may also indicate nucleation sites for unfolding or other structural transitions. The method used for solving the gammaS structure relies on the recently introduced molecular fragment replacement method, which capitalizes on the large database of protein structures previously solved by X-ray crystallography and NMR.