A quantum chemical investigation of the metal centres in cytochrome c oxidase

Quantum effects are often of key importance for the function of biological systems at molecular level. Cellular respiration, where energy is extracted from the reduction of molecular oxygen to water, is no exception. In this work, the end station of the electron transport chain in mitochondria, cytochrome c oxidase, is investigated using quantum chemical methodology. Cytochrome c oxidase contains two haems, haem a and haem a3. Haem a3, with its copper companion, CuB, is involved in the final reduction of oxygen into water. This binuclear centre receives the necessary electrons from haem a. Haem a, in turn, receives its electrons from a copper ion pair in the vicinity, called CuA. Density functional theory (DFT) has been used to clarify the charge and spin distributions of haem a, as well as changes in these during redox activity. Upon reduction, the added electron is shown to be evenly distributed over the entire haem structure, important for the accommodation of the prosthetic group within the protein. At the same time, the spin distribution of the open-shell oxidised state is more localised to the central iron. The exact spin density distribution has been disputed in the literature, however, different experiments indicating different distributions of the unpaired electron. The apparent contradiction is shown to be due to the false assumption of a unit amount of unpaired electron density; in fact, the oxidised state has about 1.3 unpaired electrons. The validity of the DFT results have been corroborated by wave function based coupled cluster calculations. Point charges, for use in classical force field based simulations, have been parameterised for the four metal centres, using a newly developed methodology. In the procedure, the subsystem for which point charges are to be obtained, is surrounded by an outer region, with the purpose of stabilising the inner region, both electronically and structurally. Finally, the possibility of vibrational promotion of the electron transfer step between haem a and a3 has been investigated. Calculating the full vibrational spectra, at DFT level, of a combined model of the two haems, revealed several normal modes that do shift electron density between the haems. The magnitude of the shift was found to be moderate, at most. The proposed mechanism could have an assisting role in the electron transfer, which still seems to be dominated by electron tunnelling.

Computationally Efficient Methods for MDL-Optimal Density Estimation and Data Clustering

The Minimum Description Length (MDL) principle is a general, well-founded theoretical formalization of statistical modeling. The most important notion of MDL is the stochastic complexity, which can be interpreted as the shortest description length of a given sample of data relative to a model class. The exact definition of the stochastic complexity has gone through several evolutionary steps. The latest instantation is based on the so-called Normalized Maximum Likelihood (NML) distribution which has been shown to possess several important theoretical properties. However, the applications of this modern version of the MDL have been quite rare because of computational complexity problems, i.e., for discrete data, the definition of NML involves an exponential sum, and in the case of continuous data, a multi-dimensional integral usually infeasible to evaluate or even approximate accurately. In this doctoral dissertation, we present mathematical techniques for computing NML efficiently for some model families involving discrete data. We also show how these techniques can be used to apply MDL in two practical applications: histogram density estimation and clustering of multi-dimensional data.