Nuclear magnetic resonance methods for high molecular weight proteins: a study involving a complex of maltose binding protein and beta-cyclodextrin

Structural basis for the association of PLEKHA7 with membrane-embedded phosphatidylinositol lipids

PLEKHA7 (pleckstrin homology domain containing family A member 7) plays key roles in intracellular signaling, cytoskeletal organization, and cell adhesion, and is associated with multiple human cancers. The interactions of its pleckstrin homology (PH) domain with membrane phosphatidyl-inositol-phosphate (PIP) lipids are critical for proper cellular localization and function, but little is known about how PLEKHA7 and other PH domains interact with membrane-embedded PIPs. Here we describe the structural basis for recognition of membrane-bound PIPs by PLEHA7. Using X-ray crystallography, nuclear magnetic resonance, molecular dynamics simulations, and isothermal titration calorimetry, we show that the interaction of PLEKHA7 with PIPs is multivalent, distinct from a discrete one-to-one interaction, and induces PIP clustering. Our findings reveal a central role of the membrane assembly in mediating protein-PIP association and provide a roadmap for understanding how the PH domain contributes to the signaling, adhesion, and nanoclustering functions of PLEKHA7.

Reconstitution and Characterization of BCL-2 Family Proteins in Lipid Bilayer Nanodiscs

The BCL-2 family proteins are key regulators of programmed cell death or apoptosis, and represent important targets for the development of anticancer drugs. Because their functions are intimately connected with intracellular membranes, it is important to perform structural and activity studies in precisely characterized samples that include phospholipids and capture the features of the native physiological environment as closely as possible. NMR studies and activity assays based on lipid bilayer nanodiscs are ideally suited for this purpose: they enable the conformations and interactions of these proteins to be probed at atomic resolution in their membrane-associated states. Here we describe detailed protocols for generating the protein components and the reconstituted nanodisc samples suitable for NMR studies and functional assays. The protocols focus on the BCL-2 family protein BCL-XL, a dominant inhibitor of programmed cell death and a major anticancer drug target. The protocols are relatively straightforward. Provided care is taken to ensure protein integrity and sample homogeneity, BCL-XL can be readily reconstituted in nanodiscs, with its hydrophobic C-terminal tail anchored through the nanodisc lipid bilayer, and its folded N-terminal head and ligand binding pocket exposed to the aqueous solution. We anticipate that BCL-2 samples prepared with these protocols will advance structural and mechanistic studies for this important protein family.

Characterization of the membrane-inserted C-terminus of cytoprotective BCL-XL

BCL-XL is a dominant inhibitor of apoptosis and a significant anti-cancer drug target. Endogenous BCL-XL is integral to the mitochondrial outer membrane (MOM). BCL-XL reconstituted in detergent-free lipid bilayer nanodiscs is anchored to the nanodisc lipid bilayer membrane by tight association of its C-terminal tail, while the N-terminal head retains the canonical structure determined for water-soluble, tail-truncated BCL-XL, with the surface groove solvent-exposed and available for BH3 ligand binding. To better understand the conformation and dynamics of this key region of BCL-XL we have developed methods for isolating the membrane-embedded C-terminal tail from its N-terminal head and for preparing protein suitable for structural and biochemical studies. Here, we outline the methods for sample preparation and characterization and describe previously unreported structural and dynamics features. We show that the C-terminal tail of BCL-XL forms a transmembrane α-helix that retains a significant degree of conformational dynamics. We also show that the presence of the intact C-terminus destabilizes the soluble state of the protein, and that the small fraction of soluble recombinant protein produced in Escherichia coli is susceptible to proteolytic degradation of C-terminal residues beyond M218. This finding impacts the numerous previous studies where recombinant soluble BCL-XL was presumed to be full-length. Nevertheless, the majority of recombinant BCL-XL produced in E. coli is insoluble and protected from proteolysis. This protein retains the complete C-terminal tail and can be reconstituted in lipid bilayers in a folded and active state.

Solution NMR studies of periplasmic binding proteins and their interaction partners

Periplasmic binding proteins (PBPs) are a crucial part of ATP-binding cassette import systems in Gram-negative bacteria. Central to their function is the ability to undergo a large-scale conformational rearrangement from open-unliganded to closed-liganded, which signals the presence of substrate and starts its translocation. Over the years, PBPs have been extensively studied not only owing to their essential role in nutrient uptake but also because they serve as excellent models for both practical applications (e.g., biosensor technology) and basic research (e.g., allosteric mechanisms). Although much of our knowledge at atomic level has been inferred from the detailed, static pictures afforded by crystallographic studies, nuclear magnetic resonance (NMR) has been able to fill certain gaps in such body of work, particularly with regard to dynamic processes. Here, we review NMR studies on PBPs, and their unique insights on conformation, dynamics, energetics, substrate binding, and interactions with related transport proteins. Based on the analysis of recent paramagnetic NMR results, as well as crystallographic and functional observations, we propose a mechanism that could explain the ability of certain PBPs to achieve a closed conformation in absence of ligand while others seem to remain open until ligand-mediated closure.

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.

¹H, ¹³C and ¹⁵N backbone resonance assignments of the monomeric human M-ficolin fibrinogen-like domain secreted by Brevibacillus choshinensis

M-ficolin, which forms trimer-based multimers, is a pathogen-recognition protein in the innate immune system, and it binds to ligands through its fibrinogen-like (FBG) domain. As the first step toward the elucidation of the molecular basis for pathogen-recognition by the M-ficolin multimers, we assigned the backbone resonances of the monomeric mutant of the M-ficolin FBG domain, recombinantly expressed by Brevibacillus choshinensis. Like the wild-type trimeric FBG domain, the monomeric FBG domain also requires His251, His284 and His297 for the ligand-binding activity, as judged by mutational analyses using zonal affinity chromatography. The secondary structure predicted by the backbone resonance assignments is similar to that of the trimeric FBG domain in the crystal, indicating that the monomeric FBG domain is folded correctly to perform its function.

Molecular structure and peptidoglycan recognition of Mycobacterium tuberculosis ArfA (Rv0899)

Mycobacterium tuberculosis ArfA (Rv0899) is a membrane protein encoded by an operon that is required for supporting bacterial growth in acidic environments. Its C-terminal domain (C domain) shares significant sequence homology with the OmpA-like family of peptidoglycan-binding domains, suggesting that its physiological function in acid stress protection may be related to its interaction with the mycobacterial cell wall. Previously, we showed that ArfA forms three independently structured modules, and we reported the structure of its central domain (B domain). Here, we describe the high-resolution structure and dynamics of the C domain, we identify ArfA as a peptidoglycan-binding protein and we elucidate the molecular basis for its specific recognition of diaminopimelate-type peptidoglycan. The C domain of ArfA adopts the characteristic fold of the OmpA-like family. It exhibits pH-dependent conformational dynamics (with significant heterogeneity at neutral pH and a more ordered structure at acidic pH), which could be related to its acid stress response. The C domain associates tightly with polymeric peptidoglycan isolated from M. tuberculosis and also associates with a soluble peptide intermediate of peptidoglycan biosynthesis. This enabled us to characterize the peptidoglycan binding site where five highly conserved ArfA residues, including two key arginines, establish the specificity for diaminopimelate- but not Lys-type peptidoglycan. ArfA is the first peptidoglycan-binding protein to be identified in M. tuberculosis. Its functions in acid stress protection and peptidoglycan binding suggest a link between the acid stress response and the physicochemical properties of the mycobacterial cell wall.

Mycobacterium tuberculosis Rv0899 adopts a mixed alpha/beta-structure and does not form a transmembrane beta-barrel

The membrane protein Rv0899 (OmpATb) from Mycobacterium tuberculosis, has been proposed to act as an outer membrane porin and to contribute to the bacterium's adaptation to the acidic environment of the phagosome during infection. The gene is restricted to pathogenic mycobacteria and, thus, is an attractive candidate for the development of anti-tuberculosis chemotherapy. The 326-residue protein contains three domains: an N-terminal domain (residues 1-72) that includes a sequence of 20 hydrophobic amino acids required for membrane translocation, a central B domain (residues 73-200) with homology to the conserved putative lipid-binding BON (bacterial OsmY and nodulation) superfamily, and a C domain (residues 201-326) with homology to the OmpA-C-like superfamily of periplasmic peptidoglycan-binding sequences, found in several types of bacterial membrane proteins, including in the C-terminus of the Escherichia coli outer membrane protein OmpA. We have characterized the structure and dynamics of the B and C domains and have determined the three-dimensional structure of the B domain. Rv0899 does not form a transmembrane beta-barrel. Residues 73-326 form a mixed alpha/beta-globular structure, encompassing two independently folded modules corresponding to the B and C domains connected by a flexible linker. The B domain folds with three parallel/antiparallel alpha-helices packed against six parallel/antiparallel beta-strands that form a flat beta-sheet. The core is hydrophobic, while the exterior is polar and predominantly acidic. The structure of a BON homology domain is revealed here for the first time. In light of this unexpected structure, it is hard to reconcile an outer membrane porin activity with the central domain of the protein. The structure of the B domain and the overall architecture of the protein suggest alternative modes of membrane association.

Simultaneous NMR study of protein structure and dynamics using conservative mutagenesis

A novel iterative procedure is described that allows both the orientation and dynamics of internuclear bond vectors to be determined from direct interpretation of NMR dipolar couplings, measured under at least three orthogonal alignment conditions. If five orthogonal alignments are available, the approach also yields information on the degree of motional anisotropy and the direction in which the largest amplitude internal motion of each bond vector takes place. The method is demonstrated for the backbone (15)N-(1)H, (13)C(alpha)-(1)H(alpha), and (13)C(alpha)-13C' interactions in the previously well-studied protein domain GB3, dissolved in a liquid crystalline suspension of filamentous phage Pf1. Alignment variation is achieved by using conservative mutations of charged surface residues. Results indicate remarkably uniform backbone dynamics, with amplitudes that agree well with those of previous (15)N relaxation studies for most residues involved in elements of secondary structure, but larger amplitude dynamics than those found by (15)N relaxation for residues in loop and turn regions. In agreement with a previous analysis of dipolar couplings, the N-H bonds in the second beta-strand, which is involved in antibody recognition, show elevated dynamics with largest amplitudes orthogonal to the chain direction.

Synthesis of 5-fluoropyrimidine nucleotides as sensitive NMR probes of RNA structure

Enzymatic synthesis methods for the fluorinated 5'-triphosphate analogues 5F-UTP and 5F-CTP have been developed to facilitate 19F-labeling of RNAs for biophysical studies. HIV-2 TAR RNAs were synthesized using these analogues by in vitro transcription reactions using T7 RNA polymerase. The uniform incorporation of 5F-U or 5F-C analogues into HIV-2 TAR RNA transcripts does not significantly alter the RNA structure or thermodynamic stability. Fluorine observed homonuclear 19F-19F and heteronuclear 19F-1H NOE experiments providing selective distance information are presented and discussed. The availability of efficient synthesis of 5F-UTP, and for the first time, 5F-CTP, will facilitate the use of 5F-labeled RNAs in structural, ligand binding, and dynamic studies of RNAs using the advantages of 19F-labeling.

NMR identification of transient complexes critical to adenylate kinase catalysis

A fundamental question in protein chemistry is how the native energy landscape of enzymes enables efficient catalysis of chemical reactions. Adenylate kinase is a small monomeric enzyme that catalyzes the reversible conversion of AMP and ATP into two ADP molecules. Previous structural studies have revealed that substrate binding is accompanied by large rate-limiting spatial displacements of both the ATP and AMP binding motifs. In this report a solution-state NMR approach was used to probe the native energy landscape of adenylate kinase in its free form, in complex with its natural substrates, and in the presence of a tight binding inhibitor. Binding of ATP induces a dynamic equilibrium in which the ATP binding motif populates both the open and the closed conformations with almost equal populations. A similar scenario is observed for AMP binding, which induces an equilibrium between open and closed conformations of the AMP binding motif. These ATP- and AMP-bound structural ensembles represent complexes that exist transiently during catalysis. Simultaneous binding of AMP and ATP is required to force both substrate binding motifs to close cooperatively. In addition, a previously unknown unidirectional energetic coupling between the ATP and AMP binding sites was discovered. On the basis of these and previous results, we propose that adenylate kinase belongs to a group of enzymes whose substrates act to shift pre-existing equilibria toward catalytically active states.

Increasing the accuracy of solution NMR structures of membrane proteins by application of residual dipolar couplings. High-resolution structure of outer membrane protein A

The structure determination of membrane proteins is one of the most challenging applications of solution NMR spectroscopy. The paucity of distance information available from the highly deuterated proteins employed requires new approaches in structure determination. Here we demonstrate that significant improvement in the structure accuracy of the membrane protein OmpA can be achieved by refinement with residual dipolar couplings (RDCs). The application of charged polyacrylamide gels allowed us to obtain two alignments and accurately measure numerous heteronuclear dipolar couplings. Furthermore, we have demonstrated that using a large set of RDCs in the refinement can yield a structure with 1 A rms deviation to the backbone of the high-resolution crystal structure. Our simulations with various data sets indicate that dipolar couplings will be critical for obtaining accurate structures of membrane proteins.

Molecular basis of Bcl-xL's target recognition versatility revealed by the structure of Bcl-xL in complex with the BH3 domain of Beclin-1

Beclin-1, originally identified as a Bcl-2 binding protein, is an evolutionarily conserved protein required for autophagy. The direct interaction between Beclin-1 and Bcl-2 or Bcl-xL provides a potential convergence point for apoptosis and autophagy, two programmed cell death processes. Given the functional significance of the interaction between Beclin-1 and Bcl-2/Bcl-xL, we performed detailed biochemical and structural characterizations of this interaction. We demonstrated that the Bcl-xL-binding domain of Beclin-1 contains a BH3 domain. Therefore, Beclin-1 is a new member of the BH3-only family proteins. The structure of Bcl-xL in complex with the Beclin-1 BH3 domain was determined at high resolution by NMR spectroscopy. Although similar to other known BH3 domains, the Beclin-1 BH3 domain displays its own distinct features in the complex with Bcl-xL. Systematic analysis of all known Bcl-xL/BH3 domain complexes helped us to identify the molecular basis underlying the capacity of Bcl-xL to recognize diverse target sequences.

Orientation of deoxyhemoglobin at high magnetic fields: structural insights from RDCs in solution

Human normal adult hemoglobin (Hb A) is a tetrameric protein molecule of ~64 kDa consisting of two identical -chains and two identical -chains of 141 and 146 amino acid residues each and four bound heme moieties. In the oxygen-free form of Hb A, also known as deoxyhemoglobin A (deoxy-Hb A), the hemes are paramagnetic with S = 2. We have measured the one-bond spin-spin couplings (1JNH + 1DNH) on (15N,2H)-labeled deoxy-Hb A in solution as a function of magnetic field strengths from 11.7 to 21.1 T and found that these couplings are linearly proportional to the square of the magnetic field. This field dependence provides an opportunity to extract the residual dipolar couplings (RDCs, 1DNH) and, thus, to compare predictions about the solution structure of deoxy-Hb A to crystal structures for this molecule. Such comparison is essential for our understanding of the structure, dynamics, and function of this allosteric protein under conditions close to the physiological state. This report illustrates the usefulness of using the magnetic-field dependent RDCs to determine the solution structure of a large paramagnetic protein. This method is especially valuable for those proteins whose structures must be determined in an oxygen-free environment.

J-GFT NMR for precise measurement of mutually correlated nuclear spin-spin couplings

G-matrix Fourier transform (GFT) NMR spectroscopy is presented for accurate and precise measurement of chemical shifts and nuclear spin-spin couplings correlated according to spin system. The new approach, named "J-GFT NMR", is based on a largely extended GFT NMR formalism and promises to have a broad impact on projection NMR spectroscopy. Specifically, constant-time J-GFT (6,2)D (HA-CA-CO)-N-HN was implemented for simultaneous measurement of five mutually correlated NMR parameters, that is, 15N backbone chemical shifts and the four one-bond spin-spin couplings 13Calpha-1Halpha, 13Calpha-13C', 15N-13C', and 15N-1HNu. The experiment was applied for measuring residual dipolar couplings (RDCs) in an 8 kDa protein Z-domain aligned with Pf1 phages. Comparison with RDC values extracted from conventional NMR experiments reveals that RDCs are measured with high precision and accuracy, which is attributable to the facts that (i) the use of constant time evolution ensures that signals do not broaden whenever multiple RDCs are jointly measured in a single dimension and (ii) RDCs are multiply encoded in the multiplets arising from the joint sampling. This corresponds to measuring the couplings multiple times in a statistically independent manner. A key feature of J-GFT NMR, i.e., the correlation of couplings according to spin systems without reference to sequential resonance assignments, promises to be particularly valuable for rapid identification of backbone conformation and classification of protein fold families on the basis of statistical analysis of dipolar couplings.

Spin-state selective carbon-detected HNCO with TROSY optimization in all dimensions and double echo-antiecho sensitivity enhancement in both indirect dimensions

A carbon-detected TROSY-optimized experiment correlating 1HN, 15N, and 13C' resonances, referred to as c-TROSY-HNCO is presented, in which the 1HN and 15N TROSY effects are maintained in both indirect dimensions, while the directly detected 13C' is doubly TROSY-optimized with respect to 1HN and 15N. A new strategy for sensitivity enhancement, the so-called double echo-antiecho (dEA), is described and implemented in the c-TROSY-HNCO experiment. dEA offers sensitivity enhancement of square root of 2 in both indirect dimensions and is generally applicable to many multidimensional experiments. A carbon-detected HNCO experiment, c-HNCO, without TROSY optimization and sensitivity enhancement is also designed for comparison purposes. Relaxation simulations show that for a protein with a rotational correlation time of 10 ns or larger, the c-TROSY-HNCO experiment displays comparable or higher signal-to-noise (S/N) ratios than the c-HNCO experiment, although the former selects only 1/4 of the initial magnetization relative to the later. The high resolution afforded in the directly detected carbon dimension allows direct measurement of the doublet splitting to extract 1JCalphaC' scalar and 1DCalphaC' residual dipolar couplings. Simulations indicate that the c-TROSY-HNCO experiment offers higher precision (lower uncertainty) compared to the c-HNCO experiment for larger proteins. The experiments are applied to 15N/13C/2H/[Leu,Val]-methyl-protonated IIBMannose, a protein of molecular mass 18.6 kDa with a correlation time of approximately 10 ns at 30 degrees C. The experimental pairwise root-mean-square deviation for the measured 1JCalphaC' couplings obtained from duplicate experiments is 0.77 Hz. By directly measuring the doublet splitting, the experiments described here are expected to be much more tolerant to nonuniform values of 1JCalphaC' (or 1JCalphaC' + 1DCalphaC' for aligned samples) and pulse imperfections due to the smaller number of applied pulses in the "out-and-stay" coherence transfer in the c-HNCO-TROSY experiment relative to conventional 1H-detected "out-and-back" quantitative J correlation experiments. A carbon-detected TROSY-optimized experiment correlating 1HN, 15N, and 13C' resonances, referred to as c-TROSY-HNCO is presented, in which the 1HN and 15N TROSY effects are maintained in both indirect dimensions, while the directly detected 13C' is doubly TROSY-optimized with respect to 1HN and 15N. A new strategy for sensitivity enhancement, the so-called double echo-antiecho (dEA), is described and implemented in the c-TROSY-HNCO experiment. dEA offers sensitivity enhancement of in both indirect dimensions and is generally applicable to many multidimensional experiments.

Geometry, energetics, and dynamics of hydrogen bonds in proteins: structural information derived from NMR scalar couplings

An accurate description of hydrogen bonds is essential to identify the determinants of protein stability and function as well as folding and misfolding behavior. We describe a method of using J couplings through hydrogen bonds as ensemble-averaged restraints in molecular dynamics simulations. Applications to the cases of ubiquitin and protein G show that these scalar couplings provide powerful structural information that, when used through the methodology that we present here, enables the description of the geometry and energetics of hydrogen bonds with an accuracy approaching that of high-resolution X-ray structures.

NMR-based Model Reveals the Structural Determinants of Mammalian Arylamine N-Acetyltransferase Substrate Specificity

Arylamine N-acetyltransferases (NATs) catalyze the acetylation of arylamines, a key step in the detoxification of many carcinogens. The determinants of NAT substrate specificity are not known, yet this knowledge is required to understand why NAT enzymes acetylate some arylamines, but not others. Here, we use NMR spectroscopy and homology modeling to reveal the structural determinants of arylamine acetylation by NATs. In particular, by using chemical shift perturbation analysis, we have identified residues that play a critical role in substrate binding and catalysis. This study reveals why human NAT1 acetylates the sunscreen additive p-aminobenzoic acid and tobacco smoke carcinogen 4-aminobiphenyl, but not o-toluidine and other arylamines linked to bladder cancer. Our results represent an important step toward predicting whether arylamines present in new products can be detoxified by mammalian NATs.

Conformational fluctuations affect protein alignment in dilute liquid crystal media

The discovery of dilute liquid crystalline media to align biological macromolecules has opened many new possibilities to study protein and nucleic acid structures by NMR spectroscopy. We inspect the basic alignment phenomenon for an ensemble of protein conformations to deduce relative contributions of each member to the residual dipolar coupling signals. We find that molecular fluctuations can affect the alignment and discover a resulting emphasis of certain conformations. However, the internal fluctuations are largely uncorrelated with those of the alignment, implying that proteins have liquidlike molecular surfaces. Furthermore, we consider the implications of a dynamic bias to structure determination using data from the weak alignment method.

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.

Solid-State 35/37Cl NMR Spectroscopy of Hydrochloride Salts of Amino Acids Implicated in Chloride Ion Transport Channel Selectivity: Opportunities at 900 MHz

The results of a detailed systematic chlorine solid-state NMR study of several hydrochloride salts of amino acids implicated in chloride ion transport channel selectivity are reported. (35)Cl and (37)Cl NMR spectra have been obtained for stationary and/or magic-angle spinning powdered samples of the following compounds on 500 and/or 900 MHz spectrometers: DL-arginine HCl monohydrate, L-lysine HCl, L-serine HCl, L-glutamic acid HCl, L-proline HCl, L-isoleucine HCl, L-valine HCl, L-phenylalanine HCl, and glycine HCl. Spectral analyses provide information on the anisotropic properties and relative orientations of the chlorine electric field gradient and chemical shift (CS) tensors, which are intimately related to the local molecular and electronic structure. Data obtained at 900 MHz provide unique examples of the effects of CS anisotropy on the NMR spectrum of a quadrupolar nucleus. The range of chlorine quadrupolar coupling constants (C(Q)) measured, -6.42 to 2.03 MHz, demonstrates the sensitivity of this parameter to the chloride ion environment and suggests the applicability of chlorine solid-state NMR as a novel experimental tool for defining chloride binding environments in larger ion channel systems. Salts of hydrophobic amino acids are observed to tend to exhibit larger values of C(Q) than salts of hydrophilic amino acids. A simple model for rationalizing the observed trend in C(Q) is proposed. For salts for which neutron diffraction structures are available, we identify a quantum chemical method which reproduces experimental values of C(Q) with a root-mean-square deviation of 0.1 MHz and a correlation coefficient of 0.9998. On the basis of this, chlorine NMR tensors are predicted for the Cl(-) binding site in ClC channels.

Structure of the p53 binding domain of HAUSP/USP7 bound to Epstein-Barr nuclear antigen 1 implications for EBV-mediated immortalization

USP7/HAUSP is a key regulator of p53 and Mdm2 and is targeted by the Epstein-Barr nuclear antigen 1 (EBNA1) protein of Epstein-Barr virus (EBV). We have determined the crystal structure of the p53 binding domain of USP7 alone and bound to an EBNA1 peptide. This domain is an eight-stranded beta sandwich similar to the TRAF-C domains of TNF-receptor associated factors, although the mode of peptide binding differs significantly from previously observed TRAF-peptide interactions in the sequence (DPGEGPS) and the conformation of the bound peptide. NMR chemical shift analyses of USP7 bound by EBNA1 and p53 indicated that p53 binds the same pocket as EBNA1 but makes less extensive contacts with USP7. Functional studies indicated that EBNA1 binding to USP7 can protect cells from apoptotic challenge by lowering p53 levels. The data provide a structural and conceptual framework for understanding how EBNA1 might contribute to the survival of Epstein-Barr virus-infected cells.

In-cell NMR spectroscopy

The role of a protein inside a cell is determined by both its location and its conformational state. Although fluorescence techniques are widely used to determine the cellular localization of proteins in vivo, these approaches cannot provide detailed information about a protein's three-dimensional state. This gap, however, can be filled by NMR spectroscopy, which can be used to investigate both the conformation as well as the dynamics of proteins inside living cells. In this chapter we describe technical aspects of these "in-cell NMR" experiments. In particular, we show that in the case of (15)N-labeling schemes the background caused by labeling all cellular components is negligible, while (13)C-based experiments suffer from high background levels and require selective labeling schemes. A correlation between the signal-to-noise ratio of in-cell NMR experiments with the overexpression level of the protein shows that the current detection limit is 150-200 muM (intracellular concentration). We also discuss experiments that demonstrate that the intracellular viscosity is not a limiting factor since the intracellular rotational correlation time is only approximately two times longer than the correlation time in water. Furthermore, we describe applications of the technique and discuss its limitations.

Characterization of Nonspecific Protein−DNA Interactions by 1H Paramagnetic Relaxation Enhancement

Nonspecific protein-DNA interactions play an important role in a variety of contexts related to DNA packaging, nucleoprotein complex formation, and gene regulation. Biophysical characterization of nonspecific protein-DNA interactions at the atomic level poses significant challenges owing to the dynamic nature of such complexes. Although NMR spectroscopy represents a powerful tool for the analysis of dynamic systems, conventional NMR techniques have provided little information on nonspecific protein-DNA interactions. We show that intermolecular (1)H paramagnetic relaxation enhancement (PRE) arising from Mn(2+) chelated to an EDTA-group covalently attached to a thymine base (dT-EDTA-Mn(2+)) in DNA provides a unique approach for probing the global dynamics and equilibrium distribution of nonspecific protein-DNA interactions. For nonspecific DNA binding, similar intermolecular (1)H-PRE profiles are observed on the (1)H resonances of the bound protein when dT-EDTA-Mn(2+) is located at either end of a DNA oligonucleotide duplex. We demonstrate the applicability of this approach to HMG-box proteins and contrast the results obtained for nonspecific DNA binding of the A-box of HMGB-1 (HMGB-1A) with sequence-specific DNA binding of the related SRY protein. Intermolecular (1)H-PRE data demonstrate unambiguously that HMGB-1A binds to multiple sites in multiple orientations even on a DNA fragment as short as 14 base pairs. Combining the (1)H-PRE data with the crystal structure of the HMGB-1 A-box/cisplatin-modified DNA complex allows one to obtain a semiquantitative estimate of the equilibrium populations at the various sites.

On the Interpretation of Residual Dipolar Couplings as Reporters of Molecular Dynamics

The analysis of residual dipolar couplings from an ensemble of conformations to extract molecular dynamics is intricate. The very mechanism that is necessary to perturb overall molecular tumbling to generate nonvanishing residual dipolar couplings gives rise to convoluted data. The measured values are essentially weighted averages over conformations. However, the weights are not simply the populations of conformations. Consequently, the observed order parameter is not exactly the true measure of motion. In the case of paramagnetic alignment, the apparent order parameter is expected to depend on the number of torsions that separate the locus of interest from the paramagnetic site. In the case of alignment due to steric obstruction, the uneven selection of conformations by their differing Saupe order matrices leads to a bias in the residual dipolar couplings-probed molecular dynamics.

Anisotropic small amplitude Peptide plane dynamics in proteins from residual dipolar couplings

The effect of small amplitude anisotropic peptide plane motion on residual dipolar couplings (RDC) measured in proteins has been investigated. RDC averaging effects in the presence of GAF (Gaussian axial fluctuation) motions are found to vary strongly depending on the peptide plane orientation. Even low amplitude dynamics can significantly affect derived alignment tensor parameters if this motion is not taken into account. An analytical description of averaged N-(N)H RDCs is introduced that includes basic GAF-like motion. The averaging depends on the orientation of the peptide plane (alpha, beta, gamma) in the alignment frame and on the motional amplitude (sigma). This expression is used to investigate the presence of anisotropic reorientational dynamics in proteins by incorporating sigma as an additional parameter into the alignment tensor analysis. Average GAF amplitudes (sigma(av)) are determined for secondary structural elements from single experimental N-(N)H RDC data sets from five different proteins, in combination with high-resolution structural models. This yields statistically significant improvement over the static description, and detects sigma(av) values ranging from 14.4 to 17.0 degrees for the different proteins. A higher value of sigma(av) = 20 degrees from loop regions was found using two independent sets of N-(N)H RDC in the protein lysozyme, for which a very high-resolution structure is available. Comparison of fitting behavior over 13 structures from lysozyme of crystal diffraction resolution ranging from 0.9 to 2.1A indicates a small spread of motional amplitudes, demonstrating that the method is robust up to this level of resolution. A combined definition of (alpha)C-C' and N-(N)H RDC under the influence of GAF motions allows simultaneous fitting of both RDC. Application to three proteins leads to similar sigma(av) values and a more significant improvement with respect to the static model. Using the GAF model to describe conformationally averaged RDC is important for two reasons: a more accurate definition of the alignment tensor magnitude can be derived, and the method can be used to detect average small amplitude motions in protein backbones from readily accessible data, on time scales not easily sampled by other NMR techniques.

Structure of the hepatocyte nuclear factor 6alpha and its interaction with DNA

Hepatocyte nuclear factor 6 (HNF-6) belongs to the family of One Cut transcription factors (also known as OC-1) and is essential for the development of the mouse pancreas, gall bladder, and the interhepatic bile ducts. HNF-6 binds to DNA as a monomer utilizing a single cut domain and a divergent homeodomain motif located at its C terminus. Here, we have used NMR methods to determine the solution structures of the 162 amino acid residue DNA-binding domain of the HNF-6alpha protein. The resulting overall structure of HNF-6alpha has two different distinct domains: the Cut domain and the Homeodomain connected by a long flexible linker. Our NMR structure shows that the Cut domain folds into a topology homologous to the POU DNA-binding domain, even though the sequences of these two protein families do not show homology. The DNA contact sequence of the HNF-6alpha was mapped with chemical shift perturbation methods. Our data also show that a proposed CREB-binding protein histone acetyltransferase protein-recruiting sequence, LSDLL, forms a helix and is involved in the hydrophobic core of the Cut domain. The structure implies that this sequence has to undergo structural changes when it interacts with CREB-binding protein.

How much backbone motion in ubiquitin is required to account for dipolar coupling data measured in multiple alignment media as assessed by independent cross-validation?

The magnitude of backbone internal motions in the small protein ubiquitin that needs to be invoked to account for dipolar coupling data measured in multiple alignment media is investigated using an intuitively straightforward approach. This involves simultaneous refinement of the coordinates (against NOE, torsion angle, and dipolar coupling restraints) and optimization of the magnitudes and orientations of the alignment tensors by means of torsion angle simulated annealing and Cartesian space minimization. We show that N-H dipolar couplings in 11 different alignment media and N-C', H(N)-C', and C alpha-C' dipolar coupling in two alignment media can be accounted for, at approximately the level of uncertainty in the experimental data, by a single structure representation. Extension to a two-member ensemble representation which provides the simplest description of anisotropic motions in the form of a two-site jump model (in which the overall calculated dipolar couplings are the averages of the calculated dipolar couplings of the individual ensemble members), results in modest, but significant, improvements in dipolar coupling R-factors for both the working set of couplings used in the refinement and for the free cross-validated set of C alpha-H alpha dipolar couplings recorded in two alignment media. Extensions to larger ensemble sizes do not result in any R-factor improvement for the cross-validated C alpha-H alpha dipolar couplings. With a few notable exceptions, the amplitudes of the anisotropic motions are small, with S(2)(jump) order parameters > or =0.8. Moreover, the structural impact of those few residues that do exhibit larger amplitude motions (S(2)(jump) ranging from 0.3 to 0.8) is minimal and can readily be accommodated by very small backbone atomic rms shifts (<0.5 A) because of compensatory changes in phi and psi backbone torsion angles. In addition, evidence for correlated motions of N-H bond vectors is observed. For most practical applications, however, refinement of NMR structures against dipolar couplings using a single structure representation is adequate and will not adversely impact coordinate accuracy within the limits of the NMR method.

On the Origin of Residual Dipolar Couplings from Denatured Proteins

Effects of steric obstruction on random flight chains are examined. Spatial probability distributions are elaborated to calculate residual dipolar couplings and residual chemical shift anisotropy, parameters that are acquired by NMR spectroscopy from solutes dissolved in dilute liquid crystals. Calculations yield chain length and residue position-dependent values in good agreement with simulations to provide understanding of recently acquired data from denatured proteins.

Tandem PDZ repeats in glutamate receptor–interacting proteins have a novel mode of PDZ domain–mediated target binding

The interaction of the glutamate receptor subunits 2 and 3 (GluR2/3) with multi-PDZ domain glutamate receptor-interacting protein (GRIP) is important for the synaptic trafficking and clustering of the receptors. Binding of GluR2/3 to GRIP requires both the fourth and fifth PDZ domains (PDZ4 and PDZ5) to be covalently linked, although only one PDZ domain is directly involved in binding to the receptor tail. To elucidate the molecular basis of this mode of PDZ domain-mediated target recognition, we solved the solution structures of the PDZ45 tandem and the isolated PDZ4 of GRIP. The two PDZ domains form a compact structure with a fixed interdomain orientation. The interdomain packing and the stable folding of both PDZ domains require a short stretch of amino acids N-terminal to PDZ4 and a conserved linker connecting PDZ4 and PDZ5. PDZ4 contains a deformed aB-bB groove that is unlikely to bind to carboxyl peptides. Instead, the domain stabilizes the structure of PDZ5.

Theoretical Analysis of Residual Dipolar Coupling Patterns in Regular Secondary Structures of Proteins

A new approach to the interpretation of residual dipolar couplings for the regular secondary structures of proteins is presented. This paper deals with the analysis of the steric and chiral requirements of protein secondary structures and establishes a quantitative correlation between structure periodicity and the experimental values of the backbone residual dipolar couplings. Building on the recent interpretation of the periodicity of residual dipolar couplings in alpha-helices (i.e., "dipolar waves"), a general parametric equation for fitting the residual dipolar couplings of any regular secondary structure is derived. This equation interprets the modulation of the residual dipolar couplings' periodicity in terms of the secondary structure orientation with respect to an arbitrary reference frame, laying the groundwork for using backbone residual dipolar couplings as a fast tool for determining protein folding by NMR spectroscopy.

EDTA-derivatized deoxythymidine as a tool for rapid determination of protein binding polarity to DNA by intermolecular paramagnetic relaxation enhancement

EDTA-derivatized deoxythymidine (dT-EDTA), incorporated into DNA and complexed to Fe2+ in the presence of dithiothreitol, is a widely used reagent for sequence-specific cleavage of duplex DNA. Using HPLC/electrospray mass spectrometry, we show that cleavage is specific to Fe2+, and no cleavage occurs when DNA-EDTA is complexed to other metal ions such as Ca2+, Mn2+, and Fe3+ even after many days. Because dT-EDTA can be incorporated at any desired position of a synthetic oligonucleotide, DNA-EDTA is ideally suited for the measurement of intermolecular paramagnetic relaxation enhancement effects between a paramagnetic ion chelated to DNA-EDTA and a bound protein. Measurements on the SRY/DNA-EDTA complex using two double-stranded oligonucleotides bearing dT-EDTA at opposite ends of the sequence indicate that intermolecular 1HN-T2 enhancement by chelated Mn2+ can be used to readily ascertain the polarity of protein binding to DNA and to derive quantitative long-range distance information for structure refinement. In the case of the SRY-DNA complex, excellent agreement between observed and calculated 1HN-T2 paramagnetic relaxation enhancement data can be achieved with insignificant shifts in the atomic coordinates ( approximately 0.25 A for all heavy atoms) while simultaneously satisfying all other experimental restraints. A unique feature of DNA-EDTA is that the relaxation enhancement effect can be tuned by judicious choice of the paramagnetic metal ion, thereby permitting a wide range of long-range intermolecular electron-proton distances, ranging from approximately 9 to 35 A, to be probed.

Self-consistency analysis of dipolar couplings in multiple alignments of ubiquitin

A self-consistency analysis of backbone N-H residual dipolar couplings of ubiquitin collected in 10 different media is described to assess the degree of structural and dynamic heterogeneous behavior across the media. The SECONDA method, which works with and without any structural or dynamic information about the molecular system, is based on a principal component analysis and is very sensitive to the presence of heterogeneities or experimental errors. It is found that the regular secondary structural elements behave highly homogeneously, while small heterogeneities are manifested in the loop region 51-63. Many residues that exhibit increased dynamics in NMR relaxation experiments are inert with respect to changes in the alignment.

Docking of protein-protein complexes on the basis of highly ambiguous intermolecular distance restraints derived from 1H/15N chemical shift mapping and backbone 15N-1H residual dipolar couplings using conjoined rigid body/torsion angle dynamics

A simple and reliable method for docking protein-protein complexes using (1)H(N)/(15)N chemical shift mapping and backbone (15)N-(1)H residual dipolar couplings is presented and illustrated with three complexes (EIN-HPr, IIA(Glc)-HPr, and IIA(Mtl)-HPr) of known structure. The (1)H(N)/(15)N chemical shift mapping data are transformed into a set of highly ambiguous, intermolecular distance restraints (comprising between 400 and 3000 individual distances) with translational and some degree of orientational information content, while the dipolar couplings provide information on relative protein-protein orientation. The optimization protocol employs conjoined rigid body/torsion angle dynamics in simulated annealing calculations. The target function also comprises three nonbonded interactions terms: a van der Waals repulsion term to prevent atomic overlap, a radius of gyration term (E(rgyr)) to avoid expansion at the protein-protein interface, and a torsion angle database potential of mean force to bias interfacial side chain conformations toward physically allowed rotamers. For the EIN-HPr and IIA(Glc)-HPr complexes, all structures satisfying the experimental restraints (i.e., both the ambiguous intermolecular distance restraints and the dipolar couplings) converge to a single cluster with mean backbone coordinate accuracies of 0.7-1.5 A. For the IIA(Mtl)-HPr complex, twofold degeneracy remains, and the structures cluster into two distinct solutions differing by a 180 degrees rotation about the z axis of the alignment tensor. The correct and incorrect solutions which have mean backbone coordinate accuracies of approximately 0.5 and approximately 10.5 A, respectively, can readily be distinguished using a variety of criteria: (a) examination of the overall (1)H(N)/(15)N chemical shift perturbation map (because the incorrect cluster predicts the presence of residues at the interface that experience only minimal chemical shift perturbations; this information is readily incorporated into the calculations in the form of ambiguous intermolecular repulsion restraints); (b) back-calculation of dipolar couplings on the basis of molecular shape; or (c) the E(rgyr) distribution which, because of its global nature, directly reflects the interfacial packing quality. This methodology should be particularly useful for high throughput, NMR-based, structural proteomics.

A novel method for the determination of stereochemistry in six-membered chairlike rings using residual dipolar couplings

A novel method for the determination of the relative stereochemistry of six-membered chairlike ring molecules by residual dipolar couplings is presented. C-H residual dipolar couplings were used to investigate the relative stereochemistry of 4,6-O-ethylidene-d-glucopyranose. For this and similar systems it is not necessary to acquire redundant dipolar couplings and to calculate the orientation order tensor. The presented methodology is a paradigmatic leap for the determination of the relative stereochemistry or remote stereochemistry in this kind of fused ring system. Residual dipolar coupling data were collected by 1D and 2D direct-measurement heteronuclear multiple quantum coherence (HMQC) spectroscopy. It was demonstrated that direct measurement of HMQC was quick and accurate for small molecules at natural abundance.

Improving the accuracy of NMR structures of RNA by means of conformational database potentials of mean force as assessed by complete dipolar coupling cross-validation

The description of the nonbonded contact terms used in simulated annealing refinement can have a major impact on nucleic acid structures generated from NMR data. Using complete dipolar coupling cross-validation, we demonstrate that substantial improvements in coordinate accuracy of NMR structures of RNA can be obtained by making use of two conformational database potentials of mean force: a nucleic acid torsion angle database potential consisting of various multidimensional torsion angle correlations; and an RNA specific base-base positioning potential that provides a simple geometric, statistically based, description of sequential and nonsequential base-base interactions. The former is based on 416 nucleic acid crystal structures solved at a resolution of </=2 A and an R-factor </=25%; the latter is based on 131 RNA crystal structures solved at a resolution of </=3 A and an R-factor of </=25%, and includes both the large and small subunits of the ribosome. The application of these two database potentials is illustrated for the structure refinement of an RNA aptamer/theophylline complex for which extensive NOE and residual dipolar coupling data have been measured in solution.

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.

Direct structure determination using residual dipolar couplings: reaction-site conformation of methionine sulfoxide reductase in solution

Residual dipolar couplings (RDC) from partially aligned molecules provide long-range structural data and are thus particularly well adapted to rapid structure validation or protein fold recognition. Extensive measurements in two alignment media can also provide precise de novo structure from RDC alone. We have applied a novel combination of these approaches to the study of methionine sulfoxide reductase (MsrA) from Erwinia chrysanthemi, a 27 kDa enzyme essential for repairing oxidative stress damage. The tertiary fold was initially validated by comparing backbone RDC to expected values from the crystal structure of the homologous MsrA from Escherichia coli. Good agreement was found throughout the chain, verifying the overall topology of the molecule, with the exception of the catalytically important peptide P196-L202, where strong and systematic RDC violation was observed. No evidence for local differential mobility in this region was detected, implying that the structure of the strand differs in the two molecules. We have therefore applied the de novo approach meccano to determine the conformation of this peptide using only RDC. A single conformation is found that is in agreement with all measured data. The aligned peptide can be docked onto the expected covalence of the rest of the template molecule while respecting its strictly defined relative orientation. In contrast to the structure of MsrA from E. coli, the reactive side chain of Cys200 is oriented toward the interior of the molecule and therefore closer to the catalytic Cys53, obviating the need for previously proposed conformational reorganization prior to formation of this disulfide intermediate. This analysis requires only backbone assignment and uses unambiguously assigned and readily measurable structural data, thereby greatly economizing investigation time compared to established nuclear Overhauser effect- (nOe-) based structure calculation methods.

A novel approach to the retrieval of structural and dynamic information from residual dipolar couplings using several oriented media in biomolecular NMR spectroscopy

The interpretation of residual dipolar couplings in terms of molecular properties of interest is complicated because of difficulties in separating structural and dynamic effects as well as the need to estimate alignment tensor parameters a priori. An approach is introduced here that allows many of these difficulties to be circumvented when data are acquired in multiple alignment media. The method allows the simultaneous extraction of both structural and dynamic information directly from the residual dipolar coupling data, in favorable cases even in the complete absence of prior structural knowledge. Application to the protein ubiquitin indicates greater amplitudes of internal motion than expected from traditional (15)N spin relaxation analysis.

Detection and characterization of xenon-binding sites in proteins by 129Xe NMR spectroscopy

Xenon-binding sites in proteins have led to a number of applications of xenon in biochemical and structural studies. Here we further develop the utility of 129Xe NMR in characterizing specific xenon-protein interactions. The sensitivity of the 129Xe chemical shift to its local environment and the intense signals attainable by optical pumping make xenon a useful NMR reporter of its own interactions with proteins. A method for detecting specific xenon-binding interactions by analysis of 129Xe chemical shift data is illustrated using the maltose binding protein (MBP) from Escherichia coli as an example. The crystal structure of MBP in the presence of 8atm of xenon confirms the binding site determined from NMR data. Changes in the structure of the xenon-binding cavity upon the binding of maltose by the protein can account for the sensitivity of the 129Xe chemical shift to MBP conformation. 129Xe NMR data for xenon in solution with a number of cavity containing phage T4 lysozyme mutants show that xenon can report on cavity structure. In particular, a correlation exists between cavity size and the binding-induced 129Xe chemical shift. Further applications of 129Xe NMR to biochemical assays, including the screening of proteins for xenon binding for crystallography are considered.