An X-ray spectrum estimation method from transmission measurement combined with scatter correction

PurposeConventional x-ray spectrum estimation methods from transmission measurement often lead to inaccurate results when extensive x-ray scatter is present in the measured projection. This study aims to apply the weighted L1-norm scatter correction algorithm in spectrum estimation for reducing residual differences between the estimated and true spectrum.MethodThe scatter correction algorithm is based on a simple radiographic scattering model where the intensity of scattered x-ray is directly estimated from a transmission measurement. Then, the scatter-corrected measurement is used for the spectrum estimation method that consists of deciding the weights of predefined spectra and representing the spectrum as a linear combination of the predefined spectra with the weights. The performances of the estimation method combined with scatter correction are evaluated on both simulated and experimental data.ResultsThe results show that the estimated spectra using the scatter-corrected projection nearly match the true spectra. The normalized-root-mean-square-error and the mean energy difference between the estimated spectra and corresponding true spectra are reduced from 5.8% and 1.33 keV without the scatter correction to 3.2% and 0.73 keV with the scatter correction for both simulation and experimental data, respectively.ConclusionsThe proposed method is more accurate for the acquisition of x-ray spectrum than the estimation method without scatter correction and the spectrum can be successfully estimated even the materials of the filters and their thicknesses are unknown. The proposed method has the potential to be used in several diagnostic x-ray imaging applications.     

Improved frequency resolution in multidimensional constant-time experiments by multidimensional Bayesian analysis

Summary The resolution of spectral frequencies in NMR data obtained from discrete Fourier transformation (DFT) along D constant-time dimensions can be improved significantly through extrapolation of the D-dimensional free induction decay (FID) by multidimensional Bayesian analysis. Starting from Bayesian probability theory for parameter estimation and model detection of one-dimensional time-domain data [Bretthorst, (1990) J. Magn. Reson., 88, 533–551; 552–570; 571–595], a theory for the D-dimensional case has been developed and implemented in an algorithm called BAMBAM (BAyesian Model Building Algorithm in Multidimensions). BAMBAM finds the most probable sinusoidal model to account for the systematic portion of any D-dimensional stationary FID. According to the parameters estimated by the algorithm, the FID is extrapolated in D dimensions prior to apodization and Fourier transformation. Multidimensional Bayesian analysis allows for the detection of signals not resolved by the DFT alone or even by sequential one-dimensional extrapolation from mirror-image linear prediction prior to the DFT. The procedure has been tested with a theoretical two-dimensional dataset and with four-dimensional HN(CO)CAHA (Kay et al. (1992) J. Magn. Reson., 98, 443–450) data from a small protein (8 kDa) where BAMBAM was applied to the ¹³C and H constant-time dimensions.To whom correspondence should be addressed.     

A 4D HCC(CO)NNH experiment for the correlation of aliphatic side-chain and backbone resonances in 13C/15N-labelled proteins

Summary We recently proposed a novel four-dimensional (4D) NMR strategy for the assignment of backbone nuclei in spectra of ¹³C/¹⁵N-labelled proteins (Boucher et al. (1992) J. Am. Chem. Soc., 114, 2262–2264 and J. Biomol. NMR, 2, 631–637). In this paper we extend this approach with a new constant time 4D HCC(CO)NNH experiment that also correlates the chemical shifts of the aliphatic sidechain (¹H and ¹³C) and backbone (¹H, ¹³C and ¹⁵N) nuclei. It separates the sidechain resonances, which may heavily overlap in spectra of proteins with large numbers of similar residues, according to the backbone nitrogen and amide proton chemical shifts. When used in conjunction with a 4D HCANNH or HNCAHA experiment it allows, in principle, complete assignment of aliphatic sidechain and backbone resonances with just two 4D NMR experiments.     

NMR Spectroscopic Study of a DNA Duplex with Mercury-Mediated T-T Base Pairs

Recently, we reported that T-T mismatches can specifically recognize Hg ^II (T-Hg ^II -T pair formation). In order to understand the properties of the T-Hg ^II -T pair, we recorded NMR spectra for a DNA duplex, d(CGCG TT GTCC) ? d(GGAC TT CGCG), with two successive T-T mismatches (Hg ^II -binding sites). We assigned ¹ H resonances for mercury-free and di-mercurated duplexes, and performed titration experiments with Hg ^II by using 1D ¹ H NMR spectra. Because of the above mentioned assignments, we could confirm the existence of mono-mercurated species, because individual components gave independent NMR signals in the titration spectra.     

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.     

Multiple quadrature detection in reduced dimensionality experiments

A new, simple procedure is proposed which enables acquisition of two or more chemical shifts encoded in a common dimension simultaneously in quadrature. For n chemical shifts projected in a single dimension, the expected effect is obtained by interleaved acquisition and appropriate combination of 2ⁿ data sets per increment of respective evolution time. The particular chemical shifts can be calculated from sums and differences of signal frequencies obtained by different combination of the acquired data sets. In comparison to the established reduced dimensionality (RD) techniques, the proposed method enhances resolution due to reduction of the number of signals and requires less evolution time increments owing to narrower spectral width in the RD-domain. We show examples of the application of the new approach to the 2D HNCA and HN(CO)CA techniques with two, and 2D HACANH with three frequencies simultaneously encoded in the t ₁ evolution period, for ¹³C,¹⁵N-labeled ubiquitin.     

In-vivo dose measurements with MOSFET dosimeters during MV portal imaging

BackgroundThe purpose of this study was to investigate the feasibility of MOSFET dosimeter in measuring eye dose during 2D MV portal imaging for setup verification in radiotherapy.Materials and methodsThe in-vivo dose measurements were performed by placing the dosimeters over the eyes of 30 brain patients during the acquisition of portal images in linear accelerator by delivering 1 MU with the field sizes of 10 × 10 cm² and 15 × 15 cm².ResultsThe mean doses received by the left and right eyes of 10 out of 30 patients when both eyes were completely inside the anterior portal field were found to be 2.56 ± 0.2 cGy and 2.75 ± 0.2, respectively. Similarly, for next 10 patients out of the same 30 patients the mean doses to left and right eyes when both eyes were completely out of the anterior portal fields were found to be 0.13 ± 0.02 cGy and 0.17 ± 0.02 cGy, respectively. The mean doses to ipsilateral and contralateral eye for the last 10 patients when one eye was inside the anterior portal field were found to be 3.28 ± 0.2 cGy and 0.36 ± 0.1 cGy, respectively.ConclusionThe promising results obtained during 2D MV portal imaging using MOSFET have shown that this dosimeter is well suitable for assessing low doses during imaging thereby enabling to optimize the imaging procedure using the dosimetric data obtained. In addition, the documentation of the dose received by the patient during imaging procedure is possible with the help of an in-built software in conjunction with the MOSFET reader module.     

A fast NMR method for resonance assignments: application to metabolomics

We present a new method for rapid NMR data acquisition and assignments applicable to unlabeled (¹²C) or ¹³C-labeled biomolecules/organic molecules in general and metabolomics in particular. The method involves the acquisition of three two dimensional (2D) NMR spectra simultaneously using a dual receiver system. The three spectra, namely: (1) G-matrix Fourier transform (GFT) (3,2)D [¹³C, ¹H] HSQC–TOCSY, (2) 2D ¹H–¹H TOCSY and (3) 2D ¹³C–¹H HETCOR are acquired in a single experiment and provide mutually complementary information to completely assign individual metabolites in a mixture. The GFT (3,2)D [¹³C, ¹H] HSQC–TOCSY provides 3D correlations in a reduced dimensionality manner facilitating high resolution and unambiguous assignments. The experiments were applied for complete ¹H and ¹³C assignments of a mixture of 21 unlabeled metabolites corresponding to a medium used in assisted reproductive technology. Taken together, the experiments provide time gain of order of magnitudes compared to the conventional data acquisition methods and can be combined with other fast NMR techniques such as non-uniform sampling and covariance spectroscopy. This provides new avenues for using multiple receivers and projection NMR techniques for high-throughput approaches in metabolomics.     

Fluorescence detected magnetic resonance (FDMR) of green sulfur photosynthetic bacteria Chlorobium sp.

Fluorescence Detected Magnetic Resonance (FDMR) spectra have been measured for whole cells and isolated chlorosomal fractions for the green photosyntheic bacteria Chlorobium phaeobacteroides (containing bacteriochlorophyll e, and isorenieratene as major carotenoid) and Chlorobium limicola (containing bacteriochlorophyll c, and chlorobactene as major carotenoid). The observed transition at 237 MHz (identical in both bacteria) and > 1100 MHz can be assigned, by analogy with published data on other carotenoids, to the 2E and D + E transitions, respectively, of Chlorobium carotenoids. Their zero field splitting (ZFS) parameters are estimated to be: |D|=0.0332 cm^–1 and |E|=0.0039 cm^–1 (chlorobactene), and |D|=0.0355 cm^–1 and |E|=0.0039 cm^–1 (isorenieratene). In the intermediate frequency range 300–1000 MHz the observed transitions can be assigned to chlorosomal bacteriochlorophylls c and e, and to bacteriochlorophyll a located in the chlorosome envelope and water-soluble protein. The bacteriochlorophyll e triplet state measured in 750 nm fluorescence (aggregated chlorosomal BChl e) is characterised by the ZFS parameters: |D|=0.0251 cm^–1 and |E|=0.0050 cm^–1.Abbreviations BChl - bacteriochlorophyll - BPh - bacteriopheophytin - Chl. - Chlorobium - F(A)(O)DMR - fluorescence (absorption) (optical) detected magnetic resonance - FF - fluorescence fading - ISC - intramolecular intersystem crossing - RC - reaction center - ZFS - zero field splitting     

Computer-assisted assignment of 2D1H NMR spectra of proteins: Basic algorithms and application to phoratoxin B

Summary A suite of computer programs (CLAIRE) is described which can be of assistance in the process of assigning 2D¹H NMR spectra of proteins. The programs embody a software implementation of the sequential assignment approach first developed by Wüthrich and co-workers (K. Wüthrich. G. Wider, G. Wagner and W. Braun (1982)J. Mol. Biol. 155, 311). After data-abstraction (peakpicking), the software can be used to detect patterns (spin systems), to find cross peaks between patterns in 2D NOE data sets and to generate assignments that are consistent with all available data and which satisfy a number of constraints imposed by the user. An interactive graphics program calledCONPAT is used to control the entire assignment process as well as to provide the essential feedback from the experimental NMR spectra. The algorithms are described in detail and the approach is demonstrated on a set of spectra from the mistletoe protein phoratoxin B, a homolog of crambin. The results obtained compare well with those reported earlier based entirely on a manual assignment process.     

Improved NMR spectra of a protein–DNA complex through rational mutagenesis and the application of a sensitivity optimized isotope-filtered NOESY experiment

The NMR spectra of the complex between the DNA-binding domain of the Dead ringer protein (DRI-DBD, Gly262-Gly398) and its DNA binding site (DRI-DBD:DNA, 26 kDa) have been optimized by biochemical and spectroscopic means. First, we demonstrate the utility of a modified 2D [F1,F2] ¹³C-filtered NOESY experiment that employs a ¹J_HC versus chemical shift optimized adiabatic ¹³C inversion pulse [Zwahlen, C. et al. (1997) J. Am. Chem. Soc., 119, 6711–6721]. The new sequence is shown to be more sensitive than previously published pulse schemes (up to 40% in favorable cases) and its utility is demonstrated using two protein–DNA complexes. Second, we demonstrate that the targeted replacement of an interfacial aromatic residue in the DRI-DBD:DNA complex substantially reduces line broadening within its NMR spectra. The spectral changes are dramatic, salvaging a protein–DNA complex that was originally ill suited for structural analysis by NMR. This biochemical approach is not a general method, but may prove useful in the spectral optimization of other protein complexes that suffer from interfacial line broadening caused by dynamic changes in proximal aromatic rings.     

A Leap-frog Algorithm for Stochastic Dynamics

Abstract

A third-order algorithm for stochastic dynamics (SD) simulations is proposed, identical to the powerful molecular dynamics leap-frog algorithm in the limit of infinitely small friction coefficient γ. It belongs to the class of SD algorithms, in which the integration time step Δt is not limited by the condition Δt ≤ γ^?1, but only by the properties of the systematic force. It is shown how constraints, such as bond length or bond angle constraints, can be incorporated in the computational scheme. It is argued that the third-order Verlet-type SD algorithm proposed earlier may be simplified without loosing its third-order accuracy. The leap-frog SD algorithm is proven to be equivalent to the verlet-type SD algorithm. Both these SD algorithms are slightly more economical on computer storage than the Beeman-type SD algorithm.     

Innovations in Generation and Analysis of 2D [C,H] COSY NMR Spectra for Metabolic Flux Analysis Purposes

2D [¹³C,¹H] COSY NMR is used by the metabolic engineering community for determining ¹³C–¹³C connectivities in intracellular compounds that contain information regarding the steady-state fluxes in cellular metabolism. This paper proposes innovations in the generation and analysis of these specific NMR spectra. These include a computer tool that allows accurate determination of the relative peak areas and their complete covariance matrices even in very complex spectra. Additionally, a method is introduced for correcting the results for isotopic non-steady-state conditions. The proposed methods were applied to measured 2D [¹³C,¹H] COSY NMR spectra. Peak intensities in a one-dimensional section of the spectrum are frequently not representative for relative peak volumes in the two-dimensional spectrum. It is shown that for some spectra a significant amount of additional information can be gained from long-range ¹³C–¹³C scalar couplings in 2D [¹³C,¹H] COSY NMR spectra. Finally, the NMR resolution enhancement by dissolving amino acid derivatives in a nonpolar solvent is demonstrated.     

Fully automated sequence-specific resonance assignments of hetero- nuclear protein spectra

Full automation of the analysis of spectra is a prerequisite for high-throughput NMR studies in structural or functional genomics. Sequence-specific assignments often form the major bottleneck. Here, we present a procedure that yields nearly complete backbone and side chain resonance assignments starting from a set of heteronuclear three-dimensional spectra. Neither manual intervention, e.g., to correct lists obtained from peak picking before feeding these to an assignment program, nor protein-specific information, e.g., structures of homologous proteins, were required. By combining two earlier published procedures, AUTOPSY [Koradi et al. (1998) J. Magn. Reson., 135, 288–297] and GARANT [Bartels et al. (1996) J. Biomol. NMR, 7, 207–213], with a new program, PICS, all necessary steps from spectra analyses to sequence-specific assignments were performed fully automatically. Characteristic features of the present approach are a flexible design allowing as input almost any combination of NMR spectra, applicability to side chains, robustness with respect to parameter choices (such as noise levels) and reproducibility. In this study, automated resonance assignments were obtained for the 14 kD blue copper protein azurin from P. aeruginosa using five spectra: HNCACB, HNHA, HCCH-TOCSY, ¹⁵N-NOESY-HSQC and ¹³C-NOESY-HSQC. Peaks from these three-dimensional spectra were filtered and calibrated with the help of two two-dimensional spectra: ¹⁵N-HSQC and ¹³C-HSQC. The rate of incorrect assignments is less than 1.5% for backbone nuclei and about 3.5% when side chain protons are also considered.     

1H NMR studies on ferricytochromec 3 fromDesulfovibrio vulgaris Miyazaki F and its interaction with ferredoxin I

Summary The¹H NMR signals of the heme methyl, propionate and related chemical groups of cytochromec ₃ fromDesulfovibrio vulgaris Miyazaki F (D.v. MF) were site-specifically assigned by means of ID NOE, 2D DQFCOSY and 2D TOCSY spectra. They were consistent with the site-specific assignments of the hemes with the highest and second-lowest redox potentials reported by Fan et al. (Biochemistry,29 (1990) 2257–2263). The site-specific heme assignments were also supported by NOE between the methyl groups of these hemes and the side chain of Val¹⁸. All the results contradicted the heme assignments forD.v. MF cytochromec ₃ made on the basis of electron spin resonance (Gayda et al. (1987)FEBS Lett.,217 57–61). Based on these assignments, the interaction of cytochromec ₃ withD.v. MF ferredoxin I was investigated by NMR. The major interaction site of cytochromec ₃ was identified as the heme with the highest redox potential, which is surrounded by the highest density of positive charges. The stoichiometry and association constant were two cytochromec ₃ molecules per monomer of ferredoxin I and 10⁸ M^–2 (at 53 mM ionic strength and 25°C), respectively.     

Simultaneous acquisition of [13C,15N]- and [15N,15N]-separated 4D gradient-enhanced NOESY spectra in proteins

Summary The simultaneous acquisition of a 4D gradient-enhanced and sensitivity-enhanced [¹³C,¹⁵N]/[¹⁵N,¹⁵N]-separated NOESY is presented for the 74-residue [¹³C,¹⁵N]-labeled N-terminal SH3 domain of mGrb2 complexed with a peptide gragment from mSOS-2 in 90% H₂O. The method readily accommodates different ¹³C and ¹⁵N spectral widths, but requires that the same number of increments be collected for both ¹³C and ¹⁵N in the simultaneous dimension (F₂). For purposes of display and analysis, the two 4D spectra can be deconvolved during the processing stage by the appropriate linear combination of separately stored FIDs. Compared to collecting each of these two 4D data sets separately, the presented method is a factor (2)^1/2 more efficient in sensitivity per unit acquisition time. The interleaved nature of this method may also lead to improved peak registration between the two 4D spectra.     

Doubly sensitivity-enhanced 3D TOCSY-HSQC

Summary Recently, strategies for double sensitivity enhancement in heteronuclear three-dimensional NMR experiments were introduced (Krishnamurthy, V.V. (1995) J. Magn. Reson., B106, 170–177; Sattler et al. (1995) J. Biomol. NMR, 6, 11–22; Sattler et al. (1995) J. Magn. Reson., B108, 235–242). Since a sensitivity enhancement of a factor 2^1/2 can be achieved for each indirect dimension, nD spectra can theoretically be enhanced up to a factor of 2^((n-1)/2). We propose and analyze a doubly enhanced three-dimensional TOCSY-HSQC sequence. The application of the doubly enhanced three-dimensional {¹⁵N, ¹H} TOCSY-HSQC sequence is shown for uniformly ¹³C-/¹⁵N- and ¹⁵N-labeled samples of the relatively large Azotobacter vinelandii flavodoxin II (179 amino acids). The main factors that contribute to the final signal-to-noise enhancement have been systematically investigated. The sensitivity enhancement obtained for the doubly enhanced TOCSY-HSQC pulse sequence as compared to the standard (unenhanced) version is close to the theoretically expected factor of two.     

New modes of data partitioning based on PARS peak alignment for improved multivariate biomarker/biopattern detection in 1H-NMR spectroscopic metabolic profiling of urine

This paper addresses the possibility of mathematically partition and process urine ¹H-NMR spectra to enhance the efficiency of the subsequent multivariate data analysis in the context of metabolic profiling of a toxicity study. We show that by processing the NMR data with the peak alignment using reduced set mapping (PARS) algorithm and the use of sparse representation of the data results in the information contained in the original NMR data being preserved with retained resolution but free of the problem of peak shifts. We can now describe a method for differential expression analysis of NMR spectra by using prior knowledge, i.e., the onset of dosing, a partitioning not possible to achieve using raw or bucketed data. In addition we also outline a scheme for soft removal of “biological noise” from the aligned data: exhaustive bio-noise subtraction (EBS). The result is a straightforward protocol for detection of peaks that appear as a consequence of the drug response. In other words, it is possible to elucidate peak origin, either from endogenous substances or from the administered drug/biomarkers. The partition of data originating from the normally regulating metabolome can, furthermore, be analyzed free of the superimposed biological noise. The proposed protocol results in enhanced interpretability of the processed data, i.e., a more refined metabolic trace, simplification of detection of consistent biomarkers, and a simplified search for metabolic end products of the administered drug.     

An automated flow injection analysis system for on-line monitoring of glucose and L-lactate during lactic acid fermentation in a recycle bioreactor

The study concerns on-line sequential analysis of glucose and L-lactate during lactic acid fermentation using a flow injection analysis (FIA) system. Enzyme electrodes containing immobilized glucose oxidase and L-lactate oxidase were used with an amperometric detection system. A 12-bit data acquisition card with 16 analog input channels and 8 digital output channels was used. The software for data acquisition was developed using Visual C++, and was devised for sampling every hour for sequential analyses of lactate and glucose. The detection range was found to be 2–100 g l^–1 for glucose and 1–60 g l^–1 for L-lactate using the biosensors. This FIA system was used for monitoring glucose utilization and L-lactate production by immobilized cells of Lactobacillus casei subsp. rhamnosus during a lactic acid fermentation process in a recycle batch reactor. After 13 h of fermentation, complete sugar utilization and maximal L-lactate production was observed. A good agreement was observed between analysis data obtained using the biosensors and data from standard analyses of reducing sugar and L-lactate. The biosensors exhibited excellent stability during continuous operation for at least 45 days.