Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

3

A Knudsen flow reactor has been used to quantify functional groups on the surface of seven different types of combustion particle samples: 3 amorphous carbons (FS 101, Printex 60, FW 2), 2 flame soots (hexane soot generated from a rich and a lean diffusion flame), and 2 Diesel particles (SRM 2975, Diesel soot recovered from a Diesel particulate filter). The technique is based on a heterogeneous titration reaction between a probe gas and a specific functional group on the particle surface. Six probe gases have been selected for the quantification of important functional groups: N(CH3)3 for the titration of acidic sites, NH2OH for carbonyl functions of aldehydes and ketones, CF3COOH and HCl for basic sites of different strength, O3 and NO2 for oxidizable groups. The limit of detection was generally well below 1% of a formal monolayer of adsorbed probe gas. Results obtained with N(CH3)3 were higher for the FW 2 amorphous carbon (post-oxidized sample, according to the manufacturer) and the Diesel particles (between 5.2·10 13 and 5.8·10 13 molecule/cm2), indicating a higher state of oxidation than for the other samples (between 1.3·10 12 and 3.7·10 12 molecule/cm2). The ratio of uptakes of CF3COOH and HCl inferred the presence of basic oxides on the particle surface, owing to the larger stability of the acetate compared to the chloride counter ion in the resulting pyrylium salt. The reactivity of the FS 101 amorphous carbon (3.7·10 15 molecule/cm2) and the hexane flame soot (between 1.9·10 15 and 2.7·10 15 molecule/cm2) towards O3 was very high, indicating the presence of a huge amount of oxidizable or reduced groups on the surface of these samples. Besides the quantification of surface functional groups, the kinetics of reactions between particles and probe gases has also been studied. The uptake coefficient γ0 was roughly correlated with the amount of probe gas taken up by the samples. Indeed, the presence of a high density of functional groups led to fast uptake of the probe gas. These different findings indicate that the particle surface appeared multi-functional, with the simultaneous presence of antagonistic functional groups which do not undergo internal chemical reactions, such as acid-base neutralization. Results also point to important differences in the surface reactivity of the samples, depending on the combustion conditions. The relative distribution of the surface functional groups may be a useful indicator for the state of oxidation and the reactivity of the particle surface.

Spontaneous mirror symmetry breaking in the limited enantioselective autocatalysis model: abyssal hydrotermal vents as scenario for the emergence of chirality in prebiotic chemistry

The emergence of chirality in enantioselective autocatalysis for compounds unable to transform according to the Frank-like reaction network is discussed with respect to the controversial limited enantioselectivity (LES) model composed of coupled enantioselective and non-enantioselective autocatalyses. The LES model cannot lead to spontaneous mirror symmetry breaking (SMSB) either in closed systems with a homogeneous temperature distribution or in closed systems with a stationary non-uniform temperature distribution. However, simulations of chemical kinetics in a two-compartment model demonstrate that SMSB may occur if both autocatalytic reactions are spatially separated at different temperatures in different compartments but coupled under the action of a continuous internal flow. In such conditions, the system can evolve, for certain reaction and system parameters, toward a chiral stationary state; that is, the system is able to reach a bifurcation point leading to SMSB. Numerical simulations in which reasonable chemical parameters have been used suggest that an ade- quate scenario for such a SMSB would be that of abyssal hydrothermal vents, by virtue of the typical temper- ature gradients found there and the role of inorganic solids mediating chemical reactions in an enzyme-like role. Key Words: Homochirality Prebiotic chemistry.

The Role of Magnesium in Concrete Deterioration, HR-355, 1994

Concretes with service lives of less than 15 years and those with lives greater than 40 years were studied with petrographic microscope, scanning electron microscope, and electron microprobe to determine why these two groups of concrete exhibit such different degrees of durability under highway conditions. Coarse aggregate used in both types of concrete were from dolomite rock, but investigation revealed that dolomite aggregate in the two groups of concretes were much different in several respects. The poorly-performing aggregate is fine-grained, has numerous euhedral and subhedral dolomite rhombohedra, and has relatively high porosity. Aggregate from durable concrete is coarse-grained, with tightly interlocked crystal fabric, anhedral dolomite boundaries, and low porosity. Aggregate in short service life concrete was found to have undergone pervasive chemical reactions with the cement which produced reaction rims on the boundaries of coarse aggregate particles and in the cement region adjacent to aggregate boundaries. Textural and porosity differences are believed to be chiefly responsible for different service lives of the two groups of concrete. The basic reaction that has occurred in the short service life concretes between coarse aggregate and cement is an alkali-dolomite reaction. In the reaction dolomite from the aggregate reacts with hydroxide ions from the cement to free magnesium ions and carbonate ions, and the magnesium ions precipitate as brucite, Mg(OH)2. Simultaneously with this reaction, a second reaction occurs in which product carbonate ions react with portlandite from the cement to form calcite and hydroxide ions. Crystal growth pressures of newly formed brucite and calcite together with other processes, e.g. hydration state changes of magnesium chloride hydrates, lead to expansion of the concretes with resultant rapid deterioration. According to this model, magnesium from any source, either from reacting dolomite or from magnesium road deicers, has a major role in highway concrete deterioration. Consequently, magnesium deicers should be used with caution, and long-term testing of the effects of magnesium deicers on highway concrete should be implemented to determine their effects on durability.

Expansive Mineral Growth and Concrete Deterioration, HR-384, 1997

A significant question is: What role does newly-formed expansive mineral growth play in the premature deterioration of concrete? These minerals (ettringite and brucite) are formed in cement paste as a result of chemical reactions involving cement and coarse/fine aggregate. Petrographic observations and SEM/EDAX analysis were conducted in order to determine chemical and mineralogical changes in the aggregate and cement paste of samples taken from Iowa concrete highways that showed premature deterioration. Mechanisms involved in deterioration were investigated. A second objective was to investigate whether deicer solutions exacerbate the formation of expansive minerals and concrete deterioration. Magnesium in deicer solutions causes the most severe paste deterioration by forming non-cementitious magnesium silicate hydrate and brucite. Chloride in deicer solutions promotes decalcification of paste and alters ettringite to chloroaluminate. Calcium magnesium acetate (CMA) and magnesium acetate (Mg-acetate) produce the most deleterious effects on concrete, with calcium acetate (Ca-acetate) being much less severe.

Analyzing peptides and proteins by mass spectrometry: principles and applications in proteomics

The study of proteins has been a key element in biomedicine and biotechnology because of their important role in cell functions or enzymatic activity. Cells are the basic unit of living organisms, which are governed by a vast range of chemical reactions. These chemical reactions must be highly regulatedin order to achieve homeostasis. Proteins are polymeric molecules that havetaken on the evolutionary process the role, along with other factors, of controlthese chemical reactions. Learning how proteins interact and control their up anddown regulations can teach us how living cells regulate their functions, as well asthe cause of certain anomalies that occur in different diseases where proteins areinvolved. Mass spectrometry (MS) is an analytical widely used technique to studythe protein content inside the cells as a biomarker point, which describesdysfunctions in diseases and increases knowledge of how proteins are working.All the methodologies involved in these descriptions are integrated in the fieldcalled Proteomics.

Investigation of Rapid Thermal Analysis Procedures for Prediction of the Service Life of PCCP Carbonate Coarse Aggregate, HR-337, 1993

The major objective of this research project was to use thermal analysis techniques in conjunction with x-ray analysis methods to identify and explain chemical reactions that promote aggregate related deterioration in portland cement concrete. Twenty-two different carbonate aggregate samples were subjected to a chemical testing scheme that included: • bulk chemistry (major, minor and selected trace elements) • bulk mineralogy (minor phases concentrated by acid extraction) • solid-solution in the major carbonate phases • crystallite size determinations for the major carbonate phases • a salt treatment study to evaluate the impact of deicer salts Test results from these different studies were then compared to information that had been obtained using thermogravimetric analysis techniques. Since many of the limestones and dolomites that were used in the study had extensive field service records it was possible to correlate many of the variables with service life. The results of this study have indicated that thermogravimetric analysis can play an important role in categorizing carbonate aggregates. In fact, with modern automated thermal analysis systems it should be possible to utilize such methods on a quality control basis. Strong correlations were found between several of the variables that were monitored in this study. In fact, several of the variables exhibited significant correlations to concrete service life. When the full data set was utilized (n = 18), the significant correlations to service life can be summarized as follows ( a = 5% level): • Correlation coefficient, r, = -0.73 for premature TG loss versus service life. • Correlation coefficient, r, = 0.74 for relative crystallite size versus service life. • Correlation coefficient, r, = 0.53 for ASTM C666 durability factor versus service life. • Correlation coefficient, r, = -0.52 for acid-insoluble residue versus service life. Separation of the carbonate aggregates into their mineralogical categories (i.e., calcites and dolomites) tended to increase the correlation coefficients for some specific variables (r sometimes approached 0.90); however, the reliability of such correlations was questionable because of the small number of samples that were present in this study.

Performance Analysis of Two Quantum Reaction Dynamics Codes: Time-Dependent and Time-Independent Strategies

The computer simulation of reaction dynamics has nowadays reached a remarkable degree of accuracy. Triatomic elementary reactions are rigorously studied with great detail on a straightforward basis using a considerable variety of Quantum Dynamics computational tools available to the scientific community. In our contribution we compare the performance of two quantum scattering codes in the computation of reaction cross sections of a triatomic benchmark reaction such as the gas phase reaction Ne + H2+ %12. NeH++ H. The computational codes are selected as representative of time-dependent (Real Wave Packet [ ]) and time-independent (ABC [ ]) methodologies. The main conclusion to be drawn from our study is that both strategies are, to a great extent, not competing but rather complementary. While time-dependent calculations advantages with respect to the energy range that can be covered in a single simulation, time-independent approaches offer much more detailed information from each single energy calculation. Further details such as the calculation of reactivity at very low collision energies or the computational effort related to account for the Coriolis couplings are analyzed in this paper.

Investigation of Rapid Thermal Analysis Procedures for Prediction of the Service Life of PCCP Carbonate Coarse Aggregate: Phase I, Progress Report, HR-337, 1992

The major objective of this research project is to utilize thermal analysis techniques in conjunction with x-ray analysis methods to identify and explain chemical reactions that promote aggregate related deterioration in Portland cement concrete. The first year of this project has been spent obtaining and analyzing limestone and dolomite samples that exhibit a wide range of field service performance. Most of the samples chosen for the study also had laboratory durability test information (ASTM C 666, method B) that was readily available. Preliminary test results indicate that a strong relationship exists between the average crystallite size of the limestone (calcite) specimens and their apparent decomposition temperatures as measured by thermogravimetric analysis. Also, premature weight loss in the thermogravimetric analysis tests appeared to be related to the apparent decomposition temperature of the various calcite test specimens.

Equilibrium of the simultaneous etherification ofisobutene and isoamylenes with ethanol inliquid-phase

The simultaneous etherification of isobutene and isoamylenes with ethanol has been studied using macroreticu-lar acid ion-exchange resins as catalyst. Most of the experiments were carried out over Amberlyst-35. In addition,Amberlyst-15 and Purolite CT-275 were also tested. Chemical equilibrium of four chemical reactions was studied:ethyl tert-butyl ether formation, tert-amyl ethyl ether formation from isoamylenes (2-methyl-1-butene and 2-methyl-2-butene) and isomerization reaction between both isoamylenes. Equilibrium data were obtained in a batchwisestirred tank reactor operated at 2.0 MPa and within the temperature range from 323 to 353 K. Experimental molarstandard enthalpy and entropy changes of reaction were determined for each reaction. From these data, the molarenthalpy change of formation of ethyl tert-butyl ether and tert-amyl ethyl ether were estimated. Besides, the chemical equilibrium between both diisobutene dimers, 2,4,4-trimethyl-1-pentene and 2,4,4-trimethyl-2-pentene, wasevaluated. A good agreement between thermodynamic results for the simultaneous etherification carried out in thiswork and those obtained for the isolated ethyl tert-butyl ether and tert-amyl ethyl ether systems was obtained.

Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

UniPathway: a resource for the exploration and annotation of metabolic pathways.

UniPathway (http://www.unipathway.org) is a fully manually curated resource for the representation and annotation of metabolic pathways. UniPathway provides explicit representations of enzyme-catalyzed and spontaneous chemical reactions, as well as a hierarchical representation of metabolic pathways. This hierarchy uses linear subpathways as the basic building block for the assembly of larger and more complex pathways, including species-specific pathway variants. All of the pathway data in UniPathway has been extensively cross-linked to existing pathway resources such as KEGG and MetaCyc, as well as sequence resources such as the UniProt KnowledgeBase (UniProtKB), for which UniPathway provides a controlled vocabulary for pathway annotation. We introduce here the basic concepts underlying the UniPathway resource, with the aim of allowing users to fully exploit the information provided by UniPathway.

Aspects on Process Engineering in the Finnish Pulp and Paper Industry : Thesis for the Degree of Licentiate in Science of Technology

Process development will be largely driven by the main equipment suppliers. The reason for this development is their ambition to supply complete plants or process systems instead of single pieces of equipment. The pulp and paper companies' interest lies in product development, as their main goal is to create winning brands and effective brand management. Design engineering companies will find their niche in detail engineering based on approved process solutions. Their development work will focus on increasing the efficiency of engineering work. Process design is a content-producing profession, which requires certain special characteristics: creativity, carefulness, the ability to work as a member of a design team according to time schedules and fluency in oral as well as written presentation. In the future, process engineers will increasingly need knowledge of chemistry as well as information and automation technology. Process engineering tools are developing rapidly. At the moment, these tools are good enough for static sizing and balancing, but dynamic simulation tools are not yet good enough for the complicated chemical reactions of pulp and paper chemistry. Dynamic simulation and virtual mill models are used as tools for training the operators. Computational fluid dynamics will certainlygain ground in process design.

Compact thermal storage: a-state-of-the-art review of experimental results

Thermal energy storage (TES) can increase the thermal energy effieresa, of a process by reusing the waste heat from industrial process, solar energy or other sources. There are different ways to store thermal energy: by sensible heat, by latest heat, by sorption process or by chemical reaction. This thesrs provides a-state-of-the-art review of the experimental performance of TES systems based on solid gas sorption process and chemical reactions. The importance of theses processes is that provides a heat loss free storage system with a high energy density.

Mitochondrial DNA damage and animal longevity: insights from comparative studies

Chemical reactions in living cells are under strict enzyme control and conform to a tightly regulated metabolic program. However, uncontrolled and potentially deleterious endogenous reactions occur, even under physiological conditions. Aging, in this chemical context, could be viewed as an entropic process, the result of chemical side reactions that chronically and cumulatively degrade the function of biological systems. Mitochondria are a main source of reactive oxygen species (ROS) and chemical sidereactions in healthy aerobic tissues and are the only known extranuclear cellular organelles in animal cells that contain their own DNA (mtDNA). ROS can modify mtDNA directly at the sugar-phosphate backbone or at the bases, producing many different oxidatively modified purines and pyrimidines, as well as single and double strand breaks and DNA mutations. In this scenario, natural selection tends to decrease the mitochondrial ROS generation, the oxidative damage to mtDNA, and the mitochondrial mutation rate in long-lived species, in agreement with the mitochondrial oxidative stress theory of aging.