337 research outputs found
Fast interior point solution of quadratic programming problems arising from PDE-constrained optimization
Interior point methods provide an attractive class of approaches for solving linear, quadratic and nonlinear programming problems, due to their excellent efficiency and wide applicability. In this paper, we consider PDE-constrained optimization problems with bound constraints on the state and control variables, and their representation on the discrete level as quadratic programming problems. To tackle complex problems and achieve high accuracy in the solution, one is required to solve matrix systems of huge scale resulting from Newton iteration, and hence fast and robust methods for these systems are required. We present preconditioned iterative techniques for solving a number of these problems using Krylov subspace methods, considering in what circumstances one may predict rapid convergence of the solvers in theory, as well as the solutions observed from practical computations
lp-Recovery of the Most Significant Subspace among Multiple Subspaces with Outliers
We assume data sampled from a mixture of d-dimensional linear subspaces with
spherically symmetric distributions within each subspace and an additional
outlier component with spherically symmetric distribution within the ambient
space (for simplicity we may assume that all distributions are uniform on their
corresponding unit spheres). We also assume mixture weights for the different
components. We say that one of the underlying subspaces of the model is most
significant if its mixture weight is higher than the sum of the mixture weights
of all other subspaces. We study the recovery of the most significant subspace
by minimizing the lp-averaged distances of data points from d-dimensional
subspaces, where p>0. Unlike other lp minimization problems, this minimization
is non-convex for all p>0 and thus requires different methods for its analysis.
We show that if 0<p<=1, then for any fraction of outliers the most significant
subspace can be recovered by lp minimization with overwhelming probability
(which depends on the generating distribution and its parameters). We show that
when adding small noise around the underlying subspaces the most significant
subspace can be nearly recovered by lp minimization for any 0<p<=1 with an
error proportional to the noise level. On the other hand, if p>1 and there is
more than one underlying subspace, then with overwhelming probability the most
significant subspace cannot be recovered or nearly recovered. This last result
does not require spherically symmetric outliers.Comment: This is a revised version of the part of 1002.1994 that deals with
single subspace recovery. V3: Improved estimates (in particular for Lemma 3.1
and for estimates relying on it), asymptotic dependence of probabilities and
constants on D and d and further clarifications; for simplicity it assumes
uniform distributions on spheres. V4: minor revision for the published
versio
Structural similarity assessment for drug sensitivity prediction in cancer
<p>Abstract</p> <p>Background</p> <p>The ability to predict drug sensitivity in cancer is one of the exciting promises of pharmacogenomic research. Several groups have demonstrated the ability to predict drug sensitivity by integrating chemo-sensitivity data and associated gene expression measurements from large anti-cancer drug screens such as NCI-60. The general approach is based on comparing gene expression measurements from sensitive and resistant cancer cell lines and deriving drug sensitivity profiles consisting of lists of genes whose expression is predictive of response to a drug. Importantly, it has been shown that such profiles are generic and can be applied to cancer cell lines that are not part of the anti-cancer screen. However, one limitation is that the profiles can not be generated for untested drugs (i.e., drugs that are not part of an anti-cancer drug screen). In this work, we propose using an existing drug sensitivity profile for drug A as a substitute for an untested drug B given high structural similarities between drugs A and B.</p> <p>Results</p> <p>We first show that structural similarity between pairs of compounds in the NCI-60 dataset highly correlates with the similarity between their activities across the cancer cell lines. This result shows that structurally similar drugs can be expected to have a similar effect on cancer cell lines. We next set out to test our hypothesis that we can use existing drug sensitivity profiles as substitute profiles for untested drugs. In a cross-validation experiment, we found that the use of substitute profiles is possible without a significant loss of prediction accuracy if the substitute profile was generated from a compound with high structural similarity to the untested compound.</p> <p>Conclusion</p> <p>Anti-cancer drug screens are a valuable resource for generating omics-based drug sensitivity profiles. We show that it is possible to extend the usefulness of existing screens to untested drugs by deriving substitute sensitivity profiles from structurally similar drugs part of the screen.</p
Increased RPA1 gene dosage affects genomic stability potentially contributing to 17p13.3 duplication syndrome
A novel microduplication syndrome involving various-sized contiguous duplications in 17p13.3 has recently been described, suggesting that increased copy number of genes in 17p13.3, particularly PAFAH1B1, is associated with clinical features including facial dysmorphism, developmental delay, and autism spectrum disorder. We have previously shown that patient-derived cell lines from individuals with haploinsufficiency of RPA1, a gene within 17p13.3, exhibit an impaired ATR-dependent DNA damage response (DDR). Here, we show that cell lines from patients with duplications specifically incorporating RPA1 exhibit a different although characteristic spectrum of DDR defects including abnormal S phase distribution, attenuated DNA double strand break (DSB)-induced RAD51 chromatin retention, elevated genomic instability, and increased sensitivity to DNA damaging agents. Using controlled conditional over-expression of RPA1 in a human model cell system, we also see attenuated DSB-induced RAD51 chromatin retention. Furthermore, we find that transient over-expression of RPA1 can impact on homologous recombination (HR) pathways following DSB formation, favouring engagement in aberrant forms of recombination and repair. Our data identifies unanticipated defects in the DDR associated with duplications in 17p13.3 in humans involving modest RPA1 over-expression
An inexact dual logarithmic barrier method for solving sparse semidefinite programs
A dual logarithmic barrier method for solving large, sparse semidefinite programs is proposed in this paper. The method avoids any explicit use of the primal variable X and therefore is well-suited to problems with a sparse dual matrix S. It relies on inexact Newton steps in dual space which are computed by the conjugate gradient method applied to the Schur complement of the reduced KKT system. The method may take advantage of low-rank representations of matrices Ai to perform implicit matrix-vector products with the Schur complement matrix and to compute only specific parts of this matrix. This allows the construction of the partial Cholesky factorization of the Schur complement matrix which serves as a good preconditioner for it and permits the method to be run in a matrix-free scheme. Convergence properties of the method are studied and a polynomial complexity result is extended to the case when inexact Newton steps are employed. A Matlab-based implementation is developed and preliminary computational results of applying the method to maximum cut and matrix completion problems are reported
Cytokines Elevated in HIV Elite Controllers Reduce HIV Replication In Vitro and Modulate HIV Restriction Factor Expression
A subset of HIV-infected individuals termed elite controllers (ECs) maintain CD4+ T cell counts and control viral replication in the absence of antiretroviral therapy (ART). Systemic cytokine responses may differentiate ECs from subjects with uncontrolled viral replication or from those who require ART to suppress viral replication. We measured 87 cytokines in four groups of women: 73 ECs, 42 with pharmacologically suppressed viremia (ART), 42 with uncontrolled viral replication (noncontrollers [NCs]), and 48 HIV-uninfected (NEG) subjects. Four cytokines were elevated in ECs but not NCs or ART subjects: CCL14, CCL21, CCL27, and XCL1. In addition, median stromal cell-derived factor-1 (SDF-1) levels were 43% higher in ECs than in NCs. The combination of the five cytokines suppressed R5 and X4 virus replication in resting CD4+ T cells, and individually SDF-1β, CCL14, and CCL27 suppressed R5 virus replication, while SDF-1β, CCL21, and CCL14 suppressed X4 virus replication. Functional studies revealed that the combination of the five cytokines upregulated CD69 and CCR5 and downregulated CXCR4 and CCR7 on CD4+ T cells. The CD69 and CXCR4 effects were driven by SDF-1, while CCL21 downregulated CCR7. The combination of the EC-associated cytokines induced expression of the anti-HIV host restriction factors IFITM1 and IFITM2 and suppressed expression of RNase L and SAMHD1. These results identify a set of cytokines that are elevated in ECs and define their effects on cellular activation, HIV coreceptor expression, and innate restriction factor expression. This cytokine pattern may be a signature characteristic of HIV-1 elite control, potentially important for HIV therapeutic and curative strategies.IMPORTANCE Approximately 1% of people infected with HIV control virus replication without taking antiviral medications. These subjects, termed elite controllers (ECs), are known to have stronger immune responses targeting HIV than the typical HIV-infected subject, but the exact mechanisms of how their immune responses control infection are not known. In this study, we identified five soluble immune signaling molecules (cytokines) in the blood that were higher in ECs than in subjects with typical chronic HIV infection. We demonstrated that these cytokines can activate CD4+ T cells, the target cells for HIV infection. Furthermore, these five EC-associated cytokines could change expression levels of intrinsic resistance factors, or molecules inside the target cell that fight HIV infection. This study is significant in that it identified cytokines elevated in subjects with a good immune response against HIV and defined potential mechanisms as to how these cytokines could induce resistance to the virus in target cells
The intellectual influence of economic journals: quality versus quantity
The evaluation of scientific output has a key role in the allocation of
research funds and academic positions. Decisions are often based on quality indicators
for academic journals, and over the years, a handful of scoring methods have
been proposed for this purpose. Discussing the most prominent methods (de facto
standards) we show that they do not distinguish quality from quantity at article level.
The systematic bias we find is analytically tractable and implies that the methods are
manipulable. We introduce modified methods that correct for this bias, and use them
to provide rankings of economic journals. Our methodology is transparent; our results
are replicable
Tumor classification: molecular analysis meets Aristotle
BACKGROUND: Traditionally, tumors have been classified by their morphologic appearances. Unfortunately, tumors with similar histologic features often follow different clinical courses or respond differently to chemotherapy. Limitations in the clinical utility of morphology-based tumor classifications have prompted a search for a new tumor classification based on molecular analysis. Gene expression array data and proteomic data from tumor samples will provide complex data that is unobtainable from morphologic examination alone. The growing question facing cancer researchers is, "How can we successfully integrate the molecular, morphologic and clinical characteristics of human cancer to produce a helpful tumor classification?" DISCUSSION: Current efforts to classify cancers based on molecular features ignore lessons learned from millennia of experience in biological classification. A tumor classification must include every type of tumor and must provide a unique place for each tumor within the classification. Groups within a classification inherit the properties of their ancestors and impart properties to their descendants. A classification was prepared grouping tumors according to their histogenetic development. The classification is simple (reducing the complexity of information received from the molecular analysis of tumors), comprehensive (providing a place for every tumor of man), and consistent with recent attempts to characterize tumors by cytogenetic and molecular features. The clinical and research value of this historical approach to tumor classification is discussed. SUMMARY: This manuscript reviews tumor classification and provides a new and comprehensive classification for neoplasia that preserves traditional nomenclature while incorporating information derived from the molecular analysis of tumors. The classification is provided as an open access XML document that can be used by cancer researchers to relate tumor classes with heterogeneous experimental and clinical tumor databases
High levels of the adhesion molecule CD44 on leukemic cells generate acute myeloid leukemia relapse after withdrawal of the initial transforming event
Multiple genetic hits are detected in patients with acute myeloid leukemia (AML). To investigate this further, we developed a tetracycline-inducible mouse model of AML, in which the initial transforming event, overexpression of HOXA10, can be eliminated. Continuous overexpression of HOXA10 is required to generate AML in primary recipient mice, but is not essential for maintenance of the leukemia. Transplantation of AML to secondary recipients showed that in established leukemias, ∼80% of the leukemia-initiating cells (LICs) in bone marrow stopped proliferating upon withdrawal of HOXA10 overexpression. However, the population of LICs in primary recipients is heterogeneous, as ∼20% of the LICs induce leukemia in secondary recipients despite elimination of HOXA10-induced overexpression. Intrinsic genetic activation of several proto-oncogenes was observed in leukemic cells resistant to inactivation of the initial transformation event. Interestingly, high levels of the adhesion molecule CD44 on leukemic cells are essential to generate leukemia after removal of the primary event. This suggests that extrinsic niche-dependent factors are also involved in the host-dependent outgrowth of leukemias after withdrawal of HOXA10 overexpression event that initiates the leukemia
Early Elevation of Matrix Metalloproteinase-8 and -9 in Pediatric ARDS Is Associated with an Increased Risk of Prolonged Mechanical Ventilation
BACKGROUND: Matrix metalloproteinases (MMP) -8 and -9 may play key roles in the modulation of neutrophilic lung inflammation seen in pediatric Acute Respiratory Distress Syndrome (ARDS). We aimed to perform a comprehensive analysis of MMP-8 and MMP-9 activity in tracheal aspirates of pediatric ARDS patients compared with non-ARDS controls, testing whether increased MMP-8 and -9 activities were associated with clinical outcomes. METHODS: Tracheal aspirates were collected from 33 pediatric ARDS patients and 21 non-ARDS controls at 48 hours of intubation, and serially for those who remained intubated greater than five days. MMPs, tissue inhibitor of metalloproteinases (TIMPs), human neutrophil elastase (HNE) and myeloperoxidase (MPO) activity were measured by ELISA, and correlated with clinical indicators of disease severity such as PRISM (Pediatric Risk of Mortality) scores, oxygen index (OI), multi-organ system failure (MOSF) and clinical outcome measures including length of intubation, ventilator-free days (VFDs) and mortality in the Pediatric Intensive Care Unit (PICU). RESULTS: Active MMP-9 was elevated early in pediatric ARDS subjects compared to non-ARDS controls. Higher MMP-8 and active MMP-9 levels at 48 hours correlated with a longer course of mechanical ventilation (r = 0.41, p = 0.018 and r = 0.75, p<0.001; respectively) and fewer number of VFDs (r = -0.43, p = 0.013 and r = -0.76, p<0.001; respectively), independent of age, gender and severity of illness. Patients with the highest number of ventilator days had the highest levels of active MMP-9. MMP-9 and to a lesser extent MMP-8 activities in tracheal aspirates from ARDS subjects were sensitive to blockade by small molecule inhibitors. CONCLUSIONS: Higher MMP-8 and active MMP-9 levels at 48 hours of disease onset are associated with a longer duration of mechanical ventilation and fewer ventilator-free days among pediatric patients with ARDS. Together, these results identify early biomarkers predictive of disease course and potential therapeutic targets for this life threatening disease
- …