Molecular motion of spin labeled side chains in the C-terminal domain of RGL2 protein: a SDSL-EPR and MD study

Five singly spin labeled side chains at surface sites in the C-terminal domain of RGL2 protein have been analyzed to investigate the general relationship between nitroxide side chain mobility and protein structure. At these sites, the structural perturbation produced by replacement of a native residue with a nitroxide side chain appears to be very slight at the level of the backbone fold. The primary determinants of the nitroxide side chain mobility are backbone dynamics and tertiary interactions. On the exposed surfaces of alpha-helices, the side chain mobility is not restricted by tertiary interactions but appears to be determined by backbone dynamics, while in loop sites, the side chain mobility is even higher. For a better understanding of the changes in the EPR spectral line shape, molecular dynamics simulations were performed and found in agreement with EPR spectral data.     

Crystal structures of spin labeled T4 lysozyme mutants: implications for the interpretation of EPR spectra in terms of structure

High resolution (1.43-1.8 A) crystal structures and the corresponding electron paramagnetic resonance (EPR) spectra were determined for T4 lysozyme derivatives with a disulfide-linked nitroxide side chain [-CH(2)-S-S-CH(2)-(3-[2,2,5,5-tetramethyl pyrroline-1-oxyl]) identical with R1] substituted at solvent-exposed helix surface sites (Lys65, Arg80, Arg119) or a tertiary contact site (Val75). In each case, electron density is clearly resolved for the disulfide group, revealing distinct rotamers of the side chain, defined by the dihedral angles X(1) and X(2). The electron density associated with the nitroxide ring in the different mutants is inversely correlated with its mobility determined from the EPR spectrum. Residue 80R1 assumes a single g(+)()g(+)() conformation (Chi(1) = 286, X(2) = 294). Residue 119R1 has two EPR spectral components, apparently corresponding to two rotamers, one similar to that for 80R1 and the other in a tg(-)() conformation (Chi(1) = 175, X(2) = 54). The latter state is apparently stabilized by interaction of the disulfide with a Gln at i + 4, a situation also observed at 65R1. R1 residues at helix surface site 65 and tertiary contact site 75 make intra- as well as intermolecular contacts in the crystal and serve to identify the kind of molecular interactions possible for the R1 side chain. A single conformation of the entire 75R1 side chain is stabilized by a variety of interactions with the nitroxide ring, including hydrophobic contacts and two unconventional C-H.O hydrogen bonds, one in which the nitroxide acts as a donor (with tyrosine) and the other in which it acts as an acceptor (with phenylalanine). The interactions revealed in these structures provide an important link between the dynamics of the R1 side chain, reflected in the EPR spectrum, and local protein structure. A library of such interactions will provide a basis for the quantitative interpretation of EPR spectra in terms of protein structure and dynamics.     

Structural features and light-dependent changes in the sequence 306-322 extending from helix VII to the palmitoylation sites in rhodopsin: a site-directed spin-labeling study.

Sixteen single-cysteine substitution mutants of rhodopsin were prepared in the sequence 306-321 which begins in transmembrane helix VII and ends at the palmitoylation sites at 322C and 323C. The substituted cysteine residues were modified with a selective reagent to generate a nitroxide side chain, and the electron paramagnetic resonance spectrum of each spin-labeled mutant was analyzed in terms of residue accessibility and mobility. The periodic behavior of these parameters along the sequence indicated that residues 306-314 were in a regular alpha-helical conformation representing the end of helix VII. This helix apparently extends about 1.5 turns above the surface of the membrane, with one face in strong tertiary interaction with the core of the protein. For the segment 315-321, substituted cysteine residues at 317, 318, 320, and 321 had low reactivity with the spin-label reagent. This segment has the most extensive tertiary interactions yet observed in the rhodopsin extra-membrane sequences at the cytoplasmic surface. Previous studies showed the spontaneous formation of a disulfide bond between cysteine residues at 65 and 316. This result indicates that at least some of the tertiary contacts made in the 315-321 segment are with the sequence connecting transmembrane helices I and II. Photoactivation of rhodopsin produces changes in structure detected by spin labels at 306, 313, and 316. The changes at 313 can be accounted for by movements in the adjacent helix VI.     

Molecular motion of spin labeled side chains in alpha-helices: analysis by variation of side chain structure

Two single cysteine substitution mutants at helix surface sites in T4 lysozyme (D72C and V131C) have been modified with a series of nitroxide methanethiosulfonate reagents to investigate the structural and dynamical origins of their electron paramagnetic resonance spectra. The novel reagents include 4-substituted derivatives of either the pyrroline or pyrrolidine series of nitroxides. The spectral line shapes were analyzed as a function of side chain structure and temperature using a simulation method with a single order parameter and diffusion rates about three orthogonal axes as parameters. Taken together, the results provide strong support for an anisotropic motional model of the side chain, which was previously proposed from qualitative features of the spectra and crystal structures of spin labeled T4 lysozyme. Site-specific differences in apparent order parameter are interpreted in terms of backbone dynamics modes with characteristic correlation times in the nanosecond or faster time scale. The saturated 4-substituted pyrrolidine nitroxides are shown to be a suitable template for novel "functionalized" side chains designed to mimic salient features of the native side chains they replace.     

Mapping backbone dynamics in solution with site-directed spin labeling: GCN4-58 bZip free and bound to DNA

In site-directed spin labeling, a nitroxide-containing side chain is introduced at selected sites in a protein. The EPR spectrum of the labeled protein encodes information about the motion of the nitroxide on the nanosecond time scale, which has contributions from the rotary diffusion of the protein, from internal motions in the side chain, and from backbone fluctuations. In the simplest model for the motion of noninteracting (surface) side chains, the contribution from the internal motion is sequence independent, as is that from protein rotary diffusion. Hence, differences in backbone motions should be revealed by comparing the sequence-dependent motions of nitroxides at structurally homologous sites. To examine this model, nitroxide side chains were introduced, one at a time, along the GCN4-58 bZip sequence, for which NMR (15)N relaxation experiments have identified a striking gradient of backbone mobility along the DNA-binding region [Bracken et al. (1999) J. Mol. Biol. 285, 2133]. Spectral simulation techniques and a simple line width measure were used to extract dynamical parameters from the EPR spectra, and the results reveal a mobility gradient similar to that observed in NMR relaxation, indicating that side chain motions mirror backbone motions. In addition, the sequence-dependent side chain dynamics were analyzed in the DNA/protein complex, which has not been previously investigated by NMR relaxation methods. As anticipated, the backbone motions are damped in the DNA-bound state, although a gradient of motion persists with residues at the DNA-binding site being the most highly ordered, similar to those of helices on globular proteins.     

Structure of the KcsA potassium channel from Streptomyces lividans: a site-directed spin labeling study of the second transmembrane segment.

KcsA is a prokaryotic potassium channel. The present study employs cysteine scanning mutagenesis and site-directed spin labeling to investigate the structure of the second transmembrane segment (residues 82-120) in functional tetrameric channels reconstituted in lipid bilayers. Spin-spin interactions are observed between nitroxide side chains at symmetry-related sites close to the 4-fold axis of symmetry. To aid in quantitative analysis of these interactions, a new diamagnetic analogue of the nitroxide side chain is used to prepare magnetically dilute samples with constant structure. Using constraints imposed by the spin-spin interactions, a packing model for this segment is deduced that is in excellent agreement with the recently reported crystal structure [Doyle, D., et al. (1998) Science 280, 69-77]. The relatively immobilized state of the nitroxide side chains suggests that the channel is rigid on the electron paramagnetic resonance time scale. Moreover, the poor sulfhydryl reactivity of the cysteine at many locations indicates that the channel is not subject to the low-frequency fluctuations that permit reaction of buried cysteines. At sites expected to be located in the pore, the accessibility of the side chains to collision with O(2) or nickel(II) ethylenediaminediacetate is low. This inaccessibility, together with the generally low mobility of the side chains throughout the sequence, makes it difficult to detect the presence of the pore based on these measurements. However, the presence of a solvated pore can be directly demonstrated using a polarity parameter deduced from the EPR spectra recorded at low temperature. These measurements also reveal the presence of a polarity gradient in the phospholipid bilayer.     

Conformation of the diphtheria toxin T domain in membranes: a site-directed spin-labeling study of the TH8 helix and TL5 loop.

The isolated T domain of diphtheria toxin was mutated by cysteine-scanning mutagenesis at 28 consecutive sites (residues 328-355) that comprise the TH8 helix and the TL5 interhelical loop in the native toxin. After derivatizing the mutant proteins with a sulfhydryl-selective nitroxide reagent, we examined the mobility of each nitroxide and its accessibility to polar and nonpolar paramagnetic reagents, before and after insertion into phospholipid bilayers. The data obtained with the proteins in solution at pH 8 are generally consistent with predictions from the crystal structure of the toxin. Upon membrane binding at pH 4.6, a major structural reorganization of the domain was seen, which dramatically reduced the accessibility of most residues in this region to the polar reagent nickel(II)-ethylenediaminediacetate complex (NiEDDA). Many of these residues also showed reduced accessibility to the nonpolar reagent O(2). Periodic accessibility of the nitroxide side chains along the sequence to these reagents shows that TH8 remains largely helical in the membrane-bound state, with one surface associated with protein and the other facing the hydrophobic interior of the bilayer. In addition, the TL5 loop also appears to become alpha-helical in the membrane, with one surface in contact with protein and the other in contact with the bilayer interior. These findings provide a structural framework for understanding how the T domain forms a transmembrane channel and mediates translocation of diphtheria toxin's enzymic moiety across a membrane.     

Structural features of the C-terminal domain of bovine rhodopsin: a site-directed spin-labeling study.

Eleven single-cysteine substitution mutants have been prepared in the sequence 325-340 of rhodopsin, corresponding to the C-terminal domain. Each of the cysteine mutants was modified with a selective nitroxide reagent to introduce a spin-labeled side chain. The electron paramagnetic resonance spectra of the labeled proteins were analyzed in terms of side chain dynamics. At all sites, the spectra reflected the presence of two populations of different mobility, although one was always dominant. The mobility of the dominant population increased in a regular fashion from the palmitoylation sites at 322C and 323C to the C-terminus, where the spectra resembled those of an unfolded protein. This apparent mobility gradient is only slightly affected in mutants lacking the palmitoyl groups, suggesting that they are not responsible for physically anchoring the C-terminal peptide at one end. Binding of a monoclonal antibody to its epitope at the C-terminus dramatically reduces the mobility of nearby residues, creating a local mobility gradient opposite that in the absence of the antibody. These results indicate that the C-terminal domain of rhodopsin, beyond the palmitoylation sites, is highly disordered and dynamic, resembling an unfolded peptide tethered at one end.     

Structure in the channel forming domain of colicin E1 bound to membranes: the 402-424 sequence

To explore the structure of the pore-forming fragment of colicin E1 in membranes, a series of 23 consecutive single cysteine substitution mutants was prepared in the sequence 402-424. Each mutant was reacted with a sulfhydryl-specific reagent to generate a nitroxide labeled side chain, and the mobility of the side chain and its accessibility to collision with paramagnetic reagents was determined from the electron paramagnetic resonance spectrum. Individual values of these quantities were used to identify tertiary contact sites and the nature of the surrounding solvent, while their periodic dependence on sequence position was used to identify secondary structure. In solution, the data revealed a regular helix of 11 residues in the region 406-416, consistent with helix IV of the crystal structure. Upon binding to negatively charged membranes at pH 4.0, helix IV apparently grows to a length of 19 residues, extending from 402-420. One face of the helix is solvated by the lipid bilayer, and the other by an environment of a polar nature. Surprisingly, a conserved charged pair, D408-R409, is located on the lipid-exposed face. Evidence is presented to suggest a transmembrane orientation of this new helix, although other topographies may exist in equilibrium.     

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.     

Side chain mobility and ligand interactions of the G strand of tear lipocalins by site-directed spin labeling.

Side chain mobility, accessibility, and backbone motion were studied by site-directed spin labeling of sequential cysteine mutants of the G strand in tear lipocalins (TL). A nitroxide scan between residues 98 and 105 revealed the alternating periodicity of mobility and accessibility to NiEDDA and oxygen, characteristic of a beta-strand. Residue 99 was the most inaccessible to NiEDDA and oxygen. EPR spectra with the fast relaxing agent, K(3)Fe(CN)(6), exhibited two nitroxide populations for most residues. The motionally constrained population was relatively less accessible to K(3)Fe(CN)(6) because of dynamic tertiary contact, probably with side chain residues of adjacent strands. With increasing concentrations of sucrose, the spectral contribution of the immobile component was greater, indicating a larger population with tertiary contact. Increased concentrations of sucrose also resulted in a restriction of mobility of spin-labeled fatty acids which were bound within the TL cavity. The data suggest that sucrose enhanced ligand affinity by slowing the backbone motion of the lipocalin. The correlation time of an MTSL derivative (I) attached to F99C resulted in the lack of side chain motion and therefore reflects the overall rotation of the TL complex. The correlation time of F99C in tears (13.5 ns) was the same as that in buffer and indicates that TL exists as a dimer under native conditions. TL-spin-labeled ligand complexes have a shorter correlation time than the protein alone, indicating that the fatty acids are not rigidly anchored in the cavity, but move within the pocket. This segmental motion of the ligand was modulated by protein backbone fluctuations. Accessibility studies with oxygen and NiEDDA were performed to determine the orientation and depth of a series of fatty acid derivatives in the cavity of TL. Fatty acids are oriented with the hydrocarbon tail buried in the cavity and the carboxyl group oriented toward the mouth. In general, the mobility of the nitroxide varied according to position such that nitroxides near the mouth had greater mobility than those located deep in the cavity. Nitroxides positioned up to 16 carbon units from the hydrocarbon tail of the ligand are motionally restricted and inaccessible, indicating the cavity extends to at least this depth. EPR spectra obtained with and without sucrose showed that the intracavitary position of lauric acid in TL is similar to that in beta-lactoglobulin. However, unlike beta-lactoglobulin, TL binds 16-doxyl stearic acid, suggesting less steric hindrance and greater promiscuity for TL.     

Accessibility of nitroxide side chains: absolute Heisenberg exchange rates from power saturation EPR

In site-directed spin labeling, the relative solvent accessibility of spin-labeled side chains is taken to be proportional to the Heisenberg exchange rate (W(ex)) of the nitroxide with a paramagnetic reagent in solution. In turn, relative values of W(ex) are determined by continuous wave power saturation methods and expressed as a proportional and dimensionless parameter Pi. In the experiments presented here, NiEDDA is characterized as a paramagnetic reagent for solvent accessibility studies, and it is shown that absolute values of W(ex) can be determined from Pi, and that the proportionality constant relating them is independent of the paramagnetic reagent and mobility of the nitroxide. Based on absolute exchange rates, an accessibility factor is defined (0 < rho < 1) that serves as a quantitative measure of side-chain solvent accessibility. The accessibility factors for a nitroxide side chain at 14 different sites in T4 lysozyme are shown to correlate with a structure-based accessibility parameter derived from the crystal structure of the protein. These results provide a useful means for relating crystallographic and site-directed spin labeling data, and hence comparing crystal and solution structures.     

Calcium- and membrane-induced changes in the structure and dynamics of three helical hairpins in annexin B12

Annexins are a family of soluble proteins that can undergo reversible Ca(2+)-dependent interaction with the interfacial region of phospholipid membranes. The helical hairpins on the convex face of the crystal structure of soluble annexins are proposed to mediate binding to membranes, but the mechanism is not defined. For this study, we used a site-directed spin labeling (SDSL) experimental approach to investigate Ca(2+) and membrane-induced structural and dynamic changes that occurred in the helical hairpins encompassing three of the four D and E helices of annexin B12. Electron paramagnetic resonance (EPR) parameters were analyzed for the soluble and Ca(2+)-dependent membrane-bound states of the following nitroxide scans of annexin B12: a continuous 24-residue scan of the D and E helices in the third repeat (residues 219-242) and short scans encompassing the D-E loop regions of the first repeat (residues 68-74) and the fourth repeat (300-305). EPR mobility and accessibility parameters of most sites were similar when the protein was in solution or in the membrane-bound state, and both sets of data were consistent with the crystal structure of the protein. However, membrane-induced changes in mobility and accessibility were observed in all three loop regions, with the most dramatic changes noted at sites corresponding to the highly conserved serine and glycine residues in the loops. EPR accessibility parameters clearly established that nitroxide side chains placed at these sites made direct contact with the bilayer. EPR mobility parameters showed that these sites were very mobile in solution, but immobilized on the EPR time scale in the membrane-bound state. Since the headgroup regions of bilayer phospholipids are relatively mobile in the absence of annexins, Ca(2+)-dependent binding of annexin B12 appears to form a complex in which the mobility of the D-E loop region of the protein and the headgroup region of the phospholipid are highly constrained. Possible biological consequences of annexin-induced restriction of membrane mobility are discussed.     

Structural determinants of nitroxide motion in spin-labeled proteins: solvent-exposed sites in helix B of T4 lysozyme

Site-directed spin labeling provides a means for exploring structure and dynamics in proteins. To interpret the complex EPR spectra that often arise, it is necessary to characterize the rotamers of the spin-labeled side chain and the interactions they make with the local environment in proteins of known structure. For this purpose, crystal structures have been determined for T4 lysozyme bearing a nitroxide side chain (R1) at the solvent-exposed helical sites 41 and 44 in the B helix. These sites are of particular interest in that the corresponding EPR spectra reveal two dynamic states of R1, one of which is relatively immobilized suggesting interactions of the nitroxide with the environment. The crystal structures together with the effect of mutagenesis of nearest neighbors on the motion of R1 suggest intrahelical interactions of 41R1 with the i + 4 residue and of 44R1 with the i + 1 residue. Such interactions appear to be specific to particular rotamers of the R1 side chain.     

Structure and dynamics of a helical hairpin and loop region in annexin 12: a site-directed spin labeling study

A 30-residue nitroxide scan encompassing a helical hairpin and an extended loop in soluble annexin 12 (helices D and E in repeat 2; residues 134-163) has been analyzed in terms of nitroxide side chain mobility and accessibility to collision with paramagnetic reagents (Pi). Values of Pi for both O(2) and a Ni(II) metal complex (NiEDDA) are remarkably well correlated with the fractional solvent accessibility of the native side chains at the corresponding positions computed from the known crystal structure. This result demonstrates the utility of Pi as an experimental measure of side chain accessibility in solution, as well as the lack of structural perturbation due to the presence of the nitroxide side chain. The pattern of side chain mobility is also in excellent agreement with predictions from the crystal structure. The results presented here extend the correlations between mobility and structure described in earlier work on other helical proteins, and suggest their generality. The periodic dependence of Pi and mobility along the sequence of annexin 12 reveals the helical segments and their orientation in the fold, as expected for a nonperturbing nitroxide side chain. However, these data do not distinguish the helix-loop-helix motif from a continuous helix, because immobilized side chains in the short loop sequence maintain the periodicity. As shown here, the ratio of Pi values for O(2) and NiEDDA clearly delineates the loop region, due to size exclusion effects between the two reagents. A new feature evident in a nitroxide scan through multiple secondary elements is a modulation of the basic Pi and mobility patterns along the sequence, apparently due to differences in helix packing and backbone motion. Thus, in the short helix D, residues are consistently more mobile and accessible throughout the sequence compared to the residues in the longer, less-solvated and more ordered helix E.     

Structural determinants of nitroxide motion in spin-labeled proteins: tertiary contact and solvent-inaccessible sites in helix G of T4 lysozyme

A nitroxide side chain (R1) has been substituted at single sites along a helix-turn-helix motif in T4 lysozyme (residues 114-135). Together with previously published data, the new sites reported complete a continuous scan through the motif. Mutants with R1 at sites 115 and 118 were selected for crystallographic analysis to identify the structural origins of the corresponding two-component EPR spectra. At 115, R1 is shown to occupy two rotamers in the room temperature crystal structure, one of which has not been previously reported. The two components in the EPR spectrum apparently arise from differential interactions of the two rotamers with the surrounding structure, the most important of which is a hydrophobic interaction of the nitroxide ring. Interestingly, the crystal structure at 100 K reveals a single rotamer, emphasizing the possibility of rotamer selection in low-temperature crystal structures. Residue 118 is at a solvent-inaccessible site in the protein core, and the structure of 118R1, the first reported for the R1 side chain at a buried site, reveals how the side chain is accommodated in an overpacked core.     

Identification of pH-sensitive regions in the mouse prion by the cysteine-scanning spin-labeling ESR technique

We analyzed the pH-induced mobility changes in moPrP(C) alpha-helix and beta-sheets by cysteine-scanning site-directed spin labeling (SDSL) with ESR. Nine amino acid residues of alpha-helix1 (H1, codon 143-151), four amino acid residues of beta-sheet1 (S1, codon 127-130), and four amino acid residues of beta-sheet2 (S2, codon 160-163) were substituted for by cysteine residues. These recombinant mouse PrP(C) (moPrP(C)) mutants were reacted with a methane thiosulfonate sulfhydryl-specific spin labeling reagent (MTSSL). The 1/deltaH of the central (14N hyperfine) component (M(I) = 0) in the ESR spectrum of spin-labeled moPrP(C) was measured as a mobility parameter of nitroxide residues (R1). The mobilities of E145R1 and Y149R1 at pH 7.4, which was identified as a tertiary contact site by a previous NMR study of moPrP, were lower than those of D143R1, R147R1, and R150R1 reported on the helix surface. Thus, the mobility in the H1 region in the neutral solution was observed with the periodicity associated with a helical structure. On the other hand, the values in the S2 region, known to be located in the buried side, were lower than those in the S1 region located in the surface side. These results indicated that the mobility parameter of the nitroxide label was well correlated with the 3D structure of moPrP. Furthermore, the present study clearly demonstrated three pH-sensitive sites in moPrP, i.e., (1) the N-terminal tertiary contact site of H1, (2) the C-terminal end of H1, and (3) the S2 region. In particular, among these pH-sensitive sites, the N-terminal tertiary contact region of H1 was found to be the most pH-sensitive one and was easily converted to a flexible structure by a slight decrease of pH in the solution. These data provided molecular evidence to explain the cellular mechanism for conversion from PrP(C) to PrP(Sc) in acidic organelles such as the endosome.     

Predicting reactivities of protein surface cysteines as part of a strategy for selective multiple labeling

A variety of biophysical methods used to study proteins requires protein modification using conjugated molecular probes. Cysteine is the main residue that can be modified without the risk of altering other residues in the protein chain. It is possible to label several cysteines in a protein using highly selective labeling reactions, if the cysteines react at very different rates. The reactivity of a cysteine residue introduced into an exposed surface site depends on the fraction of cysteine in the deprotonated state. Here, it is shown that cysteine reactivity differences can be effectively predicted by an electrostatic model that yields site-specifically the fractions of cysteinate. The model accounts for electrostatic interactions between the cysteinyl anion and side chains, the local protein backbone, and water. The energies of interaction with side chains and the main chain are calculated by using the two different dielectric constants, 40 and 22, respectively. Twenty-six mutants of Escherichia coli adenylate kinase were produced, each containing a single cysteine at the protein surface, and the rates of the reaction with 5,5'-dithiobis(2-nitrobenzoic acid) (Ellman's reagent) were measured. Cysteine residues were chosen on the basis of locations that were expected to allow modification of the protein with minimal risk of perturbing its structure. The reaction rates spanned a range of 6 orders of magnitude. The correlation between predicted fractions of cysteinate and measured reaction rates was strong (R = 92%) and especially high (R = 97%) for cysteines at the helix termini. The approach developed here allows reasonably fast, automated screening of protein surfaces to identify sites that permit efficient preparations of double- or triple-labeled protein.