Redox-linked proton transfer by cytochrome c oxidase

The respiratory chain is found in the inner mitochondrial membrane of higher organisms and in the plasma membrane of many bacteria. It consists of several membrane-spanning enzymes, which conserve the energy that is liberated from the degradation of food molecules as an electrochemical proton gradient across the membrane. The proton gradient can later be utilized by the cell for different energy requiring processes, e.g. ATP production, cellular motion or active transport of ions. The difference in proton concentration between the two sides of the membrane is a result of the translocation of protons by the enzymes of the respiratory chain, from the negatively charged (N-side) to the positively charged side (P-side) of the lipid bilayer, against the proton concentration gradient. The endergonic proton transfer is driven by the flow of electrons through the enzymes of the respiratory chain, from low redox-potential electron donors to acceptors of higher potential, and ultimately to oxygen. Cytochrome c oxidase is the last enzyme in the respiratory chain and catalyzes the reduction of dioxygen to water. The redox reaction is coupled to proton transport across the membrane by a yet unresolved mechanism. Cytochrome c oxidase has two proton-conducting pathways through which protons are taken up to the interior part of the enzyme from the N-side of the membrane. The K-pathway transfers merely substrate protons, which are consumed in the process of water formation at the catalytic site. The D-pathway transfers both substrate protons and protons that are pumped to the P-side of the membrane. This thesis focuses on the role of two conserved amino acids in proton translocation by cytochrome c oxidase, glutamate 278 and tryptophan 164. Glu278 is located at the end of the D-pathway and is thought to constitute the branching point for substrate and pumped protons. In this work, it was shown that although Glu278 has an important role in the proton transfer mechanism, its presence is not an obligatory requirement. Alternative structural solutions in the area around Glu278, much like the ones present in some distantly related heme-copper oxidases, could in the absence of Glu278 support the formation of a long hydrogen-bonded water chain through which proton transfer from the D-pathway to the catalytic site is possible. The other studied amino acid, Trp164, is hydrogen bonded to the ∆-propionate of heme a3 of the catalytic site. Mutation of this amino acid showed that it may be involved in regulation of proton access to a proton acceptor, a pump site, from which the proton later is expelled to the P-side of the membrane. The ion pair that is formed by the ∆-propionate of heme a3 and arginine 473 is likely to form a gate-like structure, which regulates proton mobility to the P-side of the membrane. The same gate may also be part of an exit path through which water molecules produced at the catalytically active site are removed towards the external side of the membrane. Time-resolved optical and electrometrical experiments with the Trp164 to phenylalanine mutant revealed a so far undetected step in the proton pumping mechanism. During the A to PR transition of the catalytic cycle, a proton is transferred from Glu278 to the pump site, located somewhere in the vicinity of the ∆-propionate of heme a3. A mechanism for proton pumping by cytochrome c oxidase is proposed on the basis of the presented results and the mechanism is discussed in relation to some relevant experimental data. A common proton pumping mechanism for all members of the heme-copper oxidase family is moreover considered.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Ecological impacts of Phlebiopsis gigantea biocontrol treatment against Heterobasidion spp. as revealed by fungal community profiling and population analyses

Wood decay fungi belonging to the species complex Heterobasidion annosum sensu lato are among the most common and economically important species causing root rot and stem decay in conifers of the northern temperate regions. New infections by these pathogens can be suppressed by tree stump treatments using chemical or biological control agents. In Finland, the corticiaceous fungus Phlebiopsis gigantea has been formulated into a commercial biocontrol agent called Rotstop (Verdera Ltd.). This thesis addresses the ecological impacts of Rotstop biocontrol treatment on the mycoflora of conifer stumps. Locally, fungal communities within Rotstop-treated and untreated stumps were analyzed using a novel method based on DGGE profiling of small subunit ribosomal DNA fragments amplified directly from wood samples. Population analyses for P. gigantea and H. annosum s.l. were conducted to evaluate possible risks associated with local and/or global distribution of the Rotstop strain. Based on molecular community profiling by DGGE, we detected a few individual wood-inhabiting fungal species (OTUs) that seemed to have suffered or benefited from the Rotstop biocontrol treatment. The DGGE analyses also revealed fungal diversity not retrieved by cultivation and some fungal sequence types untypical for decomposing conifer wood. However, statistical analysis of DGGE community profiles obtained from Rotstop-treated and untreated conifer stumps revealed that the Rotstop treatment had not caused a statistically significant reduction in the species diversity of wood-inhabiting fungi within our experimental forest plots. Locally, ISSR genotyping of cultured P. gigantea strains showed that the Rotstop biocontrol strain was capable of surviving up to six years within treated Norway spruce stumps, while in Scots pine stumps it was sooner replaced by successor fungal species. In addition, the spread of resident P. gigantea strains into Rotstop-treated forest stands seemed effective in preventing the formation of genetically monomorphic populations in the short run. On a global scale, we detected a considerable level of genetic differentiation between the interfertile European and North American populations of P. gigantea. These results strongly suggest that local biocontrol strains should be used in order to prevent global spread of P. gigantea and hybrid formation between geographically isolated populations. The population analysis for H. annosum s.l. revealed a collection of Chinese fungal strains that showed a high degree of laboratory fertility with three different allopatric H. annosum s.l. taxa. However, based on the molecular markers, the Chinese strains could be clearly affiliated with the H. parviporum taxonomical cluster, which thus appears to have a continuous distribution range from Europe through southern Siberia to northern China. Keywords: Rotstop, wood decay, DGGE, ISSR fingerprinting, ribosomal DNA

The role of productivity in the ecological and evolutionary dynamics of predator-prey interaction

Productivity is predicted to drive the ecological and evolutionary dynamics of predator-prey interaction through changes in resource allocation between different traits. However, resources are seldom constantly available and thus temporal variation in productivity could have considerable effect on the species' potential to evolve. To study this, three long-term microbial laboratory experiments were established where Serratia marcescens prey bacteria was exposed to predation of protist Tetrahymena thermophila in different prey resource environments. The consequences of prey resource availability for the ecological properties of the predator-prey system, such as trophic dynamics, stability, and virulence, were determined. The evolutionary changes in species traits and prey genetic diversity were measured. The prey defence evolved stronger in high productivity environment. Increased allocation to defence incurred cost in terms of reduced prey resource use ability, which probably constrained prey evolution by increasing the effect of resource competition. However, the magnitude of this trade-off diminished when measured in high resource concentrations. Predation selected for white, non-pigmented, highly defensive prey clones that produced predation resistant biofilm. The biofilm defence was also potentially accompanied with cytotoxicity for predators and could have been traded off with high motility. Evidence for the evolution of predators was also found in one experiment suggesting that co-evolutionary dynamics could affect the evolution and ecology of predator-prey interaction. Temporal variation in resource availability increased variation in predator densities leading to temporally fluctuating selection for prey defences and resource use ability. Temporal variation in resource availability was also able to constrain prey evolution when the allocation to defence incurred high cost. However, when the magnitude of prey trade-off was small and the resource turnover was periodically high, temporal variation facilitated the formation of predator resistant biofilm. The evolution of prey defence constrained the transfer of energy from basal to higher trophic levels, decreasing the strength of top-down regulation on prey community. Predation and temporal variation in productivity decreased the stability of populations and prey traits in general. However, predation-induced destabilization was less pronounced in the high productivity environment where the evolution of prey defence was stronger. In addition, evolution of prey defence weakened the environmental variation induced destabilization of predator population dynamics. Moreover, protozoan predation decreased the S. marcescens virulence in the insect host moth (Parasemia plantaginis) suggesting that species interactions outside the context of host-pathogen relationship could be important indirect drivers for the evolution of pathogenesis. This thesis demonstrates that rapid evolution can affect various ecological properties of predator-prey interaction. The effect of evolution on the ecological dynamics depended on the productivity of the environment, being most evident in the constant environments with high productivity.

Heat transfer to power-law fluids in thermal entrance region with viscous dissipation for constant heat flux conditions

An analytical solution of the heat transfer problem with viscous dissipation for non-Newtonian fluids with power-law model in the thermal entrance region of a circular pipe and two parallel plates under constant heat flux conditions is obtained using eigenvalue approach by suitably replacing one of the boundary conditions by total energy balance equation. Analytical expressions for the wall and the bulk temperatures and the local Nusselt number are presented. The results are in close agreement with those obtained by implicit finite-difference scheme. It is found that the role of viscous dissipation on heat transfer is completely different for heating and cooling conditions at the wall. The results for the case of cooling at the wall are of interest in the design of the oil pipe line.

Electron and proton transfer in NADH:ubiquinone oxidoreductase (Complex I) from Escherichia coli

Cells of every living organism on our planet − bacterium, plant or animal − are organized in such a way that despite differences in structure and function they utilize the same metabolic energy represented by electrochemical proton gradient across a membrane. This gradient of protons is generated by the series of membrane bound multisubunit proteins, Complex I, II, III and IV, organized in so-called respiratory or electron transport chain. In the eukaryotic cell it locates in the inner mitochondrial membrane while in the bacterial cell it locates in the cytoplasmic membrane. The function of the respiratory chain is to accept electrons from NADH and ubiquinol and transfer them to oxygen resulting in the formation of water. The free energy released upon these redox reactions is converted by respiratory enzymes into an electrochemical proton gradient, which is used for synthesis of ATP as well as for many other energy dependent processes. This thesis is focused on studies of the first member of the respiratory chain − NADH:ubiquinone oxidoreductase or Complex I. This enzyme has a boot-shape structure with hydrophilic and hydrophobic domains, the former of which has all redox groups of the protein, the flavin and eight to nine iron-sulfur clusters. Complex I serves as a proton pump coupling transfer of two electrons from NADH to ubiquinone to the translocation of four protons across the membrane. So far the mechanism of energy transduction by Complex I is unknown. In the present study we applied a set of different methods to study the electron and proton transfer reactions in Complex I from Escherichia coli. The main achievement was the experiment that showed that the electron transfer through the hydrophilic domain of Complex I is unlikely to be coupled to proton transfer directly or to conformational changes in the protein. In this work for the first time properties of all redox centers of Complex I were characterized in the intact purified bacterial enzyme. We also probed the role of several conserved amino acid residues in the electron transfer of Complex I. Finally, we found that highly conserved amino acid residues in several membrane subunits form a common pattern with a very prominent feature – the presence of a few lysines within the membrane. Based on the experimental data, we suggested a tentative principle which may govern the redox-coupled proton pumping in Complex I.

Population dynamics of blue mussels in a variable environment at the edge of their range

Spatial and temporal variation in the abundance of species can often be ascribed to spatial and temporal variation in the surrounding environment. Knowledge of how biotic and abiotic factors operate over different spatial and temporal scales in determining distribution, abundance, and structure of populations lies at the heart of ecology. The major part of the current ecological theory stems from studies carried out in central parts of the distributional range of species, whereas knowledge of how marginal populations function is inadequate. Understanding how marginal populations, living at the edge of their range, function is however in a key position to advance ecology and evolutionary biology as scientific disciplines. My thesis focuses on the factors affecting dynamics of marginal populations of blue mussels (Mytilus edulis) living close to their tolerance limits with regard to salinity. The thesis aims to highlight the dynamics at the edge of the range and contrast these with dynamics in more central parts of the range in order to understand the potential interplay between the central and the marginal part in the focal system. The objectives of the thesis are approached by studies on: (1) factors affecting regional patterns of the species, (2) long-term temporal dynamics of the focal species spaced along a regional salinity gradient, (3) selective predation by increasing populations of roach (Rutilus rutilus) when feeding on their main food item, the blue mussel, (4) the primary and secondary effects of local wave exposure gradients and (5) the role of small-scale habitat heterogeneity as determinants of large-scale pattern. The thesis shows that populations of blue mussels are largely determined by large scale changes in sea water salinity, affecting mainly recruitment success and longevity of local populations. In opposite to the traditional view, the thesis strongly indicate that vertebrate predators strongly affect abundance and size structure of blue mussel populations, and that the role of these predators increases towards the margin where populations are increasingly top-down controlled. The thesis also indicates that the positive role of biogenic habitat modifiers increases towards the marginal areas, where populations of blue mussels are largely recruitment limited. Finally, the thesis shows that local blue mussel populations are strongly dependent on high water turbulence, and therefore, dense populations are constrained to offshore habitats. Finally, the thesis suggests that ongoing sedimentation of rocky shores is detrimental for the species, affecting recruitment success and post-recruit survival, pushing stable mussel beds towards offshore areas. Ongoing large scale changes in the Baltic Sea, especially dilution processes with attendant effects, are predicted to substantially contract the distributional range of the mussel, but also affect more central populations. The thesis shows that in order to understand the functioning of marginal populations, research should (1) strive for multi-scale approaches in order to link ecosystem patterns with ecosystem processes, and (2) challenge the prevailing tenets that origin from research carried out in central areas that may not be valid at the edge.

Succinate and malate improve development of hamster eight-cell embryos in vitro: Confirmation of viability by embryo transfer

The in vitro development of hamster preimplantation embryos is supported by non-glucose energy substrates. To investigate the importance of embryonic metabolism, influence of succinate and malate on the development of hamster 8-cell embryos to blastocysts was examined using a chemically defined protein-free modified hamster embryo culture medium-2 (HECM-2m). There was a dose-dependent influence of succinate on blastocyst development; 0.5 mM succinate was optimal (85.1% ± 3.9 vs. 54.5% ± 3.5). In succinate-supplemented HECM-2m, blastocyst development was reduced by omission of lactate (68.5% ± 7.2), but not pyruvate (85.8% ± 6.2) or glutamine (84.1% ± 2.1). Succinate along with either glutamine or lactate or pyruvate poorly supported blastocyst development (28%-58%). Malate also stimulated blastocyst development; 0.01 mM malate was optimal (86.3% ± 2.8). Supplementation of both succinate and malate to HECM-2m supported maximal (100%) blastocyst development, which was inhibited 4-fold by the addition of glucose/phosphate. The mean cell numbers (MCN) of blastocysts cultured in succinate-supplemented HECM-2m was higher (28.3 ± 1.1) than it was for those cultured in the absence of glutamine or pyruvate (range 20-24). The MCN was the highest (33.4 ± 1.6) for blastocysts cultured in succinate-malate-supplemented HECM-2m followed by those in succinate (28.3 ± 1.1) or malate (24.7 ± 0.5) supplemented HECM-2m. Embryo transfer experiments showed that 29.8% (±4.5) of transferred blastocysts cultured in succinate-malate-supplemented HECM-2m produced live births, similar (P > 0.1) to the control transfers of freshly recovered 8-cells (33.5% ± 2.0) or blastocysts (28.9% ± 3.0). These data show that supplementation of succinate and malate to HECM-2m supports 100% development of hamster 8-cell embryos to high quality viable blastocysts and that non-glucose oxidizable energy substrates are the most preferred components in hamster embryo culture medium. Mol. Reprod. Dev. 47:440-447, 1997. © 1997 Wiley-Liss, Inc.

Conceptual and statistical modelling of environmental effects in population dynamics

Population dynamics are generally viewed as the result of intrinsic (purely density dependent) and extrinsic (environmental) processes. Both components, and potential interactions between those two, have to be modelled in order to understand and predict dynamics of natural populations; a topic that is of great importance in population management and conservation. This thesis focuses on modelling environmental effects in population dynamics and how effects of potentially relevant environmental variables can be statistically identified and quantified from time series data. Chapter I presents some useful models of multiplicative environmental effects for unstructured density dependent populations. The presented models can be written as standard multiple regression models that are easy to fit to data. Chapters II IV constitute empirical studies that statistically model environmental effects on population dynamics of several migratory bird species with different life history characteristics and migration strategies. In Chapter II, spruce cone crops are found to have a strong positive effect on the population growth of the great spotted woodpecker (Dendrocopos major), while cone crops of pine another important food resource for the species do not effectively explain population growth. The study compares rate- and ratio-dependent effects of cone availability, using state-space models that distinguish between process and observation error in the time series data. Chapter III shows how drought, in combination with settling behaviour during migration, produces asymmetric spatially synchronous patterns of population dynamics in North American ducks (genus Anas). Chapter IV investigates the dynamics of a Finnish population of skylark (Alauda arvensis), and point out effects of rainfall and habitat quality on population growth. Because the skylark time series and some of the environmental variables included show strong positive autocorrelation, the statistical significances are calculated using a Monte Carlo method, where random autocorrelated time series are generated. Chapter V is a simulation-based study, showing that ignoring observation error in analyses of population time series data can bias the estimated effects and measures of uncertainty, if the environmental variables are autocorrelated. It is concluded that the use of state-space models is an effective way to reach more accurate results. In summary, there are several biological assumptions and methodological issues that can affect the inferential outcome when estimating environmental effects from time series data, and that therefore need special attention. The functional form of the environmental effects and potential interactions between environment and population density are important to deal with. Other issues that should be considered are assumptions about density dependent regulation, modelling potential observation error, and when needed, accounting for spatial and/or temporal autocorrelation.

Human genetic variation in the Baltic Sea region: Features of population history and natural selection

In this thesis, the genetic variation of human populations from the Baltic Sea region was studied in order to elucidate population history as well as evolutionary adaptation in this region. The study provided novel understanding of how the complex population level processes of migration, genetic drift, and natural selection have shaped genetic variation in North European populations. Results from genome-wide, mitochondrial DNA and Y-chromosomal analyses suggested that the genetic background of the populations of the Baltic Sea region lies predominantly in Continental Europe, which is consistent with earlier studies and archaeological evidence. The late settlement of Fennoscandia after the Ice Age and the subsequent small population size have led to pronounced genetic drift, especially in Finland and Karelia but also in Sweden, evident especially in genome-wide and Y-chromosomal analyses. Consequently, these populations show striking genetic differentiation, as opposed to much more homogeneous pattern of variation in Central European populations. Additionally, the eastern side of the Baltic Sea was observed to have experienced eastern influence in the genome-wide data as well as in mitochondrial DNA and Y-chromosomal variation – consistent with linguistic connections. However, Slavic influence in the Baltic Sea populations appears minor on genetic level. While the genetic diversity of the Finnish population overall was low, genome-wide and Y-chromosomal results showed pronounced regional differences. The genetic distance between Western and Eastern Finland was larger than for many geographically distant population pairs, and provinces also showed genetic differences. This is probably mainly due to the late settlement of Eastern Finland and local isolation, although differences in ancestral migration waves may contribute to this, too. In contrast, mitochondrial DNA and Y-chromosomal analyses of the contemporary Swedish population revealed a much less pronounced population structure and a fusion of the traces of ancient admixture, genetic drift, and recent immigration. Genome-wide datasets also provide a resource for studying the adaptive evolution of human populations. This study revealed tens of loci with strong signs of recent positive selection in Northern Europe. These results provide interesting targets for future research on evolutionary adaptation, and may be important for understanding the background of disease-causing variants in human populations.

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.