Automatic learning for the classification of chemical reactions and in statistical thermodynamics

This Thesis describes the application of automatic learning methods for a) the classification of organic and metabolic reactions, and b) the mapping of Potential Energy Surfaces(PES). The classification of reactions was approached with two distinct methodologies: a representation of chemical reactions based on NMR data, and a representation of chemical reactions from the reaction equation based on the physico-chemical and topological features of chemical bonds. NMR-based classification of photochemical and enzymatic reactions. Photochemical and metabolic reactions were classified by Kohonen Self-Organizing Maps (Kohonen SOMs) and Random Forests (RFs) taking as input the difference between the 1H NMR spectra of the products and the reactants. The development of such a representation can be applied in automatic analysis of changes in the 1H NMR spectrum of a mixture and their interpretation in terms of the chemical reactions taking place. Examples of possible applications are the monitoring of reaction processes, evaluation of the stability of chemicals, or even the interpretation of metabonomic data. A Kohonen SOM trained with a data set of metabolic reactions catalysed by transferases was able to correctly classify 75% of an independent test set in terms of the EC number subclass. Random Forests improved the correct predictions to 79%. With photochemical reactions classified into 7 groups, an independent test set was classified with 86-93% accuracy. The data set of photochemical reactions was also used to simulate mixtures with two reactions occurring simultaneously. Kohonen SOMs and Feed-Forward Neural Networks (FFNNs) were trained to classify the reactions occurring in a mixture based on the 1H NMR spectra of the products and reactants. Kohonen SOMs allowed the correct assignment of 53-63% of the mixtures (in a test set). Counter-Propagation Neural Networks (CPNNs) gave origin to similar results. The use of supervised learning techniques allowed an improvement in the results. They were improved to 77% of correct assignments when an ensemble of ten FFNNs were used and to 80% when Random Forests were used. This study was performed with NMR data simulated from the molecular structure by the SPINUS program. In the design of one test set, simulated data was combined with experimental data. The results support the proposal of linking databases of chemical reactions to experimental or simulated NMR data for automatic classification of reactions and mixtures of reactions. Genome-scale classification of enzymatic reactions from their reaction equation. The MOLMAP descriptor relies on a Kohonen SOM that defines types of bonds on the basis of their physico-chemical and topological properties. The MOLMAP descriptor of a molecule represents the types of bonds available in that molecule. The MOLMAP descriptor of a reaction is defined as the difference between the MOLMAPs of the products and the reactants, and numerically encodes the pattern of bonds that are broken, changed, and made during a chemical reaction. The automatic perception of chemical similarities between metabolic reactions is required for a variety of applications ranging from the computer validation of classification systems, genome-scale reconstruction (or comparison) of metabolic pathways, to the classification of enzymatic mechanisms. Catalytic functions of proteins are generally described by the EC numbers that are simultaneously employed as identifiers of reactions, enzymes, and enzyme genes, thus linking metabolic and genomic information. Different methods should be available to automatically compare metabolic reactions and for the automatic assignment of EC numbers to reactions still not officially classified. In this study, the genome-scale data set of enzymatic reactions available in the KEGG database was encoded by the MOLMAP descriptors, and was submitted to Kohonen SOMs to compare the resulting map with the official EC number classification, to explore the possibility of predicting EC numbers from the reaction equation, and to assess the internal consistency of the EC classification at the class level. A general agreement with the EC classification was observed, i.e. a relationship between the similarity of MOLMAPs and the similarity of EC numbers. At the same time, MOLMAPs were able to discriminate between EC sub-subclasses. EC numbers could be assigned at the class, subclass, and sub-subclass levels with accuracies up to 92%, 80%, and 70% for independent test sets. The correspondence between chemical similarity of metabolic reactions and their MOLMAP descriptors was applied to the identification of a number of reactions mapped into the same neuron but belonging to different EC classes, which demonstrated the ability of the MOLMAP/SOM approach to verify the internal consistency of classifications in databases of metabolic reactions. RFs were also used to assign the four levels of the EC hierarchy from the reaction equation. EC numbers were correctly assigned in 95%, 90%, 85% and 86% of the cases (for independent test sets) at the class, subclass, sub-subclass and full EC number level,respectively. Experiments for the classification of reactions from the main reactants and products were performed with RFs - EC numbers were assigned at the class, subclass and sub-subclass level with accuracies of 78%, 74% and 63%, respectively. In the course of the experiments with metabolic reactions we suggested that the MOLMAP / SOM concept could be extended to the representation of other levels of metabolic information such as metabolic pathways. Following the MOLMAP idea, the pattern of neurons activated by the reactions of a metabolic pathway is a representation of the reactions involved in that pathway - a descriptor of the metabolic pathway. This reasoning enabled the comparison of different pathways, the automatic classification of pathways, and a classification of organisms based on their biochemical machinery. The three levels of classification (from bonds to metabolic pathways) allowed to map and perceive chemical similarities between metabolic pathways even for pathways of different types of metabolism and pathways that do not share similarities in terms of EC numbers. Mapping of PES by neural networks (NNs). In a first series of experiments, ensembles of Feed-Forward NNs (EnsFFNNs) and Associative Neural Networks (ASNNs) were trained to reproduce PES represented by the Lennard-Jones (LJ) analytical potential function. The accuracy of the method was assessed by comparing the results of molecular dynamics simulations (thermal, structural, and dynamic properties) obtained from the NNs-PES and from the LJ function. The results indicated that for LJ-type potentials, NNs can be trained to generate accurate PES to be used in molecular simulations. EnsFFNNs and ASNNs gave better results than single FFNNs. A remarkable ability of the NNs models to interpolate between distant curves and accurately reproduce potentials to be used in molecular simulations is shown. The purpose of the first study was to systematically analyse the accuracy of different NNs. Our main motivation, however, is reflected in the next study: the mapping of multidimensional PES by NNs to simulate, by Molecular Dynamics or Monte Carlo, the adsorption and self-assembly of solvated organic molecules on noble-metal electrodes. Indeed, for such complex and heterogeneous systems the development of suitable analytical functions that fit quantum mechanical interaction energies is a non-trivial or even impossible task. The data consisted of energy values, from Density Functional Theory (DFT) calculations, at different distances, for several molecular orientations and three electrode adsorption sites. The results indicate that NNs require a data set large enough to cover well the diversity of possible interaction sites, distances, and orientations. NNs trained with such data sets can perform equally well or even better than analytical functions. Therefore, they can be used in molecular simulations, particularly for the ethanol/Au (111) interface which is the case studied in the present Thesis. Once properly trained, the networks are able to produce, as output, any required number of energy points for accurate interpolations.

Classification of some penalty methods

Optimization problems arise in science, engineering, economy, etc. and we need to find the best solutions for each reality. The methods used to solve these problems depend on several factors, including the amount and type of accessible information, the available algorithms for solving them, and, obviously, the intrinsic characteristics of the problem. There are many kinds of optimization problems and, consequently, many kinds of methods to solve them. When the involved functions are nonlinear and their derivatives are not known or are very difficult to calculate, these methods are more rare. These kinds of functions are frequently called black box functions. To solve such problems without constraints (unconstrained optimization), we can use direct search methods. These methods do not require any derivatives or approximations of them. But when the problem has constraints (nonlinear programming problems) and, additionally, the constraint functions are black box functions, it is much more difficult to find the most appropriate method. Penalty methods can then be used. They transform the original problem into a sequence of other problems, derived from the initial, all without constraints. Then this sequence of problems (without constraints) can be solved using the methods available for unconstrained optimization. In this chapter, we present a classification of some of the existing penalty methods and describe some of their assumptions and limitations. These methods allow the solving of optimization problems with continuous, discrete, and mixing constraints, without requiring continuity, differentiability, or convexity. Thus, penalty methods can be used as the first step in the resolution of constrained problems, by means of methods that typically are used by unconstrained problems. We also discuss a new class of penalty methods for nonlinear optimization, which adjust the penalty parameter dynamically.

A multidimensional scaling classification of robotic sensors

This chapter analyzes the signals captured during impacts and vibrations of a mechanical manipulator. Eighteen signals are captured and several metrics are calculated between them, such as the correlation, the mutual information and the entropy. A sensor classification scheme based on the multidimensional scaling technique is presented.

An overview of the higher level classification of Pucciniomycotina based on combined analyses of nuclear large and small subunit rDNA sequences

Mycologia, Vol. 98, nº6

Rheo-NMR study of water-based cellulose liquid crystal system at high shear rates

Since long ago cellulosic lyotropic liquid crystals were thought as potential materials to produce fibers competitive with spidersilk or Kevlar, yet the processing of high modulus materials from cellulose-based precursors was hampered by their complex rheological behavior. In this work, by using the Rheo-NMR technique, which combines deuterium NMR with rheology, we investigate the high shear rate regimes that may be of interest to the industrial processing of these materials. Whereas the low shear rate regimes were already investigated by this technique in different works [1-4], the high shear rates range is still lacking a detailed study. This work focuses on the orientational order in the system both under shear and subsequent relaxation process arising after shear cessation through the analysis of deuterium spectra from the deuterated solvent water. At the analyzed shear rates the cholesteric order is suppressed and a flow-aligned nematic is observed which for the higher shear rates develops after certain time periodic perturbations that transiently annihilate the order in the system. During relaxation the flow aligned nematic starts losing order due to the onset of the cholesteric helices leading to a period of very low order where cholesteric helices with different orientations are forming from the aligned nematic, followed in the final stage by an increase in order at long relaxation times corresponding to the development of aligned cholesteric domains. This study sheds light on the complex rheological behavior of chiral nematic cellulose-based systems and opens ways to improve its processing. (C) 2015 Elsevier Ltd. All rights reserved.

Partial classification of Lorenz knots: Syllable permutations of torus knots words

We define families of aperiodic words associated to Lorenz knots that arise naturally as syllable permutations of symbolic words corresponding to torus knots. An algorithm to construct symbolic words of satellite Lorenz knots is defined. We prove, subject to the validity of a previous conjecture, that Lorenz knots coded by some of these families of words are hyperbolic, by showing that they are neither satellites nor torus knots and making use of Thurston's theorem. Infinite families of hyperbolic Lorenz knots are generated in this way, to our knowledge, for the first time. The techniques used can be generalized to study other families of Lorenz knots.

Biochemistry of grape berries: post-genomics approaches to uncover the effects of water deficits on ripening

Dissertation presented to obtain the Ph.D degree in Biochemistry, Plant Physiology

The effect of water inorganic matrix in ibuprofen adsorption onto activated carbon for water and wastewater treatment

Dissertação para obtenção do Grau de Mestre em Engenharia Química e Bioquímica

Classification of varieties of pasteurella pestis strains isolated in the interior of the state of Pernambuco-Brazil

Studies were made on the biochemical behavior of 100 strains of P.pestis isolated in Northeastern Brazil with regard to production of nitrous acid, reduction of nitrates to nitrltes, and aciáification of glycerol. Results showed that 98 strains can be classified as "orientalis variety", while the remaining two could not be included in any of the existing "varieties".

Unsupervised classification of remote sensing images combining Self Organizing Maps and segmentation techniques

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

MicroRNAs as modulators of water deficit responses in Medicago truncatula

Dissertação para a obtenção do grau de doutor em Biologia pelo Instituto de Tecnologia Química e Biológica. Universidade Nova de Lisboa

Effects of water deprivation on renal hydroelectrolytic excretion in chronically Trypanosoma cruzi-infected rats

The effect of an 8 hour-period of water deprivation on fluid and electrolyte renal excretion was investigated in male Wistar rats infected with the strain São Felipe (12SF) of Trypanosoma cruzi, in comparison with age and sex matched non-infected controls. The median percent reductions in the urinary flow (-40% v -63%) and excretion ofsodium (-57% v-79%) were smaller in chagasic than in control rats, respectively. So, chagasic rats excreted more than controls. On the other hand, the median percent decrement in the clearance of creatinine was higher in chagasic (-51%) than in controls (-39%). Thus, chagasic rats showed some disturbed renal hydroelectrolytic responses to water deprivation, expressed by smaller conservation, or higher excretion of water and sodium in association with smaller glomerularfiltration rate. This fact denoted an elevation in the fractional excretion of sodium and water.

Retrospective histopathological classification of 1,108 skin biopsies from patients clinically suspected of having leprosy from Bahia, Northeast Brazil

We report a retrospective histopathological classification carried out under laboratory conditions by the method of Ridley & Jopling of 1,108 skin biopsies from patients clinically suspected of having leprosy from Bahia, Northeast Brazil.

Storage of water molecules into biomimetic heterostructures: the role of roughness

The development of devices based on heterostructured thin films of biomolecules conveys a huge contribution on biomedical field. However, to achieve high efficiency of these devices, the storage of water molecules into these heterostructures, in order to maintain the biological molecules hydrated, is mandatory. Such hydrated environment may be achieved with lipids molecules which have the ability to rearrange spontaneously into vesicles creating a stable barrier between two aqueous compartments. Yet it is necessary to find conditions that lead to the immobilization of whole vesicles on the heterostructures. In this work, the conditions that govern the deposition of open and closed liposomes of 1.2-dipalmitoyl-sn-Glycero-3-[Phospho-rac-(1-glycerol)] (sodium Salt) (DPPG) onto polyelectrolytes cushions prepared by the layer-by-layer (LbL) method were analyzed. Electronic transitions of DPPG molecules as well as absorption coefficients were obtained by vacuum ultraviolet spectroscopy, while the elemental composition of the heterostructures was characterized by x-ray photoelectron spectroscopy (XPS). The presence of water molecules in the films was inferred by XPS and infrared spectroscopy. Quartz crystal microbalance (QCM) data analysis allowed to conclude that, in certain cases, the DPPG adsorbed amount is dependent of the bilayers number already adsorbed. Moreover, the adsorption kinetics curves of both adsorbed amount and surface roughness allowed to determine the kinetics parameters that are related with adsorption processes namely, electrostatic forces, liposomes diffusion and lipids re-organization on surface. Scaling exponents attained from atomic force microscopy images statistical analysis demonstrate that DPPG vesicles adsorption mechanism is ruled by the diffusion Villain model confirming that adsorption is governed by electrostatic forces. The power spectral density treatment enabled a thorough description of the accessible surface of the samples as well as of its inner structural properties. These outcomes proved that surface roughness influences the adsorption of DPPG liposomes onto surfaces covered by a polyelectrolyte layer. Thus, low roughness was shown to induce liposome rupture creating a lipid bilayer while high roughness allows the adsorption of whole liposomes. In addition, the fraction of open liposomes calculated from the normalized maximum adsorbed amounts decreases with the cushion roughness increase, allowing us to conclude that the surface roughness is a crucial variable that governs the adsorption of open or whole liposomes. This conclusion is fundamental for the development of well-designed sensors based on functional biomolecules incorporated in liposomes. Indeed, LbL films composed of polyelectrolytes and liposomes with and without melanin encapsulated were successfully applied to sensors of olive oil.