New techniques for the measurement of C'N and C'H(N) J coupling constants across hydrogen bonds in proteins

Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study

Hydrogen bonding is a crucial feature of biomolecules, but its characterization in glycans dissolved in aqueous solutions is challenging due to rapid hydrogen exchange between hydroxyl groups and H2O. In principle, the scalar (J) coupling constant can reveal the relative orientation of the atoms in the molecule. In contrast to J-coupling through H-bonds reported in proteins and nucleic acids, research on J-coupling through H-bonds in glycans dissolved in water is lacking. Here, we use sucrose as a model system for H-bonding studies; its structure, which consists of glucose (Glc) and fructose (Frc), is well-studied, and it is readily available. We apply the in-phase, antiphase-HSQC-TOCSY and quantify previously unreported through H-bond J-values for Frc-OH1-Glc-OH2 in H2O. While earlier reports of Brown and Levy indicate this H-bond as having only a single direction, our reported findings indicate the potential presence of two involving these same atoms, namely, G2OH ➔ F1O and F1OH ➔ G2O (where F and G stand for Frc and Glc, respectively). The calculated density functional theory J-values for the G2OH ➔ F1O agree with the experimental values. Additionally, we detected four other possible H-bonds in sucrose, which require different phi, psi (ϕ, ψ) torsion angles. The ϕ, ψ values are consistent with previous predictions of du Penhoat et al. and Venable et al. Our results will provide new insights into the molecular structure of sucrose and its interactions with proteins.

Intramolecular Hydrogen Bonding in N6-Substituted 2-Chloroadenosines: Evidence from NMR Spectroscopy

Two forms were found in the NMR spectra of N⁶-substituted 2-chloroadenosines. The proportion of the mini-form was 11-32% of the main form. It was characterized by a separate set of signals in COSY, ¹⁵N-HMBC and other NMR spectra. We assumed that the mini-form arises due to the formation of an intramolecular hydrogen bond between the N7 atom of purine and the N⁶-CH proton of the substituent. The ¹H,¹⁵N-HMBC spectrum confirmed the presence of a hydrogen bond in the mini-form of the nucleoside and its absence in the main form. Compounds incapable of forming such a hydrogen bond were synthesized. In these compounds, either the N7 atom of the purine or the N⁶-CH proton of the substituent was absent. The mini-form was not found in the NMR spectra of these nucleosides, confirming the importance of the intramolecular hydrogen bond in its formation.

Structural characterisation of natural products by means of quantum chemical calculations of NMR parameters: new insights

In this review, we focus in all aspects of NMR simulation of natural products, from the fundamentals to the new computational toolboxes available, combining advanced quantum chemical calculations with upstream data processing and machine learning.

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.

A review with comprehensive data on experimental indirect scalar NMR spin-spin coupling constants across hydrogen bonds

Scalar NMR spin-spin coupling constants across hydrogen bonds are fundamental in structural studies and as test grounds for theoretical calculations. Since they are scattered among many articles of different kinds, it seems useful to collect them in the most comprehensive way.

Characterizing the Cooperativity in H-Bonded Amino Structures

Density functional theory calculations were used to examine the effect of H-bond cooperativity on the magnitude of the NMR chemical shifts and spin-spin coupling constants in a C4h-symmetric G-quartet and in structures consisting of six cyanamide monomers. These included two ring structures (a planar C6h-symmetric structure and a nonplanar S6-symmetric structure) and two linear chain structures (a fully optimized planar Cs-symmetric chain and a planar chain structure where all intra- and intermolecular parameters were constrained to be identical). The NMR parameters were computed for the G-quartet and cyanamide structures, as well as for shorter fragments derived from these assemblies without reoptimization. In the ring structures and the chain with identical monomers, the intra- and intermolecular geometries of the cyanamides were identical, thereby allowing the study of cooperative effects in the absence of geometry changes. The magnitude of the |1JNH| coupling, 1H and 15N chemical shifts of the H-bonding amino N-H group, and the |h2JNN| H-bond coupling increased, whereas the size of the |1JNH| coupling of the non-H-bonded amino N-H bonds of the first amino group in the chain, which are roughly perpendicular to the H-bonding network, decreased in magnitude when H-bonding monomers were progressively added to extending ring or chain structures. These effects are attributed to electron redistribution induced by the presence of the nearby H-bonding guanine or cyanamide molecules.

Characterization of the monovalent ion position and hydrogen-bond network in guanine quartets by DFT calculations of NMR parameters

Conformational stability of G-quartets found in telomeric DNA quadruplex structures requires the coordination of monovalent ions. Here, an extensive Hartree-Fock and density functional theory analysis of the energetically favored position of Li+, Na+, and K+ ions is presented. The calculations show that at quartet-quartet distances observed in DNA quadruplex structures (3.3 A), the Li+ and Na+ ions favor positions of 0.55 and 0.95 A outside the plane of the G-quartet, respectively. The larger K+ ion prefers a central position between successive G-quartets. The energy barrier separating the minima in the quartet-ion-quartet model are much smaller for the Li+ and Na+ ions compared with the K+ ion; this suggests that K+ ions will not move as freely through the central channel of the DNA quadruplex. Spin-spin coupling constants and isotropic chemical shifts in G-quartets extracted from crystal structures of K+- and Na+-coordinated DNA quadruplexes were calculated with B3LYP/6-311G(d). The results show that the sizes of the trans-hydrogen-bond couplings are influenced primarily by the hydrogen bond geometry and only slightly by the presence of the ion. The calculations show that the R(N2N7) distance of the N2-H2...N7 hydrogen bond is characterized by strong correlations to both the chemical shifts of the donor group atoms and the (h2)J(N2N7) couplings. In contrast, weaker correlations between the (h3)J(N1C6') couplings and single geometric factors related to the N1-H1...O6=C6 hydrogen bond are observed. As such, deriving geometric information on the hydrogen bond through the use of trans-hydrogen-bond couplings and chemical shifts is more complex for the N1-H1...O6=C6 hydrogen bond than for the N2-H2...N7 moiety. The computed trans-hydrogen-bond couplings are shown to correlate with the experimentally determined couplings. However, the experimental values do not show such strong geometric dependencies.

The Properties of Weak and Strong Dihydrogen-Bonded DH⋅⋅⋅HA Complexes

The properties of six dihydrogen-bonded (DHB) dimers with the BeH2 molecule as a proton acceptor were calculated by MP2, CCSD(T) and B3LYP methods. The structural, energetic and spectroscopic parameters are presented and analyzed in terms of their possible correlation with the interaction energy and the intermolecular H...H separation. The symmetry-adapted perturbation theory (SAPT) calculations were performed to gain more insight into the nature of the H...H interactions. The studied complexes are divided into three groups based on the calculated intermolecular distances and the interaction energies which range from approximately -1 to -42 kJ mol(-1). The analysis of the interaction energy components indicates that, in contrast to conventional hydrogen bonds, the induction energy is the most important term in the BeH2NH4+ complex. On the other hand, there is no sharp boundary between the DHB complexes classified as hydrogen bonded and van der Waals systems. The complexation-induced changes in vibrational frequencies and in proton shielding constants show a relationship with the interaction energy. The values of the 2hJXH and 3hJBeX coupling constants correlate well with the interaction energy and with the intermolecular distance.

Molecular details of Itk activation by prolyl isomerization and phospholigand binding: the NMR structure of the Itk SH2 domain bound to a phosphopeptide

The Src homology 2 (SH2) domain of interleukin-2 tyrosine kinase (Itk) is a critical component of the regulatory apparatus controlling the activity of this immunologically important enzyme. To gain insight into the structural features associated with the activated form of Itk, we have solved the NMR structure of the SH2 domain bound to a phosphotyrosine-containing peptide (pY) and analyzed changes in trans-hydrogen bond scalar couplings ((3h)J(NC')) that result from pY binding. Isomerization of a single prolyl imide bond in this domain is responsible for simultaneous existence of two distinct SH2 conformers. Prolyl isomerization directs ligand recognition: the trans conformer preferentially binds pY. The structure of the SH2/pY complex provides insight into the ligand specificity; the BG loop in the ligand-free trans SH2 conformer is pre-arranged for optimal contacts with the pY+3 residue of the ligand. Analysis of (3h)J(NC') couplings arising from hydrogen bonds has revealed propagation of structural changes from the pY binding pocket to the CD loop containing conformationally heterogeneous proline as well as to the alphaB helix, on the opposite site of the domain. These findings offer a structural framework for understanding the roles of prolyl isomerization and pY binding in Itk regulation.

CLOPPA-IPPP Analysis of Electronic Mechanisms of Intermolecular 1hJ(A,H) and 2hJ(A,D) Spin−Spin Coupling Constants in Systems with D−H···A Hydrogen Bonds

The electronic origin of intermolecular (2h)J(A,D) and (1h)J(A,H) couplings is discussed by means of the CLOPPA-IPPP approach in several model complexes with D-H...A hydrogen bonds. It is found that the origin of these couplings is mainly due to the interaction between the acceptor sigma lone pair and vacant molecular orbitals localized in the D-H...A moiety, regardless of the donor and acceptor nuclei. The problem of the larger absolute value of (2h)J(A,D) compared to (1h)J(A,H) is also addressed.

Cofactor-Apoprotein Hydrogen Bonding in Oxidized and Fully Reduced Flavodoxin Monitored by Trans-Hydrogen-Bond Scalar Couplings

Hydrogen bonding plays a key role in the tight binding of the FMN cofactor and the regulation of its redox properties in flavodoxins. Hydrogen bonding interactions can be directly observed in solution by multidimensional heteronuclear NMR spectroscopy through the scalar couplings between donor and acceptor nuclei. Here we report on the detection of intermolecular trans-hydrogen-bond couplings ((h)J) between the flavin ring system and the backbone of Desulfovibrio vulgaris flavodoxin in the oxidized and the two-electron reduced states. For this purpose, experiments are adapted from pulse sequences previously applied to determining (h)J coupling constants in nucleic acid-base pairs and proteins. The resulting (h2)J(N,N), (h4)J(N,N), (h3)J(C,N), and (h1)J(H,N) couplings involve the (15)N(1), (13)C(2), and (15)N(3) nuclei of the pyrimidine moiety of FMN, whereas no such interactions are detectable for (13)C(4) and (15)N(5). Several long-range (15)N-(15)N, (13)C-(15)N, and (1)H-(15)N J-coupling constants within the flavin are obtained as "by-products". The magnitudes of both (h)J and regular J couplings are found to be dependent on the redox state. In general, good correlations between (h)J coupling constants and donor-group (1)H chemical shifts and also crystallographic donor-acceptor distances are observed.

Structure and ubiquitin interactions of the conserved zinc finger domain of Npl4

Ubiquitylated proteins are directed into a large number of different cellular pathways through interactions with effector proteins that contain conserved ubiquitin binding motifs. Here, we report the solution structure and ubiquitin binding properties of one such motif, the Npl4 zinc finger or RanBP2/Nup358 zinc finger (NZF) domain. Npl4 NZF forms a compact module composed of four antiparallel beta-strands linked by three ordered loops. A single zinc ion is coordinated by four conserved cysteines from the first and third loops, which form two rubredoxin knuckles. Npl4 NZF binds specifically, but weakly, to free ubiquitin using a conserved 13TF14 dipeptide to interact with the "Ile-44" surface of ubiquitin. Our studies reveal the structure of this versatile class of protein binding domains and provide a means for identifying the subset of NZF domains likely to bind ubiquitin.

Hydrogen bonds in rubredoxins from mesophilic and hyperthermophilic organisms

The extent and strength of the hydrogen bond networks in rubredoxins from the hyperthermophile Pyrococcus furiosus (PfRd), and its mesophilic analogue Clostridium pasteurianum (CpRd), are examined and compared using NMR spectroscopy. NMR parameters examined in this study include through-hydrogen bond (h3)J(NC)(') scalar couplings and (1)H, (13)C, and (15)N chemical shifts, as well as covalent (1)J(NH) and (1)J(NC)(') scalar couplings. These parameters have allowed the characterization in solution of 12 hydrogen bonds in each protein. Despite a 83% sequence homology and a low RMSD for the backbone heavy atoms (0.648 A) in the crystalline state, subtle, but definite, changes have been identified in the detailed hydrogen-bonding patterns. CpRd shows an increased number of hydrogen bonds in the triple-stranded beta-sheet and an additional hydrogen bond in the multiple-turn segment including residues 14-32. On the other hand, PfRd exhibits an overall strengthening of N-H...O=C hydrogen bonds in the loops involved at the metal binding site as well as evidence for an additional NH...S(Cys) hydrogen bond involving the alanine residue 44. These data, as well as temperature dependence of the NMR parameters, suggest that the particular NMR hydrogen bond pattern found in the hyperthermophile rubredoxin leads to an increased stabilization at the metal binding pocket. It seems to result from a subtle redistribution of hydrogen-bonding interactions between the triple-stranded beta-sheet and the actual metal binding site.

Spin-state-selective excitation in selective 1D inverse NMR experiments

A general and very simple strategy for achieving clean spin-state-selective excitation with full sensitivity in carbon-selective gradient-enhanced 1D HMQC and HSQC pulse schemes is presented. The incorporation of an additional hard 90 degrees (13)C pulse applied along a specific orthogonal axis just prior to acquisition into the conventional sequences allows us to select a simultaneous coherence transfer pathway which usually is not detected. The superimposition of this resulting antiphase magnetization to the conventional in-phase magnetization gives the exclusive excitation of the directly attached proton showing only the alpha or beta spin state of the passive (13)C nucleus. The propagation of this particular spin state to other protons can be accomplished by adding any homonuclear mixing process just after this supplementary pulse. Such an approach affords a suite of powerful selective 1D (13)C-edited NMR experiments which are helpful for resonance assignment purposes in overcrowded proton spin systems and also for the accurate determination of the magnitude and sign of long-range proton-carbon coupling constants in CH spin sytems for samples at natural abundance. Such measurements are performed by measuring the relative displacement of relayed signals in the corresponding alpha and beta 1D subspectra.

Structure determination of membrane-associated proteins from nuclear magnetic resonance data

This Review covers the delineation and optimization of protein-lipid systems for study using solution-state NMR spectroscopy. The first half presents the necessary background for a membrane protein biochemist to initiate collaboration with an NMR spectroscopist. The second half provides guidelines for the spectroscopist on data collection, analysis for obtaining conformational information, and structure generation and assessment. Although the emphasis is on the study of peptides in detergent micelles, methods are outlined for larger membrane-associated proteins and for use of other solubilizing agents.

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.

Determination of internucleotide (h)J(HN) couplings by the modified 2D J(NN)-correlated [(15)N, (1)H] TROSY

Recent developments in the direct observation of J couplings across hydrogen bonds in proteins and nucleic acids provide additional information for structure and function studies of these molecules by NMR spectroscopy. A J(NN)-correlated [(15)N, (1)H] TROSY experiment proposed by Pervushin et al. (Proc. Natl. Acad. Sci. USA 95, 14147-14151, 1998) can be applied to measure (h)J(HN) in smaller nucleic acids in an E.COSY manner. However, it cannot be effectively applied to large nucleic acids, such as tRNA(Trp), since one of the peaks corresponding to a fast relaxing component will be too weak to be observed in the spectra of large molecules. In this Communication, we proposed a modified J(NN)-correlated [(15)N, (1)H] TROSY experiment which enables direct measurement of (h)J(HN) in large nucleic acids.

3hJ coupling between C(alpha) and H(N) across hydrogen bonds in proteins

J couplings between (13)C(alpha) and (1)H(N) across hydrogen bonds in proteins are reported for the first time, and a two- or three-dimensional NMR technique for their measurement is presented. The technique exploits the TROSY effect, i.e., the degree of interference between dipolar and chemical shift anisotropy relaxation mechanisms, for sensitivity enhancement. The 2D or 3D spectra exhibit E.COSY patterns where the splittings in the (13)CO and (1)H(N) dimensions are (1)J((13)C(alpha), (13)CO) and the desired (3h)J((13)C(alpha), (1)H(N)), respectively. A demonstration of the new method is shown for the (15)N,(13)C-labeled protein chymotrypsin inhibitor 2 where 17 (3h)J((13)C(alpha), (1)H(N)) coupling constants ranging from 0 to 1.4 Hz where identified and all of positive sign.