A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

Beta-catenin Status in P?aediatric Medulloblastomas: correlation of immunohistochemical expression with mutational status, genetic profiles, and clinical characteristics.

Medulloblastoma is the most frequent malignant paediatric brain tumour. The activation of the Wnt/beta-catenin pathway occurs in 10-15% of medulloblastomas and has been recently described as a marker for favourable patient outcome. We report a series of 72 paediatric medulloblastomas evaluated for beta-catenin protein expression, CTNNB1 mutations, and comparative genomic hybridization. Gene expression profiles were also available in a subset of 40 cases. Immunostaining of beta-catenin showed extensive nuclear staining (>50% of the tumour cells) in six cases and focal nuclear staining (<10% of cells) in three cases. The other cases either exhibited a signal strictly limited to the cytoplasm (58 cases) or were negative (five cases). CTNNB1 mutations were detected in all beta-catenin extensively nucleopositive cases. The expression profiles of these cases documented strong activation of the Wnt/beta-catenin pathway. Remarkably, five out of these six tumours showed a complete loss of chromosome 6. In contrast, cases with focal nuclear beta-catenin staining, as well as tumours with negative or cytoplasmic staining, never demonstrated CTNNB1 mutation, Wnt/beta-catenin pathway activation or chromosome 6 loss. Patients with extensive nuclear staining were significantly older at diagnosis and were in continuous complete remission after a mean follow-up of 75.7 months (range 27.5-121.2 months) from diagnosis. All three patients with focal nuclear staining of beta-catenin died within 36 months from diagnosis. Altogether, these data confirm and extend previous observations that CTNNB1-mutated tumours represent a distinct molecular subgroup of medulloblastomas with favourable outcome, indicating that therapy de-escalation should be considered. International consensus on the definition criteria of this distinct medulloblastoma subgroup should be achieved.

EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

Review of CT features of saprophytic aspergillosis (Aspergilloma) and correlation with final diagnosis : P21

Purpose: 1. To review Ct features suggestive of saprophytic aspergillosis (aspergilloma) and to correlate them with the final pathological results. 2. To illustrate the wide range of differential diagnosis. Methods and materials: The electronic database of our department from 1995 to 2007 revealed CT reports of 48 patients that had been considered very suggestive of aspergilloma. Two radiologists with 6 and 12 years experience in thoracic radiology jointly reviewed the corresponding CT features including ancillary findings and the underlying lung diseases and correlated them with the final pathological diagnosis. Results: Forty patients could be included in the study (12 women, mean age 52), while in 8 patients there was no adequate clinical follow-up. In 17 patients the diagnosis "mycetoma" due to aspergillus fumigatus infection was confirmed, either by surgery, biopsy or bronchoscopy. In 23 patients, differential diagnoses were found, such as cavitating bronchial carcinoma (n = 7), bacterial abscess (n = 3), typical (n = 2) and atypical (n = 2) tuberculosis, as well as inflammatory changes due to mucoviscidosis (n = 1), Wegener's disease (n = 1) or chronic obstructive pulmonary disease (n = 3). Fibromyxoide hamartoma, lung infarction and bronchomucocele were responsible for the typical CT feature in one patient each. Conclusion: 1. The typical CT feature suggesting mycetoma is softtissue proliferation within a pre-existing wall-thickened lung cavity, oten even considered "pathognomonic". However, this diagnosis was finally confirmed by surgery or laboratory findings in less than 50% of patients only. 2. Since differential diagnoses are very large, not only including cavitating lung cancer and tuberculosis, the individual underlying lung disease needs strongly being taken into account often giving the best clue for the correct diagnosis.

Partial correlation : Network analysis

The final year project came to us as an opportunity to get involved in a topic which has appeared to be attractive during the learning process of majoring in economics: statistics and its application to the analysis of economic data, i.e. econometrics.Moreover, the combination of econometrics and computer science is a very hot topic nowadays, given the Information Technologies boom in the last decades and the consequent exponential increase in the amount of data collected and stored day by day. Data analysts able to deal with Big Data and to find useful results from it are verydemanded in these days and, according to our understanding, the work they do, although sometimes controversial in terms of ethics, is a clear source of value added both for private corporations and the public sector. For these reasons, the essence of this project is the study of a statistical instrument valid for the analysis of large datasets which is directly related to computer science: Partial Correlation Networks.The structure of the project has been determined by our objectives through the development of it. At first, the characteristics of the studied instrument are explained, from the basic ideas up to the features of the model behind it, with the final goal of presenting SPACE model as a tool for estimating interconnections in between elements in large data sets. Afterwards, an illustrated simulation is performed in order to show the power and efficiency of the model presented. And at last, the model is put into practice by analyzing a relatively large data set of real world data, with the objective of assessing whether the proposed statistical instrument is valid and useful when applied to a real multivariate time series. In short, our main goals are to present the model and evaluate if Partial Correlation Network Analysis is an effective, useful instrument and allows finding valuable results from Big Data.As a result, the findings all along this project suggest the Partial Correlation Estimation by Joint Sparse Regression Models approach presented by Peng et al. (2009) to work well under the assumption of sparsity of data. Moreover, partial correlation networks are shown to be a very valid tool to represent cross-sectional interconnections in between elements in large data sets.The scope of this project is however limited, as there are some sections in which deeper analysis would have been appropriate. Considering intertemporal connections in between elements, the choice of the tuning parameter lambda, or a deeper analysis of the results in the real data application are examples of aspects in which this project could be completed.To sum up, the analyzed statistical tool has been proved to be a very useful instrument to find relationships that connect the elements present in a large data set. And after all, partial correlation networks allow the owner of this set to observe and analyze the existing linkages that could have been omitted otherwise.

Extensor carpi ulnaris accessory tendinous slips (ECUATS): MRI features at 3.0T and correlation with tendinosis : 49

Purpose: To assess the visibility and the features of ECUATS on 3.0-T MRI studies, and evaluate their correlation with tendinosis. Methods and materials: Our retrospective study was approved by IRB, with waiver of informed consent. Fifty wrist MRI and 48 MR arthrographies from 98 patients (55 males, 43 females, mean age 42.3 years) performed between January and November 2009 on 3.0-T units were reviewed. Images (transverse T1, T2, FS Gd T1 and VIBE) were independently analyzed by two radiologists, and a consensus reached with a third reader in case of disagreement. The visibility of ECUATS was assessed on each available transverse sequence. When present, ECUATS' origins, diameters and insertions were noted. ECU tendinosis was also evaluated. Inter-rater agreement was assessed using Cohen's Kappa coefficient. Results: ECUATS observed prevalence was 23.5% (23/98). ECUATS were more frequently noted on the VIBE sequence, with a good inter-rater agreement (Kappa = 0.72). Origins were noted in 95.7% of cases: 3 were at the level of, and 20 distal to ECU subsheath. Insertions were seen in 43.5%: 2 were on 5th metacarpal bone, 8 on extensor apparatus of 5th finger. ECUATS mean shortest and longest diameters were 0.54 and 0.85 mm respectively. ECU tendinosis was statistically more frequently noted in patients with ECUATS (p <0.05). Conclusion: ECUATS are readily visible on 3.0-T MRI studies, especially on transverse GRE VIBE images. ECU tendinosis is more frequently noted in patients bearing ECUATS.

Empirical study of correlated survival times for recurrent events with proportional hazards margins and the effect of correlation and censoring.

Background: In longitudinal studies where subjects experience recurrent incidents over a period of time, such as respiratory infections, fever or diarrhea, statistical methods are required to take into account the within-subject correlation. Methods: For repeated events data with censored failure, the independent increment (AG), marginal (WLW) and conditional (PWP) models are three multiple failure models that generalize Cox"s proportional hazard model. In this paper, we revise the efficiency, accuracy and robustness of all three models under simulated scenarios with varying degrees of within-subject correlation, censoring levels, maximum number of possible recurrences and sample size. We also study the methods performance on a real dataset from a cohort study with bronchial obstruction. Results: We find substantial differences between methods and there is not an optimal method. AG and PWP seem to be preferable to WLW for low correlation levels but the situation reverts for high correlations. Conclusions: All methods are stable in front of censoring, worsen with increasing recurrence levels and share a bias problem which, among other consequences, makes asymptotic normal confidence intervals not fully reliable, although they are well developed theoretically.

Management of Pharyngeal Fistulas After Anterior Cervical Spine Surgery: A Treatment Algorithm for Severe Complications.

STUDY DESIGN:: Retrospective database- query to identify all anterior spinal approaches. OBJECTIVES:: To assess all patients with pharyngo-cutaneous fistulas after anterior cervical spine surgery. SUMMARY OF BACKGROUND DATA:: Patients treated in University of Heidelberg Spine Medical Center, Spinal Cord Injury Unit and Department of Otolaryngology (Germany), between 2005 and 2011 with the diagnosis of pharyngo-cutaneous fistulas. METHODS:: We conducted a retrospective study on 5 patients between 2005 and 2011 with PCF after ACSS, their therapy management and outcome according to radiologic data and patient charts. RESULTS:: Upon presentation 4 patients were paraplegic. 2 had PCF arising from one piriform sinus, two patients from the posterior pharyngeal wall and piriform sinus combined and one patient only from the posterior pharyngeal wall. 2 had previous unsuccessful surgical repair elsewhere and 1 had prior radiation therapy. In 3 patients speech and swallowing could be completely restored, 2 patients died. Both were paraplegic. The patients needed an average of 2-3 procedures for complete functional recovery consisting of primary closure with various vascularised regional flaps and refining laser procedures supplemented with negative pressure wound therapy where needed. CONCLUSION:: Based on our experience we are able to provide a treatment algorithm that indicates that chronic as opposed to acute fistulas require a primary surgical closure combined with a vascularised flap that should be accompanied by the immediate application of a negative pressure wound therapy. We also conclude that particularly in paraplegic patients suffering this complication the risk for a fatal outcome is substantial.

Proteomic and transcriptomic profiling of Staphylococcus aureus surface LPXTG-proteins: correlation with agr genotypes and adherence phenotypes.

Staphylococcus aureus infections involve numerous adhesins and toxins, which expression depends on complex regulatory networks. Adhesins include a family of surface proteins covalently attached to the peptidoglycan via a conserved LPXTG motif. Here we determined the protein and mRNA expression of LPXTG-proteins of S. aureus Newman in time-course experiments, and their relation to fibrinogen adherence in vitro. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa), and fibrinogen-binding protein A (ClfA), as well as during growth in iron-rich or iron-poor media. Surface proteins were recovered by trypsin-shaving of live bacteria. Released peptides were analyzed by liquid chromatography coupled to tandem mass-spectrometry. To unambiguously identify peptides unique to LPXTG-proteins, the analytical conditions were refined using a reference library of S. aureus LPXTG-proteins heterogeneously expressed in surrogate Lactococcus lactis. Transcriptomes were determined by microarrays. Sixteen of the 18 LPXTG-proteins present in S. aureus Newman were detected by proteomics. Nine LPXTG-proteins showed a bell-shape agr-like expression that was abrogated in agr-negative mutants including Spa, fibronectin-binding protein A (FnBPA), ClfA, iron-binding IsdA, and IsdB, immunomodulator SasH, functionally uncharacterized SasD, biofilm-related SasG and methicillin resistance-related FmtB. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr- mutant, whereas all other LPXTG-proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in fibrinogen-adherence tests during late growth (24 h), whereas it remained poorly detected by proteomics. On the other hand, iron-regulated IsdA-B-C increased their protein expression by >10-times in iron-poor conditions. Thus, proteomic, transcriptomic, and adherence-phenotype demonstrated differential profiles in S. aureus. Moreover, trypsin peptide signatures suggested differential protein domain exposures in various environments, which might be relevant for anti-adhesin vaccines. A comprehensive understanding of the S. aureus physiology should integrate all three approaches.

Correlation of Locally-Based Performance of Asphalts with Their Physicochemical Parameters, Progress Report 2, HR-298, 1989

This project was undertaken to study the relationships between the performance of locally available asphalts and their physicochemical properties under Iowa conditions with the ultimate objective of development of a locally and performance-based asphalt specification for durable pavements. Physical and physicochemical tests were performed on three sets of asphalt samples including: (a) twelve samples from local asphalt suppliers and their TFOT residues, (b) six core samples of known service records, and (c) a total of 79 asphalts from 10 pavement projects including original, lab aged and recovered asphalts from field mixes, as well as from lab aged mixes. Tests included standard rheological tests, HP-GPC and TMA. Some specific viscoelastic tests (at 5 deg C) were run on b samples and on some a samples. DSC and X-ray diffraction studies were performed on a and b samples. Furthermore, NMR techniques were applied to some a, b and c samples. Efforts were made to identify physicochemical properties which are correlated to physical properties known to affect field performance. The significant physicochemical parameters were used as a basis for an improved performance-based trial specification for Iowa to ensure more durable pavements.

Correlation of physical and chemical attributes of soil with sugarcane yield

The objective of this work was to evaluate the correlation between sugarcane yield and some physical and chemical attributes of soil. For this, a 42‑ha test area in Araras, SP, Brazil, was used. Soil properties were determined from samples collected at the beginning of the 2003/2004 harvest season, using a regular 100x100 m grid. Yield assessment was done with a yield monitor (Simprocana). Correlation analyses were performed between sugarcane yield and the following soil properties: pH, pH CaCl2, N, C, cone index, clay content, soil organic matter, P, K, Ca, Mg, H+AL, cation exchange capacity, and base saturation. Correlation coefficients were respectively ‑0.05, ‑0.29, 0.33, 0.41, ‑0.27, 0.22, 0.44, ‑0.24, trace, ‑0.06, 0.01, 0.32, 0.14, and 0.04. Correlations of chemical and physical attributes of soil with sugarcane yield are weak, and, per se, they are not able to explain sugarcane yield variation, which suggests that other variables, besides soil attributes, should be analysed.