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

Lessons on protein structure from interleukin-4: All disulfides are not created equal

Interleukin-4 (IL-4) is a hematopoietic cytokine composed by a four-helix bundle stabilized by an antiparallel beta-sheet and three disulfide bonds: Cys3-Cys127, Cys24-Cys65, and Cys46-Cys99. IL-4 is involved in several immune responses associated to infection, allergy, autoimmunity, and cancer. Besides its physiological relevance, IL-4 is often used as a "model" for protein design and engineering. Hence, to understand the role of each disulfide in the structure and dynamics of IL-4, we carried out several spectroscopic analyses (circular dichroism [CD], fluorescence, nuclear magnetic resonance [NMR]), and molecular dynamics (MD) simulations on wild-type IL-4 and four IL-4 disulfide mutants. All disulfide mutants showed loss of structure, altered interhelical angles, and looser core packings, showing that all disulfides are relevant for maintaining the overall fold and stability of the four-helix bundle motif, even at very low pH. In the absence of the disulfide connecting both protein termini Cys3-Cys127, C3T-IL4 showed a less packed protein core, loss of secondary structure (~9%) and fast motions on the sub-nanosecond time scale (lower S2 order parameters and larger τc correlation time), especially at the two protein termini, loops, beginning of helix A and end of helix D. In the absence of Cys24-Cys65, C24T-IL4 presented shorter alpha-helices (14% loss in helical content), altered interhelical angles, less propensity to form the small anti-parallel beta-sheet and increased dynamics. Simultaneously deprived of two disulfides (Cys3-Cys127 and Cys24-Cys65), IL-4 formed a partially folded "molten globule" with high 8-anilino-1-naphtalenesulphonic acid-binding affinity and considerable loss of secondary structure (~50%decrease), as shown by the far UV-CD, NMR, and MD data.

Synergism between x-ray crystallography and NMR residual dipolar couplings in characterizing protein dynamics

The important role of structural dynamics in protein function is widely recognized. Thermal or B-factors and their anisotropy, seen in x-ray analysis of protein structures, report on the presence of atomic coordinate heterogeneity that can be attributed to motion. However, their quantitative evaluation in terms of protein dynamics by x-ray ensemble refinement remains challenging. NMR spectroscopy provides quantitative information on the amplitudes and time scales of motional processes. Unfortunately, with a few exceptions, the NMR data do not provide direct insights into the atomic details of dynamic trajectories. Residual dipolar couplings, measured by solution NMR, are very precise parameters reporting on the time-averaged bond-vector orientations and may offer the opportunity to derive correctly weighted dynamic ensembles of structures for cases where multiple high-resolution x-ray structures are available. Applications to the SARS-CoV-2 main protease, M^pro, and ubiquitin highlight this complementarity of NMR and crystallography for quantitative assessment of internal motions.

Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, Mpro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic Rfree, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55-2.19 Å) Mpro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member "super ensemble" averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB Mpro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.

Plant Defense Elicitation by the Hydrophobin Cerato-Ulmin and Correlation with Its Structural Features

Cerato-ulmin (CU) is a 75-amino-acid-long protein that belongs to the hydrophobin family. It self-assembles at hydrophobic-hydrophilic interfaces, forming films that reverse the wettability properties of the bound surface: a capability that may confer selective advantages to the fungus in colonizing and infecting elm trees. Here, we show for the first time that CU can elicit a defense reaction (induction of phytoalexin synthesis and ROS production) in non-host plants (Arabidopsis) and exerts its eliciting capacity more efficiently when in its soluble monomeric form. We identified two hydrophobic clusters on the protein's loops endowed with dynamical and physical properties compatible with the possibility of reversibly interconverting between a disordered conformation and a β-strand-rich conformation when interacting with hydrophilic or hydrophobic surfaces. We propose that the plasticity of those loops may be part of the molecular mechanism that governs the protein defense elicitation capability.

Polymer-Nanodiscs as a Novel Alignment Medium for High-Resolution NMR-Based Structural Studies of Nucleic Acids

Residual dipolar couplings (RDCs) are increasingly used for high-throughput NMR-based structural studies and to provide long-range angular constraints to validate and refine structures of various molecules determined by X-ray crystallography and NMR spectroscopy. RDCs of a given molecule can be measured in an anisotropic environment that aligns in an external magnetic field. Here, we demonstrate the first application of polymer-based nanodiscs for the measurement of RDCs from nucleic acids. Polymer-based nanodiscs prepared using negatively charged SMA-EA polymer and zwitterionic DMPC lipids were characterized by size-exclusion chromatography, 1H NMR, dynamic light-scattering, and 2H NMR. The magnetically aligned polymer-nanodiscs were used as an alignment medium to measure RDCs from a 13C/15N-labeled fluoride riboswitch aptamer using 2D ARTSY-HSQC NMR experiments. The results showed that the alignment of nanodiscs is stable for nucleic acids and nanodisc-induced RDCs fit well with the previously determined solution structure of the riboswitch. These results demonstrate that SMA-EA-based lipid-nanodiscs can be used as a stable alignment medium for high-resolution structural and dynamical studies of nucleic acids, and they can also be applicable to study various other biomolecules and small molecules in general.

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.

Extending the applicability of P3D for structure determination of small molecules

The application of anisotropic nuclear magnetic resonance (NMR) parameters for the correct structural assignment of small molecules requires the use of partially ordered media. Previously we demonstrated that the use of P3D simulations using poly(γ -benzyl-L-glutamate) (PBLG) as an alignment medium allows for the determination of the correct diastereomer from extremely sparse NMR data. Through the analysis of the structural characteristics of small molecules in different alignment media, here we show that when steric or electrostatic factors dominate the alignment, P3D-PBLG retains its diastereomer discrimination power. We also demonstrate that P3D simulations can define the different conformations of a flexible small molecule from sparse NMR data.

Generating Five Independent Molecular Alignments for Simultaneous Protein Structure and Dynamics Determination Using Nuclear Magnetic Resonance Spectroscopy

Residual dipolar couplings (RDCs) are commonly used in NMR for protein structure and dynamics studies, but it is challenging to generate five independent RDC data sets (required for simultaneous structure and dynamics determination) for most protein molecules in the magnetic field. In this work, a reporter protein with a lanthanide tag is introduced to create five independent alignments. This reporter protein is then attached to target proteins where five independent sets of RDCs are also obtained for the target proteins. The fitting of RDCs provides important information about the structure and dynamics of the target proteins. The method is simple and effective and, in principle, can be used to generate complete sets of RDCs for different protein molecules.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

Leiomodin creates a leaky cap at the pointed end of actin-thin filaments

Improper lengths of actin-thin filaments are associated with altered contractile activity and lethal myopathies. Leiomodin, a member of the tropomodulin family of proteins, is critical in thin filament assembly and maintenance; however, its role is under dispute. Using nuclear magnetic resonance data and molecular dynamics simulations, we generated the first atomic structural model of the binding interface between the tropomyosin-binding site of cardiac leiomodin and the N-terminus of striated muscle tropomyosin. Our structural data indicate that the leiomodin/tropomyosin complex only forms at the pointed end of thin filaments, where the tropomyosin N-terminus is not blocked by an adjacent tropomyosin protomer. This discovery provides evidence supporting the debated mechanism where leiomodin and tropomodulin regulate thin filament lengths by competing for thin filament binding. Data from experiments performed in cardiomyocytes provide additional support for the competition model; specifically, expression of a leiomodin mutant that is unable to interact with tropomyosin fails to displace tropomodulin at thin filament pointed ends and fails to elongate thin filaments. Together with previous structural and biochemical data, we now propose a molecular mechanism of actin polymerization at the pointed end in the presence of bound leiomodin. In the proposed model, the N-terminal actin-binding site of leiomodin can act as a "swinging gate" allowing limited actin polymerization, thus making leiomodin a leaky pointed-end cap. Results presented in this work answer long-standing questions about the role of leiomodin in thin filament length regulation and maintenance.

Anisotropy and NMR spectroscopy

Abstract

In this paper, different aspects concerning anisotropy in Nuclear Magnetic Resonance (NMR) spectroscopy have been reviewed. In particular, the relevant theory has been presented, showing how anisotropy stems from the dependence of internal nuclear spin interactions on the molecular orientation with respect to the external magnetic field direction. The consequences of anisotropy in the use of NMR spectroscopy have been critically discussed: on one side, the availability of very detailed structural and dynamic information, and on the other side, the loss of spectral resolution. The experiments used to measure the anisotropic properties in solid and soft materials, where, in contrast to liquids, such properties are not averaged out by the molecular tumbling, have been described. Such experiments can be based either on static low-resolution techniques or on one- and two-dimensional pulse sequences exploiting Magic Angle Spinning (MAS). Examples of applications of NMR spectroscopy have been shown, which exploit anisotropy to obtain important physico-chemical information on several categories of systems, including pharmaceuticals, inorganic materials, polymers, liquid crystals, and self-assembling amphiphiles in water. Solid-state NMR spectroscopy can be considered, nowadays, one of the most powerful characterization techniques for all kinds of solid, either amorphous or crystalline, and semi-solid systems for the obtainment of both structural and dynamic properties on a molecular and supra-molecular scale.

Graphic abstract

Solution NMR spectroscopy of membrane proteins

Integral membrane proteins (IMPs) perform unique and indispensable functions in the cell, making them attractive targets for fundamental research and drug discovery. Developments in protein production, isotope labeling, sample preparation, and pulse sequences have extended the utility of solution NMR spectroscopy for studying IMPs with multiple transmembrane segments. Here we review some recent applications of solution NMR for studying structure, dynamics, and interactions of polytopic IMPs, emphasizing strategies used to overcome common technical challenges.

Determination of Complex Small-Molecule Structures Using Molecular Alignment Simulation

Correct structural assignment of small molecules and natural products is critical for drug discovery and organic chemistry. Anisotropy-based NMR spectroscopy is a powerful tool for the structural assignment of organic molecules, but it relies on the utilization of a medium that disrupts the isotropic motion of molecules in organic solvents. Here, we establish a quantitative correlation between the atomic structure of the alignment medium, the molecular structure of the small molecule, and molecule-specific anisotropic NMR parameters. The quantitative correlation uses an accurate three-dimensional molecular alignment model that predicts residual dipolar couplings of small molecules aligned by poly(γ-benzyl-l-glutamate). The technique facilitates reliable determination of the correct stereoisomer and enables unequivocal, rapid determination of complex molecular structures from extremely sparse NMR data.

Amide additives improve RDC measurements in polyacrylamide

Residual dipolar couplings (RDCs) provide valuable NMR parameters that can be used for structural calculation and verification. Measuring RDCs requires aligning macromolecules using one of various types of alignment media. Of different alignment media options, stretched or compressed polyacrylamide gels are advantageous due to their chemical stability. However, polyacrylamide interacts with proteins and significantly broadens NMR resonances. In this study, we found that the amide-containing compounds asparagine, glutamine and propionamide improve spectral quality of proteins in polyacrylamide gel without significantly reducing the magnitude of RDC values. Moreover, we showed that propionamide is an attractive additive that increases protein solubility without interfering with protein stability, ligand binding or NMR pulse width, suggesting its potential applications for our NMR methods.

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.

Paramagpy: software for fitting magnetic susceptibility tensors using paramagnetic effects measured in NMR spectra

Paramagnetic metal ions with fast-relaxing electrons generate pseudocontact shifts (PCSs), residual dipolar couplings (RDCs), paramagnetic relaxation enhancements (PREs) and cross-correlated relaxation (CCR) in the nuclear magnetic resonance (NMR) spectra of the molecules they bind to. These effects offer long-range structural information in molecules equipped with binding sites for such metal ions. Here we present the new open-source software Paramagpy, which has been written in Python 3 with a graphic user interface. Paramagpy combines the functionalities of different currently available programs to support the fitting of magnetic susceptibility tensors using PCS, RDC, PRE and CCR data and molecular coordinates in Protein Data Bank (PDB) format, including a convenient graphical user interface. Paramagpy uses efficient fitting algorithms to avoid local minima and supports corrections to back-calculated PCS and PRE data arising from cross-correlation effects with chemical shift tensors. The source code is available from 10.5281/zenodo.3594568 .

Residual Dipolar Couplings for Resolving Cysteine Bridges in Disulfide-Rich Peptides

Disulfide bridges in proteins are formed by the oxidation of pairs of cysteine residues. These cross-links play a critical role in stabilizing the 3D-structure of small disulfide rich polypeptides such as hormones and venom toxins. The arrangement of the multiple disulfide bonds directs the peptide fold into distinct structural motifs that have evolved for resistance against biochemical and physical insults. These structural scaffolds have, therefore, proven to be very attractive in bioengineering efforts to develop novel biologics with applications in health and agriculture. Structural characterization of small disulfide rich peptides (DRPs) presents unique challenges when using commonly applied biophysical methods. NMR is the most commonly used method for studying such molecules, where the relatively small size of these molecules results in highly precise structural ensembles defined by a large number of distance and dihedral angle restraints per amino acid. However, in NMR the sulfur atoms that are involved in three of the five dihedral angles in a disulfide bond cannot be readily measured. Given the central role of disulfide bonds in the structure of these molecules, it is unclear what the inherent resolution of such NMR structures is when using traditional NMR methods. Here, we use an extensive set of long-range residual dipolar couplings (RDCs) to assess the resolution of the NMR structure of a disulfide-rich peptide. We find that structures based primarily on NOEs, yield ensembles that are equivalent to a crystallographic resolution of 2-3 Å in resolution, and that incorporation of RDCs reduces this to ~1-1.5 Å resolution. At this resolution the sidechain of ordered amino acids can be defined accurately, allowing the geometry of the cysteine bridges to be better defined, and allowing for disulfide-bond connectivities to be determined with high confidence. The observed improvements in resolution when using RDCs is remarkable considering the small size of these peptides.

On the complementarity of X-ray and NMR data

X-ray crystallography and NMR contain complementary information for the structural characterization of biological macromolecules. X-ray diffraction is primarily sensitive to the overall shape of the molecule, whereas NMR is mostly sensitive to the atomic detail. Their combination can therefore provide a stronger justification for the resulting structure. For their combination we have recently proposed REFMAC-NMR, which relies on primary data from both techniques for joint refinement. This possibility raises the compelling question of how far the complementarity can be extended. In this paper, we describe an integrative approach to the refinement with NMR data of four X-ray structures of hen-egg-white lysozyme, solved at atomic resolution in four different crystal forms, and we demonstrate that the outcome critically depends on the crystal form itself, reflecting the sensitivity of NMR to fine details.

Opportunities and Challenges of Backbone, Sidechain, and RDC Experiments to Study Membrane Protein Dynamics in a Detergent-Free Lipid Environment Using Solution State NMR

Whereas solution state NMR provided a wealth of information on the dynamics landscape of soluble proteins, only few studies have investigated membrane protein dynamics in a detergent-free lipid environment. Recent developments of smaller nanodiscs and other lipid-scaffolding polymers, such as styrene maleic acid (SMA), however, open new and promising avenues to explore the function-dynamics relationship of membrane proteins as well as between membrane proteins and their surrounding lipid environment. Favorably sized lipid-bilayer nanodiscs, established membrane protein reconstitution protocols and sophisticated solution NMR relaxation methods probing dynamics over a wide range of timescales will eventually reveal unprecedented lipid-membrane protein interdependencies that allow us to explain things we have not been able to explain so far. In particular, methyl group dynamics resulting from CEST, CPMG, ZZ exchange, and RDC experiments are expected to provide new and surprising insights due to their proximity to lipids, their applicability in large 100+ kDa assemblies and their simple labeling due to the availability of commercial precursors. This review summarizes the recent developments of membrane protein dynamics with a special focus on membrane protein dynamics in lipid-bilayer nanodiscs. Opportunities and challenges of backbone, side chain and RDC dynamics applied to membrane proteins are discussed. Solution-state NMR and lipid nanodiscs bear great potential to change our molecular understanding of lipid-membrane protein interactions.

Magnetic Alignment of Polymer Macro-Nanodiscs Enables Residual-Dipolar-Coupling-Based High-Resolution Structural Studies by NMR Spectroscopy

Experimentally measured residual dipolar couplings (RDCs) are highly valuable for atomic-resolution structural and dynamic studies of molecular systems ranging from small molecules to large proteins by solution NMR spectroscopy. Here we demonstrate the first use of magnetic-alignment behavior of lyotropic liquid-crystalline polymer macro-nanodiscs (>20 nm in diameter) as a novel alignment medium for the measurement of RDCs using high-resolution NMR. The easy preparation of macro-nanodiscs, their high stability against pH changes and the presence of divalent metal ions, and their high homogeneity make them an efficient tool to investigate a wide range of molecular systems including natural products, proteins, and RNA.

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.

Partial alignment, residual dipolar couplings and molecular symmetry in solution NMR

Residual dipolar couplings (RDCs) and residual anisotropic chemical shifts (RACSs) are produced by the partial alignment of solution NMR samples. RDCs and RACSs yield high-resolution structural and dynamic information on the orientation of bonds and chemical groups in molecules. Many molecules form oligomers or have intrinsic symmetries, which may simplify the analysis of their partial alignment datasets. In this report, we explore the theory of partial alignment using an irreducible spherical representation, and we investigate the impact of molecular symmetry on the alignment of molecules. Though previous studies have reported simplified relationships on the partial alignment of molecules bearing different symmetry groups, we show that these simplified relationships may not be universal and only apply to a limited set of systems.

Accuracy of MD solvent models in RNA structure refinement assessed via liquid-crystal NMR and spin relaxation data

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA.

Solution NMR Spectroscopy for the Determination of Structures of Membrane Proteins in a Lipid Environment

NMR spectroscopy has harnessed the recent technical advances to emerge as a competitive, elegant, and eminently viable technique for determining the solution structures of membrane proteins at the level of atomic resolution. Once a good level of cell-based or cell-free expression and purification of a suitably sized membrane protein has been achieved, then NMR offers a combination of several versatile strategies, for example choice of appropriate deuterated or nondeuterated detergents, temperature, and ionic strength; isotope labeling with ²H, ¹³C, ¹⁵N, with or without protonation of Ile (δ1), Leu, and Val methyl protons; combinatorial labeling or unlabeling of specific amino acids; TROSY based-, nonuniform sampling (NUS) based-, and other NMR experiments; measurement of residual dipolar couplings using stretched polyacrylamide gels or DNA nanotubes; spin labeling and paramagnetic relaxation enhancements (PRE). Strategic combinations of these advancements together with availability of highly sensitive cryogenically cooled-probes equipped high-field NMR spectrometers (up to 1 GHz ¹H frequency) have allowed the perseverant investigator to successfully overcome several of the conventional pitfalls associated with the NMR technique and membrane proteins, viz., low sensitivity, poor sample stability, spectral crowding, and a limited number of NOEs and other constraints for structure calculations. This has resulted in an unprecedented growth in the number of successfully determined NMR structures of large and complex membrane proteins over the last two decades, and this technique now holds great promise for the structure determination of an ever larger body of membrane proteins.

13C NMR-Based Approaches for Solving Challenging Stereochemical Problems

Determining the configuration of proton-deficient molecules is challenging using conventional NMR methods including nuclear Overhauser effect (NOE) and the proton-dependent J-based configuration analysis (JBCA). The problem is exacerbated when only one stereoisomer is available. Alternative methods based on the utilization of ¹³C NMR chemical shifts, ¹³C-¹³C homonuclear couplings measured at natural abundance, and residual chemical shift anisotropy measurements in conjunction with density functional theory calculations are illustrated with a proton-deficient model compound.

Maximizing accuracy of RNA structure in refinement against residual dipolar couplings

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

Extending the Applicability of Exact Nuclear Overhauser Enhancements to Large Proteins and RNA

Distance-dependent nuclear Overhauser enhancements (NOEs) are one of the most popular and important experimental restraints for calculating NMR structures. Despite this, they are mostly employed as semiquantitative upper distance bounds, and this discards the wealth of information that is encoded in the cross-relaxation rate constant. Information that is lost includes exact distances between protons and dynamics that occur on the sub-millisecond timescale. Our recently introduced exact measurement of the NOE (eNOE) requires little additional experimental effort relative to other NMR observables. So far, we have used eNOEs to calculate multistate ensembles of proteins up to approximately 150 residues. Here, we briefly revisit eNOE methodology and present two new directions for the use of eNOEs: applications to large proteins and RNA.

Enhanced Sampling of Interdomain Motion Using Map-Restrained Langevin Dynamics and NMR: Application to Pin1

Many signaling proteins consist of globular domains connected by flexible linkers that allow for substantial domain motion. Because these domains often serve as complementary functional modules, the possibility of functionally important domain motions arises. To explore this possibility, we require knowledge of the ensemble of protein conformations sampled by interdomain motion. Measurements of NMR residual dipolar couplings (RDCs) of backbone HN bonds offer a per-residue characterization of interdomain dynamics, as the couplings are sensitive to domain orientation. A challenge in reaching this potential is the need to interpret the RDCs as averages over dynamic ensembles of domain conformations. Here, we address this challenge by introducing an efficient protocol for generating conformational ensembles appropriate for flexible, multi-domain proteins. The protocol uses map-restrained self-guided Langevin dynamics simulations to promote collective, interdomain motion while restraining the internal domain motion to near rigidity. Critically, the simulations retain an all-atom description for facile inclusion of site-specific NMR RDC restraints. The result is the rapid generation of conformational ensembles consistent with the RDC data. We illustrate this protocol on human Pin1, a two-domain peptidyl-prolyl isomerase relevant for cancer and Alzheimer's disease. The results include the ensemble of domain orientations sampled by Pin1, as well as those of a dysfunctional variant, I28A-Pin1. The differences between the ensembles corroborate our previous spin relaxation results that showed weakened interdomain contact in the I28A variant relative to wild type. Our protocol extends our abilities to explore the functional significance of protein domain motions.

NMR-based investigations into target DNA search processes of proteins

To perform their function, transcription factors and DNA-repair/modifying enzymes must first locate their targets in the vast presence of nonspecific, but structurally similar sites on genomic DNA. Before reaching their targets, these proteins stochastically scan DNA and dynamically move from one site to another on DNA. Solution NMR spectroscopy provides unique atomic-level insights into the dynamic DNA-scanning processes, which are difficult to gain by any other experimental means. In this review, we provide an introductory overview on the NMR methods for the structural, dynamic, and kinetic investigations of target DNA search by proteins. We also discuss advantages and disadvantages of these NMR methods over other methods such as single-molecule techniques and biochemical approaches.

Structural Characterization of a Heparan Sulfate Pentamer Interacting with LAR-Ig1-2

Leukocyte common antigen-related (LAR) protein is one of the type IIa receptor protein tyrosine phosphatases (RPTPs) that are important for signal transduction in biological processes, including axon growth and regeneration. Glycosaminoglycan chains, including heparan sulfate (HS) and chondroitin sulfate (CS), act as ligands that regulate LAR signaling. Here, we report the structural characterization of the first two immunoglobulin domains (Ig1-2) of LAR interacting with an HS pentasaccharide (GlcNS6S-GlcA-GlcNS3,6S-IdoA2S-GlcNS6S-OME, fondaparinux) using multiple solution-based NMR methods. In the course of the study, we extended an assignment strategy useful for sparsely labeled proteins expressed in mammalian cell culture supplemented with a single type of isotopically enriched amino acid ([¹⁵N]-Lys in this case) by including paramagnetic perturbations to NMR resonances. The folded two-domain structure for LAR-Ig1-2 seen in previous crystal structures has been validated in solution using residual dipolar coupling data, and a combination of chemical shift perturbation on titration of LAR-Ig1-2 with fondaparinux, saturation transfer difference (STD) spectra, and transferred nuclear Overhauser effects (trNOEs) have been employed in the docking program HADDOCK to generate models for the LAR-fondaparinux complex. These models are further analyzed by postprocessing energetic analysis to identify key binding interactions. In addition to providing insight into the ligand interaction mechanisms of type IIa RPTPs and the origin of opposing effects of CS and HS ligands, these results may assist in future design of therapeutic compounds for nervous system repair.

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

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

Tilted, Uninterrupted, Monomeric HIV-1 gp41 Transmembrane Helix from Residual Dipolar Couplings

Cryo-electron microscopy and X-ray crystallography have shown that the pre- and postfusion states of the HIV-1 gp41 viral coat protein, although very different from one another, each adopt C₃ symmetric structures. A stable homotrimeric structure for the transmembrane domain (TM) also was modeled and supported by experimental data. For a C₃ symmetric structure, alignment in an anisotropic medium must be axially symmetric, with the unique axis of the alignment tensor coinciding with the C₃ axis. However, NMR residual dipolar couplings (RDCs) measured under three different alignment conditions were found to be incompatible with C₃ symmetry. Subsequent measurements by paramagnetic relaxation enhancement, analytical ultracentrifugation, and DEER EPR, indicate that the transmembrane domain is monomeric. ¹⁵N NMR relaxation data and RDCs show that TM is highly ordered and uninterrupted for a total length of 32 residues, extending well into the membrane proximal external region.

Unequivocal determination of complex molecular structures using anisotropic NMR measurements

Assignment of complex molecular structures from nuclear magnetic resonance (NMR) data can be prone to interpretational mistakes. Residual dipolar couplings and residual chemical shift anisotropy provide a spatial view of the relative orientations between bonds and chemical shielding tensors, respectively, regardless of separation. Consequently, these data constitute a reliable reporter of global structural validity. Anisotropic NMR parameters can be used to evaluate investigators' structure proposals or structures generated by computer-assisted structure elucidation. Application of the method to several complex structure assignment problems shows promising results that signal a potential paradigm shift from conventional NMR data interpretation, which may be of particular utility for compounds not amenable to x-ray crystallography.

Mollib: a molecular and NMR data analysis software

Mollib is a software framework for the analysis of molecular structures, properties and data with an emphasis on data collected by NMR. It uses an open source model and a plugin framework to promote community-driven development of new and enhanced features. Mollib includes tools for the automatic retrieval and caching of protein databank (PDB) structures, the hydrogenation of biomolecules, the analysis of backbone dihedral angles and hydrogen bonds, and the fitting of residual dipolar coupling (RDC) and residual anisotropic chemical shift (RACS) data. In this article, we release version 1.0 of mollib and demonstrate its application to common molecular and NMR data analyses.

Solution NMR Spectroscopy in Target-Based Drug Discovery

Solution NMR spectroscopy is a powerful tool to study protein structures and dynamics under physiological conditions. This technique is particularly useful in target-based drug discovery projects as it provides protein-ligand binding information in solution. Accumulated studies have shown that NMR will play more and more important roles in multiple steps of the drug discovery process. In a fragment-based drug discovery process, ligand-observed and protein-observed NMR spectroscopy can be applied to screen fragments with low binding affinities. The screened fragments can be further optimized into drug-like molecules. In combination with other biophysical techniques, NMR will guide structure-based drug discovery. In this review, we describe the possible roles of NMR spectroscopy in drug discovery. We also illustrate the challenges encountered in the drug discovery process. We include several examples demonstrating the roles of NMR in target-based drug discoveries such as hit identification, ranking ligand binding affinities, and mapping the ligand binding site. We also speculate the possible roles of NMR in target engagement based on recent processes in in-cell NMR spectroscopy.

Slow conformational exchange and overall rocking motion in ubiquitin protein crystals

Proteins perform their functions in solution but their structures are most frequently studied inside crystals. Here we probe how the crystal packing alters microsecond dynamics, using solid-state NMR measurements and multi-microsecond MD simulations of different crystal forms of ubiquitin. In particular, near-rotary-resonance relaxation dispersion (NERRD) experiments probe angular backbone motion, while Bloch–McConnell relaxation dispersion data report on fluctuations of the local electronic environment. These experiments and simulations reveal that the packing of the protein can significantly alter the thermodynamics and kinetics of local conformational exchange. Moreover, we report small-amplitude reorientational motion of protein molecules in the crystal lattice with an ~3–5° amplitude on a tens-of-microseconds time scale in one of the crystals, but not in others. An intriguing possibility arises that overall motion is to some extent coupled to local dynamics. Our study highlights the importance of considering the packing when analyzing dynamics of crystalline proteins.

X-ray crystallography is the main method for protein structure determination. Here the authors combine solid-state NMR measurements and molecular dynamics simulations and show that crystal packing alters the thermodynamics and kinetics of local conformational exchange as well as overall rocking motion of protein molecules in the crystal lattice.

NMR characterization of HtpG, the E. coli Hsp90, using sparse labeling with 13C-methyl alanine

A strategy for acquiring structural information from sparsely isotopically labeled large proteins is illustrated with an application to the E. coli heat-shock protein, HtpG (high temperature protein G), a 145 kDa dimer. It uses ¹³C-alanine methyl labeling in a perdeuterated background to take advantage of the sensitivity and resolution of Methyl-TROSY spectra, as well as the backbone-centered structural information from ¹H-¹³C residual dipolar couplings (RDCs) of alanine methyl groups. In all, 40 of the 47 expected crosspeaks were resolved and 36 gave RDC data. Assignments of crosspeaks were partially achieved by transferring assignments from those made on individual domains using triple resonance methods. However, these were incomplete and in many cases the transfer was ambiguous. A genetic algorithm search for consistency between predictions based on domain structures and measurements for chemical shifts and RDCs allowed 60% of the 40 resolved crosspeaks to be assigned with confidence. Chemical shift changes of these crosspeaks on adding an ATP analog to the apo-protein are shown to be consistent with structural changes expected on comparing previous crystal structures for apo- and complex- structures. RDCs collected on the assigned alanine methyl peaks are used to generate a new solution model for the apo-protein structure.

Hybrid NMR: A Union of Solution- and Solid-State NMR

Hybrid NMR (hdNMR) is a powerful new tool that combines the strengths of solution- and solid-state NMR to measure dipolar, chemical shift, and quadrupolar tensors in aqueous solution. We introduce the theory of hdNMR and partially randomly oriented (PRO) crystalline hydrogel samples. PRO samples produce randomly oriented spectra with characteristic Pake patterns from the solid state, yet they maintain the high-resolution dispersion of solution NMR experiments. With new pulse sequences, we show how hdNMR can be used to measure with high precision the ¹H^α-¹³C^α dipolar tensor and carboxylate chemical shift anisotropy tensor of aspartate. These measurements contain detailed information on the distribution of electron density, interatomic distances, and the orientation dependence of molecular motion.

Di(ethylene glycol) methyl ether methacrylate (DEGMEMA)-derived gels align small organic molecules in methanol

Residual dipolar couplings (RDCs) constitute an important NMR parameter for structural elucidation in all areas of chemistry. In this study, di(ethylene glycol) methyl ether methacrylate (DEGMEMA)-based gels are introduced as alignment media for the measurement of RDCs of small organic molecules in polar solvents such as methanol. The low viscosity of methanol permits the execution of J-scaled BIRD HSQC experiments that yield very sharp lines in anisotropic conditions. The gels have excellent mechanical properties, and their compression and expansion in the swollen state can be reversed and performed multiple times. This process enables the easy loading and release of analytes. The excellent performance of these new aligning gels is demonstrated by analyzing the structure of the alkaloid retrorsine. Copyright © 2016 John Wiley & Sons, Ltd.

Hybrid Applications of Solution Scattering to Aid Structural Biology

Biomolecular applications of solution X-ray and neutron scattering (SAXS and SANS, respectively) started in late 1960s - early 1970s but were relatively limited in their ability to provide a detailed structural picture and lagged behind what became the two primary methods of experimental structural biology - X-ray crystallography and NMR. However, improvements in both data analysis and instrumentation led to an explosive growth in the number of studies that used small-angle scattering (SAS) for investigation of macromolecular structure, often in combination with other biophysical techniques. Such hybrid applications are nowadays quickly becoming a norm whenever scattering data are used for two reasons. First, it is generally accepted that SAS data on their own cannot lead to a uniquely defined high-resolution structural model, creating a need for supplementing them with information from complementary techniques. Second, solution scattering data are frequently applied in situations when a method such NMR or X-ray crystallography cannot provide a satisfactory structural picture, which makes these additional restraints highly desirable. Maturation of the hybrid bio-SAS approaches brings to light new questions including completeness of the conformational space sampling, model validation, and data compatibility.

Coding and decoding libraries of sequence-defined functional copolymers synthesized via photoligation

Designing artificial macromolecules with absolute sequence order represents a considerable challenge. Here we report an advanced light-induced avenue to monodisperse sequence-defined functional linear macromolecules up to decamers via a unique photochemical approach. The versatility of the synthetic strategy-combining sequential and modular concepts-enables the synthesis of perfect macromolecules varying in chemical constitution and topology. Specific functions are placed at arbitrary positions along the chain via the successive addition of monomer units and blocks, leading to a library of functional homopolymers, alternating copolymers and block copolymers. The in-depth characterization of each sequence-defined chain confirms the precision nature of the macromolecules. Decoding of the functional information contained in the molecular structure is achieved via tandem mass spectrometry without recourse to their synthetic history, showing that the sequence information can be read. We submit that the presented photochemical strategy is a viable and advanced concept for coding individual monomer units along a macromolecular chain.

Protein dynamics and function from solution state NMR spectroscopy

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

Glycosylation Alters Dimerization Properties of a Cell-surface Signaling Protein, Carcinoembryonic Antigen-related Cell Adhesion Molecule 1 (CEACAM1)

Human carcinoembryonic antigen-related cell adhesion molecule 1 (C?/Au: EACAM1) is a cell-surface signaling molecule involved in cell adhesion, proliferation, and immune response. It is also implicated in cancer angiogenesis, progression, and metastasis. This diverse set of effects likely arises as a result of the numerous homophilic and heterophilic interactions that CEACAM1 can have with itself and other molecules. Its N-terminal Ig variable (IgV) domain has been suggested to be a principal player in these interactions. Previous crystal structures of the β-sandwich-like IgV domain have been produced using Escherichia coli-expressed material, which lacks native glycosylation. These have led to distinctly different proposals for dimer interfaces, one involving interactions of ABED β-strands and the other involving GFCC'C″ β-strands, with the former burying one prominent glycosylation site. These structures raise questions as to which form may exist in solution and what the effect of glycosylation may have on this form. Here, we use NMR cross-correlation measurements to examine the effect of glycosylation on CEACAM1-IgV dimerization and use residual dipolar coupling (RDC) measurements to characterize the solution structure of the non-glycosylated form. Our findings demonstrate that even addition of a single N-linked GlcNAc at potential glycosylation sites inhibits dimer formation. Surprisingly, RDC data collected on E. coli expressed material in solution indicate that a dimer using the non-glycosylated GFCC'C″ interface is preferred even in the absence of glycosylation. The results open new questions about what other factors may facilitate dimerization of CEACAM1 in vivo, and what roles glycosylation may play in heterophylic interactions.

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

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