Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Proton translocation coupled to electron transfer reactions in terminal oxidases

Terminal oxidases are the final proteins of the respiratory chain in eukaryotes and some bacteria. They catalyze most of the biological oxygen consumption on Earth done by aerobic organisms. During the catalytic reaction terminal oxidases reduce dioxygen to water and use the energy released in this process to maintain the electrochemical proton gradient by functioning as a redox-driven proton pump. This membrane gradient of protons is extremely important for cells as it is used for many cellular processes, such as transportation of substrates and ATP synthesis. Even though the structures of several terminal oxidases are known, they are not sufficient in themselves to explain the molecular mechanism of proton pumping. In this work we have applied a complex approach using a variety of different techniques to address the properties and the mechanism of proton translocation by the terminal oxidases. The combination of direct measurements of pH changes during catalytic turnover, time-resolved potentiometric electrometry and optical spectroscopy, made it possible to obtain valuable information about various aspects of oxidase functioning. We compared oxygen binding properties of terminal oxidases from the distinct heme-copper (CcO) and cytochrome bd families and found that cytochrome bd has a high affinity for oxygen, which is 3 orders of magnitude higher than that of CcO. Interestingly, the difference between CcO and cytochrome bd is not only in higher affinity of the latter to oxygen, but also in the way that each of these enzymes traps oxygen during catalysis. CcO traps oxygen kinetically - the molecule of bound dioxygen is rapidly reduced before it can dissociate. Alternatively, cytochrome bd employs an alternative mechanism of oxygen trapping - part of the redox energy is invested into tight oxygen binding, and the price paid for this is the lack of proton pumping. A single cycle of oxygen reduction to water is characterized by translocation of four protons across the membrane. Our results make it possible to assign the pumping steps to discrete transitions of the catalytic cycle and indicate that during in vivo turnover of the oxidase these four protons are transferred, one at a time, during the P→F, F→OH, Oh→Eh, and Eh→R transitions. At the same time, each individual proton translocation step in the catalytic cycle is not just a single reaction catalyzed by CcO, but rather a complicated sequence of interdependent electron and proton transfers. We assume that each single proton translocation cycle of CcO is assured by internal proton transfer from the conserved Glu-278 to an as yet unidentified pump site above the hemes. Delivery of a proton to the pump site serves as a driving reaction that forces the proton translocation cycle to continue.

Dynamics of free-running, pump-modulated and coupled semiconductor and solid-state lasers

The output of a laser is a high frequency propagating electromagnetic field with superior coherence and brightness compared to that emitted by thermal sources. A multitude of different types of lasers exist, which also translates into large differences in the properties of their output. Moreover, the characteristics of the electromagnetic field emitted by a laser can be influenced from the outside, e.g., by injecting an external optical field or by optical feedback. In the case of free-running solitary class-B lasers, such as semiconductor and Nd:YVO4 solid-state lasers, the phase space is two-dimensional, the dynamical variables being the population inversion and the amplitude of the electromagnetic field. The two-dimensional structure of the phase space means that no complex dynamics can be found. If a class-B laser is perturbed from its steady state, then the steady state is restored after a short transient. However, as discussed in part (i) of this Thesis, the static properties of class-B lasers, as well as their artificially or noise induced dynamics around the steady state, can be experimentally studied in order to gain insight on laser behaviour, and to determine model parameters that are not known ab initio. In this Thesis particular attention is given to the linewidth enhancement factor, which describes the coupling between the gain and the refractive index in the active material. A highly desirable attribute of an oscillator is stability, both in frequency and amplitude. Nowadays, however, instabilities in coupled lasers have become an active area of research motivated not only by the interesting complex nonlinear dynamics but also by potential applications. In part (ii) of this Thesis the complex dynamics of unidirectionally coupled, i.e., optically injected, class-B lasers is investigated. An injected optical field increases the dimensionality of the phase space to three by turning the phase of the electromagnetic field into an important variable. This has a radical effect on laser behaviour, since very complex dynamics, including chaos, can be found in a nonlinear system with three degrees of freedom. The output of the injected laser can be controlled in experiments by varying the injection rate and the frequency of the injected light. In this Thesis the dynamics of unidirectionally coupled semiconductor and Nd:YVO4 solid-state lasers is studied numerically and experimentally.

Characterizing Physical Properties of Wetting Paper by Air Coupled Ultrasound and Light

This thesis reports investigations into the paper wetting process and its effects on the surface roughness and the out-of-plane (ZD) stiffness of machine-made paper. The aim of this work was to test the feasibility of employing air-borne ultrasound methods to determine surface roughness (by reflection) and ZD stiffness (by through transmission) of paper during penetration of distilled water, isopropanol and their mixtures. Air-borne ultrasound provides a non-contacting way to evaluate sample structure and mechanics during the liquid penetration event. Contrary to liquid immersion techniques, an air-borne measurement allows studying partial wetting of paper. In addition, two optical methods were developed to reveal the liquid location in paper during wetting. The laser light through transmission method was developed to monitor the liquid location in partially wetted paper. The white light reflection method was primarily used to monitor the penetration of the liquid front in the thickness direction. In the latter experiment the paper was fully wetted. The main results of the thesis were: 1) Liquid penetration induced surface roughening was quantified by monitoring the ultrasound reflection from the paper surface. 2) Liquid penetration induced stiffness alteration in the ZD of paper could be followed by measuring the change in the ultrasound ZD resonance in paper. 3) Through transmitted light revealed the liquid location in the partially wetted paper. 4) Liquid movement in the ZD of the paper could be observed by light reflection. The results imply that the presented ultrasonic means can without contact measure the alteration of paper roughness and stiffness during liquid transport. These methods can help avoiding over engineering the paper which reduces raw material and energy consumption in paper manufacturing. The presented optical means can estimate paper specific wetting properties, such as liquid penetration speed, transport mechanisms and liquid location within the paper structure. In process monitoring, these methods allow process tuning and manufacturing of paper with engineered liquid transport characteristics. With such knowledge the paper behaviour during printing can be predicted. These findings provide new methods for paper printing, surface sizing, and paper coating research.

Regeneration ecology and woody species diversity in the remnant natural forest at Wondo Genet, Ethiopia – Potential for landscape restoration

Regeneration ecology, diversity of native woody species and its potential for landscape restoration was studied in the remnant natural forest at the College of Forestry and Natural Resources at Wondo Genet, Ethiopia. The type of forest is Afromontane rainforest , with many valuable tree species like Aningeria adolfi-friederici, and it is an important provider of ecological, social and economical services for the population that lives in this area. The study contains two parts, natural regeneration studies (at the natural forest) and interviews with farmers in the nearby village of the remnant patch. The objective of the first part was to investigate the floristic composition, densitiy and regeneration profiles of native woody species in the forest, paying special attention to woody species that are considered the most relevant (socio-economic). The second part provided information on woody species preferred by the farmers and on multiple uses of the adjacent natural forest, it also provided information and analysed perceptions on forest degradation. Systematic plot sampling was used in the forest inventory. Twenty square plots of 20 x 20 m were assessed, with 38 identified woody species (the total number of species was 45), representing 26 families. Of these species 61% were trees, 13% shrubs, 11% lianas and 16% species that could have both life forms. An analysis of natural regeneration of five important tree species in the natural forest showed that Aningeria adolfi-friederici had the best regeneration results. An analysis of population structure (as determined by height classes) of two commercially important woody species in the forest, Aningeria adolfi-friederici and Podocarpus falcatus, showed a marked difference: Aningeria had a typical “reversed J” frequency distribution, while Podocarpus showed very low values in all height classes. Multi dimensional scaling (MDS) was used to map the sample plots according to their similarity in species composition, using the Sørensen quantitative index, coupled with indicator species analysis .Three groups were identified with respective indicator species: Group 1 – Adhatoda schimperiana, Group 2 – Olea hochstetteri , Group 3 – Acacia senegal and Aningeria adolfi-friederici. Thirty questionnaire interviews were conducted with farmers in the village of Gotu Onoma that use the nearby remant forest patch. Their tree preferences were exotic species such as Eucalyptus globulus for construction and fuelwood and Grevillea robusta for shade and fertility. Considering forest land degradation farmers were aware of the problem and suggested that the governmental institutions address the problem by planting more Eucalyptus globulus. The natural forest seemed to have moderate levels of disturbance and it was still floristically diverse. However, the low rate of natural regeneration of Podocarpus falcatus suggested that this species is threatened and must be a priority in conservation actions. Plantations and agroforestry seem to be possible solutions for rehabilitation of the surrounding degraded lands, thereby decreasing the existent pressure in the remnant natural forest.

Models of Dispersal and Diversification

Ecology and evolutionary biology is the study of life on this planet. One of the many methods applied to answering the great diversity of questions regarding the lives and characteristics of individual organisms, is the utilization of mathematical models. Such models are used in a wide variety of ways. Some help us to reason, functioning as aids to, or substitutes for, our own fallible logic, thus making argumentation and thinking clearer. Models which help our reasoning can lead to conceptual clarification; by expressing ideas in algebraic terms, the relationship between different concepts become clearer. Other mathematical models are used to better understand yet more complicated models, or to develop mathematical tools for their analysis. Though helping us to reason and being used as tools in the craftmanship of science, many models do not tell us much about the real biological phenomena we are, at least initially, interested in. The main reason for this is that any mathematical model is a simplification of the real world, reducing the complexity and variety of interactions and idiosynchracies of individual organisms. What such models can tell us, however, both is and has been very valuable throughout the history of ecology and evolution. Minimally, a model simplifying the complex world can tell us that in principle, the patterns produced in a model could also be produced in the real world. We can never know how different a simplified mathematical representation is from the real world, but the similarity models do strive for, gives us confidence that their results could apply. This thesis deals with a variety of different models, used for different purposes. One model deals with how one can measure and analyse invasions; the expanding phase of invasive species. Earlier analyses claims to have shown that such invasions can be a regulated phenomena, that higher invasion speeds at a given point in time will lead to a reduction in speed. Two simple mathematical models show that analysis on this particular measure of invasion speed need not be evidence of regulation. In the context of dispersal evolution, two models acting as proof-of-principle are presented. Parent-offspring conflict emerges when there are different evolutionary optima for adaptive behavior for parents and offspring. We show that the evolution of dispersal distances can entail such a conflict, and that under parental control of dispersal (as, for example, in higher plants) wider dispersal kernels are optimal. We also show that dispersal homeostasis can be optimal; in a setting where dispersal decisions (to leave or stay in a natal patch) are made, strategies that divide their seeds or eggs into fractions that disperse or not, as opposed to randomized for each seed, can prevail. We also present a model of the evolution of bet-hedging strategies; evolutionary adaptations that occur despite their fitness, on average, being lower than a competing strategy. Such strategies can win in the long run because they have a reduced variance in fitness coupled with a reduction in mean fitness, and fitness is of a multiplicative nature across generations, and therefore sensitive to variability. This model is used for conceptual clarification; by developing a population genetical model with uncertain fitness and expressing genotypic variance in fitness as a product between individual level variance and correlations between individuals of a genotype. We arrive at expressions that intuitively reflect two of the main categorizations of bet-hedging strategies; conservative vs diversifying and within- vs between-generation bet hedging. In addition, this model shows that these divisions in fact are false dichotomies.