Methods of NMR structure refinement: molecular dynamics simulations improve the agreement with measured NMR data of a C-terminal peptide of GCN4-p1

The C-terminal trigger sequence is essential in the coiled-coil formation of GCN4-p1; its conformational properties are thus of importance for understanding this process at the atomic level. A solution NMR model structure of a peptide, GCN4p16–31, encompassing the GCN4-p1 trigger sequence was proposed a few years ago. Derived using a standard single-structure refinement protocol based on 172 nuclear Overhauser effect (NOE) distance restraints, 14 hydrogen-bond and 11 ϕ torsional-angle restraints, the resulting set of 20 NMR model structures exhibits regular α-helical structure. However, the set slightly violates some measured NOE bounds and does not reproduce all 15 measured ³J(H_N-H_Cα)-coupling constants, indicating that different conformers of GCN4p16–31 might be present in solution. With the aim to resolve structures compatible with all NOE upper distance bounds and ³J-coupling constants, we executed several structure refinement protocols employing unrestrained and restrained molecular dynamics (MD) simulations with two force fields. We find that only configurational ensembles obtained by applying simultaneously time-averaged NOE distance and ³J-coupling constant restraining with either force field reproduce all the experimental data. Additionally, analyses of the simulated ensembles show that the conformational variability of GCN4p16–31 in solution admitted by the available set of 187 measured NMR data is larger than represented by the set of the NMR model structures. The conformations of GCN4p16–31 in solution differ in the orientation not only of the side-chains but also of the backbone. The inconsistencies between the NMR model structures and the measured NMR data are due to the neglect of averaging effects and the inclusion of hydrogen-bond and torsional-angle restraints that have little basis in the primary, i.e. measured NMR data.     

1H-NMR Relaxation Studies of Glycopeptides: A Dynamic Structural Investigation

Abstract

The Glycopeptide Man₅GlcNAc₄Asn (ACCB2) in water solution has been studied by means of ¹H NMR relaxation techniques in order to define molecular structure and dynamics. From the analysis of selective and non-selective proton relaxation rates of selected ACCB2 protons, a lack of internal mobility along the polysaccaride chain was observed. The presence of a conformationally well-defined molecular structure for ACCB2 is proposed.     

Applications of tritium NMR to macromolecules: A study of two nucleic acid molecules

Summary We have tritium labeled two nucleic acid molecules, an 8 kDa DNA oligomer and a 20 kDa hammerhead RNA for tritium NMR investigations. The DNA sequence studied has been previously used in homonuclear studies of DNA-bound water molecules and tritium NMR was expected to facilitate these investigations by eliminating the need to suppress the water resonance in tritium-detected ³H-¹H NOESY experiments. We observed the anticipated through-space interactions found in B-form DNA in the NOESY experiments and an unexpected antiphase cross-peak at the water frequency. T₁ measurements on the tritiated DNA molecule indicated that relaxation rates were also accelerated for tritium and protons. Tritium NMR spectra of the hammerhead RNA molecule indicated conformational dynamics in the conserved region of the molecule in the absence of Mg²⁺ and spermine, two components necessary for cleavage. The dynamics were also investigated by ¹⁵N-correlated ¹H spectroscopy and persisted after the addition of Mg²⁺ and spermine.     

Tunable paramagnetic relaxation enhancements by [Gd(DPA)<Subscript>3</Subscript>]<Superscript>3−</Superscript> for protein structure analysis

Paramagnetic relaxation enhancements (PRE) present a powerful source of structural information in nuclear magnetic resonance (NMR) studies of proteins and protein–ligand complexes. In contrast to conventional PRE reagents that are covalently attached to the protein, the complex between gadolinium and three dipicolinic acid (DPA) molecules, [Gd(DPA)₃]^3?, can bind to proteins in a non-covalent yet site-specific manner. This offers straightforward access to PREs that can be scaled by using different ratios of [Gd(DPA)₃]^3? to protein, allowing quantitative distance measurements for nuclear spins within about 15 Å of the Gd³⁺ ion. Such data accurately define the metal position relative to the protein, greatly enhancing the interpretation of pseudocontact shifts induced by [Ln(DPA)₃]^3? complexes of paramagnetic lanthanide (Ln³⁺) ions other than gadolinium. As an example we studied the quaternary structure of the homodimeric GCN4 leucine zipper.     

Refinement of protein structure against non-redundant carbonyl <Superscript>13</Superscript>C NMR relaxation

Carbonyl ¹³C′ relaxation is dominated by the contribution from the ¹³C′ chemical shift anisotropy (CSA). The relaxation rates provide useful and non-redundant structural information in addition to dynamic parameters. It is straightforward to acquire, and offers complimentary structural information to the ¹⁵N relaxation data. Furthermore, the non-axial nature of the ¹³C′ CSA tensor results in a T₁/T₂ value that depends on an additional angular variable even when the diffusion tensor of the protein molecule is axially symmetric. This dependence on an extra degree of freedom provides new geometrical information that is not available from the NH dipolar relaxation. A protocol that incorporates such structural restraints into NMR structure calculation was developed within the program Xplor-NIH. Its application was illustrated with the yeast Fis1 NMR structure. Refinement against the ¹³C′ T₁/T₂ improved the overall quality of the structure, as evaluated by cross-validation against the residual dipolar coupling as well as the ¹⁵N relaxation data. In addition, possible variations of the CSA tensor were addressed. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Sliding and target location of DNA-binding proteins:an NMR view of the lac repressor system

In non-specific lac headpiece-DNA complexes selective NMR line broadening is observed that strongly depends on length and composition of the DNA fragments. This broadening involves amide protons found in the non-specific lac-DNA structure to be interacting with the DNA phosphate backbone, and can be ascribed to DNA sliding of the protein along the DNA. This NMR exchange broadening has been used to estimate the 1D diffusion constant for sliding along non-specific DNA. The observed 1D diffusion constant of 4×10^?12 cm²/s is two orders of magnitude smaller than derived from previous kinetic experiments, but falls in the range of values determined more recently using single molecule methods. This strongly supports the notion that sliding could play at most a minor role in the association kinetics of binding of lac repressor to lac operator and that other processes such as hopping and intersegment transfer contribute to facilitate the DNA recognition process.     

Water and Backbone Dynamics in a Hydrated Protein

Rotational immobilization of proteins permits characterization of the internal peptide and water molecule dynamics by magnetic relaxation dispersion spectroscopy. Using different experimental approaches, we have extended measurements of the magnetic field dependence of the proton-spin-lattice-relaxation rate by one decade from 0.01 to 300 MHz for ¹H and showed that the underlying dynamics driving the protein ¹H spin-lattice relaxation is preserved over 4.5 decades in frequency. This extension is critical to understanding the role of ¹H₂O in the total proton-spin-relaxation process. The fact that the protein-proton-relaxation-dispersion profile is a power law in frequency with constant coefficient and exponent over nearly 5 decades indicates that the characteristics of the native protein structural fluctuations that cause proton nuclear spin-lattice relaxation are remarkably constant over this wide frequency and length-scale interval. Comparison of protein-proton-spin-lattice-relaxation rate constants in protein gels equilibrated with ²H₂O rather than ¹H₂O shows that water protons make an important contribution to the total spin-lattice relaxation in the middle of this frequency range for hydrated proteins because of water molecule dynamics in the time range of tens of ns. This water contribution is with the motion of relatively rare, long-lived, and perhaps buried water molecules constrained by the confinement. The presence of water molecule reorientational dynamics in the tens of ns range that are sufficient to affect the spin-lattice relaxation driven by ¹H dipole-dipole fluctuations should make the local dielectric properties in the protein frequency dependent in a regime relevant to catalytically important kinetic barriers to conformational rearrangements.     

Direct force measurement of single DNA–peptide interactions using atomic force microscopy

The selective interactions between DNA and miniature (39 residues) engineered peptide were directly measured at the single‐molecule level by using atomic force microscopy. This peptide (p007) contains an α‐helical recognition site similar to leucine zipper GCN4 and specifically recognizes the ATGAC sequence in the DNA with nanomolar affinity. The average rupture force was 42.1 pN, which is similar to the unbinding forces of the digoxigenin–antidigoxigenin complex, one of the strongest interactions in biological systems. The single linear fit of the rupture forces versus the logarithm of pulling rates showed a single energy barrier with a transition state located at 0.74 nm from the bound state. The smaller k_off compared with that of other similar systems was presumably due to the increased stability of the helical structure by putative folding residues in p007. This strong sequence‐specific DNA–peptide interaction has a potential to be utilized to prepare well‐defined mechanically stable DNA–protein hybrid nanostructures. Copyright © 2013 John Wiley & Sons, Ltd.     

Direct NMR detection of the "invisible" alkali metal cations tightly bound to G-quadruplex structures

We report the first direct solution NMR detection of the alkali metal cations (²³Na⁺, ³⁹K⁺, and ⁸⁷Rb⁺) residing inside G-quadruplex channel structures formed by guanosine 5′-monophosphate and a DNA oligomer, d(TG₄T). In solution, these channel alkali metal cations are tightly bound to the G-quadruplex structure and have been considered to be “invisible” to NMR spectroscopy for many years. Our finding that it is possible to directly observe these alkali metal cations by NMR spectroscopy provides a new tool for studying cation binding affinity and dynamics in G-quadruplex DNA.     

Recognition of T-rich single-stranded DNA by the cold shock protein Bs-CspB in solution

Cold shock proteins (CSP) belong to the family of single-stranded nucleic acid binding proteins with OB-fold. CSP are believed to function as ‘RNA chaperones’ and during anti-termination. We determined the solution structure of Bs-CspB bound to the single-stranded DNA (ssDNA) fragment heptathymidine (dT₇) by NMR spectroscopy. Bs-CspB reveals an almost invariant conformation when bound to dT₇ with only minor reorientations in loop β1–β2 and β3–β4 and of few aromatic side chains involved in base stacking. Binding studies of protein variants and mutated ssDNA demonstrated that Bs-CspB associates with ssDNA at almost diffusion controlled rates and low sequence specificity consistent with its biological function. A variation of the ssDNA affinity is accomplished solely by changes of the dissociation rate. ¹⁵N NMR relaxation and H/D exchange experiments revealed that binding of dT₇ increases the stability of Bs-CspB and reduces the sub-nanosecond dynamics of the entire protein and especially of loop β3–β4.     

The Sortase A Enzyme That Attaches Proteins to the Cell Wall of Bacillus anthracis Contains an Unusual Active Site Architecture

The pathogen Bacillus anthracis uses the Sortase A (SrtA) enzyme to anchor proteins to its cell wall envelope during vegetative growth. To gain insight into the mechanism of protein attachment to the cell wall in B. anthracis we investigated the structure, backbone dynamics, and function of SrtA. The NMR structure of SrtA has been determined with a backbone coordinate precision of 0.40 ± 0.07 Å. SrtA possesses several novel features not previously observed in sortase enzymes including the presence of a structurally ordered amino terminus positioned within the active site and in contact with catalytically essential histidine residue (His¹²⁶). We propose that this appendage, in combination with a unique flexible active site loop, mediates the recognition of lipid II, the second substrate to which proteins are attached during the anchoring reaction. pK_a measurements indicate that His¹²⁶ is uncharged at physiological pH compatible with the enzyme operating through a “reverse protonation” mechanism. Interestingly, NMR relaxation measurements and the results of a model building study suggest that SrtA recognizes the LPXTG sorting signal through a lock-in-key mechanism in contrast to the prototypical SrtA enzyme from Staphylococcus aureus.     

Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Structural and dynamics characterization of norovirus protease

Norovirus protease is an essential enzyme for proteolytic maturation of norovirus nonstructural proteins and has been implicated as a potential target for antiviral drug development. Although X‐ray structural studies of the protease give us wealth of structural information including interactions of the protease with its substrate and dimeric overall structure, the role of protein dynamics in the substrate recognition and the biological relevance of the protease dimer remain unclear. Here we determined the solution NMR structure of the 3C‐like protease from Norwalk virus (NV 3CLpro), a prototype strain of norovirus, and analyzed its backbone dynamics and hydrodynamic behavior in solution. ¹⁵N spin relaxation and analytical ultracentrifugation analyses demonstrate that NV 3CLpro is predominantly a monomer in solution. Solution structure of NV 3CLpro shows significant structural variation in C‐terminal domain compared with crystal structures and among lower energy structure ensembles. Also, ¹⁵N spin relaxation and Carr–Purcell–Meiboom–Gill (CPMG)‐based relaxation dispersion analyses reveal the dynamic properties of residues in the C‐terminal domain over a wide range of timescales. In particular, the long loop spanning residues T123–G133 show fast motion (ps‐ns), and the residues in the bII–cII region forming the large hydrophobic pocket (S2 site) undergo conformational exchanges on slower timescales (μs–ms), suggesting their important role in substrate recognition.     

Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR 15N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Conformational dynamics of recoverin's Ca2+-myristoyl switch probed by 15N NMR relaxation dispersion and chemical shift analysis

Recoverin, a member of the neuronal calcium sensor (NCS) branch of the calmodulin superfamily, serves as a calcium sensor in retinal rod cells. Ca²⁺‐induced conformational changes in recoverin promote extrusion of its covalently attached myristate, known as the Ca²⁺‐myristoyl switch. Here, we present nuclear magnetic resonance (NMR) relaxation dispersion and chemical shift analysis on ¹⁵N‐labeled recoverin to probe main chain conformational dynamics. ¹⁵N NMR relaxation data suggest that Ca²⁺‐free recoverin undergoes millisecond conformational dynamics at particular amide sites throughout the protein. The addition of trace Ca²⁺ levels (0.05 equivalents) increases the number of residues that show detectable relaxation dispersion. The Ca²⁺‐dependent chemical shifts and relaxation dispersion suggest that recoverin has an intermediate conformational state (I) between the sequestered apo state (T) and Ca²⁺ saturated extruded state (R): T ? I ? R. The first step is a fast conformational equilibrium ([T]/[I] < 100) on the millisecond time scale (τ_exδω < 1). The final step (I ? R) is much slower (τ_exδω > 1). The main chain structure of I is similar in part to the structure of half‐saturated E85Q recoverin with a sequestered myristoyl group. We propose that millisecond dynamics during T ? I may transiently increase the exposure of Ca²⁺‐binding sites to initiate Ca²⁺ binding that drives extrusion of the myristoyl group during I ? R. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Main-chain Dynamics of a Partially Folded Protein:15N NMR Relaxation Measurements of Hen Egg White Lysozyme Denatured in Trifluoroethanol

¹⁵N NMR relaxation measurements have been used to study the dynamic behaviour of the main-chain of hen lysozyme in a partially folded state, formed in a 70% (v/v) trifluoroethanol (TFE)/30% water mixture at 37°C and pH 2. This state is characterised by helical secondary structure in the absence of extensive tertiary interactions. The NMR relaxation data were interpreted by mapping of spectral density functions and by derivation of segmental as well as global order parameters. The results imply that the dynamics of lysozyme in TFE can, at least for the great majority of residues, be adequately described by internal motions which are superimposed on an overall isotropic tumbling of the molecule. Although the dynamic behaviour shows substantial variations along the polypeptide chain, it correlates well with the conformational preferences identified in the TFE state by other NMR parameters. Segments of the polypeptide chain which are part of persistent helical structures are highly restricted in their motion (S²> 0.8, with effective internal correlation times τ_e< 200 ps) but are also found to experience conformational exchange on a millisecond timescale. Regions which are stabilised in less persistent helical structure possess greater flexibility (0.6 <S²< 0.8, 200 ps < τ_e< 1 ns) and those which lack defined conformational preferences are highly flexible (S²< 0.6, τ_e∼1 ns). The dynamic behaviour of the main-chain was found to be correlated with other local features of the polypeptide chain, including hydrophobicity and the position of the disulphide bridges. Despite the absence of extensive tertiary interactions, preferential stabilisation of native-like secondary structure by TFE results in a pattern of main-chain dynamics which is similar to that of the native state.     

Structure and dynamic properties of the single disulfide-deficient α-amylase inhibitor [C45A/C73A]tendamistat: An NMR study

Covalent linkages such as disulfide bonds are important for the stabilization of proteins. In the present NMR study we compare the structure and the dynamics of the single disulfide-deficient variant C45A/C73A of the α-amylase inhibitor tendamistat and the wild-type protein, which contains two disulfide bonds (C11-C27 and C45-C73). Complete proton assignment was achieved by standard homonuclear 2D techniques for the variant. Chemical shift differences, intra-strand NOE effects and protected amide proton were used to compare the connectivity of the secondary structure elements of variant and wild-type. Dynamic properties of the wild-type protein were studied by ¹³C_α heteronuclear NOE experiments with carbon in natural abundance. ¹⁵N isotope labeling was necessary to obtain the relaxation parameters of the variant, because of sample degradation. The ¹⁵N resonance assignment was achieved by a ¹⁵N 3D-NOESY-HMQC. Removal of the C45-C73 bond by the C45A/C73A mutation has no influence upon the β-barrel structure of tendamistat beside very local changes at the mutation site. The relaxation data revealed only subtle differences between variant and wild-type on a subnanosecond time scale. Only the N-terminus and G62 in the connecting loop between the anti-parallel β-sheets showed an increased mobility. The results are discussed in respect to thermodynamic stability and the secretion efficiency of tendamistat. Proteins 33:285–294, 1998. © 1998 Wiley-Liss, Inc.