Unambiguous determination of the ionization state of a glycoside hydrolase active site lysine by 1H-15N heteronuclear correlation spectroscopy

Accurately Predicting Protein pKa Values Using Nonequilibrium Alchemy

The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.

Weighted AnX multiplets

That the NMR transition of a spin-1/2 nucleus is split into n evenly spaced lines by indirect dipole-dipole (J) coupling to n magnetically equivalent nuclei, whose successive amplitudes follow the n^th row of Pascal's triangle, is an elementary result in NMR. Described here are a family of less well known multiplet structures with different amplitudes for the evenly spaced lines. The amplitudes can be derived from the n^th row of Pascal's triangle by weighting the corresponding value of the row by z or z², where z is related to the summed magnetic quantum number of the magnetically equivalent spins. z¹-multiplets have been described in INEPT experiments. A z²-multiplet can be indirectly observed in HSQC experiments when the decoupling pulse during t₁ is removed, i.e., an F1-coupled HSQC. While not difficult to generate and despite some reported usefulness, to the best of our knowledge, z²-multiplets have not been rigorously described in previous literature.

Conformational features and ionization states of Lys side chains in a protein studied using the stereo-array isotope labeling (SAIL) method

Although both the hydrophobic aliphatic chain and hydrophilic ζ -amino group of the Lys side chain presumably contribute to the structures and functions of proteins, the dual nature of the Lys residue has not been fully investigated using NMR spectroscopy, due to the lack of appropriate methods to acquire comprehensive information on its long consecutive methylene chain. We describe herein a robust strategy to address the current situation, using various isotope-aided NMR technologies. The feasibility of our approach is demonstrated for the Δ + PHS/V66K variant of staphylococcal nuclease (SNase), which contains 21 Lys residues, including the engineered Lys-66 with an unusually low pK a of ∼  5.6. All of the NMR signals for the 21 Lys residues were sequentially and stereospecifically assigned using the stereo-array isotope-labeled Lys (SAIL-Lys), [U-13 C,15 N; β 2 ,γ 2 ,δ 2 ,ε 3 -D4 ]-Lys. The complete set of assigned 1 H, 13 C, and 15 N NMR signals for the Lys side-chain moieties affords useful structural information. For example, the set includes the characteristic chemical shifts for the 13 Cδ , 13 Cε , and 15 Nζ signals for Lys-66, which has the deprotonated ζ -amino group, and the large upfield shifts for the 1 H and 13 C signals for the Lys-9, Lys-28, Lys-84, Lys-110, and Lys-133 side chains, which are indicative of nearby aromatic rings. The 13 Cε and 15 Nζ chemical shifts of the SNase variant selectively labeled with either [ε -13 C;ε ,ε -D2 ]-Lys or SAIL-Lys, dissolved in H2 O and D2 O, showed that the deuterium-induced shifts for Lys-66 were substantially different from those of the other 20 Lys residues. Namely, the deuterium-induced shifts of the 13 Cε and 15 Nζ signals depend on the ionization states of the ζ -amino group, i.e., - 0.32 ppm for Δ δ 13 Cε [Nζ D3 + -Nζ H3 + ] vs. - 0.21 ppm for Δ δ 13 Cε [Nζ D2 -Nζ H2 ] and - 1.1 ppm for Δ δ 15 Nζ [Nζ D3 + -Nζ H3 + ] vs. - 1.8 ppm for Δ δ 15 Nζ [Nζ D2 -Nζ H2 ]. Since the 1D 13 C NMR spectrum of a protein selectively labeled with [ε -13 C;ε ,ε -D2 ]-Lys shows narrow (>  2 Hz) and well-dispersed 13 C signals, the deuterium-induced shift difference of 0.11 ppm for the protonated and deprotonated ζ -amino groups, which corresponds to 16.5 Hz at a field strength of 14 T (150 MHz for 13 C), could be accurately measured. Although the isotope shift difference itself may not be absolutely decisive to distinguish the ionization state of the ζ -amino group, the 13 Cδ , 13 Cε , and 15 Nζ signals for a Lys residue with a deprotonated ζ -amino group are likely to exhibit distinctive chemical shifts as compared to the normal residues with protonated ζ -amino groups. Therefore, the isotope shifts would provide a useful auxiliary index for identifying Lys residues with deprotonated ζ -amino groups at physiological pH levels.

A platform for chemical modification of mandelate racemase: characterization of the C92S/C264S and γ-thialysine 166 variants

Mandelate racemase (MR) serves as a paradigm for our understanding of enzyme-catalyzed deprotonation of a carbon acid substrate. To facilitate structure-function studies on MR using non-natural amino acid substitutions, we engineered the Cys92Ser/Cys264Ser variant (dmMR) as a platform for introducing Cys residues at specific locations for subsequent covalent modification. While the highly reactive thiol of Cys furnishes a site for chemical modification, site-specificity requires that other Cys residues be non-reactive or replaced by a non-reactive amino acid, especially if chemical modification is conducted under denaturing conditions. The catalytic efficiency of dmMR is reduced only ~2-fold relative to wild-type MR, making dmMR a viable platform for the site-specific introduction of Cys. As an example, the inactive Lys166Cys variant of dmMR was treated with ethylenimine under denaturing conditions to replace the Brønsted acid-base catalyst Lys 166 with the non-natural amino acid γ-thialysine. Comparison of the pH-activity profiles of dmMR and the active γ-thialysine variant revealed a reduction in the pKa for the side chain amino group of ~0.4 units for the latter variant. Unlike wild-type MR for which diffusion is partially rate-limiting, dmMR and the γ-thialysine variant showed no dependence on the solvent viscosity suggesting that the chemical step is fully rate-limiting.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH₃ groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

Anomalous Properties of Lys Residues Buried in the Hydrophobic Interior of a Protein Revealed with 15N-Detect NMR Spectroscopy

Ionizable residues buried in hydrophobic environments in proteins are essential for many fundamental biochemical processes. These residues titrate with anomalous pK_a values that are challenging to reproduce with structure-based calculations owing to the conformational reorganization coupled to their ionization. Detailed characterization of this conformational reorganization is of interest; unfortunately, the properties of buried Lys residues are difficult to study experimentally. Here we demonstrate the utility of ¹⁵N NMR spectroscopy to gain insight into the protonation state, state of hydration and conformational dynamics of the Nζ amino group of buried Lys residues. The experiments were applied to five variants of staphylococcal nuclease, with internal Lys residues that titrate with pK_a values ranging from 6.2 to 8.1. Direct detection of buried Lys residues with these NMR spectroscopy methods will enable correlation between thermodynamic and structural data as well as unprecedented examination of how conformational transitions coupled to their ionization affect their pK_a values.

Direct detection of lysine side chain NH3+ in protein-heparin complexes using NMR spectroscopy

Two NMR observables, the N_ζH₃⁺ peak in the HISQC spectrum and N_ζ chemical shift difference between the free and heparin-bound forms, can identify binding-interface lysines in protein-heparin complexes. Unlike backbone chemical shifts, these direct probes are stringent and are less prone to either false positives or false negatives.

Disulfide bonds and disorder in granulin-3: An unusual handshake between structural stability and plasticity

Granulins (GRNs) are a family of small (∼6 kDa) proteins generated by the proteolytic processing of their precursor, progranulin (PGRN), in many cell types. Both PGRN and GRNs are implicated in a plethora of biological functions, often in opposing roles to each other. Lately, GRNs have generated significant attention due to their implicated roles in neurodegenerative disorders. Despite their physiological and pathological significance, the structure-function relationships of GRNs are poorly defined. GRNs contain 12 conserved cysteines forming six intramolecular disulfide bonds, making them rather exceptional, even among a few proteins with high disulfide bond density. Solution NMR investigations in the past have revealed a unique structure containing putative interdigitated disulfide bonds for several GRNs, but GRN-3 was unsolvable due to its heterogeneity and disorder. In our previous report, we showed that abrogation of disulfide bonds in GRN-3 renders the protein completely disordered (Ghag et al., Prot Eng Des Sel 2016). In this study, we report the cellular expression and biophysical analysis of fully oxidized, native GRN-3. Our results indicate that both E. coli and human embryonic kidney (HEK) cells do not exclusively make GRN-3 with homogenous disulfide bonds, likely due to the high cysteine density within the protein. Biophysical analysis suggests that GRN-3 structure is dominated by irregular loops held together only by disulfide bonds, which induced remarkable thermal stability to the protein despite the lack of regular secondary structure. This unusual handshake between disulfide bonds and disorder within GRN-3 could suggest a unique adaptation of intrinsically disordered proteins towards structural stability.

A Unique and Simple Approach to Improve Sensitivity in 15N-NMR Relaxation Measurements for NH₃⁺ Groups: Application to a Protein-DNA Complex

NMR spectroscopy is a powerful tool for research on protein dynamics. In the past decade, there has been significant progress in the development of NMR methods for studying charged side chains. In particular, NMR methods for lysine side-chain NH₃⁺ groups have been proven to be powerful for investigating the dynamics of hydrogen bonds or ion pairs that play important roles in biological processes. However, relatively low sensitivity has been a major practical issue in NMR experiments on NH₃⁺ groups. In this paper, we present a unique and simple approach to improve sensitivity in ¹⁵N relaxation measurements for NH₃⁺ groups. In this approach, the efficiency of coherence transfers for the desired components are maximized, whereas undesired anti-phase or multi-spin order components are purged through pulse schemes and rapid relaxation. For lysine side-chain NH₃⁺ groups of a protein-DNA complex, we compared the data obtained with the previous and new pulse sequences under the same conditions and confirmed that the ¹⁵N relaxation parameters were consistent for these datasets. While retaining accuracy in measuring ¹⁵N relaxation, our new pulse sequences for NH₃⁺ groups allowed an 82% increase in detection sensitivity of ¹⁵N longitudinal and transverse relaxation measurements.

Unambiguous Determination of Protein Arginine Ionization States in Solution by NMR Spectroscopy

Because arginine residues in proteins are expected to be in their protonated form almost without exception, reports demonstrating that a protein arginine residue is charge-neutral are rare and potentially controversial. Herein, we present a ¹³ C-detected NMR experiment for probing individual arginine residues in proteins notwithstanding the presence of chemical and conformational exchange effects. In the experiment, the ¹⁵ N^η and ¹⁵ N^ϵ chemical shifts of an arginine head group are correlated with that of the directly attached ¹³ C^ζ . In the resulting spectrum, the number of protons in the arginine head group can be obtained directly from the ¹⁵ N-¹ H scalar coupling splitting pattern. We applied this method to unambiguously determine the ionization state of the R52 side chain in the photoactive yellow protein from Halorhodospira halophila. Although only three H^η atoms were previously identified by neutron crystallography, we show that R52 is predominantly protonated in solution.

Self-Assembled Peptide-Based System for Mitochondrial-Targeted Gene Delivery: Functional and Structural Insights

Human mitochondrial dysfunction can lead to severe and often deadly diseases, for which there are no known cures. Although the targeted delivery of therapeutic gene to mitochondria is a promising approach to alleviate these disorders, gene carrier systems for the selective delivery of functional DNA into the mitochondria of living mammalian cells are currently unavailable. Here we rationally developed dual-domain peptides containing DNA-condensing/cell-penetrating/endosome-disruptive and mitochondria-targeting sequences. Secondary structures of the dual-domain peptides were analyzed, and variations in the physicochemical properties (stability, size, and ζ potential) of peptide/DNA complexes were studied as a function of peptide-to-DNA ratio and serum addition. An optimized formulation, identified through qualitative and quantitative studies, fulfills the fundamental prerequisites for mitochondria-specific DNA delivery, successfully transfecting a high proportion (82 ± 2%) of mitochondria in a human cell line with concomitant biocompatibility. Nuclear magnetic resonance studies confirmed the effectiveness of our bipartite peptide design with segregated functions: a helical domain necessary for mitochondrial import and an unstructured region for interaction with DNA involving lysine residues. Further analyses revealed that the lysine-specific interaction assisted the self-organization of the peptide and the DNA cargo, leading to a structural arrangement within the formed complex that is crucial for its biological efficiency. Thus the reported gene vector represents a new and reliable tool to uncover the complexity of mitochondrial transfection.

NMR studies on the dynamics of hydrogen bonds and ion pairs involving lysine side chains of proteins

Hydrogen bonds and ion pairs involving side chains play vital roles in protein functions such as molecular recognition and catalysis. Despite the wealth of structural information about hydrogen bonds and ion pairs at functionally crucial sites on proteins, the dynamics of these fundamental chemical interactions are not well understood largely due to the lack of suitable experimental tools in the past. NMR spectroscopy is a powerful tool for investigations of protein dynamics, but the vast majority of NMR methods had been applicable only to the backbone or methyl groups. Recently, a substantial progress has been made in the research on the dynamics of hydrogen bonds and ion pairs involving lysine side-chain NH3+ groups. Together with computational/theoretical approaches, the new NMR methods provide unique insights into the dynamics of hydrogen bonds and ion pairs involving lysine side chains. Here, the methodology and its applications are reviewed.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

Catalytic roles of βLys87 in tryptophan synthase: (15)N solid state NMR studies

The proposed mechanism for tryptophan synthase shows βLys87 playing multiple catalytic roles: it bonds to the PLP cofactor, activates C4' for nucleophilic attack via a protonated Schiff base nitrogen, and abstracts and returns protons to PLP-bound substrates (i.e. acid-base catalysis). ε-¹⁵N-lysine TS was prepared to access the protonation state of βLys87 using ¹⁵N solid-state nuclear magnetic resonance (SSNMR) spectroscopy for three quasi-stable intermediates along the reaction pathway. These experiments establish that the protonation state of the ε-amino group switches between protonated and neutral states as the β-site undergoes conversion from one intermediate to the next during catalysis, corresponding to mechanistic steps where this lysine residue has been anticipated to play alternating acid and base catalytic roles that help steer reaction specificity in tryptophan synthase catalysis. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications. Guest Editors: Andrea Mozzarelli and Loredano Pollegioni.

pH-dependent random coil (1)H, (13)C, and (15)N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements

The pK a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the (1)H, (13)C, and (15)N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain (15)N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.

Effective strategy to assign ¹H- ¹⁵N heteronuclear correlation NMR signals from lysine side-chain NH3₃⁺ groups of proteins at low temperature

Recent studies have shown that lysine side-chain NH3(+) groups are excellent probes for NMR investigations of dynamics involving hydrogen bonds and ion pairs relevant to protein function. However, due to rapid hydrogen exchange, observation of (1)H-(15)N NMR cross peaks from lysine NH3(+) groups often requires use of a relatively low temperature, which renders difficulty in resonance assignment. Here we present an effective strategy to assign (1)H and (15)N resonances of NH3(+) groups at low temperatures. This strategy involves two new (1)H/(13)C/(15)N triple-resonance experiments for lysine side chains. Application to a protein-DNA complex is demonstrated.

Use of H/D isotope effects to gather information about hydrogen bonding and hydrogen exchange rates

Polar side-chains in proteins play important roles in forming and maintaining three-dimensional structures, and thus participate in various biological functions. Until recently, most protein NMR studies have focused on the non-exchangeable protons of amino acid residues. The exchangeable protons attached to polar groups, such as hydroxyl (OH), sulfhydryl (SH), and amino (NH2) groups, have mostly been ignored, because in many cases these hydrogen atoms exchange too quickly with water protons, making NMR observations impractical. However, in certain environments, such as deep within the hydrophobic interior of a protein, or in a strong hydrogen bond to other polar groups or interacting ligands, the protons attached to polar groups may exhibit slow hydrogen exchange rates and thus become NMR accessible. To explore the structural and biological implications of the interactions involving polar side-chains, we have developed versatile NMR methods to detect such cases by observing the line shapes of (13)C NMR signals near the polar groups, which are affected by deuterium-proton isotope shifts in a mixture of H2O and D2O. These methods allow the detection of polar side-chains with slow hydrogen-deuterium exchange rates, and therefore provide opportunities to retrieve information about the polar side-chains, which might otherwise be overlooked by conventional NMR experiments. Future prospects of applications using deuterium-proton isotope shifts to retrieve missing structural and dynamic information of proteins are discussed.

Solution and gas-phase H/D exchange of protein-small-molecule complexes: Cex and its inhibitors

The properties of noncovalent complexes of the enzyme exo-1,4-β-D-glycanase ("Cex") with three aza-sugar inhibitors, deoxynojirimycin (X(2)DNJ), isofagomine lactam (X(2)IL), and isofagomine (X(2)IF), have been studied with solution and gas-phase hydrogen deuterium exchange (H/Dx) and measurements of collision cross sections of gas-phase ions. In solution, complexes have lower H/Dx levels than free Cex because binding the inhibitors blocks some sites from H/Dx and reduces fluctuations of the protein. In mass spectra of complexes, abundant Cex ions are seen, which mostly are formed by dissociation of complexes in the ion sampling interface. Both complex ions and Cex ions formed from a solution containing complexes have lower cross sections than Cex ions from a solution of Cex alone. This suggests the Cex ions formed by dissociation "remember" their solution conformations. For a given charge, ions of the complexes have greater gas-phase H/Dx levels than ions of Cex. Unlike cross sections, H/Dx levels of the complexes do not correlate with the relative gas-phase binding strengths measured by MS/MS. Cex ions from solutions with or without inhibitors, which have different cross sections, show the same H/Dx level after 15 s, indicating the ions may fold or unfold on the seconds time scale of the H/Dx experiment. Thus, cross sections show that complexes have more compact conformations than free protein ions on the time scale of ca. 1 ms. The gas-phase H/Dx measurements show that at least some complexes retain different conformations from the Cex ions on a time scale of seconds.

Structure, dynamics, and ionization equilibria of the tyrosine residues in Bacillus circulans xylanase

We have developed NMR spectroscopic methods to investigate the tyrosines within Bacillus circulans xylanase (BcX). Four slowly exchanging buried tyrosine hydroxyl protons with chemical shifts between 7.5 and 12.5 ppm were found using a long-range (13)C-HSQC experiment that exploits the (3)J(CH) coupling between the ring (1)H(η) and (13)C(ε) nuclei. The NMR signals from these protons were assigned via (13)C-tyrosine selective labelling and a suite of scalar and (13)C,(15)N-filtered/edited NOE correlation spectra. Of the fifteen tyrosines in BcX, only the buried Tyr79 and Tyr105 showed four distinct, rather than two averaged, signals from ring (13)C-(1)H pairs, indicative of slow flipping on the chemical shift timescale. Ring flipping rate constants of ~10 and ~0.2 s(-1) were measured for the two residues, respectively, using a (13)C longitudinal exchange experiment. The hydrogen bonding properties of the Tyr79 and Tyr105 hydroxyls were also defined by complementary NOE and J-coupling measurements. The (1)H(η) hydrogen-deuterium exchange rate constants of the buried tyrosines were determined from (13)C/(15)N-filtered spectra recorded as a function of pH. These exchange rate constants correspond to estimated protection factors of ~10(4)-10(8) relative to a random coil tyrosine. The phenolic sidechain pK (a) values were also measured by monitoring their pH-dependent (13)C(ζ) chemical shifts via (1)H(ε/δ)((13)C(ε))(13)C(ζ) correlation spectra. Exposed tyrosines had unperturbed pK (a) values of ~10.2, whereas buried residues remained predominantly neutral at or even above pH 11. Combined with selective isotope labelling, these NMR experiments should prove useful for investigating the structural and electrostatic properties of tyrosines in many interesting proteins.

γ-Hydroxy-1,N2-propano-2'-deoxyguanosine DNA adduct conjugates the N-terminal amine of the KWKK peptide via a carbinolamine linkage

The γ-hydroxy-1,N(2)-propano-2'-deoxyguanosine adduct (γ-OH-PdG) was introduced into 5'-d(GCTAGCXAGTCC)-3'·5'-d(GGACTCGCTAGC)-3' (X = γ-OH-PdG). In the presence of excess peptide KWKK, (13)C isotope-edited NMR revealed the formation of two spectroscopically distinct DNA-KWKK conjugates. These involved the reaction of the KWKK N-terminal amino group with the N(2)-dG propylaldehyde tautomer of the γ-OH-PdG lesion. The guanine N1 base imino resonance at the site of conjugation was observed in isotope-edited (15)N NMR experiments, suggesting that the conjugated guanine was inserted into the duplex and that the guanine imino proton was protected from exchange with water. The conjugates could be reduced in the presence of NaCNBH(3), suggesting that they existed, in part, as imine (Schiff base) linkages. However, (13)C isotope-edited NMR failed to detect the imine linkages, suggesting that these KWKK conjugates existed predominantly as diastereomeric carbinolamines, in equilibrium with trace amounts of the imines. The structures of the diastereomeric DNA-KWKK conjugates were predicted from potential energy minimization of model structures derived from the refined structure of the fully reduced cross-link [ Huang, H., Kozekov, I. D., Kozekova, A., Rizzo, C. J., McCullough, A., Lloyd, R. S., and Stone, M. P. ( 2010 ) Biochemistry , 49 , 6155 -6164 ]. Molecular dynamics calculations carried out in explicit solvent suggested that the conjugate bearing the S-carbinolamine linkage was the major species due to its potential for intramolecular hydrogen bonding. These carbinolamine DNA-KWKK conjugates thermally stabilized duplex DNA. However, the DNA-KWKK conjugates were chemically reversible and dissociated when the DNA was denatured. In this 5'-CpX-3' sequence, the DNA-KWKK conjugates slowly converted to interstrand N(2)-dG:N(2)-dG DNA cross-links and ring-opened γ-OH-PdG derivatives over a period of weeks.

Dynamics of lysine side-chain amino groups in a protein studied by heteronuclear 1H−15N NMR spectroscopy

Despite their importance in macromolecular interactions and functions, the dynamics of lysine side-chain amino groups in proteins are not well understood. In this study, we have developed the methodology for the investigations of the dynamics of lysine NH3(+) groups by NMR spectroscopy and computation. By using 1H−15N heteronuclear correlation experiments optimized for 15NH3(+) moieties, we have analyzed the dynamic behavior of individual lysine NH3(+) groups in human ubiquitin at 2 °C and pH 5. We modified the theoretical framework developed previously for CH3 groups and used it to analyze 15N relaxation data for the NH3(+) groups. For six lysine NH3(+) groups out of seven in ubiquitin, we have determined model-free order parameters, correlation times for bond rotation, and reorientation of the symmetry axis occurring on a pico- to nanosecond time scale. From CPMG relaxation dispersion experiment for lysine NH3(+) groups, slower dynamics occurring on a millisecond time scale have also been detected for Lys27. The NH3(+) groups of Lys48, which plays a key role as the linkage site in ubiquitination for proteasomal degradation, was found to be highly mobile with the lowest order parameter among the six NH3(+) groups analyzed by NMR. We compared the experimental order parameters for the lysine NH3(+) groups with those from a 1 μs molecular dynamics simulation in explicit solvent and found good agreement between the two. Furthermore, both the computer simulation and the experimental correlation times for the bond rotations of NH3(+) groups suggest that their hydrogen bonding is highly dynamic with a subnanosecond lifetime. This study demonstrates the utility of combining NMR experiment and simulation for an in-depth characterization of the dynamics of these functionally most important side-chains of ubiquitin.

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by ¹H NMR (Bartik et al., Biophys J 1994;66:1180–1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5–1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.

Lysine and arginine side chains in glycosaminoglycan-protein complexes investigated by NMR, cross-linking, and mass spectrometry: a case study of the factor H-heparin interaction

We have used the interaction between module 7 of complement factor H (CFH approximately 7) and a fully sulfated heparin tetrasaccharide to exemplify a new approach for studying contributions of basic side chains to the formation of glycosaminoglycan (GAG)-protein complexes. We first employed HISQC and H(2)CN NMR experiments to monitor the side-chain resonances of lysines and arginines in (15)N, (13)C-labeled protein during titrations with a fully sulfated heparin tetrasaccharide under physiological conditions. Under identical conditions and using (15)N-labeled protein, we then cross-linked tetrasaccharide to CFH approximately 7 and confirmed the 1:1 stoichiometry by FT-ICR-MS. We subsequently characterized this covalent protein-GAG conjugate by NMR and further MS techniques. MALDI-TOF MS identified protein fragments obtained via trypsin digestion or chemical fragmentation, yielding information concerning the site of GAG attachment. Combining MS and NMR data allowed us to identify the side chain of K405 as the point of attachment of the cross-linked heparin oligosaccharide to CFH approximately 7. On the basis of the analysis of NMR and MS data of the noncovalent and cross-linked CFH approximately 7-tetrasaccharide complexes, we conclude that the K446 side chain is not essential for binding the tetrasaccharide, despite the large chemical shift perturbations of its backbone amide (15)N and (1)H resonances during titrations. We show that R444 provides the most important charge-charge interaction within a C-terminal heparin-binding subsite of CFH approximately 7 whereas side chains of R404, K405, and K388 are the predominant contributors to an N-terminal binding subsite located in the immediate vicinity of residue 402, which is implicated in age-related macular degeneration (AMD).

Titration_DB: storage and analysis of NMR-monitored protein pH titration curves

NMR-monitored pH titration experiments are routinely used to measure site-specific protein pKa values. Accurate experimental pKa values are essential in dissecting enzyme catalysis, in studying the pH-dependence of protein stability and ligand binding, in benchmarking pKa prediction algorithms, and ultimately in understanding electrostatic effects in proteins. However, due to the complex ways in which pH-dependent electrostatic and structural changes manifest themselves in NMR spectra, reported apparent pKa values are often dependent on the way that NMR pH-titration curves are analyzed. It is therefore important to retain the raw NMR spectroscopic data to allow for documentation and possible re-interpretation. We have constructed a database of primary NMR pH-titration data, which is accessible via a web interface. Here, we report statistics of the database contents and analyze the data with a global perspective to provide guidelines on best practice for fitting NMR titration curves. Titration_DB is available at http://enzyme.ucd.ie/Titration_DB. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Observing in-phase single-quantum 15N multiplets for NH2/NH3+ groups with two-dimensional heteronuclear correlation spectroscopy

Two-dimensional (2D) F1-(1)H-coupled HSQC experiments provide 3:1:1:3 and 1:0:1 multiplets for AX(3) and AX(2) spin systems, respectively. These multiplets occur because, in addition to the 2S(y)H(z)(a)-->2S(y)H(z)(a) process, the coherence transfers such as 2S(y)H(z)(a)-->2S(y)H(z)(b) occurring in t(1) period provide detectable magnetization during the t(2) period. Here, we present a 2D F1-(1)H-coupled (1)H-(15)N heteronuclear correlation experiment that provides a 1:3:3:1 quartet for AX(3) spin system and a 1:2:1 triplet for AX(2). The experiment is a derivative of 2D HISQC experiment [J. Iwahara, Y.S. Jung, G.M. Clore, Heteronuclear NMR spectroscopy for lysine NH(3) groups in proteins: unique effect of water exchange on (15)N transverse relaxation. J. Am. Chem. Soc. 129 (2007) 2971-2980] and contains a scheme that kills anti-phase single-quantum terms generated in the t(1) period. The purge scheme is essential to observe in-phase single-quantum multiplets. Applications to the NH(2) and NH(3)(+) groups in proteins are demonstrated.

Direct evidence for deprotonation of a lysine side chain buried in the hydrophobic core of a protein

We report direct evidence for deprotonation of a lysine side chain buried in the hydrophobic core of a protein, demonstrating heteronuclear 1H-15N NMR data on the Lys-66 side chain amine (Nzeta) group in the delta-PHS/V66K variant of staphylococcal nuclease. Previous crystallographic study has shown that the Lys-66 Nzeta group is completely buried in the hydrophobic core. On the basis of double and triple resonance experiments, we found that the 1Hzeta and 15Nzeta chemical shifts at pH 8.0 and 6 degrees C for the buried lysine are 0.81 and 23.3 ppm, respectively, which are too abnormal to correspond to the protonated (NH3+) state. Further investigations using a model system suggested that the abnormal 1H and 15N chemical shifts represent the deprotonated (NH2) state of the Lys-66 Nzeta group. More straightforward evidence for the deprotonation was obtained with 2D F1-1H-coupled 1H-15N heteronuclear correlation experiments. Observed 15N multiplets clearly indicated that the spin system for the Lys-66 Nzeta group is AX2 (NH2) rather than AX3 (NH3+). Interestingly, although the amine group is buried in the hydrophobic core, the hydrogen exchange between water and the Lys-66 Nzeta group was found to be relatively rapid (93 s(-1) at -1 degrees C), which suggests the presence of a dynamic process such as local unfolding or water penetration. The partial self-decoupling effect on 15Nzeta multiplets due to the rapid hydrogen exchange is also discussed.

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.