Probing Microsecond Time Scale Dynamics in Proteins by Methyl 1H Carr−Purcell−Meiboom−Gill Relaxation Dispersion NMR Measurements. Application to Activation of the Signaling Protein NtrCr

To study microsecond processes by relaxation dispersion NMR spectroscopy, low power deposition and short pulses are crucial and encourage the development of experiments that employ H-1 Carr-Purcell-Meiboom-Gill (CPMG) pulse trains. Herein, a method is described for the comprehensive study of microsecond to millisecond time scale dynamics of methyl groups in proteins, exploiting their high abundance and favorable relaxation properties. In our approach, protein samples are produced using [H-1, C-13]-D-glucose in similar to 100% D2O, which yields CHD2 methyl groups for alanine, valine, threonine, isoleucine, leucine, and methionine residues with high abundance, in an otherwise largely deuterated background. Methyl groups in such samples can be sequence-specifically assigned to near completion, using C-13 TOCSY NMR spectroscopy, as was recently demonstrated (Often, R.; et al. J. Am. Chem. Soc. 2010, 132, 2952-2960). In this Article, NMR pulse schemes are presented to measure H-1 CPMG relaxation dispersion profiles for CHD2 methyl groups, in a vein similar to that of backbone relaxation experiments. Because of the high deuteration level of methyl-bearing side chains, artifacts arising from proton scalar coupling during the CPMG pulse train are negligible, with the exception of Ile-delta 1 and Thr-gamma 2 methyl groups, and a pulse scheme is described to remove the artifacts for those residues. Strong C-13 scalar coupling effects, observed for several leucine residues, are removed by alternative biochemical and NMR approaches. The methodology is applied to the transcriptional activator NtrC(r), for which an inactive/active state transition was previously measured and the motions in the microsecond time range were estimated through a combination of backbone N-15 CPMG dispersion NMR spectroscopy and a collection of experiments to determine the exchange-free component to the transverse relaxation rate. Exchange contributions to the H-1 line width were detected for 21 methyl groups, and these probes were found to collectively report on a local structural rearrangement around the phosphorylation site, with a rate constant of (15.5 +/- 0.5) x 10(3) per second (i.e., tau(ex) = 64.7 +/- 1.9 mu s). The affected methyl groups indicate that, already before phosphorylation, a substantial, transient rearrangement takes place between helices 3 and 4 and strands 4 and 5. This conformational equilibrium allows the protein to gain access to the active, signaling state in the absence of covalent modification through a shift in a pre-existing dynamic equilibrium. Moreover, the conformational switching maps exactly to the regions that differ between the solution NMR structures of the fully inactive and active states. These results demonstrate that a cost-effective and quantitative study of protein methyl group dynamics by H-1 CPMG relaxation dispersion NMR spectroscopy is possible and can be applied to study functional motions on the microsecond time scale that cannot be accessed by backbone N-15 relaxation dispersion NMR. The use of methyl groups as dynamics probes extends such applications also to larger proteins

Altering dimerization specificity by changes in surface electrostatics

Arc repressor forms a homodimer in which the subunits intertwine to create a single globular domain. To obtain Arc sequences that fold preferentially as heterodimers, variants with surface patches of excess positive or negative charge were designed. Several but not all oppositely charged sequence pairs showed preferential heterodimer formation. In the most successful design pair, α helix B of one subunit contained glutamic acids at positions 43, 46, 47, 48, and 50, whereas the other subunit contained lysines or arginines at these positions. A continuum electrostatic model captures many features of the experimental results and suggests that the most successful designs include elements of both positive and negative design

Regulation of Response Regulator Autophosphorylation through Interdomain Contacts*♦

DNA-binding response regulators (RRs) of the OmpR/PhoB subfamily alternate between inactive and active conformational states, with the latter having enhanced DNA-binding affinity. Phosphorylation of an aspartate residue in the receiver domain, usually via phosphotransfer from a cognate histidine kinase, stabilizes the active conformation. Many of the available structures of inactive OmpR/PhoB family proteins exhibit extensive interfaces between the N-terminal receiver and C-terminal DNA-binding domains. These interfaces invariably involve the α4-β5-α5 face of the receiver domain, the locus of the largest differences between inactive and active conformations and the surface that mediates dimerization of receiver domains in the active state. Structures of receiver domain dimers of DrrB, DrrD, and MtrA have been determined, and phosphorylation kinetics were analyzed. Analysis of phosphotransfer from small molecule phosphodonors has revealed large differences in autophosphorylation rates among OmpR/PhoB RRs. RRs with substantial domain interfaces exhibit slow rates of phosphorylation. Rates are greatly increased in isolated receiver domain constructs. Such differences are not observed between autophosphorylation rates of full-length and isolated receiver domains of a RR that lacks interdomain interfaces, and they are not observed in histidine kinase-mediated phosphotransfer. These findings suggest that domain interfaces restrict receiver domain conformational dynamics, stabilizing an inactive conformation that is catalytically incompetent for phosphotransfer from small molecule phosphodonors. Inhibition of phosphotransfer by domain interfaces provides an explanation for the observation that some RRs cannot be phosphorylated by small molecule phosphodonors in vitro and provides a potential mechanism for insulating some RRs from small molecule-mediated phosphorylation in vivo

Energy functions for protein design I: Efficient and accurate continuum electrostatics and solvation

Electrostatics and solvation energies are important for defining protein stability, structural specificity, and molecular recognition. Because these energies are difficult to compute quickly and accurately, they are often ignored or modeled very crudely in computational protein design. To address this problem, we have developed a simple, fast, and accurate approximation for calculating Born radii in the context of protein design calculations. When these approximate Born radii are used with the generalized Born continuum dielectric model, energies calculated by the 106-fold slower finite difference Poisson-Boltzmann model are faithfully reproduced. A similar approach can be used for estimating solvent-accessible surface areas (SASAs). As an independent test, we show that these approximations can be used to accurately predict the experimentally determined pKas of >200 ionizable groups from 15 proteins