Residual dipolar coupling derived orientational constraints on ligand geometry in a 53 kDa protein-ligand complex

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.

Investigating Protein-Ligand Interactions by Solution Nuclear Magnetic Resonance Spectroscopy

Protein-ligand interactions are of fundamental importance in almost all processes in living organisms. The ligands comprise small molecules, drugs or biological macromolecules and their interaction strength varies over several orders of magnitude. Solution NMR spectroscopy offers a large repertoire of techniques to study such complexes. Here, we give an overview of the different NMR approaches available. The information they provide ranges from the simple information about the presence of binding or epitope mapping to the complete 3 D structure of the complex. NMR spectroscopy is particularly useful for the study of weak interactions and for the screening of binding ligands with atomic resolution.

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.

Modern Technologies of Solution Nuclear Magnetic Resonance Spectroscopy for Three-dimensional Structure Determination of Proteins Open Avenues for Life Scientists

Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for structural studies of chemical compounds and biomolecules such as DNA and proteins. Since the NMR signal sensitively reflects the chemical environment and the dynamics of a nuclear spin, NMR experiments provide a wealth of structural and dynamic information about the molecule of interest at atomic resolution. In general, structural biology studies using NMR spectroscopy still require a reasonable understanding of the theory behind the technique and experience on how to recorded NMR data. Owing to the remarkable progress in the past decade, we can easily access suitable and popular analytical resources for NMR structure determination of proteins with high accuracy. Here, we describe the practical aspects, workflow and key points of modern NMR techniques used for solution structure determination of proteins. This review should aid NMR specialists aiming to develop new methods that accelerate the structure determination process, and open avenues for non-specialist and life scientists interested in using NMR spectroscopy to solve protein structures.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.

Fully deuterated magnetically oriented system based on fatty acid direct hexagonal phases

There is strong demand in the field of NMR for simple oriented lipid supramolecular assemblies, the constituents of which can be fully deuterated, for specifically studying the structure of host protonated molecules (e.g., peptides, proteins...) in a lipid environment. Also, small-angle neutron scattering (SANS) in fully deuterated oriented systems is powerful for gaining information on protonated host molecules in a lipid environment by using the contrast proton/deuterium method. Here we report on a very simple system made of fatty acids (dodecanoic and tetradecanoic) and ethanolamine in water. All components of this system can be obtained commercially as perdeuterated. Depending on the molar ratio and the concentration, the system self-assembles at room temperature into a direct hexagonal phase that is oriented by moderate magnetic fields of a few tesla. The orientation occurs within the magnetic field upon cooling the system from its higher-temperature isotropic phase: the lipid cylinders of the hexagonal phase become oriented parallel to the field. This is shown by solid-state NMR using either perdeuterated fatty acids or ethanolamine. This system bears strong interest for studying host protonated molecules but also in materials chemistry for building oriented solid materials.

An introduction to biological NMR spectroscopy

NMR spectroscopy is a powerful tool for biologists interested in the structure, dynamics, and interactions of biological macromolecules. This review aims at presenting in an accessible manner the requirements and limitations of this technique. As an introduction, the history of NMR will highlight how the method evolved from physics to chemistry and finally to biology over several decades. We then introduce the NMR spectral parameters used in structural biology, namely the chemical shift, the J-coupling, nuclear Overhauser effects, and residual dipolar couplings. Resonance assignment, the required step for any further NMR study, bears a resemblance to jigsaw puzzle strategy. The NMR spectral parameters are then converted into angle and distances and used as input using restrained molecular dynamics to compute a bundle of structures. When interpreting a NMR-derived structure, the biologist has to judge its quality on the basis of the statistics provided. When the 3D structure is a priori known by other means, the molecular interaction with a partner can be mapped by NMR: information on the binding interface as well as on kinetic and thermodynamic constants can be gathered. NMR is suitable to monitor, over a wide range of frequencies, protein fluctuations that play a crucial role in their biological function. In the last section of this review, intrinsically disordered proteins, which have escaped the attention of classical structural biology, are discussed in the perspective of NMR, one of the rare available techniques able to describe structural ensembles. This Tutorial is part of the International Proteomics Tutorial Programme (IPTP 16 MCP).

Advances in nuclear magnetic resonance for drug discovery

Nuclear Magnetic Resonance (NMR) techniques are widely used in the drug discovery process. The primary feature exploited in these investigations is the large difference in mass between drugs and receptors (usually proteins) and the effect this has on the rotational or translational correlation times for drugs bound to their targets. Many NMR parameters, such as the diffusion coefficient, spin diffusion, nuclear Overhauser enhancement, and transverse and longitudinal relaxation times, are strong functions of either the overall tumbling or translation of molecules in solution. This has led to the development of a wide variety of NMR techniques applicable to the elucidation of protein and nucleic acid structure in solution, the screening of drug candidates for binding to a target of choice, and the study of the conformational changes which occur in a target upon drug binding. High-throughput screening by NMR methods has recently received a boost from the introduction of sophisticated computational techniques for reducing the time needed for the acquisition of the primary NMR data for multidimensional studies.

Characterizing weak protein-protein complexes by NMR residual dipolar couplings

Protein-protein interactions occur with a wide range of affinities from tight complexes characterized by femtomolar dissociation constants to weak, and more transient, complexes of millimolar affinity. Many of the weak and transiently formed protein-protein complexes have escaped characterization due to the difficulties in obtaining experimental parameters that report on the complexes alone without contributions from the unbound, free proteins. Here, we review recent developments for characterizing the structures of weak protein-protein complexes using nuclear magnetic resonance spectroscopy with special emphasis on the utility of residual dipolar couplings.

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.

NMR as a unique tool in assessment and complex determination of weak protein-protein interactions

Protein-protein interactions are crucial for a wide variety of biological processes. These interactions range from high affinity (K (d)<nM) to very low affinity (K (d)>mM). While much is known about the nature of high affinity protein complexes, our knowledge about structural characteristics of weak protein-protein interactions (wPPIs) remains limited: in addition to the technical difficulties associated with their investigation, historically wPPIs used to be considered physiologically irrelevant. However, emerging evidence suggests that wPPIs, either in the form of intact protein complexes or as part of large molecular machineries, are fundamentally important for promoting rapid on/off switches of signal transduction, reversible cell-cell contacts, transient assembly/disassembly of signaling complexes, and enzyme-substrate recognition. Therefore an atomic-level elucidation of wPPIs is vital to understanding a cornucopia of diverse cellular events. Nuclear magnetic resonance (NMR) is famous for its unique abilities to study wPPIs and, by utilization of the new technical developments combined with sparse data based computational analysis, it now allows rapid identification and structural characterization of wPPIs. Here we present our perspective on the NMR methods employed.

A novel oriented system made of fatty acid hexagonal phases with tuneable orientation

There is a strong demand in the field of solid state NMR for oriented lipid supramolecular assemblies. This is mainly devoted to biophysical structural studies or materials chemistry because the NMR signal depends on the orientation. Here we report a novel system made of a fatty acid hexagonal phase which self orient in the magnetic field. The orientation occurs within the magnetic field upon cooling the system from its isotropic phase. The cylinders of the hexagonal phase are then oriented parallel to the field. We take advantage that the hexagonal phase is a gel, i.e., the orientation is maintained fixed within the sample tube to investigate the orientational dependence of the deuterium solid state NMR signal using deuterated fatty acids and D(2)O by manually rotating the sample tube within the coil probe. As expected, the oriented signal follows the low |3cos(2)theta-1| where theta is the angle between the long cylindrical axis and the field. We expect this system to be of interest in materials chemistry and structural biology.

Accurate characterization of weak macromolecular interactions by titration of NMR residual dipolar couplings: application to the CD2AP SH3-C:ubiquitin complex

The description of the interactome represents one of key challenges remaining for structural biology. Physiologically important weak interactions, with dissociation constants above 100 μM, are remarkably common, but remain beyond the reach of most of structural biology. NMR spectroscopy, and in particular, residual dipolar couplings (RDCs) provide crucial conformational constraints on intermolecular orientation in molecular complexes, but the combination of free and bound contributions to the measured RDC seriously complicates their exploitation for weakly interacting partners. We develop a robust approach for the determination of weak complexes based on: (i) differential isotopic labeling of the partner proteins facilitating RDC measurement in both partners; (ii) measurement of RDC changes upon titration into different equilibrium mixtures of partially aligned free and complex forms of the proteins; (iii) novel analytical approaches to determine the effective alignment in all equilibrium mixtures; and (iv) extraction of precise RDCs for bound forms of both partner proteins. The approach is demonstrated for the determination of the three-dimensional structure of the weakly interacting CD2AP SH3-C:Ubiquitin complex (K_d = 132 ± 13 μM) and is shown, using cross-validation, to be highly precise. We expect this methodology to extend the remarkable and unique ability of NMR to study weak protein–protein complexes.

Structure determination of protein-ligand complexes by transferred paramagnetic shifts

Rational drug design depends on the knowledge of the three-dimensional (3D) structure of complexes between proteins and lead compounds of low molecular weight. A novel nuclear magnetic resonance (NMR) spectroscopy strategy based on the paramagnetic effects from lanthanide ions allows the rapid determination of the 3D structure of a small ligand molecule bound to its protein target in solution and, simultaneously, its location and orientation with respect to the protein. The method relies on the presence of a lanthanide ion in the protein target and on fast exchange between bound and free ligand. The binding affinity of the ligand and the paramagnetic effects experienced in the bound state are derived from concentration-dependent (1)H and (13)C spectra of the ligand at natural isotopic abundance. Combined with prior knowledge of the crystal or solution structure of the protein and of the magnetic susceptibility tensor of the lanthanide ion, the paramagnetic data define the location and orientation of the bound ligand molecule with respect to the protein from simple 1D NMR spectra. The method was verified with the ternary 30 kDa complex between the lanthanide-labeled N-terminal domain of the epsilon exonuclease subunit from the Escherichia coli DNA polymerase III, the subunit theta, and thymidine. The binding mode of thymidine was found to be very similar to that of thymidine monophosphate present in the crystal structure.

Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition

Many regulatory RNAs undergo large changes in structure upon recognition of proteins and ligands, but the mechanism by which this occurs remains poorly understood. Using NMR residual dipolar coupling (RDCs), we characterized Na+-induced changes in the structure and dynamics of the bulge-containing HIV-1 transactivation response element (TAR) RNA that mirrors changes induced by small molecules bearing a different number of cationic groups. Increasing the Na+ concentration from 25 to 320 mM led to a continuous reduction in the average inter-helical bend angle (from 46 degrees to 22 degrees ), inter-helical twist angle (from 66 degrees to -18 degrees ), and inter-helix flexibility (as measured by an increase in the internal generalized degree of order from 0.56 to 0.74). Similar conformational changes were observed with Mg2+, indicating that nonspecific electrostatic interactions drive the conformational transition, although results also suggest that Na+ and Mg2+ may associate with TAR in distinct modes. The transition can be rationalized on the basis of a population-weighted average of two ensembles comprising an electrostatically relaxed bent and flexible TAR conformation that is weakly associated with counterions and a globally rigid coaxial conformation that has stronger electrostatic potential and association with counterions. The TAR inter-helical orientations that are stabilized by small molecules fall around the metal-induced conformational pathway, indicating that counterions may help predispose the TAR conformation for target recognition. Our results underscore the intricate sensitivity of RNA conformational dynamics to environmental conditions and demonstrate the ability to detect subtle conformational changes using NMR RDCs.

Solution NMR of large molecules and assemblies

Solution NMR spectroscopy represents a powerful tool for examining the structure and function of biological macromolecules. The advent of multidimensional (2D-4D) NMR, together with the widespread use of uniform isotopic labeling of proteins and RNA with the NMR-active isotopes, 15N and 13C, opened the door to detailed analyses of macromolecular structure, dynamics, and interactions of smaller macromolecules (< approximately 25 kDa). Over the past 10 years, advances in NMR and isotope labeling methods have expanded the range of NMR-tractable targets by at least an order of magnitude. Here we briefly describe the methodological advances that allow NMR spectroscopy of large macromolecules and their complexes and provide a perspective on the wide range of applications of NMR to biochemical problems.

Structure-based drug design: NMR-based approach for ligand-protein interactions

The realization of the powerfulness in analyzing ligand-protein interactions at the atomic resolution has made NMR techniques increasingly attractive in drug discovery and development. With some significant new method developments during the past few years, NMR-based approaches will undoubtedly be helpful in high throughput screening and in providing structural and interaction information beneficial for new drug developments. Here in this review, instead of providing an exhaustive account of applications of NMR, we will seek to highlight some key points related to the solution state NMR methods for extracting information from both ligands and proteins.:

NMR spectroscopic characterization of the membrane affinity of polyols

Residual dipolar couplings (RDCs) are applied here for the analysis of weak, transient binding events between phosphatidylcholine bilayers and polyols. Large signal responses are observed even for low percentages of 'ligand-receptor complexes', making RDCs a sensitive tool for the analysis of molecular recognition events. The different degree of alignment in solution can be compared as a result of the calculation of the alignment tensor elements. By varying polarity and/or charge of the molecules under investigation, nonspecific hydrophobic effects can be excluded.

Partial 13C isotopic enrichment of nucleoside monophosphates: useful reporters for NMR structural studies

Analysis of the ¹³C isotopic labeling patterns of nucleoside monophosphates (NMPs) extracted from Escherichia coli grown in a mixture of C-1 and C-2 glucose is presented. By comparing our results to previous observations on amino acids grown in similar media, we have been able to rationalize the labeling pattern based on the well-known biochemistry of nucleotide biosynthesis. Except for a few notable absences of label (C4 in purines and C3′ in ribose) and one highly enriched site (C1′ in ribose), most carbons are randomly enriched at a low level (an average of 13%). These sparsely labeled NMPs give less complex NMR spectra than their fully isotopically labeled analogs due to the elimination of most ¹³C–¹³C scalar couplings. The spectral simplicity is particularly advantageous when working in ordered systems, as illustrated with guanosine diphosphate (GDP) bound to ADP ribosylation factor 1 (ARF1) aligned in a liquid crystalline medium. In this system, the absence of scalar couplings and additional long-range dipolar couplings significantly enhances signal to noise and resolution.

Weak protein-protein interactions as probed by NMR spectroscopy

Weak protein-protein interactions (PPIs) are fundamental to many cellular processes, such as reversible cell-cell contact, rapid enzyme turnover and transient assembly and/or reassembly of large signaling complexes. However, structural and functional characterizations of weak PPIs have been technically challenging and lagged behind those for strong PPIs. Here, we describe nuclear magnetic resonance (NMR) spectroscopy as a highly effective tool for unraveling the atomic details of weak PPIs. We highlight the recent advances of how NMR can be used to rapidly detect and structurally determine extremely weak PPIs (K(d)>10(-4)M). Coupled with functional approaches, NMR has the potential to look into a wide variety of biologically important weak PPIs at the detailed molecular level, thereby facilitating a thorough view of how proteins function in living cells.

Application of residual dipolar couplings in organic compounds

Residual dipolar couplings (RDCs) induced by anisotropic media are a powerful tool for the structure determination of biomolecules through NMR spectroscopy. Recent advances have proven it to be a valuable tool for determination of the stereochemistry of organic molecules. By simple inspection or order matrix calculations, RDCs provide unambiguous information about the relative configurations or complete stereochemistry of organic compounds.

Rapid analysis of large protein-protein complexes using NMR-derived orientational constraints: the 95 kDa complex of LpxA with acyl carrier protein

Characterization of protein-protein interactions that are critical to the specific function of many biological systems has become a primary goal of structural biology research. Analysis of these interactions by structural techniques is, however, challenging due to inherent limitations of the techniques and because many of the interactions are transient, and suitable complexes are difficult to isolate. In particular, structural studies of large protein complexes by traditional solution NMR methods are difficult due to a priori requirement of extensive assignments and a large number of intermolecular restraints for the complex. An approach overcoming some of these challenges by utilizing orientational restraints from residual dipolar couplings collected on solution NMR samples is presented. The approach exploits existing structures of individual components, including the symmetry properties of some of these structures, to assemble rapidly models for relatively large protein-protein complexes. An application is illustrated with a 95 kDa homotrimeric complex of the acyltransferase protein, LpxA (UDP-N-acetylglucosamine acyltransferase), and acyl carrier protein. LpxA catalyzes the first step in the biosynthesis of the lipid A component of lipopolysaccharide in Gram-negative bacteria. The structural model generated for this complex can be useful in the design of new anti-bacterial agents that inhibit the biosynthesis of lipid A.

Weak alignment of membrane proteins in stressed polyacrylamide gels

Residual dipolar couplings are important as angular constraints for the structure determination of membrane proteins in micelles. Strained polyacrylamide gels are one of the few available mechanisms available for inducing the requisite weak alignment for these samples. However, their use is frequently limited by the ability to incorporate proteins and buffer solutions into the gel matrix. The implementation of several methods of incorporating membrane proteins into gels are described. Conditions for copolymerizing the protein in the absence of a change in pH are detailed. Electrophoresis is also shown to be a useful method to incorporate proteins. Weak alignment of the protein-micelle complex in the gel matrix is subsequently achieved using either vertical or radial compression. The magnitude of alignment can be controlled by altering the gel concentration, the acrylamide/bisacrylamide ratio, and the compression ratio. The alignment tensor can be altered relative to uncharged polyacrylamide gels by copolymerizing samples with acrylamide/acrylic acid to incorporate negative charges in the strained polyacrylamide gel to provide an alternate orientation.

Residual dipolar couplings in NMR structure analysis

Residual dipolar couplings (RDCs) have recently emerged as a new tool in nuclear magnetic resonance (NMR) with which to study macromolecular structure and function in a solution environment. RDCs are complementary to the more conventional use of NOEs to provide structural information. While NOEs are local-distance restraints, RDCs provide long-range orientational information. RDCs are now widely utilized in structure calculations. Increasingly, they are being used in novel applications to address complex issues in structural biology such as the accurate determination of the global structure of oligonucleotides and the relative orientation of protein domains. This review briefly describes the theory and methods for obtaining RDCs and then describes the range of biological applications where RDCs have been used.

REDCAT: a residual dipolar coupling analysis tool

Recent advancements in the utilization of residual dipolar couplings (RDCs) as a means of structure validation and elucidation have demonstrated the need for, not only a more user friendly, but also a more powerful RDC analysis tool. In this paper, we introduce a software package named REsidual Dipolar Coupling Analysis Tool (REDCAT) designed to address the above issues. REDCAT is a user-friendly program with its graphical-user-interface developed in Tcl/Tk, which is highly portable. Furthermore, the computational engine behind this GUI is written in C/C++ and its computational performance is therefore excellent. The modular implementation of REDCAT's algorithms, with separation of the computational engine from the graphical engine allows for flexible and easy command line interaction. This feature can be utilized for the design of automated data analysis sessions. Furthermore, this software package is portable to Linux clusters for high throughput applications. In addition to basic utilities to solve for order tensors and back calculate couplings from a given order tensor and proposed structure, a number of improved algorithms have been incorporated. These include the proper sampling of the Null-space (when the system of linear equations is under-determined), more sophisticated filters for invalid order-tensor identification, error analysis for the identification of the problematic measurements and simulation of the effects of dynamic averaging processes.

Use of residual dipolar couplings as restraints inab initio protein structure prediction

NMR residual dipolar couplings (RDCs), in the form of the projection angles between the respective internuclear bond vectors, are used as structural restraints in the ab initio structure prediction of a test set of six proteins. The restraints are applied using a recently developed SICHO (SIde-CHain-Only) lattice protein model that employs a replica exchange Monte Carlo (MC) algorithm to search conformational space. Using a small number of RDC restraints, the quality of the predicted structures is improved as reflected by lower RMSD/dRMSD (root mean square deviation/distance root mean square deviation) values from the corresponding native structures and by the higher correlation of the most cooperative mode of motion of each predicted structure with that of the native structure. The latter, in particular, has possible implications for the structure-based functional analysis of predicted structures.

Structural characterization of a mannose-binding protein-trimannoside complex using residual dipolar couplings

The ligand-binding properties of a 53 kDa homomultimeric trimer from mannose-binding protein (MBP) have been investigated using residual dipolar couplings (RDCs) that are easily measured from NMR spectra of the ligand and isotopically labeled protein. Using a limited set of 1H-15N backbone amide NMR assignments for MBP and orientational information derived from the RDC measurements in aligned media, an order tensor for MBP has been determined that is consistent with symmetry-based predictions of an axially symmetric system. 13C-1H couplings for a bound trisaccharide ligand, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside (trimannoside) have been determined at natural abundance and used as orientational constraints. The bound ligand geometry and orientational constraints allowed docking of the trimannoside ligand in the binding site of MBP to produce a structural model for MBP-oligosaccharide interactions.

Conformational differences in liganded and unliganded states of Galectin-3

The conformation of the carbohydrate recognition domain of Galectin-3, a lectin known to bind galactose containing oligosaccharides in mammalian systems, has been investigated in the absence of ligand and in the presence of N-acetylactosamine. A new methodology based on the measurement of residual dipolar couplings from NMR spectra has been used to characterize differences in protein structure along the backbone in the presence and absence of ligand, as well as the binding geometry of the ligand itself. The data on the ligand are consistent with the ligand binding geometry found in a crystal structure of the complexed state. However, a significant rearrangement of backbone loops near the binding site appears to occur in the absence of ligand. The implications for ligand specificity and protein functionality are discussed.

Refined solution structure of the iron-responsive element RNA using residual dipolar couplings

The iron-responsive element (IRE) is a 30nt RNA motif located in the non-coding regions of mRNAs of proteins involved in iron regulation. In humans, the IRE plays a direct role in the control of iron levels by post-transcriptional regulation of the ferritin and transferrin receptor proteins through highly specific recognition by IRE-binding proteins. The IRE fold is representative of many RNA motifs that contain helical domains separated by a bulge or internal loop. The global structures of such extended multi-domain RNAs are not well defined by conventional NMR-distance and torsion angle structural restraints. Residual dipolar couplings (RDCs) are employed here to better define the global structure of the IRE RNA in solution. RDCs contain valuable long-range structural information that compliments the short-range structural data derived from standard NOE-distance and torsion angle restraints. Several approaches for estimating alignment tensor parameters and incorporating RDCs into RNA structure determinations are compared. Both the local and global structure of the IRE are improved significantly by refinement with RDCs. These RDC refinements provide insight on the conformational dynamics of the IRE. These studies highlight some issues that need to be addressed when incorporating RDCs in solution structure determinations of nucleic acids. The approach used here should prove valuable for structure determinations of various multi-domain systems.

Conformation of glycomimetics in the free and protein-bound state: structural and binding features of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man

The conformational properties of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man in solution have been carefully analyzed by a combination of NMR spectroscopy and time-averaged restrained molecular dynamics. It has been found that both the alpha-1,3- and the alpha-1,6-glycosidic linkages show a major conformational averaging. Unusual Phi ca. 60 degrees orientations for both Phi torsion angles are found. Moreover, a major conformational distinction between the natural compound and the glycomimetic affects to the behavior of the omega(16) torsion angle around the alpha-1 --> 6-linkage. Despite this increased flexibility, the C-glycosyl analogue is recognized by three mannose binding lectins, as shown by NMR (line broadening, TR-NOE, and STD) and surface plasmon resonance (SPR) methods. Moreover, a process of conformational selection takes place, so that these lectins probably bind the glycomimetic similarly to the way they recognize the natural analogue. Depending upon the architecture and extension of the binding site of the lectin, loss or gain of binding affinity with respect to the natural analogue is found.

Detecting protein kinase recognition modes of calmodulin by residual dipolar couplings in solution NMR

Calmodulin-regulated serine/threonine kinases (CaM kinases) play crucial roles in Ca2+-dependent signaling transduction pathways in eukaryotes. Despite having a similar overall molecular architecture of catalytic and regulatory domains, CaM kinases employ different binding modes for Ca2+/CaM recruitment which is required for their activation. Here we present a residual dipolar coupling (RDC)-based NMR approach to characterizing the molecular recognition of CaM with five different CaM kinases. Our analyses indicate that CaM kinase I and likely IV use the same CaM binding mode as myosin light chain kinase (1-14 motif), distinct from those of CaM kinase II (1-10 motif) and CaM kinase kinase (1-16- motif). This NMR approach provides an efficient experimental guide for homology modeling and structural characterization of CaM-target complexes.

Residual dipolar couplings in nucleic acid structure determination

Solution NMR spectroscopy of nucleic acids has been limited by the short-range nature of the nuclear Overhauser effect and scalar coupling restraints normally used in structure determination. The addition of residual dipolar couplings, obtained from slightly oriented mixtures, provides bond vector angles relative to a universal alignment tensor. The accurate determination of helix curvature, domain orientation and the stoichiometry of homomultimeric nucleic acid complexes is now possible.

Conformational studies of blood group A and blood group B oligosaccharides using NMR residual dipolar couplings

The conformations of two synthetic trisaccharides of blood group A and B (alpha-L-Fucp-(1-->2)-[alpha-D-GalpNAc-(1-->3)]-alpha-D-Galp and alpha-L-Fucp-(1-->2)-[alpha-D-Galp-(1-->3)]-alpha-D-Galp, respectively) and of a type A tetrasaccharide alditol, Fucp-(1-->2)-[alpha-D-GalpNAc-(1-->3)]-beta-D-Galp-(1-->3)-GalNAc-ol, were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. The conformations of the three oligosaccharides were analyzed by generating thousands of structures using a Monte-Carlo method. Two different strategies were applied to calculate theoretical dipolar couplings for these structures. In the first method, the orientation of the molecule was calculated from the optimal fit of the molecular model to the experimental data, while in the second method the orientation tensor was calculated directly from the moment of inertia of the molecular model. Both methods of analysis give similar results but with slightly better agreement with experiment for the former one. The analysis of the results implies a single unique conformation for both blood group epitopes in solution in disagreement with theoretical models suggesting the existence of two conformers in solution.

Conformation of a trimannoside bound to mannose-binding protein by nuclear magnetic resonance and molecular dynamics simulations

A model of the carbohydrate recognition domain of the serum form of mannose-binding protein (MBP) from rat complexed with methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside is presented. Allowed conformations for the bound sugar were derived from simulated annealing protocols incorporating distance restraints computed from transferred NOESY spectra. The resulting sugar conformations were then modeled into the MBP binding site, and these models of the complex were refined using molecular dynamics (MD) simulations in the presence of solvent water. These studies indicate that only one of the two major conformations of the alpha(1-->6) linkage found in solution is significantly populated in the bound state (omega = 60 degrees ), whereas the alpha(1-->3) linkage samples at least two states, similar to its behavior in free solution. The bound conformation allows direct hydrogen bonds to form between the sugar and K182 of MBP, in addition to other water-mediated hydrogen bonds. Estimates of binding constants of candidate complexes based on changes in solvent-accessible surface areas upon binding support the NMR and MD results. These estimates further suggest that the enthalpic gains of the additional sugar-MBP interactions in a trisaccharide as opposed to a monosaccharide are offset by entropic penalties, offering an explanation for previous binding data.

Conformational studies of Lewis X and Lewis A trisaccharides using NMR residual dipolar couplings

The conformations of the histo-blood group carbohydrate antigens Lewis X (Le(x)) and Lewis A (Le(a)) were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. A strategy for rapid calculation of the difference between theoretical and experimental dipolar couplings of a large number of model structures generated by computer simulations was developed, resulting in an accurate model structure for the compounds. Monte Carlo simulations were used to generate models for the trisaccharides, and orientations of each model were sought that could reproduce the experimental residual dipolar coupling values. For both, Le(a) and Le(x), single low energy models giving excellent agreement with experiment were found, implying a compact rigidly folded conformation for both trisaccharides. The new approach was also applied to the pentasaccharides lacto-N-fucopentaose 2 (LNF-2) and lacto-N-fucopentaose 3 (LNF-3) proving its consistency and robustness. For describing the conformation of tightly folded oligosaccharides, a definition for characterization of ring planes in pyranoside chairs is proposed and applied to the analysis of the relation between the fucose and galactose residues in the epitopes, revealing the structural similarity between them.

Refinement of local and long-range structural order in theophylline-binding RNA using (13)C-(1)H residual dipolar couplings and restrained molecular dynamics

13C-(1)H residual dipolar couplings (RDC) have been measured for the bases and sugars in the theophylline-binding RNA aptamer, dissolved in filamentous phage medium, and used to investigate the long-range structural and dynamic behavior of the molecule in the solution state. The orientation dependent RDC provide additional restraints to further refine the overall structure of the RNA-theophylline complex, whose long-range order was poorly defined in the NOE-based structural ensemble. Structure refinement using RDC normally assumes that molecular alignment can be characterized by a single tensor and that the molecule is essentially rigid. To address the validity of this assumption for the complex of interest, we have analyzed distinct domains of the RNA molecule separately, so that local structure and alignment tensors experienced by each region are independently determined. Alignment tensors for the stem regions of the molecule were allowed to float freely during a restrained molecular dynamics structure refinement protocol and found to converge to similar magnitudes. During the second stage of the calculation, a single alignment tensor was thus applied for the whole molecule and an average molecular conformation satisfying all experimental data was determined. Semirigid-body molecular dynamics calculations were used to reorient the refined helical regions to a relative orientation consistent with this alignment tensor, allowing determination of the global conformation of the molecule. Simultaneously, the local structure of the theophylline-binding core of the molecule was refined under the influence of this common tensor. The final ensemble has an average pairwise root mean square deviation of 1.50 +/- 0.19 A taken over all heavy atoms, compared to 3.5 +/- 1.1 A for the ensemble determined without residual dipolar coupling. This study illustrates the importance of considering both the local and long-range nature of RDC when applying these restraints to structure refinements of nucleic acids.

A novel approach for assessing macromolecular complexes combining soft-docking calculations with NMR data

We present a novel and efficient approach for assessing protein-protein complex formation, which combines ab initio docking calculations performed with the protein docking algorithm BiGGER and chemical shift perturbation data collected with heteronuclear single quantum coherence (HSQC) or TROSY nuclear magnetic resonance (NMR) spectroscopy. This method, termed "restrained soft-docking," is validated for several known protein complexes. These data demonstrate that restrained soft-docking extends the size limitations of NMR spectroscopy and provides an alternative method for investigating macromolecular protein complexes that requires less experimental time, effort, and resources. The potential utility of this novel NMR and simulated docking approach in current structural genomic initiatives is discussed.

Dipolar couplings as a probe of molecular dynamics and structure in solution

The introduction of residual dipolar coupling methodology has increased the scope of structural biological problems that can be addressed by NMR spectroscopy. Conformational changes, the relative orientation of domains, and intermolecular complexes can now be characterized accurately and rapidly using NMR. The development of residual dipolar coupling methodology for the rapid recognition of homologous protein folds and for studies of submillisecond timescale dynamics has also seen considerable progress.

Partial alignment of biomolecules: an aid to NMR characterization

Partial alignment of biomolecules in solution has added a new dimension to structural investigation by high-resolution NMR methods. Applications to proteins, nucleic acids and carbohydrates now abound. Limitations initially associated with compatibility of biomolecules with the liquid-crystal media commonly used to achieve alignment have begun to disappear. This is, in part, a result of the introduction of a wide variety of new media. Future applications to biologically important problems such as the structural organization of multi-domain proteins and multi-protein assemblies look very promising.

Field- and phage-induced dipolar couplings in a homodimeric DNA quadruplex: relative orientation of G.(C-A) triad and G-tetrad motifs and direct determination of C2 symmetry axis orientation

We present a new NMR procedure for determining the three-dimensional fold of C2-symmetric nucleic acid homodimers that relies on long-range orientational constraints derived from the measurement of two independent sets of residual dipolar couplings under two alignment conditions. The application is demonstrated on an (15)N/(13)C-enriched deoxyoligonucleotide sequence, d(G-G-G-T-T-C-A-G-G), shown previously to dimerize into a quadruplex in solution and form a pair of G.(C-A) triads and G-G-G-G tetrads (G-tetrad) motifs. One-bond (1)H-(15)N ((1)D(NH)) and (1)H-(13)C ((1)D(CH)) residual dipolar couplings have been measured between nuclei in the bases of these motifs using bacteriophage as an ordering medium, and under direct magnetic field alignment (800 MHz). By combining the two dipolar data sets in an order matrix analysis, the orientation of the G.(C-A) triad relative to the G-tetrad within a contiguous monomeric unit can directly be determined, even in the presence of interstrand/intrastrand NOE ambiguity. We further demonstrate that the orientation of the C2-axis of molecular symmetry in the homodimer relative to the G.(C-A) triad and G-tetrad motifs can unambiguously be determined using the two sets of independent dipolar coupling measurements. The three-dimensional fold of the homodimer determined using this procedure is very regular and in excellent agreement with a previously determined high-resolution NOE-based NMR structure, where interstrand/intrastrand NOEs were treated as ambiguous and where noncrystallographic symmetry constraints were implicitly imposed during the structure calculation.

Conformational analysis of a flexible oligosaccharide using residual dipolar couplings

We present a new approach to the analysis of the conformational and the motional properties of an oligosaccharide, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside. The approach relies on an order matrix analysis of residual dipolar couplings in the solution state. By combining a number of different types of couplings, (1)D(CH), (2)D(CH), and D(HH), an order matrix is solved for each ring of the trimannoside. The resulting order parameters indicate the internal motion at the alpha (1,3) linkage to be limited, while significant motion is suggested at the alpha (1,6) linkage. Two structures for the trimannoside were determined by aligning the order tensor principal axes obtained from two different orienting media, bicelles and phage. The very similar conformations at the alpha (1,3) linkage of these two structures confirm that the internal motion at the alpha (1,3) linkage is small and the conformation is a good representation of a single preferred structure. The different conformations at the alpha (1,6) linkage suggest that the motional amplitudes are large and the conformations must be viewed as virtual conformers. Compared with traditional NMR methods, data acquisition is easy and data analysis is straightforward.

Solution conformations of a trimannoside from nuclear magnetic resonance and molecular dynamics simulations

N-linked oligosaccharides often act as ligands for receptor proteins in a variety of cell recognition processes. Knowledge of the solution conformations, as well as protein-bound conformations, of these oligosaccharides is required to understand these important interactions. In this paper we present a model for the solution conformations sampled by a simple trimannoside, methyl 3, 6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside, which contains two of the most commonly found glycosidic linkages in N-linked oligosaccharides. This model was derived from simulated annealing protocols incorporating distance restraints extracted from NOESY spectra along with torsional restraints computed from three-bond (1)H-(13)C coupling constants measured across the glycosidic bonds. The model was refined in light of unrestrained molecular dynamics simulations conducted in the presence of solvent water. The resulting model depicts a molecule undergoing conformational averaging in solution, adopting four major and two minor conformations. The four major conformations arise from a pair of two-state transitions, one each at the alpha(1-->3) and alpha(1-->6) linkages, whereas the minor conformations result from an additional transition of the alpha(1-->6) linkage. Our data also suggest that the alpha(1-->3) transition is fast and changes the molecular shape slightly, whereas the alpha(1-->6) is much slower and alters the molecular shape dramatically.