Protein dynamics using frequency-dependent order parameters from analysis of NMR relaxation data

Modeling of Cu(II)-based protein spin labels using rotamer libraries

The bifunctional spin label double-histidine copper-(II) capped with nitrilotriacetate [dHis-Cu(II)-NTA], used in conjunction with electron paramagnetic resonance (EPR) methods can provide high-resolution distance data for investigating protein structure and backbone conformational diversity. Quantitative utilization of this data is limited due to a lack of rapid and accurate dHis-Cu(II)-NTA modeling methods that can be used to translate experimental data into modeling restraints. Here, we develop two dHis-Cu(II)-NTA rotamer libraries using a set of recently published molecular dynamics simulations and a semi-empirical meta-dynamics-based conformational ensemble sampling tool for use with the recently developed chiLife bifunctional spin label modeling method. The accuracy of both the libraries and the modeling method are tested by comparing model predictions to experimentally determined distance distributions. We show that this method is accurate with absolute deviation between the predicted and experimental modes between 0.0-1.2 Å with an average of 0.6 Å over the test data used. In doing so, we also validate the generality of the chiLife bifunctional label modeling method. Taken together, the increased structural resolution and modeling accuracy of dHis-Cu(II)-NTA over other spin labels promise improvements in the accuracy and resolution of protein models by EPR.

Beyond structure: Deciphering site-specific dynamics in proteins from double histidine-based EPR measurements

Site-specific dynamics in proteins are at the heart of protein function. While electron paramagnetic resonance (EPR) has potential to measure dynamics in large protein complexes, the reliance on flexible nitroxide labels is limitating especially for the accurate measurement of site-specific β-sheet dynamics. Here, we employed EPR spectroscopy to measure site-specific dynamics across the surface of a protein, GB1. Through the use of the double Histidine (dHis) motif, which enables labeling with a Cu(II) - nitrilotriacetic acid (NTA) complex, dynamics information was obtained for both α-helical and β-sheet sites. Spectral simulations of the resulting CW-EPR report unique site-specific fluctuations across the surface of GB1. Additionally, we performed molecular dynamics (MD) simulations to complement the EPR data. The dynamics observed from MD agree with the EPR results. Furthermore, we observe small changes in gǁ values for different sites, which may be due to small differences in coordination geometry and/or local electrostatics of the site. Taken together, this work expands the utility of Cu(II)NTA-based EPR measurements to probe information beyond distance constraints.

Backbone Motions of Free and Pheromone-Bound Major Urinary Protein I Studied by Molecular Dynamics Simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.

An investigation of the dynamics of spermine bound to duplex and quadruplex DNA by (13)C NMR spectroscopy

A detailed analysis of the (13)C relaxation of (13)C-labelled spermine bound to duplex and quadruplex DNA is presented. T(1), T(2) and heteronuclear NOE data were collected at four (13)C frequencies (75.4, 125.7, 150.9 and 201.2 MHz). The data were analyzed in terms of a frequency-dependent order parameter, S (2)(omega), to estimate the generalized order parameter and the contributions to the relaxation from different motional frequencies in the picosecond-nanosecond timescale and from any exchange processes that may be occurring on the microsecond-millisecond timescale. The relaxation data was surprisingly similar for spermine bound to two different duplexes and a linear parallel quadruplex. Analysis of the relaxation data from these complexes confirmed the conclusions of previous studies that the dominant motion of spermine is independent of the macroscopic tumbling of the DNA and has an effective correlation time of approximately 50 ps. In contrast, spermine bound to a folded antiparallel quadruplex had faster relaxation rates, especially R (2). As with the other complexes, a fast internal motion of the order of 50 ps makes a substantial contribution to the relaxation. The generalized order parameter for spermine bound to duplex DNA and the linear quadruplex is small but is larger for spermine bound to the folded quadruplex. In the latter case, there is evidence for exchange between at least two populations of spermine occurring on the microsecond-millisecond timescale.

Dynamics of the Escherichia coli O91 O-Antigen Polysaccharide in Solution as Studied by Carbon-13 NMR Relaxation

The dynamics of the O-antigen part of the lipopolysaccharide from the enterohemorrhagic Escherichia coli O91 has been determined in solution using (13)C NMR relaxation measurements at two magnetic field strengths, 9.4 and 14.1 T, thereby facilitating the testing of several dynamical models. The biological repeating unit, consisting of five sugar residues and substituents, could be determined by spectral analysis of different (1)H,(13)C correlations and corroborated by the relaxation data. The site specifically (13)C-labeled material was shown to have approximately 10 repeating units with a narrow distribution. A model-free analysis of the relaxation data revealed a complex dynamical behavior where the sugar residues could be described by a global correlation time (tau(m) = 5.4 ns), generalized order parameters (S(2) approximately 0.63), and different correlation times for internal motions related to their position in the repeating unit along the polymer (tau(e) approximately 360-520 ps). One of the sugar residues showed, in addition, a chemical exchange contribution. Furthermore, a substituent on another sugar residue was described by two order parameters (S(f)(2) = 0.51 and S(s)(2) = 0.21). The solution dynamics of the polysaccharide are thus described by highly intricate motions, both in amplitude and time scales. These results are of significance in the general description of polysaccharides surrounding bacterial cell surfaces and in the presentation of antigenic epitopes to the immune system of an invaded host.

Heat capacities and a snapshot of the energy landscape in protein GB1 from the pre-denaturation temperature dependence of backbone NH nanosecond fluctuations

Protein stability is usually characterized calorimetrically by a melting temperature and related thermodynamic parameters. Despite its importance, the microscopic origin of the melting transition and the relationship between thermodynamic stability and dynamics remains a mystery. Here, NMR relaxation parameters were acquired for backbone 15NH groups of the 56 residue immunoglobulin-binding domain of streptococcal protein G over a pre-denaturation temperature range of 5-50 degrees C. Relaxation data were analyzed using three methods: the standard three-Lorentzian model free approach; the F(omega)=2omegaJ(omega) spectral density approach that yields motional correlation time distributions, and a new approach that determines frequency-dependent order parameters. Regardless of the method of analysis, the temperature dependence of internal motional correlation times and order parameters is essentially the same. Nanosecond time-scale internal motions are found for all NHs in the protein, and their temperature dependence yields activation energies ranging up to about 33kJ/mol residue. NH motional barrier heights are structurally correlated, with the largest energy barriers being found for residues in the most "rigid" segments of the fold: beta-strands 1 and 4 and the alpha-helix. Trends in this landscape also parallel the free energy of folding-unfolding derived from hydrogen-deuterium (H-D) exchange measurements, indicating that the energetics for internal motions occurring on the nanosecond time-scale mirror those occurring on the much slower time-scale of H-D exchange. Residual heat capacities, derived from the temperature dependence of order parameters, range from near zero to near 100J/mol K residue and correlate with this energy landscape. These results provide a unique picture of this protein's energy landscape and a relationship between thermodynamic stability and dynamics that suggests thermosensitive regions in the fold that could initiate the melting process.