The accuracy of absolute differential abundance analysis from relative count data

Concerns have been raised about the use of relative abundance data derived from next generation sequencing as a proxy for absolute abundances. For example, in the differential abundance setting, compositional effects in relative abundance data may give rise to spurious differences (false positives) when considered from the absolute perspective. In practice however, relative abundances are often transformed by renormalization strategies intended to compensate for these effects and the scope of the practical problem remains unclear. We used simulated data to explore the consistency of differential abundance calling on renormalized relative abundances versus absolute abundances and find that, while overall consistency is high, with a median sensitivity (true positive rates) of 0.91 and specificity (1—false positive rates) of 0.89, consistency can be much lower where there is widespread change in the abundance of features across conditions. We confirm these findings on a large number of real data sets drawn from 16S metabarcoding, expression array, bulk RNA-seq, and single-cell RNA-seq experiments, where data sets with the greatest change between experimental conditions are also those with the highest false positive rates. Finally, we evaluate the predictive utility of summary features of relative abundance data themselves. Estimates of sparsity and the prevalence of feature-level change in relative abundance data give reasonable predictions of discrepancy in differential abundance calling in simulated data and can provide useful bounds for worst-case outcomes in real data.     

Pulse radiolytic oxidation of beta-carotene with halogenated alkylperoxyl radicals in a quaternary microemulsion: formation of retinol

Pulse radiolytically generated halogenated alkylperoxyl radicals, namely CCl3OO*, CBr3OO*, CHCl2OO*, CHBr2OO*, etc. have been employed for a study of the oxidation of beta-carotene in a quaternary microemulsion. All these halogenated alkylperoxyl radicals produce a radical cation (absorption maximum at 840 nm), either via the initial formation of an adduct (absorption maximum at 740 nm), or by direct reaction. There was a considerable blue shift of both absorption peaks in the present system as compared to those earlier reported in non-polar as well as in micellar media. An additional intense absorption peak at approximately 345 nm was noted at a longer time scale and was stable up to 5 ms. On irradiation with a large number of electron pulses, a stable product with an absorption maximum at 315 nm appeared. On excitation at 315 nm, this product gave a fluorescence spectrum with an emission maximum at 525 nm. The absorption and fluorescence spectra of the stable product compared with those of retinol; this was further confirmed by HPLC analysis. A suitable mechanism has been proposed.     

Microsecond simulations of mdm2 and its complex with p53 yield insight into force field accuracy and conformational dynamics

The p53‐MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein‐ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit‐solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb‐ildn, ff99sb‐ildn‐nmr, ff99sb‐ildn‐phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200‐ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo‐like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands. Proteins 2015; 83:1665–1676. © 2015 Wiley Periodicals, Inc.     

E,E,E-1-(4-Arylamino-4-oxo-2-butenoyl)-3,5-bis(arylidene)-4-piperidones: a topographical study of some novel potent cytotoxins

A series of E,E,E-3,5-bis(arylidene)-1-(4-arylamino-4-oxo-2-butenoyl)-4-piperidones 4 (phenylidene) and 5 (4-nitrophenylidene) were prepared in order to explore the structural features of the N-acyl group which affects the cytotoxic potency. Evaluation toward human Molt 4/C8 and CEM T-lymphocytes revealed that many of the IC(50) figures were submicromolar and lower than melphalan. Marked inhibitory potencies toward murine leukemia L1210 cells were also noted. When evaluated against a panel of human tumor cell lines, three representative compounds in series 4 displayed selective toxicity to leukemia and colon cancer cell lines and were significantly more potent than the reference drug melphalan. Molecular modeling of representative compounds in both series 4 and the analogs, in which the configuration of the olefinic double bond was changed from E to Z (series 3), revealed that the torsion angles of the arylidene aryl rings and locations of the terminal arylaminocarbonyl groups may have contributed to the greater cytotoxic properties displayed in 3. Compounds 4c (3,4-dichlorophenylamino), d (4-methylphenylamino) and 5c (3,4-dichlorophenylamino), d (4-methylphenylamino) inhibited the activity of human N-myristoyltransferase by approximately 50% at concentrations of 50-100 microM. The compounds in series 4 and 5 were well tolerated in a short-term toxicity study in mice.