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.

Production and characterization of Chitin-Glucan Complex by Komagataella pastoris

This thesis was focused on the production, extraction and characterization of chitin:β-glucan complex (CGC). In this process, glycerol byproduct from the biodiesel industry was used as carbon source. The selected CGC producing yeast was Komagataella pastoris (formerly known as Pichia pastoris), due the fact that to achieved high cell densities using as carbon source glycerol from the biodiesel industry. Firstly, a screening of K. pastoris strains was performed in shake flask assays, in order to select the strain of K. pastoris with better performance, in terms of growth, using glycerol as a carbon source. K. pastoris strain DSM 70877 achieved higher final cell densities (92-97 g/l), using pure glycerol (99%, w/v) and in glycerol from the biodiesel industry (86%, w/v), respectively, compared to DSM 70382 strain (74-82 g/l). Based on these shake flask assays results, the wild type DSM 70877 strain was selected to proceed for cultivation in a 2 l bioreactor, using glycerol byproduct (40 g/l), as sole carbon source. Biomass production by K. pastoris was performed under controlled temperature and pH (30.0 ºC and 5.0, respectively). More than 100 g/l biomass was obtained in less than 48 h. The yield of biomass on a glycerol basis was 0.55 g/g during the batch phase and 0.63 g/g during the fed-batch phase. In order to optimize the downstream process, by increasing extraction and purification efficiency of CGC from K. pastoris biomass, several assays were performed. It was found that extraction with 5 M NaOH at 65 ºC, during 2 hours, associated to neutralization with HCl, followed by successive washing steps with deionised water until conductivity of ≤20μS/cm, increased CGC purity. The obtained copolymer, CGCpure, had a chitin:glucan molar ratio of 25:75 mol% close to commercial CGC samples extracted from A. niger mycelium, kiOsmetine from Kitozyme (30:70 mol%). CGCpure was characterized by solid-state Nuclear Magnetic Resonance (NMR) spectroscopy and Differential Scanning Calorimetry (DCS), revealing a CGC with higher purity than a CGC commercial (kiOsmetine). In order to optimize CGC production, a set of batch cultivation experiments was performed to evaluate the effect of pH (3.5–6.5) and temperature (20–40 ºC) on the specific cell growth rate, CGC production and polymer composition. Statistical tools (response surface methodology and central composite design) were used. The CGC content in the biomass and the volumetric productivity (rp) were not significantly affected within the tested pH and temperature ranges. In contrast, the effect of pH and temperature on the CGC molar ratio was more pronounced. The highest chitin: β-glucan molar ratio (> 14:86) was obtained for the mid-range pH (4.5-5.8) and temperatures (26–33 ºC). The ability of K. pastoris to synthesize CGC with different molar ratios as a function of pH and temperature is a feature that can be exploited to obtain tailored polymer compositions.(...)

Membrane design and characterisation by incorporation of ionic liquids

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia Química, especialidade de Operações Unitárias e Fenómenos de Transferência, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

The making of saudade. National identity and ethnic psychology in Portugal

Roots and rituals.The construction of ethnic identities, Ton Dekker, John Helsloot Carla Wijers editors, p. 267-268; Selected papers of the 6TH SIEF conference on 'Roots & rituals', Amsterdam 20-25 April 1998.

Purification and Preliminary Characterization of Tetraheme Cytochrome and Adenylylsulf te Reductase from the Peptidolytic Sulfate-Reducing Bacterium Desulfovibrio aminophilus DSM 12254

Bioinorganic Chemistry and Applications Volume 3 (2005), Issue 1-2, Pages 81-91

Structural and functional studies of PpcA: a key protein in the electron transfer pathways of Geobacter sulfurreducens

Dissertação para obtenção do Grau de Doutor em Bioquímica, ramo de Biotecnologia

Preparation of gas-filled porous microparticles (GPPs) and microbubbles (MBs) by PGSS method

Dissertation to obtain the Master Degree in Biotechnology

Metabolic flux analysis of neural cell metabolism in primary cultures

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

Eye insight

Dissertação para Obtenção do grau de Mestrado em Arte e Ciência do Vidro

Perspectives of Identity in Julio Cortázar’s Rayuela

My aim was to produce a dissertation, based on Rayuela, which focuses on Cortázar’s questioning of identity. With this objective in mind, I have studied some of the salient elements in the novel that relate to this topic and the subsequent, interrelated, areas of study that arose in doing so. The cities of Paris and Buenos Aires are placed in contrast within the novel and reflect a dichotomy that reflects Oliveira’s condition as a “foreigner,” (more specifically as a South American in Europe). This duality is further reflected in Cortázar’s use of gender, and the development of the notions of active and passive, and an investigation into the traditional modes of thought, symbols, and stereotypes, and an open-ended questioning of their validity. These topics are framed by a notion of Judeo-Christian History that is in many ways flawed and, as such, contrasts with a more intuitive (or “oriental”) perception of reality, which is centred in figures such as la Maga. I found many explicit references to Zen philosophy, and related oriental references, that led me to believe that this area was worthy of further investigation. Rayuela is considered a classic novel within the canon of Spanish language literature. It’s famous “tabla,” like the rules for a game between the writer and the receptor that produce alternative readings, has led to many discussions regarding the novel’s structure and form, and also created a certain amount of polemic with the use of concepts such as the “lector hembra.” Many consider Cortázar a greater short story writer than a novelist, but nevertheless this novel had a profound effect on young readers upon its publication, much to Cortázar’s surprise, and continues to attract readers, dealing as it does with issues that continue to be relevant to many people.

Design of bio-inspired ionic liquids for protein stabilisation

Ionic Liquids (ILs) consist in organic salts that are liquid at/or near room temperature. Since ILs are entirely composed of ions, the formation of ion pairs is expected to be one essential feature for describing solvation in ILs. In recent years, protein - ionic liquid (P-IL) interactions have been the subject of intensive studies mainly because of their capability to promote folding/unfolding of proteins. However, the ion pairs and their lifetimes in ILs in P-IL thematic is dismissed, since the action of ILs is therefore the result of a subtle equilibrium between anion-cation interaction, ion-solvent and ion-protein interaction. The work developed in this thesis innovates in this thematic, once the design of ILs for protein stabilisation was bio-inspired in the high concentration of organic charged metabolites found in cell milieu. Although this perception is overlooked, those combined concentrations have been estimated to be ~300 mM among the macromolecules at concentrations exceeding 300 g/L (macromolecular crowding) and transient ion-pair can naturally occur with a potential specific biological role. Hence the main objective of this work is to develop new bio-ILs with a detectable ion-pair and understand its effects on protein structure and stability, under crowding environment, using advanced NMR techniques and calorimetric techniques. The choline-glutamate ([Ch][Glu]) IL was synthesized and characterized. The ion-pair was detected in water solutions using mainly the selective NOE NMR technique. Through the same technique, it was possible to detect a similar ion-pair promotion under synthetic and natural crowding environments. Using NMR spectroscopy (protein diffusion, HSQC experiments, and hydrogen-deuterium exchange) and differential scanning calorimetry (DSC), the model protein GB1 (production and purification in isotopic enrichment media) it was studied in the presence of [Ch][Glu] under macromolecular crowding conditions (PEG, BSA, lysozyme). Under dilute condition, it is possible to assert that the [Ch][Glu] induces a preferential hydration by weak and non-specific interactions, which leads to a significant stabilisation. On the other hand, under crowding environment, the [Ch][Glu] ion pair is promoted, destabilising the protein by favourable weak hydrophobic interactions , which disrupt the hydration layer of the protein. However, this capability can mitigates the effect of protein crowders. Overall, this work explored the ion-pair existence and its consequences on proteins in conditions similar to cell milieu. In this way, the charged metabolites found in cell can be understood as key for protein stabilisation.

Application of NMR RDC’s, relaxation and self-diffusion for the study of dynamic processes of small molecules in solution

Spin-lattice Relaxation, self-Diffusion coefficients and Residual Dipolar Couplings (RDC’s) are the basis of well established Nuclear Magnetic Resonance techniques for the physicochemical study of small molecules (typically organic compounds and natural products with MW < 1000 Da), as they proved to be a powerful and complementary source of information about structural dynamic processes in solution. The work developed in this thesis consists in the application of the earlier-mentioned NMR techniques to explore, analyze and systematize patterns of the molecular dynamic behavior of selected small molecules in particular experimental conditions. Two systems were chosen to investigate molecular dynamic behavior by these techniques: the dynamics of ion-pair formation and ion interaction in ionic liquids (IL) and the dynamics of molecular reorientation when molecules are placed in oriented phases (alignment media). The application of NMR spin-lattice relaxation and self-diffusion measurements was applied to study the rotational and translational molecular dynamics of the IL: 1-butyl-3-methylimidazolium tetrafluoroborate [BMIM][BF4]. The study of the cation-anion dynamics in neat and IL-water mixtures was systematically investigated by a combination of multinuclear NMR relaxation techniques with diffusion data (using by H1, C13 and F19 NMR spectroscopy). Spin-lattice relaxation time (T1), self-diffusion coefficients and nuclear Overhauser effect experiments were combined to determine the conditions that favor the formation of long lived [BMIM][BF4] ion-pairs in water. For this purpose and using the self-diffusion coefficients of cation and anion as a probe, different IL-water compositions were screened (from neat IL to infinite dilution) to find the conditions where both cation and anion present equal diffusion coefficients (8% water fraction at 25 ºC). This condition as well as the neat IL and the infinite dilution were then further studied by 13C NMR relaxation in order to determine correlation times (c) for the molecular reorientational motion using a mathematical iterative procedure and experimental data obtained in a temperature range between 273 and 353 K. The behavior of self-diffusion and relaxation data obtained in our experiments point at the combining parameters of molar fraction 8 % and temperature 298 K as the most favorable condition for the formation of long lived ion-pairs. When molecules are subjected to soft anisotropic motion by being placed in some special media, Residual Dipolar Couplings (RDCs), can be measured, because of the partial alignment induced by this media. RDCs are emerging as a powerful routine tool employed in conformational analysis, as it complements and even outperforms the approaches based on the classical NMR NOE or J3 couplings. In this work, three different alignment media have been characterized and evaluated in terms of integrity using 2H and 1H 1D-NMR spectroscopy, namely the stretched and compressed gel PMMA, and the lyotropic liquid crystals CpCl/n-hexanol/brine and cromolyn/water. The influence that different media and degrees of alignment have on the dynamic properties of several molecules was explored. Different sized sugars were used and their self-diffusion was determined as well as conformation features using RDCs. The results obtained indicate that no influence is felt by the small molecules diffusion and conformational features studied within the alignment degree range studied, which was the 3, 5 and 6 % CpCl/n-hexanol/brine for diffusion, and 5 and 7.5 % CpCl/n-hexanol/brine for conformation. It was also possible to determine that the small molecules diffusion verified in the alignment media presented close values to the ones observed in water, reinforcing the idea of no conditioning of molecular properties in such media.

NMR studies of tetrahaem cytochomes from D. desulfuricans and S. frigidimarina

Dissertation Presented to obtain the Ph.D. in Biochemistry

Metabolism of primary astrocytes studied by NMR: metabolic trafficking and neuroprotection

Dissertation presented to obtain a Ph. D. degree in Biochemistry by Universidade Nova de Lisboa, Instituto de Tecnologia Química e Biológica.

Gd(III) chelates as NMR probes of protein-protein interactions. Case study: rubredoxin and cytochrome c3

Inorganic Chemistry 50(21):10600-7