Peak fitting in 2D <Superscript>1</Superscript>H–<Superscript>13</Superscript>C HSQC NMR spectra for metabolomic studies

A modified Lorentzian distribution function is used to model peaks in two-dimensional (2D) ¹H–¹³C heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra. The model fit is used to determine accurate chemical shifts from genuine signals in complex metabolite mixtures such as blood. The algorithm can be used to extract features from a set of spectra from different samples for exploratory metabolomics. First a reference spectrum is created in which the peak intensities are given by the median value over all samples at each point in the 2D spectra so that ¹H–¹³C correlations in any spectra are accounted for. The mathematical model provides a footprint for each peak in the reference spectrum, which can be used to bin the ¹H–¹³C correlations in each HSQC spectrum. The binned intensities are then used as variables in multivariate analyses and those found to be discriminatory are rapidly identified by cross referencing the chemical shifts of the bins with a database of ¹³C and ¹H chemical shift correlations from known metabolites.     

Different effects of SNP and GSNO on mitochondrial O<Stack><Subscript>2</Subscript><Superscript>.−</Superscript></Stack>/H<Subscript>2</Subscript>O<Subscript>2</Subscript> production

Sodium Nitroprusside (SNP) and S-Nitrosoglutathione (GSNO) differently affect mitochondrial H₂O₂ release at Complex-I. mM SNP increases while GSNO decreases the release induced by succinate alone or added on top of NAD-linked substrates. Stimulation likely depends on Nitric Oxide ( ^. NO) (released by SNP but not by GSNO) inhibiting cytochrome oxidase and mitochondrial respiration. Preincubations with SNP or high GSNO (10 mM plus DTE to increases its ^. NO release) induces an inhibition of the succinate dependent H₂O₂ production consistent with a ^. NO dependent covalent modification. However maximal inhibition of the succinate dependent H₂O₂ release is obtained in the presence of low GSNO (20–100 μM), but not with SNP. This inhibition appears independent of ^. NO release since μM GSNO does not affect mitochondrial respiration, or the H₂O₂ detection systems and its effect is very rapid. Inhibition may be partly due to an increased removal of O₂^.− since GSNO chemically competes with NBT and cytochrome C in O₂^.− detection.     

Negative molecular ion of deuterium D<Stack><Subscript>2</Subscript><Superscript>−</Superscript></Stack>

D₂⁻ ions produced in collisions of D⁻ ions with relative energies of 2.5–9.2 eV were detected for the first time. It is shown that the effective cross section for this reaction is no less than 1.5 × 10⁻¹⁴ cm². Along with the theoretically predicted short-lived state of negative molecular deuterium ions, a state existing for more than 1 μs was observed.     

<Superscript>2</Superscript>H–<Superscript>13</Superscript>C correlation solid-state NMR for investigating dynamics and water accessibilities of proteins and carbohydrates

Site-specific determination of molecular motion and water accessibility by indirect detection of ²H NMR spectra has advantages over dipolar-coupling based techniques due to the large quadrupolar couplings and the ensuing high angular resolution. Recently, a Rotor Echo Short Pulse IRrAdiaTION mediated cross polarization (^RESPIRATIONCP) technique was developed, which allowed efficient transfer of ²H magnetization to ¹³C at moderate ²H radiofrequency field strengths available on most commercial MAS probes. In this work, we investigate the ²H–¹³C magnetization transfer characteristics of one-bond perdeuterated CD _n spin systems and two-bond H/D exchanged C–(O)–D and C–(N)–D spin systems in carbohydrates and proteins. Our results show that multi-bond, broadband ²H–¹³C polarization transfer can be achieved using ²H radiofrequency fields of ~50 kHz, relatively short contact times of 1.3–1.7 ms, and with sufficiently high sensitivity to enable 2D ²H–¹³C correlation experiments with undistorted ²H spectra in the indirect dimension. To demonstrate the utility of this ²H–¹³C technique for studying molecular motion, we show ²H–¹³C correlation spectra of perdeuterated bacterial cellulose, whose surface glucan chains exhibit a motionally averaged C6 ²H quadrupolar coupling that indicates fast trans-gauche isomerization about the C5–C6 bond. In comparison, the interior chains in the microfibril core are fully immobilized. Application of the ²H–¹³C correlation experiment to H/D exchanged Arabidopsis primary cell walls show that the O–D quadrupolar spectra of the highest polysaccharide peaks can be fit to a two-component model, in which 74% of the spectral intensity, assigned to cellulose, has a near-rigid-limit coupling, while 26% of the intensity, assigned to matrix polysaccharides, has a weakened coupling of 50 kHz. The latter O–D quadrupolar order parameter of 0.22 is significantly smaller than previously reported C–D dipolar order parameters of 0.46–0.55 for pectins, suggesting that additional motions exist at the C–O bonds in the wall polysaccharides. ²H–¹³C polarization transfer profiles are also compared between statistically deuterated and H/D exchanged GB1.     

Selective suppression of excipient signals in 2D <Superscript>1</Superscript>H–<Superscript>13</Superscript>C methyl spectra of biopharmaceutical products

While the use of ¹H–¹³C methyl correlated NMR spectroscopy at natural isotopic abundance has been demonstrated as feasible on protein therapeutics as large as monoclonal antibodies, spectral interference from aliphatic excipients remains a significant obstacle to its widespread application. These signals can cause large baseline artifacts, obscure protein resonances, and cause dynamic range suppression of weak peaks in non-uniform sampling applications, thus hampering both traditional peak-based spectral analyses as well as emerging chemometric methods of analysis. Here we detail modifications to the 2D ¹H–¹³C gradient-selected HSQC experiment that make use of selective pulsing techniques for targeted removal of interfering excipient signals in spectra of the NISTmAb prepared in several different formulations. This approach is demonstrated to selectively reduce interfering excipient signals while still yielding 2D spectra with only modest losses in protein signal. Furthermore, it is shown that spectral modeling based on the SMILE algorithm can be used to simulate and subtract any residual excipient signals and their attendant artifacts from the resulting 2D NMR spectra.     

Bomb <Superscript>14</Superscript>C enrichment indicates decadal C pool in deep soil?

Studies of changes in soil organic carbon (SOC) stocks normally limit their focus to the upper 20–30 cm of soil, yet 0–20 cm SOC stocks are only ∼40% of 0–1 m SOC. Accounting for only the upper 20–30 cm of SOC has been justifiable assuming that deeper SOC is unreactive since it displays ¹⁴C-derived mean residence times of hundreds or thousands of years. The dramatic increase in the ¹⁴C content of the atmosphere resulting from thermonuclear testing circa 1963 allows the unreactivity of deep SOC to be tested by examining whether deep soils show evidence of ‘bomb-¹⁴C’ incorporation. At depths of 40–100 cm, a well-studied New Zealand soil under stable pastoral management displays progressive enrichment of over 200‰ across samplings in 1959, 1974 and 2002, indicating substantial incorporation of bomb ¹⁴C. This pattern of deep ¹⁴C enrichment—previously observed in 2 well-drained California grassland soils—leads to the hypothesis that roots and/or dissolved organic C transport contribute to a decadally-reactive SOC pool comprising ∼10–40% of SOC below 50 cm. Deep reactive SOC may be important in the global C cycle because it can react to land-use or vegetation change and may respond to different processes than the reactive SOC in the upper 20–30 cm of soil.     

Sensitive <Superscript>13</Superscript>C–<Superscript>13</Superscript>C correlation spectra of amyloid fibrils at very high spinning frequencies and magnetic fields

Sensitive 2D solid-state ¹³C–¹³C correlation spectra of amyloid β fibrils have been recorded at very fast spinning frequencies and very high magnetic fields. It is demonstrated that PARIS-xy recoupling using moderate rf amplitudes can provide structural information by promoting efficient magnetization transfer even under such challenging experimental conditions. Furthermore, it has been shown both experimentally and by numerical simulations that the method is not very sensitive to dipolar truncation effects and can reveal direct transfer across distances of about 3.5–4Å.     

The neural γ<Subscript>2</Subscript>α<Subscript>1</Subscript>β<Subscript>2</Subscript>α<Subscript>1</Subscript>β<Subscript>2</Subscript> gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models—one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2–M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.     

Insight into π-hole interactions containing the inorganic heterocyclic compounds S<Subscript>2</Subscript>N<Subscript>2</Subscript>/SN<Subscript>2</Subscript>P<Subscript>2</Subscript>

Similar to σ-hole interactions, the π-hole interaction has attracted much attention in recent years. According to the positive electrostatic potentials above and below the surface of inorganic heterocyclic compounds S₂N₂ and three SN₂P₂ isomers (heterocyclic compounds 1–4), and the negative electrostatic potential outside the X atom of XH₃ (X = N, P, As), S₂N₂/SN₂P₂?XH₃ (X = N, P, As) complexes were constructed and optimized at the MP2/aug-cc-pVTZ level. The X atom of XH₃ (X = N, P, As) is almost perpendicular to the ring of the heterocyclic compounds. The π-hole interaction energy becomes greater as the trend goes from 1?XH₃ to 4?XH₃. These π-hole interactions are weak and belong to “closed-shell” noncovalent interactions. According to the energy decomposition analysis, of the three attractive terms, the dispersion energy contributes more than the electrostatic energy. The polarization effect also plays an important role in the formation of π-hole complexes, with the contrasting phenomena of decreasing electronic density in the π-hole region and increasing electric density outside the X atom of XH₃ (X = N, P, As).

Open image in new window

Graphical abstract Computed density difference plots for the complexes 3?NH ₃ (a ₁), 3?PH ₃ (b ₁), 3?AsH ₃ (c ₁) and electron density shifts for the complexes 3?NH ₃ (a ₂), 3?PH ₃ (b ₂),3?AsH ₃ (c ₂) on the 0.001 a.u. contour

     

Determination of the lifetime of D<Stack><Subscript>2</Subscript><Superscript>−</Superscript></Stack> and HD<Superscript>−</Superscript> ions

A method for determining the lifetime of unstable ions is described. The method is based on measuring the decrease in the ion beam current onto a fixed detector with increasing path length of the ion beam from the ion source to the detector. The measurements performed for D^? ₂ and HD^? molecular ions have shown that their lifetimes are 3.5 ± 0.1 and 4.4 ± 0.1 μs, respectively.     

Succinate is the controller of O<Stack><Subscript>2</Subscript><Superscript>−</Superscript></Stack>/H<Subscript>2</Subscript>O<Subscript>2</Subscript> release at mitochondrial complex I : negative modulation by malate,positive by cyanide

Mitochondrial production of H₂O₂ is low with NAD substrates (glutamate/pyruvate, 3 and 2 mM) (G/P) and increases over ten times upon further addition of succinate, with the formation of a sigmoidal curve (semimaximal value at 290 μM, maximal H₂O₂ production at 600 μM succinate). Malate counteracts rapidly the succinate induced increased H₂O₂ release and moves the succinate dependent H₂O₂ production curve to the right. Nitric oxide (NO) and carbon monoxide (CO) are cytochrome c oxidase inhibitors which increase mitochondrial ROS production. Cyanide (CN⁻) was used to mimic NO and CO. In the presence of G/P and succinate (300 μM), CN⁻ progressively increased the H₂O₂ release rate, starting at 1.5 μM. The succinate dependent H₂O₂ production curve was moved to the left by 30 μM CN⁻. The V_max was little modified. We conclude that succinate is the controller of mitochondrial H₂O₂ production, modulated by malate and CN⁻. We propose that succinate promotes an interaction between Complex II and Complex I, which activates O₂⁻ production.     

Presence of Ornithine in the Urate-binding α<Subscript>1</Subscript>—<Subscript>2</Subscript>Globulin

THE urate-binding α₁–α₂ globulin has been isolated from human plasma in a highly purified state¹. The protein was purified by DEAE-‘Sephadex’, ammonium sulphate precipitation and semi-preparative Polyacrylamide gel electrophoresis. The urate-binding α₁–α₂ globulin is a rod-shaped glycoprotein, containing 12.1% carbohydrate, with an isoelectric point of 4.6 and a molecular weight of 67,000 ± 4,000. Amino-acid analysis indicated an unknown basic compound which appeared as an extra peak just in front of lysine¹. To identify this compound, high voltage paper electrophoresis has been carried out on a plate electrophoresis apparatus in pyridine-acetate buffer pH 3.5. A spot separated out corresponding to ornithine. Amino-acid analysis on a BC-200 automatic analyser (Bio-Cal Instruments Co., West Germany), with a 54 cm column at 55° C and with 0.35 M sodium citrate buffer, pH 5.28, as elution buffer at a flow-rate of 150 ml./h, showed that ornithine was present. The presence of ornithine in the protein hydrolysate was also verified by gas chromatography/mass spectrometry².     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C,and <Superscript>15</Superscript>N resonance assignments for the pro-inflammatory cytokine interleukin-36α

Interleukin-36α (IL-36α) is a recently characterised member of the interleukin-1 superfamily. It is involved in the pathogenesis of inflammatory arthritis in one third of psoriasis patients. By binding of IL-36α to its receptor IL-36R via the NF-κB pathway other cytokines involved in inflammatory and apoptotic cascade are activated. The efficacy of complex formation is controlled by N-terminal processing. To obtain a more detailed view on the structure function relationship we performed a heteronuclear multidimensional NMR investigation and here report the ¹H, ¹³C, and ¹⁵N resonance assignments for the backbone and side chain nuclei of the pro-inflammatory cytokine interleukin-36α.     

A simple strategy for <Superscript>13</Superscript>C,<Superscript>1</Superscript>H labeling at the Ile-γ2 methyl position in highly deuterated proteins

A straightforward approach for the production of highly deuterated proteins labeled with ¹³C and ¹H at Ile-γ2 methyl positions is described. The utility of the methodology is illustrated with an application involving the half proteasome (360 kDa). High quality 2D Ile ¹³C^γ2,¹H^γ2 HMQC data sets, exploiting the methyl-TROSY principle, are recorded with excellent sensitivity and resolution, that compare favorably with Ile ¹³C^δ1,¹H^δ1 spectra. This labeling scheme adds to a growing list of different approaches that are significantly impacting the utility of solution NMR spectroscopy in studies of supra-molecular systems.     

Accurate measurement of <Superscript>15</Superscript>N–<Superscript>13</Superscript>C residual dipolar couplings in nucleic acids

New 3D HCN quantitative J (QJ) pulse schemes are presented for the precise and accurate measurement of one-bond ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 residual dipolar couplings (RDCs) in weakly aligned nucleic acids. The methods employ ¹H–¹³C multiple quantum (MQ) coherence or TROSY-type pulse sequences for optimal resolution and sensitivity. RDCs are obtained from the intensity ratio of H₁–C₁–N_1/9 (MQ-HCN-QJ) or H_6/8–C_6/8–N_1/9 (TROSY-HCN-QJ) correlations in two interleaved 3D NMR spectra, with dephasing intervals of zero (reference spectrum) and 1/(2J_NC) (attenuated spectrum). The different types of ¹⁵N–¹³C couplings can be obtained by using either the 3D MQ-HCN-QJ or TROSY-HCN-QJ pulse scheme, with the appropriate setting of the duration of the constant-time ¹⁵N evolution period and the offset of two frequency-selective ¹³C pulses. The methods are demonstrated for a uniformly ¹³C, ¹⁵N-enriched 24-nucleotide stem-loop RNA sequence, helix-35, aligned in the magnetic field using phage Pf1. For measurements of RDCs systematic errors are found to be negligible, and experiments performed on a 1.5 mM helix-35 sample result in an estimated precision of ca. 0.07 Hz for ¹D_NC, indicating the utility of the measured RDCs in structure validation and refinement. Indeed, for a complete set of ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 dipolar couplings obtained for the stem nucleotides, the measured RDCs are in excellent agreement with those predicted for an NMR structure of helix-35, refined using independently measured observables, including ¹³C–¹H, ¹³C–¹³C and ¹H–¹H dipolar couplings.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-0646-2.     

β<Subscript>2</Subscript>-Integrin-Mediated Adhesion and Intracellular Ca<Superscript>2+</Superscript> Release in Human Eosinophils

Human eosinophils spontaneously adhere to various substrates in the absence of exogenously added activators. In the present study a method was developed for characterizing eosinophil adhesion by measuring changes in impedance. Impedance measurements were performed in HCO₃-buffered HybriCare medium maintained in a humidified 5% CO₂ incubator at 37°C. Impedance increased by more than 1 kΩ within minutes after eosinophils made contact with the substrate, reaching a peak within 20 min. Blocking mobilization of intracellular [Ca²⁺] that precedes adhesion with BAPTA-AM (10 μM) completely inhibited the rise in impedance as well as the changes in cell shape typically observed in adherent cells. However, lowering the extracellular [Ca²⁺] with 2.5 mM EGTA did not inhibit the increase in impedance. Pretreatment with anti-CD18 antibody to block substrate interactions with β₂-integrins, or jasplakinolide (2 μM) to block actin reorganization, abolished the increase in impedance and adherent morphology of the cells. Exposure of eosinophils to the phosphatidylinositol 3 kinase inhibitor LY294002 (5 μM) or treatment with protein kinase C zeta pseudosubstrate to competitively inhibit activity of the enzyme significantly reduced the increase in impedance and inhibited the cell spreading associated with adhesion. These results demonstrate a novel method for measuring eosinophil adhesion and showed that, following formation of a tethered attachment, a rapid increase in intracellular [Ca²⁺] precedes the cytoskeletal rearrangements required for cell shape changes and plasma membrane-substrate interactions associated with adhesion.     

CACA-TOCSY with alternate <Superscript>13</Superscript>C–<Superscript>12</Superscript>C labeling: a <Superscript>13</Superscript>C<Superscript>α</Superscript> direct detection experiment for mainchain resonance assignment,dihedral angle information,and amino acid type identification

We present a ¹³C direct detection CACA-TOCSY experiment for samples with alternate ¹³C–¹²C labeling. It provides inter-residue correlations between ¹³C^α resonances of residue i and adjacent C^αs at positions i − 1 and i + 1. Furthermore, longer mixing times yield correlations to C^α nuclei separated by more than one residue. The experiment also provides C^α-to-sidechain correlations, some amino acid type identifications and estimates for ψ dihedral angles. The power of the experiment derives from the alternate ¹³C–¹²C labeling with [1,3-¹³C] glycerol or [2-¹³C] glycerol, which allows utilizing the small scalar ³J_CC couplings that are masked by strong ¹J_CC couplings in uniformly ¹³C labeled samples.     

Comprehensive determination of <Superscript>3</Superscript>J<Subscript>HNHα</Subscript> for unfolded proteins using <Superscript>13</Superscript>C′-resolved spin-echo difference spectroscopy

An experiment is presented to determine ³J_HNHα coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D ¹H, or 2D and 3D ¹H-¹⁵N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording ¹H-¹³C′ correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.     

The effect of metal ions on the binding of ethanol to human alcohol dehydrogenase β<Subscript>2</Subscript>β<Subscript>2</Subscript>

Molecular docking simulations were performed in this study to investigate the importance of both structural and catalytic zinc ions in the human alcohol dehydrogenase beta(2)beta(2) on substrate binding. The structural zinc ion is not only important in maintaining the structural integrity of the enzyme, but also plays an important role in determining substrate binding. The replacement of the catalytic zinc ion or both catalytic and structural zinc ions with Cu(2+) results in better substrate binding affinity than with the wild-type enzyme. The width of the bottleneck formed by L116 and V294 in the substrate binding pocket plays an important role for substrate entrance. In addition, unfavorable contacts between the substrate and T48 and F93 prevent the substrate from moving too close to the metal ion. The optimal binding position occurs between 1.9 and 2.4 A from the catalytic metal ion.