Reductive methylation and pKa determination of the lysine side chains in calbindin D9k

Rotamer Jumps, Proton Exchange, and Amine Inversion Dynamics of Dimethylated Lysine Residues in Proteins Resolved by pH-Dependent 1H and 13C NMR Relaxation Dispersion

Post-translational methylation of lysine side chains is of great importance for protein regulation, including epigenetic control. Here, we present specific ¹³CHD₂ labeling of dimethylated lysines as a sensitive probe of the structure, interactions, and dynamics of these groups, and outline a theoretical and experimental framework for analyzing their conformational dynamics using ¹H and ¹³C CPMG relaxation dispersion experiments. Dimethylated lysine side chains in calcium-loaded calmodulin show a marked pH dependence of their Carr-Purcell-Meiboom-Gill (CPMG) dispersion profiles, indicating complex exchange behavior. Combined analysis of ¹H and ¹³C CPMG relaxation dispersions requires consideration of 12-state correlated exchange of the two methyl groups due to circular three-state rotamer jumps around the Cε-Nζ axis combined with proton exchange and amine inversion. Taking into account a number of fundamental constraints, the exchange model can be reduced to include only three fitted parameters, namely, the geometric average of the rotamer-jump rate constants, the rate constant of deprotonation of Nζ, and the chemical shift difference between the trans and gauge positions of the ¹³C or ¹H nuclei. The pH dependence indicates that protonation of the end group dramatically slows down rotamer exchange for some lysine residues, whereas deprotonation leads to rapid amine inversion coupled with rotamer scrambling. The observed variation among residues in their exchange behavior appears to depend on the structural environment of the side chain. Understanding this type of exchange process is critical to correctly interpreting NMR spectra of methylated lysine side chains. The exchange model presented here forms the basis for studying the structure and dynamics of epigenetically modified lysine side chains and perturbations caused by changes in pH or interactions with target proteins.

Cysteine-ethylation of tissue-extracted membrane proteins as a tool to detect conformational states by solid-state NMR spectroscopy

Solid-state NMR (ssNMR) is an ideal tool to study structure and dynamics of membrane proteins in their native lipid environment. In principle, ssNMR has no size limitations. However, this feature is rarely exploited as large membrane proteins display severe resonance overlap. In addition, dismal yields from recombinant bacterial expression systems limit severely spectroscopic characterization of membrane proteins. For very large mammalian membrane proteins, extraction from the original organism remains the most viable approach. In this case, NMR-observable nuclei must be introduced post-translationally, but the approaches developed so far are rather scarce. Here, we detail the synthesis and engineering of a reactive ¹³C-ethylmethanethiosulfonate (¹³C-EMTS) reagent for the post-translational alkylation of cysteine sidechains of a 110kDa sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) extracted from rabbit skeletal muscle tissue. When reconstituted into liposomes, it is possible to resolve the resonances of the engineered ethyl groups by magic-angle spinning (MAS) 2D [¹³C,¹³C]-DARR experiments. Notably, the ethyl-group modification does not perturb the function of SERCA, yielding well-resolved ¹³C-¹³C fingerprints that are used to image its structural states in the catalytic cycle and filtering out overwhelming naturally-abundant ¹³C nuclei signals arising from the enzyme and lipids. We anticipate that this approach will be used together with ¹⁹F NMR to monitor conformational transitions of enzymes and proteins that are difficult to express recombinantly.

Salvage of failed protein targets by reductive alkylation

The growth of diffraction-quality single crystals is of primary importance in protein X-ray crystallography. Chemical modification of proteins can alter their surface properties and crystallization behavior. The Midwest Center for Structural Genomics (MCSG) has previously reported how reductive methylation of lysine residues in proteins can improve crystallization of unique proteins that initially failed to produce diffraction-quality crystals. Recently, this approach has been expanded to include ethylation and isopropylation in the MCSG protein crystallization pipeline. Applying standard methods, 180 unique proteins were alkylated and screened using standard crystallization procedures. Crystal structures of 12 new proteins were determined, including the first ethylated and the first isopropylated protein structures. In a few cases, the structures of native and methylated or ethylated states were obtained and the impact of reductive alkylation of lysine residues was assessed. Reductive methylation tends to be more efficient and produces the most alkylated protein structures. Structures of methylated proteins typically have higher resolution limits. A number of well-ordered alkylated lysine residues have been identified, which make both intermolecular and intramolecular contacts. The previous report is updated and complemented with the following new data; a description of a detailed alkylation protocol with results, structural features, and roles of alkylated lysine residues in protein crystals. These contribute to improved crystallization properties of some proteins.

Lysine methylation strategies for characterizing protein conformations by NMR

In the presence of formaldehyde and a mild reducing agent, reductive methylation is known to achieve near complete dimethylation of protein amino groups under non-denaturing conditions, in aqueous media. Amino methylation of proteins is employed in mass spectrometric, crystallographic, and NMR studies. Where biosynthetic labeling is prohibitive, amino (13)C-methylation provides an attractive option for monitoring folding, kinetics, protein-protein and protein-DNA interactions by NMR. Here, we demonstrate two improvements over traditional (13)C-reductive methylation schemes: (1) By judicious choice of stoichiometry and pH, ε-aminos can be preferentially monomethylated. Monomethyl tags are less perturbing and generally exhibit improved resolution over dimethyllysines, and (2) By use of deuterated reducing agents and (13)C-formaldehyde, amino groups can be labeled with (13)CH(2)D tags. Use of deutero-(13)C-formaldehyde affords either (13)CHD(2), or (13)CD(3) probes depending on choice of reducing agent. Making use of (13)C-(2)H scalar couplings, we demonstrate a filtering scheme that eliminates natural abundance (13)C signal.

New surface contacts formed upon reductive lysine methylation: improving the probability of protein crystallization

Surface lysine methylation (SLM) is a technique for improving the rate of success of protein crystallization by chemically methylating lysine residues. The exact mechanism by which SLM enhances crystallization is still not clear. To study these mechanisms, and to analyze the conditions where SLM will provide the optimal benefits for rescuing failed crystallization experiments, we compared 40 protein structures containing N,N-dimethyl-lysine (dmLys) to a nonredundant set of 18,972 nonmethylated structures from the PDB. By measuring the relative frequency of intermolecular contacts (where contacts are defined as interactions between the residues in proximity with a distance of 3.5 A or less) of basic residues in the methylated versus nonmethylated sets, dmLys-Glu contacts are seen more frequently than Lys-Glu contacts. Based on observation of the 10 proteins with both native and methylated structures, we propose that the increased rate of contact for dmLys-Glu is due to both a slight increase in the number of amine-carboxyl H-bonds and to the formation of methyl C--H...O interactions. By comparing the relative contact frequencies of dmLys with other residues, the mechanism by which methylation of lysines improves the formation of crystal contacts appears to be similar to that of Lys to Arg mutation. Moreover, analysis of methylated structures with the surface entropy reduction (SER) prediction server suggests that in many cases SLM of predicted SER sites may contribute to improved crystallization. Thus, tools that analyze protein sequences and mark residues for SER mutation may identify proteins with good candidate sites for SLM.

Residue-specific pKa determination of lysine and arginine side chains by indirect 15N and 13C NMR spectroscopy: application to apo calmodulin

Electrostatic interactions in proteins can be probed experimentally through determination of residue-specific acidity constants. We describe here triple-resonance NMR techniques for direct determination of lysine and arginine side-chain protonation states in proteins. The experiments are based on detection of nonexchangeable protons over the full range of pH and temperature and therefore are well suited for pKa determination of individual amino acid side chains. The experiments follow the side-chain 15Nzeta (lysine) and 15Nepsilon or 13Czeta (arginine) chemical shift, which changes due to sizable changes in the heteronuclear electron distribution upon (de)protonation. Since heteronuclear chemical shifts are overwhelmed by the charge state of the amino acid side chain itself, these methods supersede 1H-based NMR in terms of accuracy, sensitivity, and selectivity. Moreover, the 15Nzeta and 15Nepsilon nuclei may be used to probe changes in the local electrostatic environment. Applications to three proteins are described: apo calmodulin, calbindin D9k, and FKBP12. For apo calmodulin, residue-specific pKa values of lysine side chains were determined to fall between 10.7 and 11.2 as a result of the high net negative charge on the protein surface. Ideal two-state titration behavior observed for all lysines indicates the absence of significant direct charge interactions between the basic residues. These results are compared with earlier studies based on chemical modification.

Pharmacokinetic characterization of 14C-vascular endothelial growth factor controlled release microspheres using a rat model

The objectives of this study were to characterize the pharmacokinetics of vascular endothelial growth factor (VEGF) in poly(lactic-co-glycolic) acid (PLGA) microspheres using a rat model, and to develop a pharmacokinetic model for this controlled release formulation. 14C-VEGF was encapsulated using a solid-in-oil-in-water emulsification method. The microspheres were administered subcutaneously to rats and the pharmacokinetic parameters were compared with those of protein solutions. Intravenous administration of protein solutions resulted in short half-lives and subcutaneous administration resulted in rapid clearance from the subcutaneous tissue, with high plasma concentrations as expressed by rapid absorption and elimination. The subcutaneous administration of the VEGF microspheres produced low plasma concentrations and high subcutaneous concentrations over a period of 7 weeks. The area under the curve (AUC), the time required to achieve the maximum concentration (tmax), the maximum concentration (Cmax) in blood samples and the elimination rate constant (kel) values at the subcutaneous tissue site were selected to compare the pharmacokinetic characterization of VEGF microspheres with that of protein solutions. The in-vivo release profiles of the proteins were slower than the in-vitro release profiles and they followed the same trend as the in-vitro and in-vivo PLGA degradation rates. The PLGA microsphere degradation was the determinant step for VEGF release from the microspheres and its absorption at the subcutaneous site. Microspheres appear to be an attractive system for the localized rate-controlled delivery of VEGF. 14C-Methylation via reductive alkylation of VEGF did not affect its mitogenic activity, however approximately 25% activity was lost following release from PLGA microspheres. This loss of activity may be due to degradation in an acidic environment as a result of PLGA degradation.

Chemical modification of a variant of human MIP-1alpha; implications for dimer structure

A sequence variant of human MIP-1alpha, in which Asp26 has been replaced by Al alpha, has been chemically modified by the addition of 13C-labeled methyl groups at each of the lysine residues and the N-terminus. The sites of methylation have been verified by a combination of MALDI-TOF mass spectrometric experiments and tryptic digestion followed by N-terminal mapping. The effect of the modification on the structure and activity of the protein have been determined by analytical ultra-centrifugation, 13C NMR spectroscopy and receptor binding studies. The results of these experiments suggest that huMIP-alpha D26A (BB10010), when present as a dimer, adopts a globular structure, like MCP-3, rather than the elongated or cylindrical structure determined for dimers of huMIP-1beta and RANTES.

NMR study of the sites of human hemoglobin acetylated by aspirin

Acetylation of hemoglobin by aspirin and other acetylating agents has been used to generate hemoglobin analogs with altered structural and functional properties, and may prove useful in the treatment of sickle cell disease. We have studied the acetylation of human hemoglobin using [1'-(13)C]acetylsalicylic acid in combination with two-dimensional HMQC and HSQC NMR analysis. The spectra of the acetylated hemoglobin exhibit a number of well resolved resonances. Several spectral assignment strategies were used: blocking the 2, 3-DPG binding site non-covalently with inositol hexaphosphate or covalently with a cross-linking agent, selective carbamylation of the N-terminal valine amino groups with cyanate, spin-labeling the hemoglobin at betaCys93, and analysis of a hemoglobin triple mutant: betaV1MH2DeltaK144R, in which betaLys144 is replaced by an arginine residue. These studies support the conclusion that the most rapidly acetylated residue is betaLys82 rather than betaLys144, as previously reported. Further, it is apparent that acetyl betaLys82 can give rise to several resonances due to additional acetylation of betaLys82' or other nearby residues. An additional assignment strategy involving comparison of the chemical shifts of the acetyl resonances observed for adducts of diamagnetic carbonmonoxyhemoglobin with the shifts observed in paramagnetic cyanomethemoglobin provides information about the location of the acetyl derivatives relative to the heme irons. This approach is limited, however, by the lack of well defined structural information for the lysine residues on the protein surface. Additional tentative assignments have also been made, using the above approaches.

An NMR analysis of the reaction of ubiquitin with [acetyl-1-13C]aspirin

The acetylation of ubiquitin by [acetyl-1-13C]aspirin has been studied using 2D NMR methods. Studies performed in a 50:50 H2O:D2O medium show doubling of the acetyl carbonyl resonances, indicating that all of the stable adducts formed involved amide linkages. Assignment of the heteronuclear multiple quantum coherence (HMQC) resonances was accomplished based on comparison of resonance intensities with the results of an Edman degradation analysis, pH titration studies of acetylated ubiquitin, and analysis of two ubiquitin mutants, K33R and K63R. The presence of a single tyrosine residue in close proximity to lysine-48 suggested another assignment strategy. Nitration of tyrosine-59 resulted in a small, pH-dependent shift of the resonance assigned to lysine-48, with a pK of 7.0, close to that expected for the nitrotyrosyl hydroxyl group. An additional adduct resonance with very low intensity also was observed and tentatively assigned to the acetylated N-terminal methionine residue. The relative rates of acetylation of the various lysine residues were obtained from time-dependent HMQC studies. Since no sample preparation artifacts were introduced, the levels of modification of the various residues could be determined with relatively high accuracy. Based on the time-dependent intensity data, the relative rate constants for modification of K6, K48, K63, K11, K33, and M1 were 1.0, 0.59, 0.43, 0.26, 0.23, and 0.03, respectively. These results were in much better agreement with amino accessibility predictions based on the crystal structure of the ubiquitin monomer than with predictions based on the ubiquitin structure in the crystallized dimeric and tetrameric forms. This approach provides a useful basis for understanding how local environmental factors can influence protein adduct formation, as well as for comparing the extent and specificity of various acetylation reagents.

N epsilon,N epsilon-dimethyl-lysine cytochrome c as an NMR probe for lysine involvement in protein-protein complex formation

The reductively dimethylated derivatives of horse and yeast iso-1-ferricytochromes c have been prepared and characterized for use as NMR probes of the complexes formed by cytochrome c with bovine liver cytochrome b5 and yeast cytochrome c peroxidase. The electrostatic properties and structures of the derivatized cytochromes are not significantly perturbed by the modifications; neither are the electrostatics of protein-protein complex formation or rates of interprotein electron transfer. Two-dimensional 1H-13C NMR spectroscopy of the complexes formed by the derivatized cytochromes with cytochrome b5 and cytochrome c peroxidase has been used to investigate the number and identity of lysine residues of cytochrome c that are involved in interprotein interactions of cytochrome c. The NMR data are incompatible with simple static models proposed previously for the complexes formed by these proteins, but are consistent with the presence of multiple, interconverting complexes of comparable stability, consistent with studies employing Brownian dynamics to model the complexes. The NMR characteristics of the Nepsilon,Nepsilon-dimethyl-lysine groups, their chemical shift dispersion, oxidation state and temperature dependences and the possibility of chemical exchange phenomena are discussed with relevance to the utility of Nepsilon, Nepsilon-dimethyl-lysine's being a generally useful derivative for characterizing protein-protein complexes.

Chemical modification of lyophilized proteins in nonaqueous environments

Lyophilized proteins were reacted in vacuo with a volatile reagent or dispersed in octane and reacted with dissolved reagent. Three novel derivatives were formed with iodomethane: (a) quaternized trimethyl amino groups, (b) N1,N3-dimethylimidazolium cation, and (c) phenolic O-methyl ether. Acid anhydrides acylated amino groups and formed mixed anhydrides with side-chain carboxyl groups. Under nonaqueous conditions it was observed that: (i) The same derivatives are formed as under aqueous conditions. (ii) Hydrolytic breakdown of protein is prevented. (iii) Less reagent is required. (iv) Unreacted reagent can be recovered. (v) Water-labile derivatives can be isolated as stable intermediates. (vi) The yield of a derivatized functional group was directly related to its pK(a), its surface exposure, and the pH of the solution from which the protein was lyophilized. (vii) The physicochemical factors governing the reactivity of protein functional groups in nonaqueous environments appear to reflect the protein solution structure prior to lyophilization.