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 Spectroscopic Overview of Intramolecular Hydrogen Bonds of NH…O,S,N Type

Intramolecular NH…O,S,N interactions in non-tautomeric systems are reviewed in a broad range of compounds covering a variety of NH donors and hydrogen bond acceptors. 1H chemical shifts of NH donors are good tools to study intramolecular hydrogen bonding. However in some cases they have to be corrected for ring current effects. Deuterium isotope effects on 13C and 15N chemical shifts and primary isotope effects are usually used to judge the strength of hydrogen bonds. Primary isotope effects are investigated in a new range of magnitudes. Isotope ratios of NH stretching frequencies, νNH/ND, are revisited. Hydrogen bond energies are reviewed and two-bond deuterium isotope effects on 13C chemical shifts are investigated as a possible means of estimating hydrogen bond energies.

Phosphoserine inhibits neighboring arginine methylation in the RKS motif of histone H3

The effects of phosphorylation of histone H3 at serine 10 have been studied in the context of other posttranslational modifications such as lysine methylation. We set out to investigate the impact of phosphoserine-10 on arginine-8 methylation. We performed methylation reactions using peptides based on histone H3 that contain a phosphorylated serine and compared the extent of arginine methylation with unmodified peptides. Results obtained via fluorography indicate that peptides containing a phosphorylated serine-10 inhibit deposition of methyl groups to arginine-8 residues. To further explore the effects of phosphoserine on neighboring arginine residues, we physically characterized the non-covalent interactions between histone H3 phosphoserine-10 and arginine-8 using 31P NMR spectroscopy. A salt bridge was detected between the negatively charged phosphoserine-10 and the positively charged unmodified arginine-8 residue. This salt bridge was not detected when arginine-8 was symmetrically dimethylated. Finally, molecular simulations not only confirm the presence of a salt bridge but also identify a subset of electrostatic interactions present when arginine is replaced with alanine. Taken together, our work suggests that the negatively charged phosphoserine maximizes its interactions. By limiting its exposure and creating new contacts with neighboring residues, it will inhibit deposition of neighboring methyl groups, not through steric hindrance, but by forming intrapeptide interactions that may mask substrate recognition. Our work provides a mechanistic framework for understanding the role of phosphoserine on nearby amino acid residues and arginine methylation.

Isotope effects on chemical shifts in the study of hydrogen bonded biological systems

This review deals with biological systems and with deuterium isotope effects on chemical shifts caused by the replacement of OH, NH or SH protons by deuterons. Hydrogen bonding is clearly of central importance. Isotope effects on chemical shifts seems very suitable for use in studies of structures and reactions in the interior of proteins, as exchange of the label can be expected to be slow. One-bond deuterium isotope effects on ¹⁵N chemical shifts, and two-bond effects on ¹H chemical shifts for N(D)H_x systems can be used to gauge hydrogen bond strength in proteins as well as in salt bridges. Solvent isotope effects on ¹⁹F chemical shifts show promise in monitoring solvent access. Equilibrium isotope effects need in some cases to be taken into account. Schemes for calculation of deuterium isotope effects on chemical shifts are discussed and it is demonstrated how calculations may be used in the study of complex biological systems.

The Structure of the "Vibration Hole" around an Isotopic Substitution-Implications for the Calculation of Nuclear Magnetic Resonance (NMR) Isotopic Shifts

Calculations of nuclear magnetic resonance (NMR) isotopic shifts often rest on the unverified assumption that the “vibration hole”, that is, the change of the vibration motif upon an isotopic substitution, is strongly localized around the substitution site. Using our recently developed difference-dedicated (DD) second-order vibrational perturbation theory (VPT2) method, we test this assumption for a variety of molecules. The vibration hole turns out to be well localized in many cases but not in the interesting case where the H/D substitution site is involved in an intra-molecular hydrogen bond. For a series of salicylaldehyde derivatives recently studied by Hansen and co-workers (Molecules2019, 24, 4533), the vibrational hole was found to stretch over the whole hydrogen-bond moiety, including the bonds to the neighbouring C atoms, and to be sensitive to substituent effects. We discuss consequences of this finding for the accurate calculation of NMR isotopic shifts and point out directions for the further improvement of our DD-VPT2 method.

Supramolecular assemblies involving salt bridges: DFT and X-ray evidence of bipolarity

Three new aminopyridinium/4,4′-oxydibenzoate salts have been synthesized and structurally characterized. A common feature of these compounds is the formation of antiparallel π-stacked salt bridges.

Efficient calculation of NMR isotopic shifts: Difference-dedicated vibrational perturbation theory

We present difference-dedicated second-order vibrational perturbation theory (VPT2) as an efficient method for the computation of nuclear magnetic resonance (NMR) isotopic shifts, which reflect the geometry dependence of the NMR property in combination with different vibration patterns of two isotopologues. Conventional calculations of isotopic shifts, e.g., by standard VPT2, require scanning the geometry dependence over the whole molecule, which becomes expensive rapidly as the molecule size increases. In DD-VPT2, this scan can be restricted to a small region around the substitution site. At the heart of DD-VPT2 is a set of localized vibration modes common to the two isotopologues and designed such that the difference between the vibration patterns is caught by a small subset of them (usually fewer than 10). We tested the DD-VPT2 method for a series of molecules with increasing size and found that this method provides results with the same quality as VPT2 and in good agreement with the experiment, with computational savings up to 95% and less numerical instabilities. The method is easy to automatize and straightforward to generalize to other molecular properties.

Quantifying the Strength of a Salt Bridge by Neutron Scattering and Molecular Dynamics

The molecular structure and strength of a model salt bridge between a guanidinium cation (side chain group of arginine) and the acetate carboxylic group in an aqueous solution is characterized by a combination of neutron diffraction with isotopic substitution and molecular dynamics simulations. The present neutron scattering experiments provide direct information about ion pairing in the solution. At the same time, these measurements are used to assess the quality of the force field employed in the simulation. We show that a standard nonpolarizable force field overestimates the strength of salt bridges. In contrast, accounting for electronic polarization effects via charge scaling allows to quantitatively reproduce the experiment. Such simulations are used to quantify the weak character of a fully hydrated salt bridge. Finally, on top of the canonical hydrogen-bonding binding mode, we uncover another interaction motif involving an out-of-plane hydrophobic contact of the acetate methyl group with the guanidinium cation.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

How well do force fields capture the strength of salt bridges in proteins?

Salt bridges form between pairs of ionisable residues in close proximity and are important interactions in proteins. While salt bridges are known to be important both for protein stability, recognition and regulation, we still do not have fully accurate predictive models to assess the energetic contributions of salt bridges. Molecular dynamics simulation is one technique that may be used study the complex relationship between structure, solvation and energetics of salt bridges, but the accuracy of such simulations depends on the force field used. We have used NMR data on the B1 domain of protein G (GB1) to benchmark molecular dynamics simulations. Using enhanced sampling simulations, we calculated the free energy of forming a salt bridge for three possible lysine-carboxylate ionic interactions in GB1. The NMR experiments showed that these interactions are either not formed, or only very weakly formed, in solution. In contrast, we show that the stability of the salt bridges is overestimated, to different extents, in simulations of GB1 using seven out of eight commonly used combinations of fixed charge force fields and water models. We also find that the Amber ff15ipq force field gives rise to weaker salt bridges in good agreement with the NMR experiments. We conclude that many force fields appear to overstabilize these ionic interactions, and that further work may be needed to refine our ability to model quantitatively the stability of salt bridges through simulations. We also suggest that comparisons between NMR experiments and simulations will play a crucial role in furthering our understanding of this important interaction.

A 13C-detected 15N double-quantum NMR experiment to probe arginine side-chain guanidinium 15Nη chemical shifts

Arginine side-chains are often key for enzyme catalysis, protein-ligand and protein-protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine ¹³C^ζ-detected NMR experiment in which a double-quantum coherence involving the two ¹⁵N^η nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two ¹⁵N^η; this new approach is shown to eliminate the previously deleterious line broadenings of ¹⁵N^η resonances caused by the partially restricted rotation about the C^ζ-N^ε bond. Consequently, sharp and well-resolved ¹⁵N^η resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine ¹⁵N^η nuclei becomes possible using the double-quantum experiment.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Hydrogen bond mediated stabilization of the salt bridge structure for the glycine dimer anion

The formation of a salt bridge in deprotonated glycine dimer anions in a solvent-free environment is investigated using both infrared multiple photon dissociation spectroscopy between 600 and 1800 cm(-1) and theory. The zwitterionic and nonzwitterionic forms of glycine in this complex are computed to be nearly iso-energetic, yet predominantly the zwitterionic form is observed experimentally. The zwitterion stability is attributed to both the Coulombic attraction and the high stabilization from intramolecular hydrogen bonding that drives the energetic cost of proton transfer in a solvent free environment. These results show that there is a fine balance between the stabilities of these two forms of the anion. Elucidating the role of intrinsic factors, such as hydrogen bonding, can lead to a better understanding of the stabilities of salt bridges in the interiors of large proteins or at protein interfaces.

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.

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.