On the use of partial area under the ROC curve for comparison of two diagnostic tests

Evaluation of diagnostic performance is typically based on the receiver operating characteristic (ROC) curve and the area under the curve (AUC) as its summary index. The partial area under the curve (pAUC) is an alternative index focusing on the range of practical/clinical relevance. One of the problems preventing more frequent use of the pAUC is the perceived loss of efficiency in cases of noncrossing ROC curves. In this paper, we investigated statistical properties of comparisons of two correlated pAUCs. We demonstrated that outside of the classic model there are practically reasonable ROC types for which comparisons of noncrossing concave curves would be more powerful when based on a part of the curve rather than the entire curve. We argue that this phenomenon stems in part from the exclusion of noninformative parts of the ROC curves that resemble straight‐lines. We conducted extensive simulation studies in families of binormal, straight‐line, and bigamma ROC curves. We demonstrated that comparison of pAUCs is statistically more powerful than comparison of full AUCs when ROC curves are close to a “straight line”. For less flat binormal ROC curves an increase in the integration range often leads to a disproportional increase in pAUCs’ difference, thereby contributing to an increase in statistical power. Thus, efficiency of differences in pAUCs of noncrossing ROC curves depends on the shape of the curves, and for families of ROC curves that are nearly straight‐line shaped, such as bigamma ROC curves, there are multiple practical scenarios in which comparisons of pAUCs are preferable.     

In vivo random mutagenesis of streptomycetes using mariner-based transposon Himar1

We report here the in vivo expression of the synthetic transposase gene himar1(a) in Streptomyces coelicolor M145 and Streptomyces albus. Using the synthetic himar1(a) gene adapted for Streptomyces codon usage, we showed random insertion of the transposon into the streptomycetes genome. The insertion frequency for the Himar1-derived minitransposons is nearly 100 % of transformed Streptomyces cells, and insertions are stably inherited in the absence of an antibiotic selection. The minitransposons contain different antibiotic resistance selection markers (apramycin, hygromycin, and spectinomycin), site-specific recombinase target sites (rox and/or loxP), I-SceI meganuclease target sites, and an R6Kγ origin of replication for transposon rescue. We identified transposon insertion loci by random sequencing of more than 100 rescue plasmids. The majority of insertions were mapped to putative open-reading frames on the S. coelicolor M145 and S. albus chromosomes. These insertions included several new regulatory genes affecting S. coelicolor M145 growth and actinorhodin biosynthesis.     

The maize lipoxygenase,ZmLOX10, mediates green leaf volatile,jasmonate and herbivore‐induced plant volatile production for defense against insect attack

Fatty acid derivatives are of central importance for plant immunity against insect herbivores; however, major regulatory genes and the signals that modulate these defense metabolites are vastly understudied, especially in important agro‐economic monocot species. Here we show that products and signals derived from a single Zea mays (maize) lipoxygenase (LOX), ZmLOX10, are critical for both direct and indirect defenses to herbivory. We provide genetic evidence that two 13‐LOXs, ZmLOX10 and ZmLOX8, specialize in providing substrate for the green leaf volatile (GLV) and jasmonate (JA) biosynthesis pathways, respectively. Supporting the specialization of these LOX isoforms, LOX8 and LOX10 are localized to two distinct cellular compartments, indicating that the JA and GLV biosynthesis pathways are physically separated in maize. Reduced expression of JA biosynthesis genes and diminished levels of JA in lox10 mutants indicate that LOX10‐derived signaling is required for LOX8‐mediated JA. The possible role of GLVs in JA signaling is supported by their ability to partially restore wound‐induced JA levels in lox10 mutants. The impaired ability of lox10 mutants to produce GLVs and JA led to dramatic reductions in herbivore‐induced plant volatiles (HIPVs) and attractiveness to parasitoid wasps. Because LOX10 is under circadian rhythm regulation, this study provides a mechanistic link to the diurnal regulation of GLVs and HIPVs. GLV‐, JA‐ and HIPV‐deficient lox10 mutants display compromised resistance to insect feeding, both under laboratory and field conditions, which is strong evidence that LOX10‐dependent metabolites confer immunity against insect attack. Hence, this comprehensive gene to agro‐ecosystem study reveals the broad implications of a single LOX isoform in herbivore defense.     

Assessment of chemical-crosslink-assisted protein structure modeling in CASP13

With the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest-to-date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows.     

Protein structure prediction assisted with sparse NMR data in CASP13

CASP13 has investigated the impact of sparse NMR data on the accuracy of protein structure prediction. NOESY and ¹⁵N-¹H residual dipolar coupling data, typical of that obtained for ¹⁵N,¹³C-enriched, perdeuterated proteins up to about 40 kDa, were simulated for 11 CASP13 targets ranging in size from 80 to 326 residues. For several targets, two prediction groups generated models that are more accurate than those produced using baseline methods. Real NMR data collected for a de novo designed protein were also provided to predictors, including one data set in which only backbone resonance assignments were available. Some NMR-assisted prediction groups also did very well with these data. CASP13 also assessed whether incorporation of sparse NMR data improves the accuracy of protein structure prediction relative to nonassisted regular methods. In most cases, incorporation of sparse, noisy NMR data results in models with higher accuracy. The best NMR-assisted models were also compared with the best regular predictions of any CASP13 group for the same target. For six of 13 targets, the most accurate model provided by any NMR-assisted prediction group was more accurate than the most accurate model provided by any regular prediction group; however, for the remaining seven targets, one or more regular prediction method provided a more accurate model than even the best NMR-assisted model. These results suggest a novel approach for protein structure determination, in which advanced prediction methods are first used to generate structural models, and sparse NMR data is then used to validate and/or refine these models.     

High dimensional mediation analysis with latent variables

We propose a model for high dimensional mediation analysis that includes latent variables. We describe our model in the context of an epidemiologic study for incident breast cancer with one exposure and a large number of biomarkers (i.e., potential mediators). We assume that the exposure directly influences a group of latent, or unmeasured, factors which are associated with both the outcome and a subset of the biomarkers. The biomarkers associated with the latent factors linking the exposure to the outcome are considered “mediators.” We derive the likelihood for this model and develop an expectation‐maximization algorithm to maximize an L1‐penalized version of this likelihood to limit the number of factors and associated biomarkers. We show that the resulting estimates are consistent and that the estimates of the nonzero parameters have an asymptotically normal distribution. In simulations, procedures based on this new model can have significantly higher power for detecting the mediating biomarkers compared with the simpler approaches. We apply our method to a study that evaluates the relationship between body mass index, 481 metabolic measurements, and estrogen‐receptor positive breast cancer.     

Molecular Dynamic Simulations Reveal the Structural Determinants of Fatty Acid Binding to Oxy-Myoglobin

The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a “U” shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.     

Molecular Biophysics of Orai Store-Operated Ca2+ Channels

Upon endoplasmic reticulum Ca²⁺ store depletion, Orai channels in the plasma membrane are activated directly by endoplasmic reticulum-resident STIM proteins to generate the Ca²⁺-selective, Ca²⁺ release-activated Ca²⁺ (CRAC) current. After the molecular identification of Orai, a plethora of functional and biochemical studies sought to compare Orai homologs, determine their stoichiometry, identify structural domains responsible for the biophysical fingerprint of the CRAC current, identify the physiological functions, and investigate Orai homologs as potential therapeutic targets. Subsequently, the solved crystal structure of Drosophila Orai (dOrai) substantiated many findings from structure-function studies, but also revealed an unexpected hexameric structure. In this review, we explore Orai channels as elucidated by functional and biochemical studies, analyze the dOrai crystal structure and its implications for Orai channel function, and present newly available information from molecular dynamics simulations that shed light on Orai channel gating and permeation.     

Inhibition of Yersinia protein tyrosine phosphatase by phosphonate derivatives of calixarenes

Inhibition of Yersinia protein tyrosine phosphatase by calix[4]arene mono-, bis-, and tetrakis(methylenebisphosphonic) acids as well as calix[4]arene and thiacalix[4]arene tetrakis(methylphosphonic) acids have been investigated. The kinetic studies revealed that some compounds in this class are potent competitive inhibitors of Yersinia PTP with inhibition constants in the low micromolar range. The binding modes of macrocyclic phosphonate derivatives in the enzyme active center have been explained using computational docking approach. The results obtained indicate that calix[4]arenes are promising scaffolds for the development of inhibitors of Yersinia PTP.     

Metabolic engineering for high glycerol production by the anaerobic cultures of <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

Glycerol is used by the cosmetic, paint, automotive, food, and pharmaceutical industries and for production of explosives. Currently, glycerol is available in commercial quantities as a by-product from biodiesel production, but the purity and the cost of its purification are prohibitive. The industrial production of glycerol by glucose aerobic fermentation using osmotolerant strains of the yeasts Candida sp. and Saccharomyces cerevisiae has been described. A major drawback of the aerobic process is the high cost of production. For this reason, the development of yeast strains that effectively convert glucose to glycerol anaerobically is of great importance. Due to its ability to grow under anaerobic conditions, the yeast S. cerevisiae is an ideal system for the development of this new biotechnological platform. To increase glycerol production and accumulation from glucose, we lowered the expression of TPI1 gene coding for triose phosphate isomerase; overexpressed the fused gene consisting the GPD1 and GPP2 parts coding for glycerol-3-phosphate dehydrogenase and glycerol-3-phosphate phosphatase, respectively; overexpressed the engineered FPS1 gene that codes for aquaglyceroporin; and overexpressed the truncated gene ILV2 that codes for acetolactate synthase. The best constructed strain produced more than 20 g of glycerol/L from glucose under micro-aerobic conditions and 16 g of glycerol/L under anaerobic conditions. The increase in glycerol production led to a drop in ethanol and biomass accumulation.