Nuclear Magnetic Resonance Restraints Research Articles

Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions.

In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the β and the α+β structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.

Aromatic Interactions Drive the Coupled Folding and Binding of the Intrinsically Disordered Sesbania mosaic Virus VPg Protein.

The plant Sesbania mosaic virus [a (+)-ssRNA sobemovirus] VPg protein is intrinsically disordered in solution. For the virus life cycle, the VPg protein is essential for replication and for polyprotein processing that is carried out by a virus-encoded protease. The nuclear magnetic resonance (NMR)-derived tertiary structure of the protease-bound VPg shows it to have a novel tertiary structure with an α-β-β-β topology. The quaternary structure of the high-affinity protease-VPg complex (≈27 kDa) has been determined using HADDOCK protocols with NMR (residual dipolar coupling, dihedral angle, and nuclear Overhauser enhancement) restraints and mutagenesis data as inputs. The geometry of the complex is in excellent agreement with long-range orientational restraints such as residual dipolar couplings and ring-current shifts. A "vein" of aromatic residues on the protease surface is pivotal for the folding of VPg via intermolecular edge-to-face π···π stacking between Trp271 and Trp368 of the protease and VPg, respectively, and for the CH···π interactions between Leu361 of VPg and Trp271 of the protease. The structure of the protease-VPg complex provides a molecular framework for predicting sites of important posttranslational modifications such as RNA linkage and phosphorylation and a better understanding of the coupled folding upon binding of intrinsically disordered proteins. The structural data presented here augment the limited structural data available on viral proteins, given their propensity for structural disorder.

Conformational ensembles of an RNA hairpin using molecular dynamics and sparse NMR data.

Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.

Integrative Protein Modeling in RosettaNMR from Sparse Paramagnetic Restraints

Computational methods to predict protein structure from nuclear magnetic resonance (NMR) restraints that only require assignment of backbone signals, hold great potential to study larger proteins. Ideally, computational methods designed to work with sparse data need to add atomic detail that is missing in the experimental restraints. We introduce a comprehensive framework into the Rosetta suite that uses NMR restraints derived from paramagnetic labeling. Specifically, RosettaNMR incorporates pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. It continues to use backbone chemical shifts and nuclear Overhauser effect distance restraints. We assess RosettaNMR for protein structure prediction by folding 28 monomeric proteins and 8 homo-oligomeric proteins. Furthermore, the general applicability of RosettaNMR is demonstrated on two protein-protein and three protein-ligand docking examples. Paramagnetic restraints generated more accurate models for 85% of the benchmark proteins and, when combined with chemical shifts, sampled high-accuracy models (≤2Å) in 50% of the cases.

Integrative Approaches in Structural Biology: A More Complete Picture from the Combination of Individual Techniques.

With the recent technological and computational advancements, structural biology has begun to tackle more and more difficult questions, including complex biochemical pathways and transient interactions among macromolecules. This has demonstrated that, to approach the complexity of biology, one single technique is largely insufficient and unable to yield thorough answers, whereas integrated approaches have been more and more adopted with successful results. Traditional structural techniques (X-ray crystallography and Nuclear Magnetic Resonance (NMR)) and the emerging ones (cryo-electron microscopy (cryo-EM), Small Angle X-ray Scattering (SAXS)), together with molecular modeling, have pros and cons which very nicely complement one another. In this review, three examples of synergistic approaches chosen from our previous research will be revisited. The first shows how the joint use of both solution and solid-state NMR (SSNMR), X-ray crystallography, and cryo-EM is crucial to elucidate the structure of polyethylene glycol (PEG)ylated asparaginase, which would not be obtainable through any of the techniques taken alone. The second deals with the integrated use of NMR, X-ray crystallography, and SAXS in order to elucidate the catalytic mechanism of an enzyme that is based on the flexibility of the enzyme itself. The third one shows how it is possible to put together experimental data from X-ray crystallography and NMR restraints in order to refine a protein model in order to obtain a structure which simultaneously satisfies both experimental datasets and is therefore closer to the ‘real structure’.

Nuclear magnetic resonance structure-based drug design.

Nuclear magnetic resonance structure-based drug design.

Characterisation of the membrane topology and molecular structure of LITAF to provide insights into the molecular pathogenesis of Charcot-Marie-Tooth disease type 1C

BackgroundThe Charcot-Marie-Tooth (CMT) diseases are the commonest inherited neuromuscular disorders. These disorders are characterised by progressive degenerative peripheral neuropathies and are associated with mutations in over 80 genes. Mutations in LITAF cause CMT type 1C (CMT1C), an autosomal dominant demyelinating subtype of the disease. Why mutations in LITAF cause disease remains unclear, so we aimed to gather data about its molecular structure, which are lacking despite being highly conserved throughout evolution. MethodsWe applied biochemical and cell biological techniques (differential membrane preparation, semi-permeabilisation and immunofluorescence microscopy, and cell-free protein expression) to determine the membrane topology of LITAF. Micro particle-induced X-ray emission (microPIXE) was used to establish whether LITAF associates with metal ions. We then analysed recombinantly expressed LITAF contructs by solution nuclear magnetic resonance (NMR) spectroscopy to gain structural insights to the LITAF protein, thereby allowing us to use the Rosetta software to generate a structural molecular model. FindingsWe found that LITAF behaves as an integral membrane protein and that membrane association is dependent on six leucine residues found within the C-terminal LITAF domain. The membrane topology of the LITAF domain was such that both N-terminals and C-terminals were found on the cytosolic surface of endosomes and that these regions coordinated a single zinc atom. NMR spectroscopy was then used to characterise the secondary structural motifs of the LITAF domain, resulting in a correlation spectrum which was over 90% assigned. These data enabled us to create a structural model using NMR restraints and Rosetta. A LITAF construct harbouring a CMT1C-associated pathogenic mutation was then analysed by NMR, revealing marked structural shifts across the LITAF domain. InterpretationWe present the first structural model, to our knowledge, of the highly-conserved LITAF domain, which is a membrane-anchoring motif forming a zinc-binding structure on the surface of endosomes. We also show that a CMT1C-associated mutation causes large structural shifts across the LITAF domain, suggesting that structural instability on the surface of Scwann cell endosomes underpins the neuropathy seen in CMT1C. These insights will provide a platform for the development of future treatments. FundingWellcome Trust.

Using Pseudocontact Shifts and Residual Dipolar Couplings as Exact NMR Restraints for the Determination of Protein Structural Ensembles.

Nuclear magnetic resonance (NMR) spectroscopy provides detailed information about the structure and dynamics of proteins by exploiting the conformational dependence of the magnetic properties of certain atomic nuclei. The mapping between NMR measurements and molecular structures, however, often requires approximated descriptions based on the fitting of a number of parameters, thus reducing the quality of the information available from the experiments. To improve on this limitation, we show here that it is possible to use pseudocontact shifts and residual dipolar couplings as "exact" NMR restraints. We implement this strategy by using a replica-averaging method and illustrate its application by calculating an ensemble of structures representing the dynamics of the two-domain protein calmodulin.

Conjoined use of EM and NMR in RNA structure refinement.

More than 40% of the RNA structures have been determined using nuclear magnetic resonance (NMR) technique. NMR mainly provides local structural information of protons and works most effectively on relatively small biomacromolecules. Hence structural characterization of large RNAs can be difficult for NMR alone. Electron microscopy (EM) provides global shape information of macromolecules at nanometer resolution, which should be complementary to NMR for RNA structure determination. Here we developed a new energy term in Xplor-NIH against the density map obtained by EM. We conjointly used NMR and map restraints for the structure refinement of three RNA systems — U2/U6 small-nuclear RNA, genome-packing motif (ΨCD)2 from Moloney murine leukemia virus, and ribosome-binding element from turnip crinkle virus. In all three systems, we showed that the incorporation of a map restraint, either experimental or generated from known PDB structure, greatly improves structural precision and accuracy. Importantly, our method does not rely on an initial model assembled from RNA duplexes, and allows full torsional freedom for each nucleotide in the torsion angle simulated annealing refinement. As increasing number of macromolecules can be characterized by both NMR and EM, the marriage between the two techniques would enable better characterization of RNA three-dimensional structures.

Small-molecule ligand docking into comparative models with Rosetta

Structure-based drug design is frequently used to accelerate the development of small-molecule therapeutics. Although substantial progress has been made in X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy, the availability of high-resolution structures is limited owing to the frequent inability to crystallize or obtain sufficient NMR restraints for large or flexible proteins. Computational methods can be used to both predict unknown protein structures and model ligand interactions when experimental data are unavailable. This paper describes a comprehensive and detailed protocol using the Rosetta modeling suite to dock small-molecule ligands into comparative models. In the protocol presented here, we review the comparative modeling process, including sequence alignment, threading and loop building. Next, we cover docking a small-molecule ligand into the protein comparative model. In addition, we discuss criteria that can improve ligand docking into comparative models. Finally, and importantly, we present a strategy for assessing model quality. The entire protocol is presented on a single example selected solely for didactic purposes. The results are therefore not representative and do not replace benchmarks published elsewhere. We also provide an additional tutorial so that the user can gain hands-on experience in using Rosetta. The protocol should take 5-7 h, with additional time allocated for computer generation of models.

Three-Dimensional Protein Fold Determination from Backbone Amide Pseudocontact Shifts Generated by Lanthanide Tags at Multiple Sites

Site-specific attachment of paramagnetic lanthanide ions to a protein generates pseudocontact shifts (PCS) in the nuclear magnetic resonance (NMR) spectra of the protein that are easily measured as changes in chemical shifts. By labeling the protein with lanthanide tags at four different sites, PCSs are observed for most amide protons and accurate information is obtained about their coordinates in three-dimensional space. The approach is demonstrated with the chaperone ERp29, for which large differences have been reported between X-ray andNMR structures of the C-terminal domain, ERp29-C. The results unambiguously show that the structure of rat ERp29-C in solution is similar to the crystal structure of human ERp29-C. PCSs of backbone amides were the only structural restraints required. Because these can be measured for more dilute protein solutions than other NMR restraints, the approach greatly widens the range of proteins amenable to structural studies in solution.

Structural features of the apelin receptor N-terminal tail and first transmembrane segment implicated in ligand binding and receptor trafficking

G-protein coupled receptors (GPCRs) comprise a large family of membrane proteins with rich functional diversity. Signaling through the apelin receptor (AR or APJ) influences the cardiovascular system, central nervous system and glucose regulation. Pathophysiological involvement of apelin has been shown in atherosclerosis, cancer, human immunodeficiency virus-1 (HIV-1) infection and obesity. Here, we present the high-resolution nuclear magnetic resonance (NMR) spectroscopy-based structure of the N-terminus and first transmembrane (TM) segment of AR (residues 1–55, AR55) in dodecylphosphocholine micelles. AR55 consists of two disrupted helices, spanning residues D14-K25 and A29-R551.59. Molecular dynamics (MD) simulations of AR built from a hybrid of experimental NMR and homology model-based restraints allowed validation of the AR55 structure in the context of the full-length receptor in a hydrated bilayer. AR55 structural features were functionally probed using mutagenesis in full-length AR through monitoring of apelin-induced extracellular signal-regulated kinase (ERK) phosphorylation in transiently transfected human embryonic kidney (HEK) 293A cells. Residues E20 and D23 form an extracellular anionic face and interact with lipid headgroups during MD simulations in the absence of ligand, producing an ideal binding site for a cationic apelin ligand proximal to the membrane–water interface, lending credence to membrane-catalyzed apelin-AR binding. In the TM region of AR55, N461.50 is central to a disruption in helical character. G421.46, G451.49 and N461.50, which are all involved in the TM helical disruption, are essential for proper trafficking of AR. In summary, we introduce a new correlative NMR spectroscopy and computational biochemistry methodology and demonstrate its utility in providing some of the first high-resolution structural information for a peptide-activated GPCR TM domain.

Requirements on Paramagnetic Relaxation Enhancement Data for Membrane Protein Structure Determination by NMR

Nuclear magnetic resonance (NMR) structure calculations of the α-helical integral membrane proteins DsbB, GlpG, and halorhodopsin show that distance restraints from paramagnetic relaxation enhancement (PRE) can provide sufficient structural information to determine their structure with an accuracy of about 1.5Å in the absence of other long-range conformational restraints. Our systematic study with simulated NMR data shows that about one spin label per transmembrane helix is necessary for obtaining enough PRE distance restraints to exclude wrong topologies, such as pseudo mirror images, if only limited other NMR restraints are available. Consequently, an experimentally realistic amount of PRE data enables α-helical membrane protein structure determinations that would not be feasible with the very limited amount of conventional NOESY data normally available for these systems. These findings are in line with our recent first de novo NMR structure determination of a heptahelical integral membrane protein, proteorhodopsin, that relied extensively on PRE data.

Structural Basis of the RNase H1 Activity on Stereo Regular Borano Phosphonate DNA/RNA Hybrids

Numerous DNA chemistries for improving oligodeoxynucleotide (ODN)-based RNA targeting have been explored. The majority of the modifications render the ODN/RNA target insensitive to RNase H1. Borano phosphonate ODN's are among the few modifications that are tolerated by RNase H1. To understand the effect of the stereochemistry of the BH(3) modification on the nucleic acid structure and RNase H1 enzyme activity, we have investigated two DNA/RNA hybrids containing either a R(P) or S(P) BH(3) modification by nuclear magnetic resonance (NMR) spectroscopy. T(M) studies show that the stabilities of R(P) and S(P) modified DNA/RNA hybrids are essentially identical (313.8 K) and similar to that of an unmodified control (312.9 K). The similarity is also reflected in the imino proton spectra. To characterize such similar structures, we used a large number of NMR restraints (including dipolar couplings and backbone torsion angles) to determine structural features that were important for RNase H1 activity. The final NMR structures exhibit excellent agreement with the data (total R(x) values of <6%) with helical properties between those of an A and B helix. Subtle backbone variations are observed in the DNA near the modification, while the RNA strands are relatively unperturbed. In the case of the S(P) modification, for which more perturbations are recorded, a slightly narrower minor groove is also obtained. Unique NOE base contacts localize the S(P) BH(3) group in the major groove while the R(P) BH(3) group points away from the DNA. However, this creates a potential clash of the R(P) BH(3) groups with important RNase H1 residues in a complex, while the S(P) BH(3) groups could be tolerated. We therefore predict that on the basis of our NMR structures a fully R(P) BH(3) DNA/RNA hybrid would not be a substrate for RNase H1.

Combining NMR and EPR methods for homodimer protein structure determination.

There is a general need to develop more powerful and more robust methods for structural characterization of homodimers, homo-oligomers, and multiprotein complexes using solution-state NMR methods. In recent years, there has been increasing emphasis on integrating distinct and complementary methodologies for structure determination of multiprotein complexes. One approach not yet widely used is to obtain intermediate and long-range distance constraints from paramagnetic relaxation enhancements (PRE) and electron paramagnetic resonance (EPR)-based techniques such as double electron electron resonance (DEER), which, when used together, can provide supplemental distance constraints spanning to 10-70 A. In this Communication, we describe integration of PRE and DEER data with conventional solution-state nuclear magnetic resonance (NMR) methods for structure determination of Dsy0195, a homodimer (62 amino acids per monomer) from Desulfitobacterium hafniense. Our results indicate that combination of conventional NMR restraints with only one or a few DEER distance constraints and a small number of PRE constraints is sufficient for the automatic NMR-based structure determination program CYANA to build a network of interchain nuclear Overhauser effect constraints that can be used to accurately define both the homodimer interface and the global homodimer structure. The use of DEER distances as a source of supplemental constraints as described here has virtually no upper molecular weight limit, and utilization of the PRE constraints is limited only by the ability to make accurate assignments of the protein amide proton and nitrogen chemical shifts.

Non-Native Interactions Are Critical for Mechanical Strength in PKD Domains

SummaryExperimental observation has led to the commonly held view that native state protein topology is the principle determinant of mechanical strength. However, the PKD domains of polycystin-1 challenge this assumption: they are stronger than predicted from their native structure. Molecular dynamics simulations suggest that force induces rearrangement to an intermediate structure, with nonnative hydrogen bonds, that resists unfolding. Here we test this hypothesis directly by introducing mutations designed to prevent formation of these nonnative interactions. We find that these mutations, which only moderately destabilize the native state, reduce the mechanical stability dramatically. The results demonstrate that nonnative interactions impart significant mechanical stability, necessary for the mechanosensor function of polycystin-1. Remarkably, such nonnative interactions result from force-induced conformational change: the PKD domain is strengthened by the application of force.

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2′-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study.

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.

Structural insights into the GTPase domain of Escherichia coli MnmE protein

The Escherichia coli MnmE protein is a 50-kDa multidomain GTPase involved in tRNA modification. Its homologues in eukaryotes are crucial for mitochondrial respiration and, thus, it is thought that the human protein might be involved in mitochondrial diseases. Unlike Ras, MnmE shows a high intrinsic GTPase activity and requires effective GTP hydrolysis, and not simply GTP binding, to be functionally active. The isolated MnmE G-domain (165 residues) conserves the GTPase activity of the entire protein, suggesting that it contains the catalytic residues for GTP hydrolysis. To explore the GTP hydrolysis mechanism of MnmE, we analyzed the effect of low pH on binding and hydrolysis of GTP, as well as on the formation of a MnmE transition state mimic. GTP hydrolysis by MnmE, but not GTP binding or formation of a complex with mant-GDP and aluminium fluoride, is impaired at acidic pH, suggesting that the chemistry of the transition state mimic is different to that of the true transition state, and that some residue(s), critical for GTP hydrolysis, is severely affected by low pH. We use a nuclear magnetic resonance (NMR)-based approach to get insights into the MnmE structure and properties. The combined use of NMR restraints and homology structural information allowed the determination of the MnmE G-domain structure in its free form. Chemical shift structure-based prediction provided a good basis for structure refinement and validation. Our data support that MnmE, unlike other GTPases, does not use an arginine finger to drive catalysis, although Arg252 may play a role in stabilization of the transition state.

Amplitudes and directions of internal protein motions from a JAM analysis of15N relaxation data

A method has been developed for characterizing dynamic structures of proteins in solution by using nuclear magnetic resonance (NMR) restraints and 15N relaxation data. This method is based on the concept of the jumping-among-minima (JAM) model. In this model we assume that protein dynamics can be described on the basis of conformational substates, and involves intra- and inter-substate motion. A set of substates is created by picking energy-minimized conformations from the conformational space consistent with the geometric NMR restraints. Intra-substate motions, which occur on the timescale of approximately 10 ps, are simulated with molecular dynamics (MD) calculations with force-field energy terms. Statistical weights of the conformational substates are determined to reproduce the NMR relaxation parameters. The refinement procedure consists of four stages: (i) determination of the ensemble of structures that satisfy NMR restraints, (ii) determination of intra-substate fluctuation, (iii) determination of statistical weights of conformational substates to reproduce model-free relaxation parameters, and (iv) analysis of the resulting dynamic structure to determine amplitudes and directions of internal protein motions. This method was employed to investigate structure and dynamics of the adhesion domain of human CD2 (hCD2) in solution. Two major collective modes, whose contributions to atomic mean-square fluctuations are 77.1% in total, are identified by the refinement. The first mode is interpreted as a rigid-body motion of a protein segment consisting of a part of the B--C loop, a part of the F strand, and the F--G loop. Another type of smaller-amplitude mode is indicated for the C'--C'' loop. The motions affect primarily the curvature of the slightly concave counterreceptor-binding site and represent transitions between a concave (closed) and flat (open) binding face. By comparing the ensemble of structures in solution to the complex structure with counterreceptor CD58, we found that these two types of motions resemble the change upon counterreceptor binding.