Predicting <Superscript>13</Superscript>C<Superscript>α</Superscript> chemical shifts for validation of protein structures

The ¹³C^α chemical shifts for 16,299 residues from 213 conformations of four proteins (experimentally determined by X-ray crystallography and Nuclear Magnetic Resonance methods) were computed by using a combination of approaches that includes, but is not limited to, the use of density functional theory. Initially, a validation test of this methodology was carried out by a detailed examination of the correlation between computed and observed ¹³C^α chemical shifts of 10,564 (of the 16,299) residues from 139 conformations of the human protein ubiquitin. The results of this validation test on ubiquitin show agreement with conclusions derived from computation of the chemical shifts at the ab initio Hartree–Fock level. Further, application of this methodology to 5,735 residues from 74 conformations of the three remaining proteins that differ in their number of amino acid residues, sequence and three-dimensional structure, together with a new scoring function, namely the conformationally averaged root-mean-square-deviation, enables us to: (a) offer a criterion for an accurate assessment of the quality of NMR-derived protein conformations; (b) examine whether X-ray or NMR-solved structures are better representations of the observed ¹³C^α chemical shifts in solution; (c) provide evidence indicating that the proposed methodology is more accurate than automated predictors for validation of protein structures; (d) shed light as to whether the agreement between computed and observed ¹³C^α chemical shifts is influenced by the identity of an amino acid residue or its location in the sequence; and (e) provide evidence confirming the presence of dynamics for proteins in solution, and hence showing that an ensemble of conformations is a better representation of the structure in solution than any single conformation. Electronic Supplementary Material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

<Superscript>1</Superscript>H–<Superscript>15</Superscript>N correlation spectroscopy of nanocrystalline proteins

The limits of resolution that can be obtained in ¹H–¹⁵N 2D NMR spectroscopy of isotopically enriched nanocrystalline proteins are explored. Combinations of frequency switched Lee–Goldburg (FSLG) decoupling, fast magic angle sample spinning (MAS), and isotopic dilution via deuteration are investigated as methods for narrowing the amide ¹H resonances. Heteronuclear decoupling of ¹⁵N from the ¹H resonances is also studied. Using human ubiquitin as a model system, the best resolution is most easily obtained with uniformly ²H and ¹⁵N enriched protein where the amides have been exchanged in normal water, MAS at 20 kHz, and WALTZ-16 decoupling of the ¹⁵N nuclei. The combination of these techniques results in average ¹H lines of only 0.26 ppm full width at half maximum. Techniques for optimizing instrument stability and ¹⁵N decoupling are described for achieving the best possible performance in these experiments.     

Synthesis of the β-amyloid fragment <Superscript>5</Superscript>RHDSGY<Superscript>10</Superscript> and its isomers

The peptide RHDSGY, a fragment of the human β-amyloid Zn-binding site, and its isomers RH(D-Asp)SGY and RH(β-Asp)SGY have been obtained as amides by means of solid-phase synthesis and analyzed by HPLC and various mass spectrometric methods. The problem of low yield of the RHDSGY peptide and its isomers attributed to 9-fluorenylmethoxycarbonyl (Fmoc)-amino acids and/or formation of such side-products as RH(β-Asp)SGY (or RHDSGY during synthesis of RH(β-Asp)SGY) and RH(Asp-imide) SGY was solved via selection of individual reagents for removal of Fmoc groups from α-amino groups of the growing peptide chain.     

Rapid measurement of <Superscript>3</Superscript>J(H<Superscript>N</Superscript>–H<Superscript>α</Superscript>) and <Superscript>3</Superscript>J(N–H<Superscript>β</Superscript>) coupling constants in polypeptides

We present two NMR experiments, (3,2)D HNHA and (3,2)D HNHB, for rapid and accurate measurement of 3J(H N-H alpha) and 3J(N-H beta) coupling constants in polypeptides based on the principle of G-matrix Fourier transform NMR spectroscopy and quantitative J-correlation. These experiments, which facilitate fast acquisition of three-dimensional data with high spectral/digital resolution and chemical shift dispersion, will provide renewed opportunities to utilize them for sequence specific resonance assignments, estimation/characterization of secondary structure with/without prior knowledge of resonance assignments, stereospecific assignment of prochiral groups and 3D structure determination, refinement and validation. Taken together, these experiments have a wide range of applications from structural genomics projects to studying structure and folding in polypeptides.     

Sequence specific <Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N resonance assignments of a cataract-related variant G57W of human γS-crystallin

γS-crystallin is a major structural component of the human eye lens, which maintains its stability over the lifetime of an organism with negligible turnover. The G57W mutant of human γS-crystallin (abbreviated hereafter as γS-G57W) is associated with dominant congenital cataracts. In order to provide a structural basis for the ability of γS-G57W causing cataract, we have cloned, overexpressed, isolated and purified the protein. The 2D [¹⁵N–¹H]-HSQC spectrum recorded with uniformly ¹³C/¹⁵N-labelled γS-G57W was highly dispersed indicating the protein to adopt an ordered conformation. In this paper, we report almost complete sequence-specific ¹H, ¹³C and ¹⁵N resonance assignments of γS-G57W using a suite of heteronuclear 3D NMR experiments.     

Impact of two-bond <Superscript>15</Superscript>N–<Superscript>15</Superscript>N scalar couplings on <Superscript>15</Superscript>N transverse relaxation measurements for arginine side chains of proteins

NMR relaxation of arginine (Arg) ¹⁵N_ε nuclei is useful for studying side-chain dynamics of proteins. In this work, we studied the impact of two geminal ¹⁵N–¹⁵N scalar couplings on measurements of transverse relaxation rates (R ₂ ) for Arg side-chain ¹⁵N_ε nuclei. For 12 Arg side chains of the DNA-binding domain of the Antp protein, we measured the geminal ¹⁵N–¹⁵N couplings ( ² J _NN ) of the ¹⁵N_ε nuclei and found that the magnitudes of the ² J _NN coupling constants were virtually uniform with an average of 1.2 Hz. Our simulations, assuming ideal 180° rotations for all ¹⁵N nuclei, suggested that the two ² J _NN couplings of this magnitude could in principle cause significant modulation in signal intensities during the Carr–Purcell-Meiboom–Gill (CPMG) scheme for Arg ¹⁵N_ε R ₂ measurements. However, our experimental data show that the expected modulation via two ² J _NN couplings vanishes during the ¹⁵N CPMG scheme. This quenching of J modulation can be explained by the mechanism described in Dittmer and Bodenhausen (Chemphyschem 7:831–836, 2006). This effect allows for accurate measurements of R ₂ relaxation rates for Arg side-chain ¹⁵N_ε nuclei despite the presence of two ² J _NN couplings. Although the so-called recoupling conditions may cause overestimate of R ₂ rates for very mobile Arg side chains, such conditions can readily be avoided through appropriate experimental settings.     

The role of Ca<Superscript>2+</Superscript>, H<Superscript>+</Superscript>, and Cl<Superscript>−</Superscript> ions in generation of variation potential in pumpkin plants

The contributions of Ca²⁺, H⁺, and Cl⁻ in generation of variation potentials (VP) in 3- to 4-week-old pumpkin (Cucurbita pepo L., cv. Mozoleevskaya) plants were assessed. During VP generation, transient alkalinization of the medium around the stem was recorded with a potentiometric method. The pH changes were kinetically similar to the electric potential changes and were apparently due to temporal suppression of the plasma-membrane electrogenic H⁺ pump. These data and the observed inhibition of VP in the stem zone treated locally with a metabolic inhibitor (NaN₃) indicate that the VP generation is related to the reversible suppression of the H⁺-pump. The anion channel blocker (ethacrynic acid) decelerated significantly the front slope of VP and reduced the VP amplitude. A short-term increase in external Cl⁻ concentration around the stem was observed during potential transients representing the VP front slope and the pulses integrated into VP. The removal of Ca²⁺ from extracellular medium inhibited the VP generation. It is proposed that Ca²⁺ plays a role in activation of anion channels and in the H⁺-pump inactivation. The VP generation is probably determined by a complex mechanism, with contributions from passive ion fluxes (Ca²⁺, Cl⁻) moving along the electrochemical gradients and from changes in the electrogenic pump activity.     

H<Superscript>N</Superscript>, N,C<Superscript>α</Superscript> and C<Superscript>β</Superscript> assignments of the two periplasmic domains of <Emphasis Type="Italic">Neisseria meningitidis</Emphasis> DsbD

DsbD is a disulfide bond reductase present in the inner membrane of many Gamma-Proteobacteria. In the human pathogen Neisseria meningitidis, DsbD is required for viability and represents a potential target for the development of antibiotics. Here we report the chemical shift assignments (H^N, N, C^α and C^β) for the reduced and oxidized forms of the two periplasmic domains of N. meningitidis DsbD, n-NmDsbD and c-NmDsbD. The backbone amide resonances in all four forms were completely assigned, and the secondary structures for the core regions of the proteins were calculated using ¹³C^αβ shifts. The reduced and oxidized forms of each domain have similar secondary shifts suggesting they retain the same fold. We anticipate that these data will provide an important basis for studying the interaction between n-NmDsbD and c-NmDsbD, which is required for electron transfer across the bacterial cytoplasmic membrane.     

Comparison of Alexa Fluor<Superscript>®</Superscript> and CyDye<Superscript>™</Superscript> for practical DNA microarray use

Microarrays are a powerful tool for comparison and understanding of gene expression levels in healthy and diseased states. The method relies upon the assumption that signals from microarray features are a reflection of relative gene expression levels of the cell types under investigation. It has previously been reported that the classical fluorescent dyes used for microarray technology, Cy3 and Cy5, are not ideal due to the decreased stability and fluorescence intensity of the Cy5 dye relative to the Cy3, such that dye bias is an accepted phenomena necessitating dye swap experimental protocols and analysis of differential dye affects. The incentive to find new fluorophores is based on alleviating the problem of dye bias through synonymous performance between counterpart dyes. Alexa Fluor 555 and Alexa Fluor 647 are increasingly promoted as replacements for CyDye in microarray experiments. Performance relates to the molecular and steric similarities, which will vary for each new pair of dyes as well as the spectral integrity for the specific application required. Comparative analysis of the performance of these two competitive dye pairs in practical microarray applications is warranted towards this end. The findings of our study showed that both dye pairs were comparable but that conventional CyDye resulted in significantly higher signal intensities (P < 0.05) and signal minus background levels (P < 0.05) with no significant difference in background values (P > 0.05). This translated to greater levels of differential gene expression with CyDye than with the Alexa Fluor counterparts. However, CyDye fluorophores and in particular Cy5, were found to be less photostable over time and following repeated scans in microarray experiments. These results suggest that precautions against potential dye affects will continue to be necessary and that no one dye pair negates this need.     

A Voltage-Dependent Ca<Superscript>2+</Superscript> Influx Pathway Regulates the Ca<Superscript>2+</Superscript>-Dependent Cl<Superscript>−</Superscript> Conductance of Renal IMCD-3 Cells

We have previously shown that the membrane conductance of mIMCD-3 cells at a holding potential of 0 mV is dominated by a Ca²⁺-dependent Cl⁻ current (I_CLCA). Here we report that I_CLCA activity is also voltage dependent and that this dependence on voltage is linked to the opening of a novel Al³⁺-sensitive, voltage-dependent, Ca²⁺ influx pathway. Using whole-cell patch-clamp recordings at a physiological holding potential (−60 mV), I_CLCA was found to be inactive and resting currents were predominantly K⁺ selective. However, membrane depolarization to 0 mV resulted in a slow, sigmoidal, activation of I_CLCA (T _0.5 ~ 500 s), while repolarization in turn resulted in a monoexponential decay in I_CLCA (T _0.5 ~ 100 s). The activation of I_CLCA by depolarization was reduced by lowering extracellular Ca²⁺ and completely inhibited by buffering cytosolic Ca²⁺ with EGTA, suggesting a role for Ca²⁺ influx in the activation of I_CLCA. However, raising bulk cytosolic Ca²⁺ at −60 mV did not produce sustained I_CLCA activity. Therefore I_CLCA is dependent on both an increase in intracellular Ca²⁺ and depolarization to be active. We further show that membrane depolarization is coupled to opening of a Ca²⁺ influx pathway that displays equal permeability to Ca²⁺ and Ba²⁺ ions and that is blocked by extracellular Al³⁺ and La³⁺. Furthermore, Al³⁺ completely and reversibly inhibited depolarization-induced activation of I_CLCA, thereby directly linking Ca²⁺ influx to activation of I_CLCA. We speculate that during sustained membrane depolarization, calcium influx activates I_CLCA which functions to modulate NaCl transport across the apical membrane of IMCD cells.     

Measurement of imino <Superscript>1</Superscript>H–<Superscript>1</Superscript>H residual dipolar couplings in RNA

Imino ¹H–¹⁵N residual dipolar couplings (RDCs) provide additional structural information that complements standard ¹H–¹H NOEs leading to improvements in both the local and global structure of RNAs. Here, we report measurement of imino ¹H–¹H RDCs for the Iron Responsive Element (IRE) RNA and native E. coli tRNA^Val using a BEST-Jcomp-HMQC2 experiment. ¹H–¹H RDCs are observed between the imino protons in G–U wobble base pairs and between imino protons on neighboring base pairs in both RNAs. These imino ¹H–¹H RDCs complement standard ¹H–¹⁵N RDCs because the ¹H–¹H vectors generally point along the helical axis, roughly perpendicular to ¹H–¹⁵N RDCs. The use of longitudinal relaxation enhancement increased the signal-to-noise of the spectra by ~3.5-fold over the standard experiment. The ability to measure imino ¹H–¹H RDCs offers a new restraint, which can be used in NMR domain orientation and structural studies of RNAs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

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.     

Estimation of denitrification potential in a karst aquifer using the <Superscript>15</Superscript>N and <Superscript>18</Superscript>O isotopes of NO<Stack><Subscript>3</Subscript><Superscript>−</Superscript></Stack>

A confined aquifer in the Malm Karst of the Franconian Alb, South Germany was investigated in order to understand the role of the vadose zone in denitrifiaction processes. The concentrations of chemical tracers Sr²⁺ and Cl^– and concentrations of stable isotope ¹⁸O were measured in spring water and precipitation during storm events. Based on these measurements a conceptual model for runoff was constructed. The results indicate that pre-event water, already stored in the system at the beginning of the event, flows downslope on vertical and lateral preferential flow paths. Chemical tracers used in a mixing model for hydrograph separation have shown that the pre-event water contribution is up to 30%. Applying this information to a conceptual runoff generation model, the values of ¹⁵N and ¹⁸O in nitrate could be calculated. Field observations showed the occurence of significant microbial denitrification processes above the soil/bedrock interface before nitrate percolates through to the deeper horizon of the vadose zone. The source of nitrate could be determined and denitrification processes were calculated. Assuming that the nitrate reduction follows a Rayleigh process one could approximate a nitrate input concentration of about 170 mg/l and a residual nitrate concentration of only about 15%. The results of the chemical and isotopic tracers postulate fertilizers as nitrate source with some influence of atmospheric nitrate. The combined application of hydrograph separation and determination of isotope values in ¹⁵N and ¹⁸O of nitrate lead to an improved understanding of microbial processes (nitrification, denitrification) in dynamic 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.     

Variability and directionality of temporal changes in δ<Superscript>13</Superscript>C and δ<Superscript>15</Superscript>N of aquatic invertebrate primary consumers

Seasonal oscillations in the carbon (δ¹³C) and nitrogen (δ¹⁵N) isotope signatures of aquatic algae can cause seasonal enrichment–depletion cycles in the isotopic composition of planktonic invertebrates (e.g., copepods). Yet, there is growing evidence that seasonal enrichment–depletion cycles also occur in the isotope signatures of larger invertebrate consumers, taxa used to define reference points in isotope-based trophic models (e.g., trophic baselines). To evaluate the general assumption of temporal stability in non-zooplankton aquatic invertebrates, δ¹³C and δ¹⁵N time series data from the literature were analyzed for seasonality and the influence of biotic (feeding group) and abiotic (trophic state, climate regime) factors on isotope temporal patterns. The amplitude of δ¹³C and δ¹⁵N enrichment–depletion cycles was negatively related to body size, although all size-classes of invertebrates displayed a winter-to-summer enrichment in δ¹³C and depletion in δ¹⁵N. Among feeding groups, periphytic grazers were more variable and displayed larger temporal changes in δ¹³C than detritivores. For nitrogen, temporal variability and magnitude of directional change of δ¹⁵N was most strongly related to ecosystem trophic state (eutrophic > mesotrophic, oligotrophic). This study provides evidence of seasonality in the isotopic composition of aquatic invertebrates across very broad geographical and ecological gradients as well as identifying factors that are likely to modulate the strength and variability of seasonality. These results emphasize the need for researchers to recognize the likelihood of temporal changes in non-zooplankton aquatic invertebrate consumers at time scales relevant to seasonal studies and, if present, to account for temporal dynamics in isotope trophic models.     

Study of Electric Explosion of Flat Micron-Thick Foils at Current Densities of (5−50)×10<Superscript>8</Superscript> A/cm<Superscript>2</Superscript>

Electric explosions of flat Al, Тi, Ni, Cu, and Та foils with thicknesses of 1?16 μm, widths of 1?8 mm, and lengths of 5?11 mm were studied experimentally on the BIN, XP, and COBRA high-current generators at currents of 40?1000 kA and current densities of (5–50) × 10⁸ A/cm². The images of the exploded foils were taken at different angles to the foil surface by using point projection radiography with an X-pinch hot spot as the radiation source, the spatial resolution and exposure time being 3 μm and 50 ps, respectively, as well by the laser probing method with a spatial resolution of 20 μm and an exposure time of 180 ps. In the course of foil explosion, rapidly expanding objects resembling the core and corona of an exploded wire were observed. It is shown that the core of the exploded foil has a complicated time-varying structure.     

A <Superscript>13</Superscript>C-detected <Superscript>15</Superscript>N double-quantum NMR experiment to probe arginine side-chain guanidinium <Superscript>15</Superscript>N<Superscript>η</Superscript> 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.     

Accuracy of the StressEraser<Superscript>®</Superscript> in the Detection of Cardiac Rhythms

StressEraser is a commercially marketed biofeedback device designed to enhance heart rate variability. StressEraser makes its internal calculations on beat-to-beat measures of finger pulse intervals. However, the accuracy and precision of StressEraser in quantifying interbeat intervals using finger pulse intervals has not been evaluated against standard laboratory equipment using R-R intervals. Accuracy was assessed by simultaneously recording interbeat intervals using StressEraser and a standard laboratory ECG system. The interbeat intervals were highly correlated between the systems. The average deviation in interbeat interval recordings between the systems was approximately 6 ms. Moreover, correlations approached unity between the systems on estimates of heart period, heart rate, and heart rate variability. Feedback from StressEraser is based on an interbeat time series that provides sufficient information to provide an excellent estimate of the dynamic changes in heart rate and heart rate variability. The slight variations between StressEraser and the laboratory equipment in quantifying heart rate and heart rate variability are due to features related to monitoring heart rate with finger pulse: (1) a lack in precision in the peak of the finger pulse relative to the clearly defined inflection point in the R-wave, and (2) contribution of variations in pulse transit time.