An Efficient Local Search for Partial Latin Square Extension Problem

A partial Latin square (PLS) is a partial assignment of n symbols to an nxn grid such that, in each row and in each column, each symbol appears at most once. The partial Latin square extension problem is an NP-hard problem that asks for a largest extension of a given PLS. In this paper we propose an efficient local search for this problem. We focus on the local search such that the neighborhood is defined by (p,q)-swap, i.e., removing exactly p symbols and then assigning symbols to at most q empty cells. For p in {1,2,3}, our neighborhood search algorithm finds an improved solution or concludes that no such solution exists in O(n^{p+1}) time. We also propose a novel swap operation, Trellis-swap, which is a generalization of (1,q)-swap and (2,q)-swap. Our Trellis-neighborhood search algorithm takes O(n^{3.5}) time to do the same thing. Using these neighborhood search algorithms, we design a prototype iterated local search algorithm and show its effectiveness in comparison with state-of-the-art optimization solvers such as IBM ILOG CPLEX and LocalSolver.Comment: 17 pages, 2 figure

Avoiding catastrophic failure in correlated networks of networks

Networks in nature do not act in isolation but instead exchange information, and depend on each other to function properly. An incipient theory of Networks of Networks have shown that connected random networks may very easily result in abrupt failures. This theoretical finding bares an intrinsic paradox: If natural systems organize in interconnected networks, how can they be so stable? Here we provide a solution to this conundrum, showing that the stability of a system of networks relies on the relation between the internal structure of a network and its pattern of connections to other networks. Specifically, we demonstrate that if network inter-connections are provided by hubs of the network and if there is a moderate degree of convergence of inter-network connection the systems of network are stable and robust to failure. We test this theoretical prediction in two independent experiments of functional brain networks (in task- and resting states) which show that brain networks are connected with a topology that maximizes stability according to the theory.Comment: 40 pages, 7 figure

Design and rationale of a multi-center, pragmatic, open-label randomized trial of antimicrobial therapy - the study of clinical efficacy of antimicrobial therapy strategy using pragmatic design in Idiopathic Pulmonary Fibrosis (CleanUP-IPF) clinical trial

Compelling data have linked disease progression in patients with idiopathic pulmonary fibrosis (IPF) with lung dysbiosis and the resulting dysregulated local and systemic immune response. Moreover, prior therapeutic trials have suggested improved outcomes in these patients treated with either sulfamethoxazole/ trimethoprim or doxycycline. These trials have been limited by methodological concerns. This trial addresses the primary hypothesis that long-term treatment with antimicrobial therapy increases the time-to-event endpoint of respiratory hospitalization or all-cause mortality compared to usual care treatment in patients with IPF. We invoke numerous innovative features to achieve this goal, including: 1) utilizing a pragmatic randomized trial design; 2) collecting targeted biological samples to allow future exploration of 'personalized' therapy; and 3) developing a strong partnership between the NHLBI, a broad range of investigators, industry, and philanthropic organizations. The trial will randomize approximately 500 individuals in a 1:1 ratio to either antimicrobial therapy or usual care. The site principal investigator will declare their preferred initial antimicrobial treatment strategy (trimethoprim 160 mg/ sulfamethoxazole 800 mg twice a day plus folic acid 5 mg daily or doxycycline 100 mg once daily if body weight is < 50 kg or 100 mg twice daily if ≥50 kg) for the participant prior to randomization. Participants randomized to antimicrobial therapy will receive a voucher to help cover the additional prescription drug costs. Additionally, those participants will have 4-5 scheduled blood draws over the initial 24 months of therapy for safety monitoring. Blood sampling for DNA sequencing and genome wide transcriptomics will be collected before therapy. Blood sampling for transcriptomics and oral and fecal swabs for determination of the microbiome communities will be collected before and after study completion. As a pragmatic study, participants in both treatment arms will have limited in-person visits with the enrolling clinical center. Visits are limited to assessments of lung function and other clinical parameters at time points prior to randomization and at months 12, 24, and 36. All participants will be followed until the study completion for the assessment of clinical endpoints related to hospitalization and mortality events. TRIAL REGISTRATION: ClinicalTrials.gov identifier NCT02759120

An exact expression to calculate the derivatives of position-dependent observables in molecular simulations with flexible constraints

In this work, we introduce an algorithm to compute the derivatives of physical observables along the constrained subspace when flexible constraints are imposed on the system (i.e., constraints in which the hard coordinates are fixed to configuration-dependent values). The presented scheme is exact, it does not contain any tunable parameter, and it only requires the calculation and inversion of a sub-block of the Hessian matrix of second derivatives of the function through which the constraints are defined. We also present a practical application to the case in which the sought observables are the Euclidean coordinates of complex molecular systems, and the function whose minimization defines the constraints is the potential energy. Finally, and in order to validate the method, which, as far as we are aware, is the first of its kind in the literature, we compare it to the natural and straightforward finite-differences approach in three molecules of biological relevance: methanol, N-methyl-acetamide and a tri-glycine peptideComment: 13 pages, 8 figures, published versio

Linkage analysis of longitudinal data and design consideration

BACKGROUND: Statistical methods have been proposed recently to analyze longitudinal data in genetic studies. So far, little attention has been paid to examine the relationship among key factors in genetic longitudinal studies including power, the number of families or sibships, and the number of repeated measures per individual subjects. RESULTS: We proposed a variance component model that extends classic variance component models for a single quantitative trait to mapping longitudinal traits. Our model includes covariate effects and allows genetic effects to vary over time. Using our proposed model, we examined the power, pedigree structures, and sample size through simulation experiments. CONCLUSION: Our simulation results provide useful insights into the study design for genetic, longitudinal studies. For example, collecting a small number of large sibships is much more powerful than collecting a large number of small sibships or increasing the number of repeated measures, when the total number of measurements is comparable

The use of a bayesian hierarchy to develop and validate a co-morbidity score to predict mortality for linked primary and secondary care data from the NHS in England

Background: We have assessed whether the linkage between routine primary and secondary care records provided an opportunity to develop an improved population based co-morbidity score with the combined information on co-morbidities from both health care settings. Methods: We extracted all people older than 20 years at the start of 2005 within the linkage between the Hospital Episodes Statistics, Clinical Practice Research Datalink, and Office for National Statistics death register in England. A random 50% sample was used to identify relevant diagnostic codes using a Bayesian hierarchy to share information between similar Read and ICD 10 code groupings. Internal validation of the score was performed in the remaining 50% and discrimination was assessed using Harrell’s C statistic. Comparisons were made over time, age, and consultation rate with the Charlson and Elixhauser indexes. Results: 657,264 people were followed up from the 1st January 2005. 98 groupings of codes were derived from the Bayesian hierarchy, and 37 had an adjusted weighting of greater than zero in the Cox proportional hazards model. 11 of these groupings had a different weighting dependent on whether they were coded from hospital or primary care. The C statistic reduced from 0.88 (95% confidence interval 0.88–0.88) in the first year of follow up, to 0.85 (0.85–0.85) including all 5 years. When we stratified the linked score by consultation rate the association with mortality remained consistent, but there was a significant interaction with age, with improved discrimination and fit in those under 50 years old (C=0.85, 0.83–0.87) compared to the Charlson (C=0.79, 0.77–0.82) or Elixhauser index (C=0.81, 0.79–0.83). Conclusions: The use of linked population based primary and secondary care data developed a co-morbidity score that had improved discrimination, particularly in younger age groups, and had a greater effect when adjusting for co-morbidity than existing scores

Holocene evolution of a barrier island system, Ria Formosa, South Portugal

Holocene evolution of the Ria Formosa barrier island system was studied through the examination of a large subsurface dataset acquired from 191 boreholes and five seismic refraction profiles. Two boreholes with total depths of 26 and 16.5 m were selected for a multi-proxy detailed laboratory analysis, including mean grain size distribution, organic matter (OM) content, color variation, shell identification, and benthic foraminifera assemblages. Selected cores are thought to be representative of the identified depositional sub-basins. Subsurface age data from 16 AMS C-14 dated samples were plotted against depth and resulted in a coherent age model of sedimentary infill. The system evolution was largely controlled by sediment availability, accommodation space, and Holocene sea level rise, first at a rapid rate of 7 mm/yr from 10 kcal yr BP to 7.25 kcal yr BP, followed by a slowdown to 1.1 mm/yr until present. A conceptual model for the origin and Holocene evolution of the Ria Formosa barrier island system implies three main steps, leading to the present system geomorphology: (1) marine flooding of incised palaeovalleys by the rapid transgression of palaeovalleys in the early Holocene(2) development of a proto-barrier island chain perched on Pleistocene detritic headlands and steeper interfluve areas during the early to middle Holoceneand (3) full development of the barrier islands chain and enclosing of the coastal lagoon, followed by the maturation of the system with subsequent siltation and salt marsh expansion from the middle Holocene until present. The onset of barrier system formation dates back to ca. 8 kcal yr BP, predating previously proposed age.SIHER project [PTDC/CTE-GIX112236/2009]EU Erasmus Mundus Joint Doctorate in Marine and Coastal Management (MACOMA) fellowship grant, under University of AlgarveEU Erasmus Mundus Joint Doctorate in Marine and Coastal Management (MACOMA) fellowship grant, under University of Cadi

Benznidazole biotransformation and multiple targets in <i>Trypanosoma</i> cruzi revealed by metabolomics

&lt;b&gt;Background&lt;/b&gt;&lt;p&gt;&lt;/p&gt; The first line treatment for Chagas disease, a neglected tropical disease caused by the protozoan parasite Trypanosoma cruzi, involves administration of benznidazole (Bzn). Bzn is a 2-nitroimidazole pro-drug which requires nitroreduction to become active, although its mode of action is not fully understood. In the present work we used a non-targeted MS-based metabolomics approach to study the metabolic response of T. cruzi to Bzn.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Methodology/Principal findings&lt;/b&gt;&lt;p&gt;&lt;/p&gt; Parasites treated with Bzn were minimally altered compared to untreated trypanosomes, although the redox active thiols trypanothione, homotrypanothione and cysteine were significantly diminished in abundance post-treatment. In addition, multiple Bzn-derived metabolites were detected after treatment. These metabolites included reduction products, fragments and covalent adducts of reduced Bzn linked to each of the major low molecular weight thiols: trypanothione, glutathione, γ-glutamylcysteine, glutathionylspermidine, cysteine and ovothiol A. Bzn products known to be generated in vitro by the unusual trypanosomal nitroreductase, TcNTRI, were found within the parasites, but low molecular weight adducts of glyoxal, a proposed toxic end-product of NTRI Bzn metabolism, were not detected.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Conclusions/significance&lt;/b&gt;&lt;p&gt;&lt;/p&gt; Our data is indicative of a major role of the thiol binding capacity of Bzn reduction products in the mechanism of Bzn toxicity against T. cruzi