Partial alignment of biomolecules: an aid to NMR characterization

Concurrent Identification and Characterization of Protein Structure and Continuous Internal Dynamics with REDCRAFT

Internal dynamics of proteins can play a critical role in the biological function of some proteins. Several well documented instances have been reported such as MBP, DHFR, hTS, DGCR8, and NSP1 of the SARS-CoV family of viruses. Despite the importance of internal dynamics of proteins, there currently are very few approaches that allow for meaningful separation of internal dynamics from structural aspects using experimental data. Here we present a computational approach named REDCRAFT that allows for concurrent characterization of protein structure and dynamics. Here, we have subjected DHFR (PDB-ID 1RX2), a 159-residue protein, to a fictitious, mixed mode model of internal dynamics. In this simulation, DHFR was segmented into 7 regions where 4 of the fragments were fixed with respect to each other, two regions underwent rigid-body dynamics, and one region experienced uncorrelated and melting event. The two dynamical and rigid-body segments experienced an average orientational modification of 7° and 12° respectively. Observable RDC data for backbone C′-N, N-HN, and C′-HN were generated from 102 uniformly sampled frames that described the molecular trajectory. The structure calculation of DHFR with REDCRAFT by using traditional Ramachandran restraint produced a structure with 29 Å of structural difference measured over the backbone atoms (bb-rmsd) over the entire length of the protein and an average bb-rmsd of more than 4.7 Å over each of the dynamical fragments. The same exercise repeated with context-specific dihedral restraints generated by PDBMine produced a structure with bb-rmsd of 21 Å over the entire length of the protein but with bb-rmsd of less than 3 Å over each of the fragments. Finally, utilization of the Dynamic Profile generated by REDCRAFT allowed for the identification of different dynamical regions of the protein and the recovery of individual fragments with bb-rmsd of less than 1 Å. Following the recovery of the fragments, our assembly procedure of domains (larger segments consisting of multiple fragments with a common dynamical profile) correctly assembled the four fragments that are rigid with respect to each other, categorized the two domains that underwent rigid-body dynamics, and identified one dynamical region for which no conserved structure could be defined. In conclusion, our approach was successful in identifying the dynamical domains, recovery of structure where it is meaningful, and relative assembly of the domains when possible.

REDCRAFT: A computational platform using residual dipolar coupling NMR data for determining structures of perdeuterated proteins in solution

Nuclear Magnetic Resonance (NMR) spectroscopy is one of the three primary experimental means of characterizing macromolecular structures, including protein structures. Structure determination by solution NMR spectroscopy has traditionally relied heavily on distance restraints derived from nuclear Overhauser effect (NOE) measurements. While structure determination of proteins from NOE-based restraints is well understood and broadly used, structure determination from Residual Dipolar Couplings (RDCs) is relatively less well developed. Here, we describe the new features of the protein structure modeling program REDCRAFT and focus on the new Adaptive Decimation (AD) feature. The AD plays a critical role in improving the robustness of REDCRAFT to missing or noisy data, while allowing structure determination of larger proteins from less data. In this report we demonstrate the successful application of REDCRAFT in structure determination of proteins ranging in size from 50 to 145 residues using experimentally collected data, and of larger proteins (145 to 573 residues) using simulated RDC data. In both cases, REDCRAFT uses only RDC data that can be collected from perdeuterated proteins. Finally, we compare the accuracy of structure determination from RDCs alone with traditional NOE-based methods for the structurally novel PF.2048.1 protein. The RDC-based structure of PF.2048.1 exhibited 1.0 Å BB-RMSD with respect to a high-quality NOE-based structure. Although optimal strategies would include using RDC data together with chemical shift, NOE, and other NMR data, these studies provide proof-of-principle for robust structure determination of largely-perdeuterated proteins from RDC data alone using REDCRAFT.

Residual Dipolar Couplings have the potential to improve the accuracy and reduce the time needed to characterize protein structures. In addition, RDC data have been demonstrated to concurrently elucidate structure of proteins, provide assignment of resonances, and characterize the internal dynamics of proteins. Given all the advantages associated with the study of proteins from RDC data, based on the statistics provided by the Protein Databank (PDB), surprisingly only 124 proteins (out of nearly 150,000 proteins) have utilized RDCs as part of their structure determination. Even a smaller subset of these proteins (approximately 7) have utilized RDCs as the primary source of data for structure determination. One key factor in the use of RDCs is the challenging computational and analytical aspects of this source of data. In this report, we demonstrate the success of the REDCRAFT software package in structure determination of proteins using RDC data that can be collected from small and large proteins in a routine fashion. REDCRAFT accomplishes the challenging task of structure determination from RDCs by introducing a unique search and optimization technique that is both robust and computationally tractable. Structure determination from routinely collectable RDC data using REDCRAFT can complement existing methods to provide faster and more accurate studies of larger and more complex protein structures by NMR spectroscopy in solution state.

Magnetically Induced Alignment of Natural Products for Stereochemical Structure Determination via NMR

Anisotropic NMR has gained increasing popularity to determine the structure and specifically the configuration of small, flexible, non‐crystallizable molecules. However, it suffers from the necessity to dissolve the analyte in special media such as liquid crystals or polymer gels. Generally, small degrees of alignment are also caused by an anisotropic magnetic susceptibility of the molecule, for example, induced by aromatic moieties. For this mechanism, the alignment can be predicted via density functional theory. Here we show that both residual dipolar couplings and residual chemical shift anisotropies can be acquired from natural products without special sample preparation using magnetically induced alignment. On the two examples of the novel natural product gymnochrome G and the alkaloid strychnine, these data, together with the predicted alignment, yield the correct configuration with high certainty.

Structure of Monomeric Transthyretin Carrying the Clinically Important T119M Mutation

Mutations in the protein transthyretin can cause as well as protect individuals from transthyretin amyloidosis, an incurable fatal inherited disease. Little is known, however, about the structural basis of pathogenic and clinically protective transthyretin mutants. Here we determined the solution structure of a transthyretin monomer that carries the clinically important T119M mutation. The structure displays a non-native arrangement that is distinct from all known structures of transthyretin and highlights the importance of high-resolution studies in solution for understanding molecular processes that lead to amyloid diseases.

Hybrid Approaches to Structural Characterization of Conformational Ensembles of Complex Macromolecular Systems Combining NMR Residual Dipolar Couplings and Solution X-ray Scattering

Solving structures or structural ensembles of large macromolecular systems in solution poses a challenging problem. While NMR provides structural information at atomic resolution, increased spectral complexity, chemical shift overlap, and short transverse relaxation times (associated with slow tumbling) render application of the usual techniques that have been so successful for medium sized systems (<50 kDa) difficult. Solution X-ray scattering, on the other hand, is not limited by molecular weight but only provides low resolution structural information related to the overall shape and size of the system under investigation. Here we review how combining atomic resolution structures of smaller domains with sparse experimental data afforded by NMR residual dipolar couplings (which yield both orientational and shape information) and solution X-ray scattering data in rigid-body simulated annealing calculations provides a powerful approach for investigating the structural aspects of conformational dynamics in large multidomain proteins. The application of this hybrid methodology is illustrated for the 128 kDa dimer of bacterial Enzyme I which exists in a variety of open and closed states that are sampled at various points in the catalytic cycles, and for the capsid protein of the human immunodeficiency virus.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

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.

REDOR NMR for drug discovery

Rotational-echo double-resonance (REDOR) NMR is a powerful and versatile solid-state NMR measurement that has been recruited to elucidate drug modes of action and to drive the design of new therapeutics. REDOR has been implemented to examine composition, structure, and dynamics in diverse macromolecular and whole-cell systems, including taxol-bound microtubules, enzyme-cofactor-inhibitor ternary complexes, and antibiotic-whole-cell complexes. The REDOR approach involves the integrated design of specific isotopic labeling strategies and the selection of appropriate REDOR experiments. By way of example, this digest illustrates the versatility of the REDOR approach, with an emphasis on the practical considerations of experimental design and data interpretation.

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.

Protein structure validation and identification from unassigned residual dipolar coupling data using 2D-PDPA

More than 90% of protein structures submitted to the PDB each year are homologous to some previously characterized protein structure. The extensive resources that are required for structural characterization of proteins can be justified for the 10% of the novel structures, but not for the remaining 90%. This report presents the 2D-PDPA method, which utilizes unassigned residual dipolar coupling in order to address the economics of structure determination of routine proteins by reducing the data acquisition and processing time. 2D-PDPA has been demonstrated to successfully identify the correct structure of an array of proteins that range from 46 to 445 residues in size from a library of 619 decoy structures by using unassigned simulated RDC data. When using experimental data, 2D-PDPA successfully identified the correct NMR structures from the same library of decoy structures. In addition, the most homologous X-ray structure was also identified as the second best structural candidate. Finally, success of 2D-PDPA in identifying and evaluating the most appropriate structure from a set of computationally predicted structures in the case of a previously uncharacterized protein Pf2048.1 has been demonstrated. This protein exhibits less than 20% sequence identity to any protein with known structure and therefore presents a compelling and practical application of our proposed work.

Partial alignment and measurement of residual dipolar couplings of proteins under high hydrostatic pressure

High-pressure NMR spectroscopy has emerged as a complementary approach for investigating various structural and thermodynamic properties of macromolecules. Noticeably absent from the array of experimental restraints that have been employed to characterize protein structures at high hydrostatic pressure is the residual dipolar coupling, which requires the partial alignment of the macromolecule of interest. Here we examine five alignment media that are commonly used at ambient pressure for this purpose. We find that the spontaneous alignment of Pf1 phage, d(GpG) and a C12E5/n-hexnanol mixture in a magnetic field is preserved under high hydrostatic pressure. However, DMPC/DHPC bicelles and collagen gel are found to be unsuitable. Evidence is presented to demonstrate that pressure-induced structural changes can be identified using the residual dipolar coupling.

Very large residual dipolar couplings from deuterated ubiquitin

Main-chain (1)H(N)-(15)N residual dipolar couplings (RDCs) ranging from approximately -200 to 200 Hz have been measured for ubiquitin under strong alignment conditions in Pf1 phage. This represents a ten-fold increase in the degree of alignment over the typical weakly aligned samples. The measurements are made possible by extensive proton-dilution of the sample, achieved by deuteration of the protein with partial back-substitution of labile protons from 25 % H(2)O / 75 % D(2)O buffer. The spectral quality is further improved by application of deuterium decoupling. Since standard experiments using fixed-delay INEPT elements cannot accommodate a broad range of couplings, the measurements were conducted using J-resolved and J-modulated versions of the HSQC and TROSY sequences. Due to unusually large variations in dipolar couplings, the trosy (sharp) and anti-trosy (broad) signals are often found to be interchanged in the TROSY spectra. To distinguish between the two, we have relied on their respective (15)N linewidths. This strategy ultimately allowed us to determine the signs of RDCs. The fitting of the measured RDC values to the crystallographic coordinates of ubiquitin yields the quality factor Q = 0.16, which confirms the perturbation-free character of the Pf1 alignment. Our results demonstrate that RDC data can be successfully acquired not only in dilute liquid crystals, but also in more concentrated ones. As a general rule, the increase in liquid crystal concentration improves the stability of alignment media and makes them more tolerant to variations in sample conditions. The technical ability to measure RDCs under moderately strong alignment conditions may open the door for development of alternative alignment media, including new types of media that mimic biologically relevant systems.

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.

A solution model of the complex formed by adrenodoxin and adrenodoxin reductase determined by paramagnetic NMR spectroscopy

Lanthanide tags offer the opportunity to retrieve long-range distance information from NMR experiments that can be used to guide protein docking. To determine whether sufficient restraints can be retrieved for proteins with low solubility and availability, Ln tags were applied in the study of the 65 kDa membrane-associated protein complex formed by the electron carrier adrenodoxin and its electron donor, adrenodoxin reductase. The reductase is only monomeric at low concentration, and the paramagnetic iron-sulfur cluster of adrenodoxin broadens many of the resonances of nuclei in the interface. Guided by the paramagnetic restraints obtained using two Ln-tag attachment sites, protein docking yields a cluster of solutions with an rmsd of 3.2 A. The mean structure is close to the crystal structure of the cross-linked complex, with an rmsd of 4.0 A. It is concluded that with the application of Ln tags paramagnetic NMR restraints for structure determination can be retrieved even for difficult, low-concentration protein complexes.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

The conformational properties of the Glc3Man unit suggest conformational biasing within the chaperone-assisted glycoprotein folding pathway

A major puzzle is: are all glycoproteins routed through the ER calnexin pathway irrespective of whether this is required for their correct folding? Calnexin recognizes the terminal Glcalpha1-3Manalpha linkage, formed by trimming of the Glcalpha1-2Glcalpha1-3Glcalpha1-3Manalpha (Glc3Man) unit in Glc3Man9GlcNAc2. Different conformations of this unit have been reported. We have addressed this problem by studying the conformation of a series of N-glycans; i.e. Glc3ManOMe, Glc3Man(4,5,7)GlcNAc2 and Glc1Man9GlcNAc2 using 2D NMR NOESY, ROESY, T-ROESY and residual dipolar coupling experiments in a range of solvents, along with solution molecular dynamics simulations of Glc3ManOMe. Our results show a single conformation for the Glcalpha1-2Glcalpha and Glcalpha1-3Glcalpha linkages, and a major (65%) and a minor (30%) conformer for the Glcalpha1-3Manalpha linkage. Modeling of the binding of Glc1Man9GlcNAc2 to calnexin suggests that it is the minor conformer that is recognized by calnexin. This may be one of the mechanisms for controlling the rate of recruitment of proteins into the calnexin/calreticulin chaperone system and enabling proteins that do not require such assistance for folding to bypass the system. This is the first time evidence has been presented on glycoprotein folding that suggests the process may be optimized to balance the chaperone-assisted and chaperone-independent pathways.

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.

Estimation of relative order tensors, and reconstruction of vectors in space using unassigned RDC data and its application

Advances in NMR instrumentation and pulse sequence design have resulted in easier acquisition of Residual Dipolar Coupling (RDC) data. However, computational and theoretical analysis of this type of data has continued to challenge the international community of investigators because of their complexity and rich information content. Contemporary use of RDC data has required a-priori assignment, which significantly increases the overall cost of structural analysis. This article introduces a novel algorithm that utilizes unassigned RDC data acquired from multiple alignment media (nD-RDC, n3) for simultaneous extraction of the relative order tensor matrices and reconstruction of the interacting vectors in space. Estimation of the relative order tensors and reconstruction of the interacting vectors can be invaluable in a number of endeavors. An example application has been presented where the reconstructed vectors have been used to quantify the fitness of a template protein structure to the unknown protein structure. This work has other important direct applications such as verification of the novelty of an unknown protein and validation of the accuracy of an available protein structure model in drug design. More importantly, the presented work has the potential to bridge the gap between experimental and computational methods of structure determination.

Rapid classification of protein structure models using unassigned backbone RDCs and probability density profile analysis (PDPA)

A method of identifying the best structural model for a protein of unknown structure from a list of structural candidates using unassigned 15N1H residual dipolar coupling (RDC) data and probability density profile analysis (PDPA) is described. Ten candidate structures have been obtained for the structural genomics target protein PF2048.1 using ROBETTA. 15N1H residual dipolar couplings have been measured from NMR spectra of the protein in two alignment media and these data have been analyzed using PDPA to rank the models in terms of their ability to represent the actual structure. A number of advantages in using this method to characterize a protein structure become apparent. RDCs can easily and rapidly be acquired, and without the need for assignment, the cost and duration of data acquisition is greatly reduced. The approach is quite robust with respect to imprecise and missing data. In the case of PF2048.1, a 79 residue protein, only 58 and 55 of the total RDC data were observed. The method can accelerate structure determination at higher resolution using traditional NMR spectroscopy by providing a starting point for the addition of NOEs and other NMR structural data.

NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Assessment of solvent effects: do weak alignment media affect the structure of the solute?

Alignment media used for measuring residual dipolar couplings, such as solutions of filamentous phages, phospholipid mixtures, polyacrylamide gels and various lyotropic liquid crystalline systems were investigated with respect to solvent effects on molecular structure. Structural parameters of the small rigid model compound 13C-acetonitrile were calculated from dipolar couplings and variations from expectation values were used for assessment of solvent effects. Only minor solvent effects were observed for most of the media employed and the measured structural data are in good agreement with microwave data and theoretical predictions.

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.

Morphology of magnetically aligning DMPC/DHPC aggregates-perforated sheets, not disks

The morphology of DMPC/DHPC mixtures at total lipid concentration cL = 5% (w/w) and DMPC/DHPC ratio q approximately 3, doped with small amounts of DMPG or CTAB, was investigated. 31P NMR was used to identify the magnetically aligning phase, and cryo-transmission electron microscopy (cryo-TEM) was employed for structural characterization. Magnetic alignment was found to occur between approximately 30 and approximately 45 degrees C, and cryo-TEM showed that the magnetically aligning phase consisted of extended sheets with a lacelike structure. The aggregates are best described as intermediates between two-dimensional networks of flattened, highly branched, cylindrical micelles and lamellar sheets perforated by large irregular holes. DHPC most likely covers the edges of the holes, while DMPC makes up the bilayer bulk of the aggregates. However, 20-43% of the DHPC takes part in the bilayer, corresponding to 6-12% of the bilayer being made up of DHPC. This fraction increases with increasing temperature. At temperatures above 45 degrees C, the aligning phase collapses.

Composite alignment media for the measurement of independent sets of NMR residual dipolar couplings

The measurement of independent sets of NMR residual dipolar couplings (RDCs) in multiple alignment media can provide a detailed view of biomolecular structure and dynamics, yet remains experimentally challenging. It is demonstrated here that independent sets of RDCs can be measured for ubiquitin using just a single alignment medium composed of aligned bacteriophage Pf1 particles embedded in a strained polyacrylamide gel matrix. Using this composite medium, molecular alignment can be modulated by varying the angle between the directors of ordering for the Pf1 and strained gel matrix, or by varying the ionic strength or concentration of the Pf1 particles. This approach offers significant advantages in that greater experimental control can be exercised over the acquisition of multi-alignment RDC data while a homogeneous chemical environment is maintained across all of the measured RDC data.

NMR analysis of protein interactions

Recent technological advances in NMR spectroscopy have alleviated the size limitations for the determination of biomolecular structures in solution. At the same time, novel NMR parameters such as residual dipolar couplings are providing greater accuracy. As this review shows, the structures of protein-protein and protein-nucleic acid complexes up to 50 kDa can now be accurately determined. Although de novo structure determination still requires considerable effort, information on interaction surfaces from chemical shift perturbations is much easier to obtain. Advances in modelling and data-driven docking procedures allow this information to be used for determining approximate structures of biomolecular complexes. As a result, a wealth of information has become available on the way in which proteins interact with other biomolecules. Of particular interest is the fact that these NMR-based methods can be applied to weak and transient protein-protein complexes that are difficult to study by other structural methods.

Correlation of gene and protein structures in the FXYD family proteins

The FXYD family proteins are auxiliary subunits of the Na,K-ATPase, expressed primarily in tissues that specialize in fluid or solute transport, or that are electrically excitable. These proteins range in size from about 60 to 160 amino acid residues, and share a core homology of 35 amino acid residues in and around a single transmembrane segment. Despite their relatively small sizes, they are all encoded by genes with six to nine small exons. We show that the helical secondary structures of three FXYD family members, FXYD1, FXYD3, and FXYD4, determined in micelles by NMR spectroscopy, reflect the structures of their corresponding genes. The coincidence of helical regions, and connecting segments, with the positions of intron-exon junctions in the genes, support the hypothesis that the FXYD proteins may have been assembled from discrete structural modules through exon shuffling.

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.

Toward direct determination of conformations of protein building units from multidimensional NMR experiments VI. Chemical shift analysis of his to gain 3D structure and protonation state information

NMR--chemical shift structure correlations were investigated by using GIAO-RB3LYP/6-311++G(2d,2p) formalism. Geometries and chemical shifts (CSI values) of 103 different conformers of N'-formyl-L-histidinamide were determined including both neutral and charged protonation forms. Correlations between amino acid torsional angle values and chemical shifts were investigated for the first time for an aromatic and polar amino acid residue whose side chain may carry different charges. Linear correlation coefficients of a significant level were determined between chemical shifts and dihedral angles for CSI[(1)H(alpha)]/phi, CSI[(13)C(alpha)]/phi, and CSI[(13)C(alpha)]/psi. Protonation of the imidazole ring induces the upfield shift of CSI[(13)C(alpha)] for positively charged histidines and an opposite effect for the negative residue. We investigated the correspondence of theoretical and experimental (13)C(alpha), (13)C(beta), and (1)H(alpha) chemical shifts and the nine basic conformational building units characteristic for proteins. These three chemical shift values allow the identification of conformational building units at 80% accuracy. These results enable the prediction of additional regular secondary structural elements (e.g., polyProlineII, inverse gamma-turns) and loops beyond the assignment of chemical shifts to alpha-helices and beta-pleated sheets. Moreover, the location of the His residue can be further specified in a beta-sheet. It is possible to determine whether the appropriate residue is located at the middle or in a first/last beta-strand within a beta-sheet based on calculated CSI values. Thus, the attractive idea of establishing local residue specific backbone folding parameters in peptides and proteins by employing chemical shift information (e.g., (1)H(alpha) and (13)C(alpha)) obtained from selected heteronuclear correlation NMR experiments (e.g., 2D-HSQC) is reinforced.

Backbone solution structures of proteins using residual dipolar couplings: application to a novel structural genomics target

Structural genomics (or proteomics) activities are critically dependent on the availability of high-throughput structure determination methodology. Development of such methodology has been a particular challenge for NMR based structure determination because of the demands for isotopic labeling of proteins and the requirements for very long data acquisition times. We present here a methodology that gains efficiency from a focus on determination of backbone structures of proteins as opposed to full structures with all sidechains in place. This focus is appropriate given the presumption that many protein structures in the future will be built using computational methods that start from representative fold family structures and replace as many as 70% of the sidechains in the course of structure determination. The methodology we present is based primarily on residual dipolar couplings (RDCs), readily accessible NMR observables that constrain the orientation of backbone fragments irrespective of separation in space. A new software tool is described for the assembly of backbone fragments under RDC constraints and an application to a structural genomics target is presented. The target is an 8.7 kDa protein from Pyrococcus furiosus, PF1061, that was previously not well annotated, and had a nearest structurally characterized neighbor with only 33% sequence identity. The structure produced shows structural similarity to this sequence homologue, but also shows similarity to other proteins, which suggests a functional role in sulfur transfer. Given the backbone structure and a possible functional link this should be an ideal target for development of modeling methods.

Prediction of charge-induced molecular alignment of biomolecules dissolved in dilute liquid-crystalline phases

Alignment of macromolecules in nearly neutral aqueous lyotropic liquid-crystalline media such as bicelles, commonly used in macromolecular NMR studies, can be predicted accurately by a steric obstruction model (Zweckstetter and Bax, 2000). A simple extension of this model is described that results in improved predictions for both the alignment orientation and magnitude of protein and DNA solutes in charged nematic media, such as the widely used medium of filamentous phage Pf1. The extended model approximates the electrostatic interaction between a solute and an ordered phage particle as that between the solute's surface charges and the electric field of the phage. The model is evaluated for four different proteins and a DNA oligomer. Results indicate that alignment in charged nematic media is a function not only of the solute's shape, but also of its electric multipole moments of net charge, dipole, and quadrupole. The relative importance of these terms varies greatly from one macromolecule to another, and evaluation of the experimental data indicates that these terms scale differently with ionic strength. For several of the proteins, the calculated alignment is sensitive to the precise position of the charged groups on the protein surface. This suggests that NMR alignment measurements can potentially be used to probe protein electrostatics. Inclusion of electrostatic interactions in addition to steric effects makes the extended model applicable to all liquid crystals used in biological NMR to date.

Determination of protein backbone structures from residual dipolar couplings

There are a number of circumstances in which a focus on determination of the backbone structure of a protein, as opposed to a complete all-atom structure, may be appropriate. This is particularly the case for structures determined as a part of a structural genomics initiative in which computational modeling of many sequentially related structures from the backbone of a single family representative is anticipated. It is, however, also the case when the backbone may be a stepping-stone to more targeted studies of ligand interaction or protein-protein interaction. Here an NMR protocol is described that can produce a backbone structure of a protein without the need for extensive experiments directed at side chain resonance assignment or the collection of structural information on side chains. The procedure relies primarily on orientational constraints from residual dipolar couplings as opposed to distance constraints from NOEs. Procedures for sample preparation, data acquisition, and data analysis are described, along with examples from application to small target proteins of a structural genomics project.

Molecular-orientation analysis based on alignment-induced TROSY chemical shift changes

We present a new NMR technique for determining the alignment tensor of a weakly aligned protein using only alignment-induced 15N transverse relaxation optimized spectroscopy (TROSY) chemical shift changes. Alignment-induced TROSY chemical shift changes reflect the combined contributions from two different anisotropic spin interactions including the residual dipolar couplings (RDCs) and the residual chemical shift anisotropy effects (RCSAs). We show here that these two residual anisotropic spin interactions' values, encoded in the TROSY chemical shift changes, can be used to determine a weakly aligned protein's alignment tensor. To prove the significance of this method, we show that our TROSY-based analysis gives the consistent alignment angles with those determined using RDCs for 15N-labeled ubiquitin (8.6 kDa) in an aligned medium, within an uncertainty range estimated by considering experimental and structural noises, being 5 degrees at most. Because our approach requires a pre-determined 15N CSA tensor value, we also estimated the uncertainties associated with the resultant alignment tensor values caused by variation in 15N CSA tensors. In spite of the significant variations in literature-reported 15N CSA tensors, they gave consistent orientation angles within an uncertainty range. These results ensure that our TROSY-based approach is a useful alternative to the RDC-based method to determine the alignment angles especially for large proteins in a weakly aligned state.

Direct observation and characterization of DMPC/DHPC aggregates under conditions relevant for biological solution NMR

We have used cryo-transmission electron microscopy (cryo-TEM) for inspection of aggregates formed by dimyristoylphosphatidylcholine (DMPC) and dihexanoylphosphatidylcholine (DHPC) in aqueous solution at total phospholipid concentrations cL < or = 5% and DMPC/DHPC ratios q < or = 4.0. In combination with ocular inspections, we are able to sketch out this part of phase-diagram at T = 14-80 degrees C. The temperature and the ratio q are the dominating variables for changing sample morphology, while cL to a lesser extent affects the aggregate structure. At q = 0.5, small, possibly disc-shaped, aggregates with a diameter of approximately 6 nm are formed. At higher q-values, distorted discoidal micelles that tend to short cylindrical micelles are observed. The more well-shaped discs have a diameter of around 20 nm. Upon increasing q or the temperature, long slightly flattened cylindrical micelles that eventually branch are formed. A holey lamellar phase finally appears upon further elevation of q or temperature. The implications for biological NMR work are two. First, discs prepared as membrane mimics are frequently much smaller than predicted by current "ideal bicelle" models. Second, the q approximately 3 preparations used for aligning water-soluble biomolecules in magnetic fields consist of perforated lamellar sheets. Furthermore, the discovered sequence of morphological transitions may have important implications for the development of bicelle-based membrane protein crystallization methods.

Simultaneous Assignment of All Diastereotopic Protons in Strychnine Using RDCs: PELG as Alignment Medium for Organic Molecules

The concept of using residual dipolar couplings (RDCs) for the structure determination of organic molecules is applied to the simultaneous assignment of all diastereotopic protons in strychnine. To use this important NMR parameter the molecule has to be aligned in the magnetic field. Here we present a new alignment medium for organic substrates. The optimization of the alignment properties of mixtures of poly-gamma-ethyl-L-glutamate (PELG) and CDCl(3) are described and the alignment properties of PELG at different concentrations are evaluated. A comparison of PELG with poly-gamma-benzyl-L-glutamate (PBLG) shows considerable differences in the magnitude of alignment for strychnine in the two alignment media. PELG induces a lower degree of order and makes the measurement of residual dipolar couplings (RDCs) in strychnine possible. All one-bond C-H RDCs of strychnine in PELG were determined by using 2D heteronuclear single quantum coherence (HSQC) spectroscopy. The strategy for the extraction of RDCs for methylene groups is described in detail. The RDCs and order parameters are used to assign pairs of diastereotopic protons. This methodology can distinguish not only one pair of diastereotopic protons but it can be used to assign all pairs of diastereotopic protons simultaneously. Two different calculation approaches to achieve this task are described in detail.