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

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

Solution NMR Experiment for Measurement of (15)N-(1)H Residual Dipolar Couplings in Large Proteins and Supramolecular Complexes

NMR residual dipolar couplings (RDCs) are exquisite probes of protein structure and dynamics. A new solution NMR experiment named 2D SE2 J-TROSY is presented to measure N-H RDCs for proteins and supramolecular complexes in excess of 200 kDa. This enables validation and refinement of their X-ray crystal and solution NMR structures and the characterization of structural and dynamic changes occurring upon complex formation. Accurate N-H RDCs were measured at 750 MHz (1)H resonance frequency for 11-mer 93 kDa (2)H,(15)N-labeled Trp RNA-binding attenuator protein tumbling with a correlation time τc of 120 ns. This is about twice as long as that for the most slowly tumbling system, for which N-H RDCs could be measured, so far, and corresponds to molecular weights of ∼200 kDa at 25 °C. Furthermore, due to the robustness of SE2 J-TROSY with respect to residual (1)H density from exchangeable protons, increased sensitivity at (1)H resonance frequencies around 1 GHz promises to enable N-H RDC measurement for even larger systems.

Solution structure of the TatB component of the twin-arginine translocation system

The twin-arginine protein transport (Tat) system translocates fully folded proteins across lipid membranes. In Escherichia coli, the Tat system comprises three essential components: TatA, TatB and TatC. The protein translocation process is proposed to initiate by signal peptide recognition and substrate binding to the TatBC complex. Upon formation of the TatBC-substrate protein complex, the TatA subunits are recruited and form the protein translocation pore. Experimental evidences suggest that TatB forms a tight complex with TatC at 1:1 molar ratio and the TatBC complex contains multiple copies of both proteins. Cross-linking experiments demonstrate that TatB functions in tetrameric units and interacts with both TatC and substrate proteins. However, structural information of the TatB protein is still lacking, and its functional mechanism remains elusive. Herein, we report the solution structure of TatB in DPC micelles determined by Nuclear Magnetic Resonance (NMR) spectroscopy. Overall, the structure shows an extended 'L-shape' conformation comprising four helices: a transmembrane helix (TMH) α1, an amphipathic helix (APH) α2, and two solvent exposed helices α3 and α4. The packing of TMH and APH is relatively rigid, whereas helices α3 and α4 display notably higher mobility. The observed floppiness of helices α3 and α4 allows TatB to sample a large conformational space, thus providing high structural plasticity to interact with substrate proteins of different sizes and shapes.

Molecular mechanisms of inhibition of influenza by surfactant protein D revealed by large-scale molecular dynamics simulation

Surfactant protein D (SP-D), a mammalian C-type lectin, is the primary innate inhibitor of influenza A virus (IAV) in the lung. Interactions of SP-D with highly branched viral N-linked glycans on hemagglutinin (HA), an abundant IAV envelope protein and critical virulence factor, promote viral aggregation and neutralization through as yet unknown molecular mechanisms. Two truncated human SP-D forms, wild-type (WT) and double mutant D325A+R343V, representing neck and carbohydrate recognition domains are compared in this study. Whereas both WT and D325A+R343V bind to isolated glycosylated HA, WT does not inhibit IAV in neutralization assays; in contrast, D325A+R343V neutralization compares well with that of full-length native SP-D. To elucidate the mechanism for these biochemical observations, we have determined crystal structures of D325A+R343V in the presence and absence of a viral nonamannoside (Man9). On the basis of the D325A+R343V-Man9 structure and other crystallographic data, models of complexes between HA and WT or D325A+R343V were produced and subjected to molecular dynamics. Simulations reveal that whereas WT and D325A+R343V both block the sialic acid receptor site of HA, the D325A+R343V complex is more stable, with stronger binding caused by additional hydrogen bonds and hydrophobic interactions with HA residues. Furthermore, the blocking mechanism of HA differs for WT and D325A+R343V because of alternate glycan binding modes. The combined results suggest a mechanism through which the mode of SP-D-HA interaction could significantly influence viral aggregation and neutralization. These studies provide the first atomic-level molecular view of an innate host defense lectin inhibiting its viral glycoprotein target.

Review: Use of residual dipolar couplings to determine the structure of carbohydrates

Solution nuclear magnetic resonance spectroscopy is especially useful in the carbohydrate field. The measurement of residual dipolar couplings provides long-range structural information, a valuable complement for the structural study of carbohydrates either in its free form or in the bound state to proteins. They permit to deduce the geometry and the flexibility of the glycosidic linkages, which have a major influence on the conformation of carbohydrates and their overall shape. This article reviews the current application of the residual dipolar couplings methodology to carbohydrates.

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.

Three-dimensional structure of the weakly associated protein homodimer SeR13 using RDCs and paramagnetic surface mapping

The traditional NMR-based method for determining oligomeric protein structure usually involves distinguishing and assigning intra- and intersubunit NOEs. This task becomes challenging when determining symmetric homo-dimer structures because NOE cross-peaks from a given pair of protons occur at the same position whether intra- or intersubunit in origin. While there are isotope-filtering strategies for distinguishing intra from intermolecular NOE interactions in these cases, they are laborious and often prove ineffectual in cases of weak dimers, where observation of intermolecular NOEs is rare. Here, we present an efficient procedure for weak dimer structure determination based on residual dipolar couplings (RDCs), chemical shift changes upon dilution, and paramagnetic surface perturbations. This procedure is applied to the Northeast Structural Genomics Consortium protein target, SeR13, a negatively charged Staphylococcus epidermidis dimeric protein (K(d) 3.4 ± 1.4 mM) composed of 86 amino acids. A structure determination for the monomeric form using traditional NMR methods is presented, followed by a dimer structure determination using docking under orientation constraints from RDCs data, and scoring under residue pair potentials and shape-based predictions of RDCs. Validation using paramagnetic surface perturbation and chemical shift perturbation data acquired on sample dilution is also presented. The general utility of the dimer structure determination procedure and the possible relevance of SeR13 dimer formation are discussed.

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.

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.

Use of residual dipolar couplings in structural analysis of protein-ligand complexes by solution NMR spectroscopy

Investigation of structure-function relationships in protein complexes, specifically protein-ligand interactions, carry great significance in elucidating the structural and mechanistic bases of molecular recognition events and their role in regulating cell processes. Nuclear magnetic resonance (NMR) spectroscopy is one of the leading structural and analytical techniques in in-depth studies of protein-ligand interactions. Recent advances in NMR methodology such as transverse relaxation-optimized spectroscopy (TROSY) and residual dipolar couplings (RDCs) measured in liquid crystalline alignment medium, offer a viable alternative to traditional nuclear Overhauser enhancement (NOE)-based approaches for structure determination of large protein complexes. RDCs provide a way to constrain the relative orientation of two molecules in complex with each other by aligning their independently determined order tensors. The potential for utilization of RDCs can be extended to proteins with multiple ligands or even multimeric protein-ligand complexes, where symmetry properties of the protein can be taken advantage of. Availability of effective RDC data collection and analysis protocols can certainly aid this process by their incorporation into structure calculation protocols using intramolecular and intermolecular orientational restraints. This chapter discusses in detail some of these protocols including methods for sample preparation in liquid crystalline media, NMR experiments for RDC data collection, as well as software tools for RDC data analysis and protein-ligand complex structure determination.

Solution NMR structure of putidaredoxin-cytochrome P450cam complex via a combined residual dipolar coupling-spin labeling approach suggests a role for Trp106 of putidaredoxin in complex formation

The 58-kDa complex formed between the [2Fe-2S] ferredoxin, putidaredoxin (Pdx), and cytochrome P450cam (CYP101) from the bacterium Pseudomonas putida has been investigated by high-resolution solution NMR spectroscopy. Pdx serves as both the physiological reductant and effector for CYP101 in the enzymatic reaction involving conversion of substrate camphor to 5-exo-hydroxycamphor. In order to obtain an experimental structure for the oxidized Pdx-CYP101 complex, a combined approach using orientational data on the two proteins derived from residual dipolar couplings and distance restraints from site-specific spin labeling of Pdx has been applied. Spectral changes for residues in and near the paramagnetic metal cluster region of Pdx in complex with CYP101 have also been mapped for the first time using (15)N and (13)C NMR spectroscopy, leading to direct identification of the residues strongly affected by CYP101 binding. The new NMR structure of the Pdx-CYP101 complex agrees well with results from previous mutagenesis and biophysical studies involving residues at the binding interface such as formation of a salt bridge between Asp38 of Pdx and Arg112 of CYP101, while at the same time identifying key features different from those of earlier modeling studies. Analysis of the binding interface of the complex reveals that the side chain of Trp106, the C-terminal residue of Pdx and critical for binding to CYP101, is located across from the heme-binding loop of CYP101 and forms non-polar contacts with several residues in the vicinity of the heme group on CYP101, pointing to a potentially important role in complex formation.

Structure determination of a Galectin-3-carbohydrate complex using paramagnetism-based NMR constraints

The determination of the location and conformation of a natural ligand bound to a protein receptor is often a first step in the rational design of molecules that can modulate receptor function. NMR observables, including NOEs, often provide the basis for these determinations. However, when ligands are carbohydrates, interactions mediated by extensive hydrogen-bonding networks often reduce or eliminate NOEs between ligand and protein protons. In these cases, it is useful to look to other distance- and orientation-dependent observables that can constrain the geometry of ligand-protein complexes. Here we illustrate the use of paramagnetism-based NMR constraints, including pseudo-contact shifts (PCS) and field-induced residual dipolar couplings (RDCs). When a paramagnetic center can be attached to the protein, field-induced RDCs and PCS reflect only bound-state properties of the ligand, even when averages over small fractions of bound states and large fractions of free states are observed. The effects can also be observed over a long range, making it possible to attach a paramagnetic center to a remote part of the protein. The system studied here is a Galectin-3-lactose complex. A lanthanide-binding peptide showing minimal flexibility with respect to the protein was integrated into the C terminus of an expression construct for the Galectin-3-carbohydrate-binding domain. Dysprosium ion, which has a large magnetic susceptibility anisotropy, was complexed to the peptide, making it possible to observe both PCSs and field-induced RDCs for the protein and the ligand. The structure determined from these constraints shows agreement with a crystal structure of a Galectin-3-N-acetyllactosamine complex.

920 MHz ultra-high field NMR approaches to structural glycobiology

Although NMR spectroscopy has great potential to provide us with detailed structural information on oligosaccharides and glycoconjugates, the carbohydrate NMR analyses have been hampered by the severe spectral overlapping and the insufficiency of the conformational restraints. Recently, ultra-high field NMR spectrometers have become available for applications to structural analyses of biological macromolecules. Here we demonstrate that ultra-high fields offer not only increases in sensitivity and chemical shift dispersion but also potential benefits for providing unique information on chemical exchange and relaxation, by displaying NMR spectral data of oligosaccharide, glycoprotein, and glycolipid systems recorded at a 21.6 T magnetic field (corresponding to 920 MHz (1)H observation frequency). The ultra-high field NMR spectroscopy combined with sugar library and stable-isotope labeling approaches will open new horizons in structural glycobiology.

Bound-state residual dipolar couplings for rapidly exchanging ligands of His-tagged proteins

The study of bound-state conformations of ligands interacting with proteins is important to the understanding of protein function and the design of drugs that alter function. Traditionally, transferred nuclear Overhauser effects (trNOEs), measured from NMR spectra of ligands in rapid exchange between bound and free states, have been used in these studies, owing to the inherent heavy weighting of bound-state data in the averaged ligand signals. In principle, residual dipolar couplings (RDCs) provide a useful complement to NOE data in that they provide orientational constraints as opposed to distance constraints, but use in ligand-binding applications has been limited due to the absence of heavy weighting of bound-state data. A widely applicable approach to increasing the weighting of bound-state data in averaged RDCs measured on ligands is presented. The approach rests on association of a His-tagged protein with a nickel-chelate-carrying lipid inserted into the lipid bilayer-like alignment media used in the acquisition of RDCs. The approach is validated through the observation of bound-state RDCs for the disaccharide, lactose, bound to the carbohydrate recognition domain of the mammalian lectin, galectin-3.

Paramagnetic tagging of diamagnetic proteins for solution NMR

In this article, approaches towards the paramagnetic tagging of diamagnetic proteins are reviewed. Alignment can be achieved by adding paramagnetic fusion proteins or peptides to the C- or the N-terminus or by attaching paramagnetic tags to Cystein residues. Applications for the study of homodimer structures and protein/ligand interactions, as well as protein domain dynamics, are reviewed.

Enhancement of bound-state residual dipolar couplings: conformational analysis of lactose bound to Galectin-3

Residual dipolar couplings (RDCs) have proven to be a valuable NMR tool that can provide long-range constraints for molecular structure determination. The constraints are orientational in nature and are, thus, highly complementary to conventional distance constraints from NOE data. This complementarity would seem to extend to the study of the geometry of ligands bound to proteins. However, unlike transferred NOEs, where collection, even with a large excess of free ligand, results in measurements dominated by bound contributions, RDCs of exchanging ligands can be dominated by free-state contributions. Here we present a strategy for enhancement of RDCs from bound states that is based on specifically enhancing the alignment of the protein to which a ligand will bind. The protein is modified by addition of a hydrophobic alkyl tail that anchors it to the bicelles that are a part of the ordering medium needed for RDC measurement. As an illustration, we have added a propyl chain to the C terminus of the carbohydrate recognition domain of the protein, Galectin-3, and report enhanced RDCs that prove consistent with known bound-ligand geometries for this protein.

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.

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.

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.