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

Tandem "One-Shot" Measurement of Residual Chemical Shift Anisotropy and Residual Dipolar Coupling using Biphasic Supramolecular Peptide Liquid Crystals

Both solitary and tandem applications of residual chemical shift anisotropy (RCSA) and residual dipolar coupling (RDC) show great potential for the structural and configurational determination of organic molecules. A critical component of both RDC and RCSA methodologies is the alignment medium, whose availability is limited, especially for RCSA measurement. Moreover, reported RDC and RCSA acquisitions mainly rely on two experiments conducted under two different conditions, which are relatively time-consuming and easily cause experimental errors. Herein, a biphasic supramolecular lyotropic liquid crystalline (LLC) system was developed through the self-assembly of C21H43-CONH-V4K3-CONH2, which could act as an alignment medium for not only RDC but also RCSA extraction in DMSO-d6. Notably, the RCSA extraction was easily achieved via one-shot measurement from a single one-dimensional 13C NMR experiment, with no need for special instruments, devices, and correction. Relying on the biphasic LLC medium, meanwhile, RDC data were simply extracted from a single F1-coupled HSQC experiment, different from the standard protocol that requires two spectral acquisitions corresponding to the isotropic and anisotropic conditions. Collectively, the biphasic LLC medium is applicable for tandem RCSA and RDC measurements in one single sample, advancing the stereochemical elucidation of molecules of interest.

Probing Molecular Packing of Lipid Nanoparticles from 31P Solution and Solid-State NMR

Lipid nanoparticles (LNPs) are intricate multicomponent systems widely recognized for their efficient delivery of oligonucleotide cargo to host cells. Gaining insights into the molecular properties of LNPs is crucial for their effective design and characterization. However, analysis of their internal structure at the molecular level presents a significant challenge. This study introduces 31P nuclear magnetic resonance (NMR) methods to acquire structural and dynamic information about the phospholipid envelope of LNPs. Specifically, we demonstrate that the 31P chemical shift anisotropy (CSA) parameters serve as a sensitive indicator of the molecular assembly of distearoylphosphatidylcholine (DSPC) lipids within the particles. An analytical protocol for measuring 31P CSA is developed, which can be implemented using either solution NMR or solid-state NMR, offering wide accessibility and adaptability. The capability of this method is demonstrated using both model DSPC liposomes and real-world pharmaceutical LNP formulations. Furthermore, our method can be employed to investigate the impact of formulation processes and composition on the assembly of specifically LNP particles or, more generally, phospholipid-based delivery systems. This makes it an indispensable tool for evaluating critical pharmaceutical properties such as structural homogeneity, batch-to-batch reproducibility, and the stability of the particles.

Minimizing the risk of deducing wrong natural product structures from NMR data

There continues to be a disturbing number of natural products reported in the literature whose structures are incorrect. At least in part, this reflects the fact that many natural product chemists have limited formal nuclear magnetic resonance training. Gaps in training and lack of awareness regarding the challenges and ambiguities associated with two-dimensional nuclear magnetic resonance data interpretation can easily lead to errors in structure elucidation. The purpose of this tutorial is to point out some of these issues, highlight the kinds of errors that have been made and provide specific advice on how to avoid these missteps such that the risk of reporting a wrong structure is minimized.

Elucidating natural product structures using a robust measurement of carbon residual chemical shift anisotropy combined with DFT

Determination of configurations and conformations is an important step in the structural characterization of small molecules. Apart from utilizing isotropic J-couplings and nuclear overhauser effect (NOEs) measured in isotropic solution, anisotropic Nuclear Magnetic resonance (NMR) data such as residual dipolar couplings and residual chemical shift anisotropies (RCSAs) were also used to elucidate complex small molecule structures. Measuring RCSA has always been historically difficult due to the isotropic shift effect accompanied by molecular alignment and therefore only occasionally applied in a few examples. Here, we present a robust measurement of carbon RCSAs using a smaller gel-stretching device to determine the structures of a few small molecules. A systematic study on how different density functional theory computed anisotropies of the chemical shift anisotropy tensors impact RCSA data interpretation has also been discussed. We also discuss the effect of utilizing various carbons as reference nuclei for RCSA data extraction as well as the orientation behavior of estrone in orthogonal alignment media.

Relative configuration of micrograms of natural compounds using proton residual chemical shift anisotropy

3D molecular structure determination is a challenge for organic compounds or natural products available in minute amounts. Proton/proton and proton/carbon correlations yield the constitution. J couplings and NOEs oftentimes supported by one-bond ¹H,¹³C residual dipolar couplings (RDCs) or by ¹³C residual chemical shift anisotropies (RCSAs) provide the relative configuration. However, these RDCs or carbon RCSAs rely on 1% natural abundance of ¹³C preventing their use for compounds available only in quantities of a few 10’s of µgs. By contrast, ¹H RCSAs provide similar information on spatial orientation of structural moieties within a molecule, while using the abundant ¹H spin. Herein, ¹H RCSAs are accurately measured using constrained aligning gels or liquid crystals and applied to the 3D structural determination of molecules with varying complexities. Even more, deuterated alignment media allow the elucidation of the relative configuration of around 35 µg of a briarane compound isolated from Briareum asbestinum.

Determination of 3D molecular structures remains challenging for natural products or organic compounds available in minute amounts. Here, the authors determine the structure of complex molecules, including few micrograms of briarane B-3 isolated from Briareum asbestinums, through measurement of ¹H residual chemical shift anisotropy.

Structural basis of UCUU RNA motif recognition by splicing factor RBM20

The vertebrate splicing factor RBM20 (RNA binding motif protein 20) regulates protein isoforms important for heart development and function, with mutations in the gene linked to cardiomyopathy. Previous studies have identified the four nucleotide RNA motif UCUU as a common element in pre-mRNA targeted by RBM20. Here, we have determined the structure of the RNA Recognition Motif (RRM) domain from mouse RBM20 bound to RNA containing a UCUU sequence. The atomic details show that the RRM domain spans a larger region than initially proposed in order to interact with the complete UCUU motif, with a well-folded C-terminal helix encoded by exon 8 critical for high affinity binding. This helix only forms upon binding RNA with the final uracil, and removing the helix reduces affinity as well as specificity. We therefore find that RBM20 uses a coupled folding-binding mechanism by the C-terminal helix to specifically recognize the UCUU RNA motif.

An order-to-disorder structural switch activates the FoxM1 transcription factor

Intrinsically disordered transcription factor transactivation domains (TADs) function through structural plasticity, adopting ordered conformations when bound to transcriptional co-regulators. Many transcription factors contain a negative regulatory domain (NRD) that suppresses recruitment of transcriptional machinery through autoregulation of the TAD. We report the solution structure of an autoinhibited NRD-TAD complex within FoxM1, a critical activator of mitotic gene expression. We observe that while both the FoxM1 NRD and TAD are primarily intrinsically disordered domains, they associate and adopt a structured conformation. We identify how Plk1 and Cdk kinases cooperate to phosphorylate FoxM1, which releases the TAD into a disordered conformation that then associates with the TAZ2 or KIX domains of the transcriptional co-activator CBP. Our results support a mechanism of FoxM1 regulation in which the TAD undergoes switching between disordered and different ordered structures.

Complex macrocycle exploration: parallel, heuristic, and constraint-based conformer generation using ForceGen

ForceGen is a template-free, non-stochastic approach for 2D to 3D structure generation and conformational elaboration for small molecules, including both non-macrocycles and macrocycles. For conformational search of non-macrocycles, ForceGen is both faster and more accurate than the best of all tested methods on a very large, independently curated benchmark of 2859 PDB ligands. In this study, the primary results are on macrocycles, including results for 431 unique examples from four separate benchmarks. These include complex peptide and peptide-like cases that can form networks of internal hydrogen bonds. By making use of new physical movements (“flips” of near-linear sub-cycles and explicit formation of hydrogen bonds), ForceGen exhibited statistically significantly better performance for overall RMS deviation from experimental coordinates than all other approaches. The algorithmic approach offers natural parallelization across multiple computing-cores. On a modest multi-core workstation, for all but the most complex macrocycles, median wall-clock times were generally under a minute in fast search mode and under 2 min using thorough search. On the most complex cases (roughly cyclic decapeptides and larger) explicit exploration of likely hydrogen bonding networks yielded marked improvements, but with calculation times increasing to several minutes and in some cases to roughly an hour for fast search. In complex cases, utilization of NMR data to constrain conformational search produces accurate conformational ensembles representative of solution state macrocycle behavior. On macrocycles of typical complexity (up to 21 rotatable macrocyclic and exocyclic bonds), design-focused macrocycle optimization can be practically supported by computational chemistry at interactive time-scales, with conformational ensemble accuracy equaling what is seen with non-macrocyclic ligands. For more complex macrocycles, inclusion of sparse biophysical data is a helpful adjunct to computation.

Chemically cross-linked polyacrylonitrile. A DMSO compatible NMR alignment medium for measurement of residual dipolar couplings and residual chemical shift anisotropies

Chemically cross-linked polyacrylontrile polymer gels, have been prepared as an alignment medium compatible with DMSO-d₆. These gels allow measurement of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) with good accuracy as tested with brucine and α-santonin natural compounds. The gels can be reversibly compressed allowing easy measurement of RCSAs. They also present good physical homogeneity, clean HSQC spectra with little background ¹H signals, and allow unambiguous referencing of ¹³C spectra for RCSA extraction.

Synthesis and O-Glycosidic Linkage Conformational Analysis of 13C-Labeled Oligosaccharide Fragments of an Antifreeze Glycolipid

NMR studies of two ¹³C-labeled disaccharides and a tetrasaccharide were undertaken that comprise the backbone of a novel thermal hysteresis glycolipid containing a linear glycan sequence of alternating [βXyl p-(1→4)-βMan p-(1→4)] _n dimers. Experimental trans-glycoside NMR J-couplings, parameterized equations obtained from density functional theory (DFT) calculations, and an in-house circular statistics package ( MA'AT) were used to derive conformational models of linkage torsion angles ϕ and ψ in solution, which were compared to those obtained from molecular dynamics simulations. Modeling using different probability distribution functions showed that MA'AT models of ϕ in βMan(1→4)βXyl and βXyl(1→4)βMan linkages are very similar in the disaccharide building blocks, whereas MA'AT models of ψ differ. This pattern is conserved in the tetrasaccharide, showing that linkage context does not influence linkage geometry in this linear system. Good agreement was observed between the MA'AT and MD models of ψ with respect to mean values and circular standard deviations. Significant differences were observed for ϕ, indicating that revision of the force-field employed by GLYCAM is probably needed. Incorporation of the experimental models of ϕ and ψ into the backbone of an octasaccharide fragment leads to a helical amphipathic topography that may affect the thermal hysteresis properties of the glycolipid.

Use of Circular Statistics To Model αMan-(1→2)-αMan and αMan-(1→3)-α/βMan O-Glycosidic Linkage Conformation in 13C-Labeled Disaccharides and High-Mannose Oligosaccharides

A new experimental method, MA' AT analysis, has been applied to investigate the conformational properties of O-glycosidic linkages in several biologically important mannose-containing di- and oligosaccharides. Methyl α-d-mannopyranosyl-(1→2)-α-d-mannopyranoside (2), methyl α-d-mannopyranosyl-(1→3)-α-d-mannopyranoside (3), and methyl α-d-mannopyranosyl-(1→3)-β-d-mannopyranoside (4) were prepared with selective ¹³C-enrichment to enable the measurement of NMR scalar couplings across their internal O-glycosidic linkages. Density functional theory (DFT) was used to parameterize equations for J_CH and J_CC values in 2-4 that are sensitive to phi (ϕ) and psi (ψ). The experimental J-couplings and parameterized equations were treated using a circular statistics algorithm encoded in the MA' AT program. Conformations about ϕ and ψ treated using single-state von Mises models gave excellent fits to the ensembles of redundant J-couplings. Mean values and circular standard deviations (CSDs) for each linkage torsion angle ϕ (CSD) and ψ (CSD) in 2, -29° (25°) and 20° (22°); in 3, -36° (36°) and 8° (27°); in 4, -37° (34°) and 10° (26°); ϕ = H1'-C1'-O1'-CX and ψ = C1'-O1'-CX-HX (CX = aglycone carbon) were compared to histograms obtained from 1 μs aqueous molecular dynamics (MD) simulations and X-ray database statistical analysis. MA' AT-derived models of ψ were in very good agreement with the MD and X-ray data, but not those of ϕ, suggesting a need for force field revision. The effect of structural context on linkage conformation was also investigated in four selectively ¹³C-labeled homomannose tri- and tetrasaccharides using the MA' AT method. In the cases examined, context effects were found to be small.

The role of computer-assisted structure elucidation (CASE) programs in the structure elucidation of complex natural products

Computer-assisted structure elucidation can help to determine the structures of complex natural products while minimizing the risk of structure errors.

Insight into old and new pure shift nuclear magnetic resonance methods for enantiodiscrimination

Enantiodiscrimination and their quantification using nuclear magnetic resonance (NMR) spectroscopy has always been a subject of great interest. Proton is the nucleus of choice for enantiodiscrimination due to its high sensitivity and ubiquitous presence in nature. Despite its advantages, enantiodiscrimination suffers from extensive signal splitting by the proton-proton scalar couplings, which give complex multiplets that spread over a frequency range of some tens of hertz. These multiplets often overlap, further complicating interpretation of the spectra and quantifications. In the present review, we discuss some of the recent developments in the pure shift ¹ H NMR based methods for enantiomer resolution and enantiodiscrimination. We also compare various pure shift methods used for enantiodiscrimination and measurement of enantiomeric excess, considering the fact that conventional ¹ H NMR fails to provide any detailed insight.

Measurement of residual chemical shift anisotropies in compressed polymethylmethacrylate gels. Automatic compensation of gel isotropic shift contribution

Mechanical compression of polymer gels provides a simple way for the measurement of residual chemical shift anisotropies, which then can be employed, on its own, or in combination with residual dipolar couplings, for structural elucidation purposes. Residual chemical shift anisotropies measured using compression devices needed a posteriori correction to account for the increase of the polymer to solvent ratio inside the swollen gel. This correction has been cast before in terms of a single-free parameter which, as shown here, can be simultaneously optimized along with the components of the alignment tensor while still retaining discriminating power of the different relative configurations as illustrated in the stereochemical analysis of α-santonin and 10-epi-8-deoxycumambrin B.

Enhanced measurement of residual chemical shift anisotropy for small molecule structure elucidation

A method is introduced to measure residual chemical shift anisotropies conveniently and accurately in the mesophase of poly-γ-(benzyl-l-glutamate). The alignment amplitude is substantially enhanced over common methods which greatly benefits measurements particularly on sp3 carbons. The approach offers significant improvements in data accuracy and utility for small molecule structure determination.

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.

Conformational Populations of β-(1→4) O-Glycosidic Linkages Using Redundant NMR J-Couplings and Circular Statistics

Twelve disaccharides containing β-(1→4) linkages and displaying systematic structural variations in the vicinity of these linkages were selectively labeled with ¹³C to facilitate measurements of multiple NMR spin-spin (scalar; J) coupling constants (J_CH and J_CC values) across their O-glycosidic linkages. Ensembles of spin-couplings (²J_COC, ³J_COCH, ³J_COCC) sensitive to the two linkage torsion angles, phi (ϕ) and psi (ψ), were analyzed by using parametrized equations obtained from density functional theory (DFT) calculations, Fredholm theory, and circular statistics to calculate experiment-based rotamer populations for ϕ and ψ in each disaccharide. With the statistical program MA'AT, torsion angles ϕ and ψ were modeled as a single von Mises distribution, which yielded two parameters, the mean position and the circular standard deviation (CSD) for each angle. The NMR-derived rotamer populations were compared to those obtained from 1 μs aqueous molecular dynamics (MD) simulations and crystallographic database statistical analyses. Conformer populations obtained exclusively from the MA'AT treatment of redundant J-couplings were in very good agreement with those obtained from the MD simulations, providing evidence that conformational populations can be determined by NMR for mobile molecular elements such as O-glycosidic linkages with minimal input from theory. The approach also provides an experimental means to validate the conformational preferences predicted from MD simulations. The conformational behaviors of ϕ in the 12 disaccharides were very similar, but those of ψ varied significantly, allowing a classification of the 12 disaccharides based on preferred linkage conformation in solution.

Nuclear Magnetic Resonance Insight into the Multiple Glycosaminoglycan Binding Modes of the Link Module from Human TSG-6

Tumor necrosis factor-stimulated gene-6 (TSG-6) is a hyaluronan (HA)-binding protein that is essential for stabilizing and remodeling the extracellular matrix (ECM) during ovulation and inflammatory disease processes such as arthritis. The Link module, one of the domains of TSG-6, is responsible for binding hyaluronan and other glycosaminoglycans found in the ECM. In this study, we used a well-defined chondroitin sulfate (CS) hexasaccharide (ΔC444S) to determine the structure of the Link module, in solution, in its chondroitin sulfate-bound state. A variety of nuclear magnetic resonance techniques were employed, including chemical shift perturbation, residual dipolar couplings (RDCs), nuclear Overhauser effects, spin relaxation measurements, and paramagnetic relaxation enhancements from a spin-labeled analogue of ΔC444S. The binding site for ΔC444S on the Link module overlapped with that of HA. Surprisingly, ΔC444S binding induced dimerization of the Link module (as confirmed by analytical ultracentrifugation), and a second weak binding site that partially overlapped with a previously identified heparin site was detected. A dimer model was generated using chemical shift perturbations and RDCs as restraints in the docking program HADDOCK. We postulate that the molecular cross-linking enhanced by the multiple binding modes of the Link module might be critical for remodeling the ECM during inflammation/ovulation and might contribute to other functions of TSG-6.

Ensembles of a small number of conformations with relative populations

In our previous work, we proposed a new way to represent protein native states, using ensembles of a small number of conformations with relative Populations, or ESP in short. Using Ubiquitin as an example, we showed that using a small number of conformations could greatly reduce the potential of overfitting and assigning relative populations to protein ensembles could significantly improve their quality. To demonstrate that ESP indeed is an excellent alternative to represent protein native states, in this work we compare the quality of two ESP ensembles of Ubiquitin with several well-known regular ensembles or average structure representations. Extensive amount of significant experimental data are employed to achieve a thorough assessment. Our results demonstrate that ESP ensembles, though much smaller in size comparing to regular ensembles, perform equally or even better sometimes in all four different types of experimental data used in the assessment, namely, the residual dipolar couplings, residual chemical shift anisotropy, hydrogen exchange rates, and solution scattering profiles. This work further underlines the significance of having relative populations in describing the native states.

Spatial reorientation experiments for NMR of solids and partially oriented liquids

Motional reorientation experiments are extensions of Magic Angle Spinning (MAS) where the rotor axis is changed in order to average out, reintroduce, or scale anisotropic interactions (e.g. dipolar couplings, quadrupolar interactions or chemical shift anisotropies). This review focuses on Variable Angle Spinning (VAS), Switched Angle Spinning (SAS), and Dynamic Angle Spinning (DAS), all of which involve spinning at two or more different angles sequentially, either in successive experiments or during a multidimensional experiment. In all of these experiments, anisotropic terms in the Hamiltonian are scaled by changing the orientation of the spinning sample relative to the static magnetic field. These experiments vary in experimental complexity and instrumentation requirements. In VAS, many one-dimensional spectra are collected as a function of spinning angle. In SAS, dipolar couplings and/or chemical shift anisotropies are reintroduced by switching the sample between two different angles, often 0° or 90° and the magic angle, yielding a two-dimensional isotropic-anisotropic correlation spectrum. Dynamic Angle Spinning (DAS) is a related experiment that is used to simultaneously average out the first- and second-order quadrupolar interactions, which cannot be accomplished by spinning at any unique rotor angle in physical space. Although motional reorientation experiments generally require specialized instrumentation and data analysis schemes, some are accessible with only minor modification of standard MAS probes. In this review, the mechanics of each type of experiment are described, with representative examples. Current and historical probe and coil designs are discussed from the standpoint of how each one accomplishes the particular objectives of the experiment(s) it was designed to perform. Finally, applications to inorganic materials and liquid crystals, which present very different experimental challenges, are discussed. The review concludes with perspectives on how motional reorientation experiments can be applied to current problems in chemistry, molecular biology, and materials science, given the many advances in high-field NMR magnets, fast spinning, and sample preparation realized in recent years.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Structural characterization of a flexible two-domain protein in solution using small angle X-ray scattering and NMR data

Multidomain proteins in which individual domains are connected by linkers often possess inherent interdomain flexibility that significantly complicates their structural characterization in solution using either nuclear magnetic resonance (NMR) spectroscopy or small-angle X-ray scattering (SAXS) alone. Here, we report a protocol for joint refinement of flexible multidomain protein structures against NMR distance and angular restraints, residual dipolar couplings, and SAXS data. The protocol is based on the ensemble optimization method principle (Bernadó et al., 2007) and is compared with different refinement strategies for the structural characterization of the flexible two-domain protein sf3636 from Shigella flexneri 2a. The results of our refinement suggest the existence of a dominant population of configurational states in solution possessing an overall elongated shape and restricted relative twisting of the two domains.

Improvements to REDCRAFT: a software tool for simultaneous characterization of protein backbone structure and dynamics from residual dipolar couplings

Within the past two decades, there has been an increase in the acquisition of residual dipolar couplings (RDC) for investigations of biomolecular structures. Their use however is still not as widely adopted as the traditional methods of structure determination by NMR, despite their potential for extending the limits in studies that examine both the structure and dynamics of biomolecules. This is in part due to the difficulties associated with the analysis of this information-rich data type. The software analysis tool REDCRAFT was previously introduced to address some of these challenges. Here we describe and evaluate a number of additional features that have been incorporated in order to extend its computational and analytical capabilities. REDCRAFT's more traditional enhancements integrate a modified steric collision term, as well as structural refinement in the rotamer space. Other, non-traditional improvements include: the filtering of viable structures based on relative order tensor estimates, decimation of the conformational space based on structural similarity, and forward/reverse folding of proteins. Utilizing REDCRAFT's newest features we demonstrate de-novo folding of proteins 1D3Z and 1P7E to within less than 1.6 Å of the corresponding X-ray structures, using as many as four RDCs per residue and as little as two RDCs per residue, in two alignment media. We also show the successful folding of a structure to less than 1.6 Å of the X-ray structure using {C(i-1)-N(i), N(i)-H(i), and C(i-1)-H(i)} RDCs in one alignment medium, and only {N(i)-H(i)} in the second alignment medium (a set of data which can be collected on deuterated samples). The program is available for download from our website at http://ifestos.cse.sc.edu .

Advances in the REDCAT software package

Background

Residual Dipolar Couplings (RDCs) have emerged in the past two decades as an informative source of experimental restraints for the study of structure and dynamics of biological macromolecules and complexes. The REDCAT software package was previously introduced for the analysis of molecular structures using RDC data. Here we report additional features that have been included in this software package in order to expand the scope of its analyses. We first discuss the features that enhance REDCATs user-friendly nature, such as the integration of a number of analyses into one single operation and enabling convenient examination of a structural ensemble in order to identify the most suitable structure. We then describe the new features which expand the scope of RDC analyses, performing exercises that utilize both synthetic and experimental data to illustrate and evaluate different features with regard to structure refinement and structure validation.

Results

We establish the seamless interaction that takes place between REDCAT, VMD, and Xplor-NIH in demonstrations that utilize our newly developed REDCAT-VMD and XplorGUI interfaces. These modules enable visualization of RDC analysis results on the molecular structure displayed in VMD and refinement of structures with Xplor-NIH, respectively. We also highlight REDCAT’s Error-Analysis feature in reporting the localized fitness of a structure to RDC data, which provides a more effective means of recognizing local structural anomalies. This allows for structurally sound regions of a molecule to be identified, and for any refinement efforts to be focused solely on locally distorted regions.

Conclusions

The newly engineered REDCAT software package, which is available for download via the WWW from http://ifestos.cse.sc.edu, has been developed in the Object Oriented C++ environment. Our most recent enhancements to REDCAT serve to provide a more complete RDC analysis suite, while also accommodating a more user-friendly experience, and will be of great interest to the community of researchers and developers since it hides the complications of software development.

Chemokine oligomerization in cell signaling and migration

Chemokines are small proteins best known for their role in controlling the migration of diverse cells, particularly leukocytes. Upon binding to their G-protein-coupled receptors on the leukocytes, chemokines stimulate the signaling events that cause cytoskeletal rearrangements involved in cell movement, and migration of the cells along chemokine gradients. Depending on the cell type, chemokines also induce many other types of cellular responses including those related to defense mechanisms, cell proliferation, survival, and development. Historically, most research efforts have focused on the interaction of chemokines with their receptors, where monomeric forms of the ligands are the functionally relevant state. More recently, however, the importance of chemokine interactions with cell surface glycosaminoglycans has come to light, and in most cases appears to involve oligomeric chemokine structures. This review summarizes existing knowledge relating to the structure and function of chemokine oligomers, and emerging methodology for determining structures of complex chemokine assemblies in the future.

NMR detection and characterization of sialylated glycoproteins and cell surface polysaccharides

Few solution NMR pulse sequences exist that are explicitly designed to characterize carbohydrates (glycans). This is despite the essential role carbohydrate motifs play in cell-cell communication, microbial pathogenesis, autoimmune disease progression and cancer metastasis, and despite that fact that glycans, often shed to extra-cellular fluids, can be diagnostic of disease. Here we present a suite of two dimensional coherence experiments to measure three different correlations (H3-C2, H3-C1, and C1-C2) on sialic acids, a group of nine-carbon carbohydrates found on eukaryotic cell surfaces that often play a key role in disease processes. The chemical shifts of the H3, C2, and C1 nuclei of sialic acids are sensitive to carbohydrate linkage, linkage conformation, and ionization state of the C1 carboxylate. The experiments reported include rigorous filter elements to enable detection and characterization of isotopically labeled sialic acids with high sensitivity in living cells and crude isolates with minimal interference from unwanted signals arising from the ~1% (13)C-natural abundance of cellular metabolites. Application is illustrated with detection of sialic acids on living cells, in unpurified mixtures, and at the terminus of the N-glycan on the 55 kDa immunoglobulin G Fc.

Oligomeric structure of the chemokine CCL5/RANTES from NMR, MS, and SAXS data

CCL5 (RANTES) is a proinflammatory chemokine known to activate leukocytes through its receptor, CCR5. Although the monomeric form of CCL5 is sufficient to cause cell migration in vitro, CCL5's propensity for aggregation is essential for migration in vivo, T cell activation and apoptosis, and HIV entry into cells. However, there is currently no structural information on CCL5 oligomers larger than the canonical CC chemokine dimer. In this study the solution structure of a CCL5 oligomer was investigated using an integrated approach, including NMR residual dipolar couplings to determine allowed relative orientations of the component monomers, SAXS to restrict overall shape, and hydroxyl radical footprinting and NMR cross-saturation experiments to identify interface residues. The resulting model of the CCL5 oligomer provides a basis for explaining the disaggregating effect of E66 and E26 mutations and suggests mechanisms by which glycosaminoglycan binding may promote oligomer formation and facilitate cell migration in vivo.

Solution structure of 4'-phosphopantetheine - GmACP3 from Geobacter metallireducens: a specialized acyl carrier protein with atypical structural features and a putative role in lipopolysaccharide biosynthesis

GmACP3 from Geobacter metallireducens is a specialized acyl carrier protein (ACP) whose gene, gmet_2339, is located near genes encoding many proteins involved in lipopolysaccharide (LPS) biosynthesis, indicating a likely function for GmACP3 in LPS production. By overexpression in Escherichia coli, about 50% holo-GmACP3 and 50% apo-GmACP3 were obtained. Apo-GmACP3 exhibited slow precipitation and non-monomeric behavior by (15)N NMR relaxation measurements. Addition of 4'-phosphopantetheine (4'-PP) via enzymatic conversion by E. coli holo-ACP synthase resulted in stable >95% holo-GmACP3 that was characterized as monomeric by (15)N relaxation measurements and had no indication of conformational exchange. We have determined a high-resolution solution structure of holo-GmACP3 by standard NMR methods, including refinement with two sets of NH residual dipolar couplings, allowing for a detailed structural analysis of the interactions between 4'-PP and GmACP3. Whereas the overall four helix bundle topology is similar to previously solved ACP structures, this structure has unique characteristics, including an ordered 4'-PP conformation that places the thiol at the entrance to a central hydrophobic cavity near a conserved hydrogen-bonded Trp-His pair. These residues are part of a conserved WDSLxH/N motif found in GmACP3 and its orthologs. The helix locations and the large hydrophobic cavity are more similar to medium- and long-chain acyl-ACPs than to other apo- and holo-ACP structures. Taken together, structural characterization along with bioinformatic analysis of nearby genes suggests that GmACP3 is involved in lipid A acylation, possibly by atypical long-chain hydroxy fatty acids, and potentially is involved in synthesis of secondary metabolites.

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.

Contemporary methods in structure determination of membrane proteins by solution NMR

Integral membrane proteins are vital to life, being responsible for information and material exchange between a cell and its environment. Although high-resolution structural information is needed to understand how these functions are achieved, membrane proteins remain an under-represented subset of the protein structure databank. Solution NMR is increasingly demonstrating its ability to help address this knowledge shortfall, with the development of a diverse array of techniques to counter the challenges presented by membrane proteins. Here we document the advances that are helping to define solution NMR as an effective tool for membrane protein structure determination. Developments introduced over the last decade in the production of isotope-labeled samples, reconstitution of these samples into the growing selection of NMR-compatible membrane-mimetic systems, and the approaches used for the acquisition and application of structural restraints from these complexes are reviewed.

The use of residual dipolar coupling in studying proteins by NMR

The development of residual dipolar coupling (RDC) in protein NMR spectroscopy, over a decade ago, has become a useful and almost routine tool for accurate protein solution structure determination. RDCs provide orientation information of magnetic dipole-dipole interaction vectors within a common reference frame. Its measurement requires a nonisotropic orientation, through a direct or indirect magnetic field alignment, of the protein in solution. There has been recent progress in the developments of alignment methods to allow the measurement of RDC and of methods to analyze the resulting data. In this chapter we briefly go through the mathematical expressions for the RDC and common descriptions of the alignment tensor, which may be represented using either Saupe order or the principal order matrix. Then we review the latest developments in alignment media. In particular we looked at the lipid-compatible media that allow the measurement of RDCs for membrane proteins. Other methods including conservative surface residue mutation have been invented to obtain up to five orthogonal alignment tensors that provide a potential for de novo structure and dynamics study using RDCs exclusively. We then discuss approximations assumed in RDC interpretations and different views on dynamics uncovered from the RDC method. In addition to routine usage of RDCs in refining a single structure, novel applications such as ensemble refinement against RDCs have been implemented to represent protein structure and dynamics in solution. The RDC application also extends to the study of protein-substrate interaction as well as to solving quaternary structure of oligomer in equilibrium with a monomer, opening an avenue for RDCs in high-order protein structure determination.