Fourier Acceleration of Langevin Molecular Dynamics

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD), or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.Comment: 11 pages, two illustrations included using graphic

Low-cost rapid detection of rifampicin resistant tuberculosis using bacteriophage in Kampala, Uganda

BACKGROUND: Resistance to anti-tuberculosis drugs is a serious public health problem. Multi-drug resistant tuberculosis (MDR-TB), defined as resistance to at least rifampicin and isoniazid, has been reported in all regions of the world. Current phenotypic methods of assessing drug susceptibility of M. tuberculosis are slow. Rapid molecular methods to detect resistance to rifampicin have been developed but they are not affordable in some high prevalence countries such as those in sub Saharan Africa. A simple multi-well plate assay using mycobacteriophage D29 has been developed to test M. tuberculosis isolates for resistance to rifampicin. The purpose of this study was to investigate the performance of this technology in Kampala, Uganda. METHODS: In a blinded study 149 M. tuberculosis isolates were tested for resistance to rifampicin by the phage assay and results compared to those from routine phenotypic testing in BACTEC 460. Three concentrations of drug were used 2, 4 and 10 μg/ml. Isolates found resistant by either assay were subjected to sequence analysis of a 81 bp fragment of the rpoB gene to identify mutations predictive of resistance. Four isolates with discrepant phage and BACTEC results were tested in a second phenotypic assay to determine minimal inhibitory concentrations. RESULTS: Initial analysis suggested a sensitivity and specificity of 100% and 96.5% respectively for the phage assay used at 4 and 10 μg/ml when compared to the BACTEC 460. However, further analysis revealed 4 false negative results from the BACTEC 460 and the phage assay proved the more sensitive and specific of the two tests. Of the 39 isolates found resistant by the phage assay 38 (97.4%) were found to have mutations predictive of resistance in the 81 bp region of the rpoB gene. When used at 2 μg/ml false resistant results were observed from the phage assay. The cost of reagents for testing each isolate was estimated to be 1.3US$ when testing a batch of 20 isolates on a single 96 well plate. Results were obtained in 48 hours. CONCLUSION: The phage assay can be used for screening of isolates for resistance to rifampicin, with high sensitivity and specificity in Uganda. The test may be useful in poorly resourced laboratories as a rapid screen to differentiate between rifampicin susceptible and potential MDR-TB cases

Comparison of rapid tests for detection of rifampicin-resistant Mycobacterium tuberculosis in Kampala, Uganda

BACKGROUND: Drug resistant tuberculosis (TB) is a growing concern worldwide. Rapid detection of resistance expedites appropriate intervention to control the disease. Several technologies have recently been reported to detect rifampicin resistant Mycobacterium tuberculosis directly in sputum samples. These include phenotypic culture based methods, tests for gene mutations and tests based on bacteriophage replication. The aim of the present study was to assess the feasibility of implementing technology for rapid detection of rifampicin resistance in a high disease burden setting in Africa. METHODS: Sputum specimens from re-treatment TB patients presenting to the Mulago Hospital National TB Treatment Centre in Kampala, Uganda, were examined by conventional methods and simultaneously used in one of the four direct susceptibility tests, namely direct BACTEC 460, Etest, "in-house" phage test, and INNO- Rif.TB. The reference method was the BACTEC 460 indirect culture drug susceptibility testing. Test performance, cost and turn around times were assessed. RESULTS: In comparison with indirect BACTEC 460, the respective sensitivities and specificities for detecting rifampicin resistance were 100% and 100% for direct BACTEC and the Etest, 94% and 95% for the phage test, and 87% and 87% for the Inno-LiPA assay. Turn around times ranged from an average of 3 days for the INNO-LiPA and phage tests, 8 days for the direct BACTEC 460 and 20 days for the Etest. All methods were faster than the indirect BACTEC 460 which had a mean turn around time of 24 days. The cost per test, including labour ranged from $18.60 to$41.92 (USD). CONCLUSION: All four rapid technologies were shown capable of detecting rifampicin resistance directly from sputum. The LiPA proved rapid, but was the most expensive. It was noted, however, that the LiPA test allows sterilization of samples prior to testing thereby reducing the risk of accidental laboratory transmission. In contrast the Etest was low cost, but slow and would be of limited assistance when treating patients. The phage test was the least reproducible test studied with failure rate of 27%. The test preferred by the laboratory personnel, direct BACTEC 460, requires further study to determine its accuracy in real-time treatment decisions in Uganda

Across-arc geochemical variations in the Southern Volcanic Zone, Chile (34.5- 38.0°S): Constraints on Mantle Wedge and Input Compositions

Crustal assimilation (e.g. Hildreth and Moorbath, 1988) and/or subduction erosion (e.g. Stern, 1991; Kay et al., 2005) are believed to control the geochemical variations along the northern portion of the Chilean Southern Volcanic Zone. In order to evaluate these hypotheses, we present a comprehensive geochemical data set (major and trace elements and O-Sr-Nd-Hf-Pb isotopes) from Holocene primarily olivine-bearing volcanic rocks across the arc between 34.5-38.0°S, including volcanic front centers from Tinguiririca to Callaqui, the rear arc centers of Infernillo Volcanic Field, Laguna del Maule and Copahue, and extending 300 km into the backarc. We also present an equivalent data set for Chile Trench sediments outboard of this profile. The volcanic arc (including volcanic front and rear arc) samples primarily range from basalt to andesite/trachyandesite, whereas the backarc rocks are low-silica alkali basalts and trachybasalts. All samples show some characteristic subduction zone trace element enrichments and depletions, but the backarc samples show the least. Backarc basalts have higher Ce/Pb, Nb/U, Nb/Zr, and Ta/Hf, and lower Ba/Nb and Ba/La, consistent with less of a slab-derived component in the backarc and, consequently, lower degrees of mantle melting. The mantle-like δ18O in olivine and plagioclase phenocrysts (volcanic arc = 4.9-5.6 and backarc = 5.0-5.4 per mil) and lack of correlation between δ18O and indices of differentiation and other isotope ratios, argue against significant crustal assimilation. Volcanic arc and backarc samples almost completely overlap in Sr and Nd isotopic composition. High precision (double-spike) Pb isotope ratios are tightly correlated, precluding significant assimilation of older sialic crust but indicating mixing between a South Atlantic Mid Ocean-Ridge Basalt (MORB) source and a slab component derived from subducted sediments and altered oceanic crust. Hf-Nd isotope ratios define separate linear arrays for the volcanic arc and backarc, neither of which trend toward subducting sediment, possibly reflecting a primarily asthenospheric mantle array for the volcanic arc and involvement of enriched Proterozoic lithospheric mantle in the backarc. We propose a quantitative mixing model between a mixed-source, slab-derived melt and a heterogeneous mantle beneath the volcanic arc. The model is consistent with local geodynamic parameters, assuming water-saturated conditions within the slab

Treatment outcomes of new tuberculosis patients hospitalized in Kampala, Uganda: a prospective cohort study.

BACKGROUND: In most resource limited settings, new tuberculosis (TB) patients are usually treated as outpatients. We sought to investigate the reasons for hospitalisation and the predictors of poor treatment outcomes and mortality in a cohort of hospitalized new TB patients in Kampala, Uganda. METHODS AND FINDINGS: Ninety-six new TB patients hospitalised between 2003 and 2006 were enrolled and followed for two years. Thirty two were HIV-uninfected and 64 were HIV-infected. Among the HIV-uninfected, the commonest reasons for hospitalization were low Karnofsky score (47%) and need for diagnostic evaluation (25%). HIV-infected patients were commonly hospitalized due to low Karnofsky score (72%), concurrent illness (16%) and diagnostic evaluation (14%). Eleven HIV uninfected patients died (mortality rate 19.7 per 100 person-years) while 41 deaths occurred among the HIV-infected patients (mortality rate 46.9 per 100 person years). In all patients an unsuccessful treatment outcome (treatment failure, death during the treatment period or an unknown outcome) was associated with duration of TB symptoms, with the odds of an unsuccessful outcome decreasing with increasing duration. Among HIV-infected patients, an unsuccessful treatment outcome was also associated with male sex (P = 0.004) and age (P = 0.034). Low Karnofsky score (aHR = 8.93, 95% CI 1.88 - 42.40, P = 0.001) was the only factor significantly associated with mortality among the HIV-uninfected. Mortality among the HIV-infected was associated with the composite variable of CD4 and ART use, with patients with baseline CD4 below 200 cells/µL who were not on ART at a greater risk of death than those who were on ART, and low Karnofsky score (aHR = 2.02, 95% CI 1.02 - 4.01, P = 0.045). CONCLUSION: Poor health status is a common cause of hospitalisation for new TB patients. Mortality in this study was very high and associated with advanced HIV Disease and no use of ART

Fluorescence-Tagged Transgenic Lines Reveal Genetic Defects in Pollen Growth—Application to the Eif3 Complex

BACKGROUND: Mutations in several subunits of eukaryotic translation initiation factor 3 (eIF3) cause male transmission defects in Arabidopsis thaliana. To identify the stage of pollen development at which eIF3 becomes essential it is desirable to examine viable pollen and distinguish mutant from wild type. To accomplish this we have developed a broadly applicable method to track mutant alleles that are not already tagged by a visible marker gene through the male lineage of Arabidopsis. METHODOLOGY/PRINCIPAL FINDINGS: Fluorescence tagged lines (FTLs) harbor a transgenic fluorescent protein gene (XFP) expressed by the pollen-specific LAT52 promoter at a defined chromosomal position. In the existing collection of FTLs there are enough XFP marker genes to track nearly every nuclear gene by virtue of its genetic linkage to a transgenic marker gene. Using FTLs in a quartet mutant, which yields mature pollen tetrads, we determined that the pollen transmission defect of the eif3h-1 allele is due to a combination of reduced pollen germination and reduced pollen tube elongation. We also detected reduced pollen germination for eif3e. However, neither eif3h nor eif3e, unlike other known gametophytic mutations, measurably disrupted the early stages of pollen maturation. CONCLUSION/SIGNIFICANCE: eIF3h and eIF3e both become essential during pollen germination, a stage of vigorous translation of newly transcribed mRNAs. These data delimit the end of the developmental window during which paternal rescue is still possible. Moreover, the FTL collection of mapped fluorescent protein transgenes represents an attractive resource for elucidating the pollen development phenotypes of any fine-mapped mutation in Arabidopsis

Shower development of particles with momenta from 15 GeV to 150 GeV in the CALICE scintillator-tungsten hadronic calorimeter

We present a study of showers initiated by electrons, pions, kaons, and protons with momenta from 15 GeV to 150 GeV in the highly granular CALICE scintillator-tungsten analogue hadronic calorimeter. The data were recorded at the CERN Super Proton Synchrotron in 2011. The analysis includes measurements of the calorimeter response to each particle type as well as measurements of the energy resolution and studies of the longitudinal and radial shower development for selected particles. The results are compared to Geant4 simulations (version 9.6.p02). In the study of the energy resolution we include previously published data with beam momenta from 1 GeV to 10 GeV recorded at the CERN Proton Synchrotron in 2010.Comment: 35 pages, 21 figures, 8 table