With the advent of ultra-long MD simulations it becomes possible to model microsecond time-scale protein dynamics and, in particular, the exchange broadening effects (R(ex)) as probed by NMR relaxation dispersion measurements. This new approach allows one to identify the exchanging species, including the elusive "excited states". It further helps to map out the exchange network, which is potentially far more complex than the commonly assumed 2- or 3-site schemes. Under fast exchange conditions, this method can be useful for separating the populations of exchanging species from their respective chemical shift differences, thus paving the way for structural analyses. In this study, recent millisecond-long MD trajectory of protein BPTI (Shaw et al. Science 2010, 330, 341) is employed to simulate the time variation of amide (15)N chemical shifts. The results are used to predict the exchange broadening of (15)N lines and, more generally, the outcome of the relaxation dispersion measurements using Carr-Purcell-Meiboom-Gill sequence. The simulated R(ex) effect stems from the fast (~10-100 μs) isomerization of the C14-C38 disulfide bond, in agreement with the prior experimental findings (Grey et al. J. Am. Chem. Soc. 2003, 125, 14324).

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as 1 H N is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the 1 H N residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between 1 H N spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.Additional Supporting Information may be found in the online version of this article.Abbreviations: drkN SH3, N-terminal SH3 domain of the Drosophila adapter protein drk; EDTA, ethylenediaminetetraacetic acid; ESR, electron spin resonance; FRET, fluorescence resonance energy transfer; MTSL, methanethiosulfonate spin label; NOE, nuclear Overhauser effect; NAG, N-acetyl-glycine.

In this study, we have shown that substitution of chloride ligand for imidazole (Im) ring in the cyclometalated platinum complex Pt(phpy)(PPh)Cl (1; phpy, 2-phenylpyridine; PPh, triphenylphosphine), which is nonemissive in solution, switches on phosphorescence of the resulting compound. Crystallographic and nuclear magnetic resonance (NMR) spectroscopic studies of the substitution product showed that the luminescence ignition is a result of Im coordination to give the [Pt(phpy)(Im)(PPh)]Cl complex. The other imidazole-containing biomolecules, such as histidine and histidine-containing peptides and proteins, also trigger luminescence of the substitution products. The complex 1 proved to be highly selective toward the imidazole ring coordination that allows site-specific labeling of peptides and proteins with 1 using the route, which is orthogonal to the common bioconjugation schemes via lysine, aspartic and glutamic acids, or cysteine and does not require any preliminary modification of a biomolecule. The utility of this approach was demonstrated on (i) site-specific modification of the ubiquitin, a small protein that contains only one His residue in its sequence, and (ii) preparation of nonaggregated HSA-based Pt phosphorescent probe. The latter particles easily internalize into the live HeLa cells and display a high potential for live-cell phosphorescence lifetime imaging (PLIM) as well as for advanced correlation PLIM and FLIM experiments.

Backbone ( 15 N) NMR relaxation is one of the main sources of information on dynamics of disordered proteins. Yet, we do not know very well what drives 15 N relaxation in such systems, i.e., how different forms of motion contribute to the measurable relaxation rates. To address this problem, we have investigated, both experimentally and via molecular dynamics simulations, the dynamics of a 26-residue peptide imitating the N-terminal portion of the histone protein H4. One part of the peptide was found to be fully flexible, whereas the other part features some transient structure (a hairpin stabilized by hydrogen bonds). The following motional modes proved relevant for 15 N relaxation. 1) Sub-picosecond librations attenuate relaxation rates according to S 2 $0.85-0.90. 2) Axial peptide-plane fluctuations along a stretch of the peptide chain contribute to relaxation-active dynamics on a fast timescale (from tens to hundreds of picoseconds). 3) 4/j backbone jumps contribute to relaxation-active dynamics on both fast (from tens to hundreds of picoseconds) and slow (from hundreds of picoseconds to a nanosecond) timescales. The major contribution is from polyproline II (PPII) 4 b transitions in the Ramachandran space; in the case of glycine residues, the major contribution is from PPII 4 (b) 4 rPPII transitions, in which rPPII is the mirror-image (right-handed) version of the PPII geometry, whereas b geometry plays the role of an intermediate state. 4) Reorientational motion of certain (sufficiently long-lived) elements of transient structure, i.e., rotational tumbling, contributes to slow relaxation-active dynamics on $1-ns timescale (however, it is difficult to isolate this contribution). In conclusion, recent advances in the area of force-field development have made it possible to obtain viable Molecular Dynamics models of protein disorder. After careful validation against the experimental relaxation data, these models can provide a valuable insight into mechanistic origins of spin relaxation in disordered peptides and proteins.

No abstract