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.

Higher-resolution structure of the human insulin receptor ectodomain: Multi-modal inclusion of the insert domain

Insulin receptor (IR) signaling is critical to controlling nutrient uptake and metabolism. However, only a low-resolution (3.8 Å) structure currently exists for the IR ectodomain, with some segments ill-defined or unmodeled due to disorder. Here, we revise this structure using new diffraction data to 3.3 Å resolution that allow improved modeling of the N-linked glycans, the first and third fibronectin type III domains, and the insert domain. A novel haptic interactive molecular dynamics strategy was used to aid fitting to low-resolution electron density maps. The resulting model provides a foundation for investigation of structural transitions in IR upon ligand binding.

Missing-In-Metastasis (MIM) regulates cell morphology by promoting plasma membrane and actin cytoskeleton dynamics

The cells of multicellular organisms have differentiated to carry out specific functions that are often accompanied by distinct cell morphology. The actin cytoskeleton is one of the key regulators of cell shape subsequently controlling multiple cellular events including cell migration, cell division, endo- and exocytosis. A large set of actin regulating proteins has evolved to achieve and tightly coordinate this wide range of functions. Some actin regulator proteins have so-called house keeping roles and are essential for all eukaryotic cells, but some have evolved to meet the requirements of more specialized cell-types found in higher organisms enabling complex functions of differentiated organs, such as liver, kidney and brain. Often processes mediated by the actin cytoskeleton, like formation of cellular protrusions during cell migration, are intimately linked to plasma membrane remodeling. Thus, a close cooperation between these two cellular compartments is necessary, yet not much is known about the underlying molecular mechanisms. This study focused on a vertebrate-specific protein called missing-in-metastasis (MIM), which was originally characterized as a metastasis suppressor of bladder cancer. We demonstrated that MIM regulates the dynamics of actin cytoskeleton via its WH2 domain, and is expressed in a cell-type specific manner. Interestingly, further examination showed that the IM-domain of MIM displays a novel membrane tubulation activity, which induces formation of filopodia in cells. Following studies demonstrated that this membrane deformation activity is crucial for cell protrusions driven by MIM. In mammals, there are five members of IM-domain protein family. Functions and expression patterns of these family members have remained poorly characterized. To understand the physiological functions of MIM, we generated MIM knockout mice. MIM-deficient mice display no apparent developmental defects, but instead suffer from progressive renal disease and increased susceptibility to tumors. This indicates that MIM plays a role in the maintenance of specific physiological functions associated with distinct cell morphologies. Taken together, these studies implicate MIM both in the regulation of the actin cytoskeleton and the plasma membrane. Our results thus suggest that members of MIM/IRSp53 protein family coordinate the actin cytoskeleton:plasma membrane interface to control cell and tissue morphogenesis in multicellular organisms.

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.

Dissipative dynamics of a harmonic oscillator: A nonperturbative approach

Starting from a microscopic theory, we derive a master equation for a harmonic oscillator coupled to a bath of noninteracting oscillators. We follow a nonperturbative approach, proposed earlier by us for the free Brownian particle. The diffusion constants are calculated analytically and the positivity of the master equation is shown to hold above a critical temperature. We compare the long time behavior of the average kinetic and potential energies with known thermodynamic results. In the limit of vanishing oscillator frequency of the system, we recover the results of the free Brownian particle.

Entanglement production due to quench dynamics of an anisotropic XY chain in a transverse field

We compute concurrence and negativity as measures of two-spin entanglement generated by a power-law quench (characterized by a rate tau(-1) and an exponent alpha) which takes an anisotropic XY chain in a transverse field through a quantum critical point (QCP). We show that only spins separated by an even number of lattice spacings get entangled in such a process. Moreover, there is a critical rate of quench, tau(-1)(c), above which no two-spin entanglement is generated; the entire entanglement is multipartite. The ratio of the entanglements between consecutive even neighbors can be tuned by changing the quench rate. We also show that for large tau, the concurrence (negativity) scales as root alpha/tau(alpha/tau), and we relate this scaling behavior to defect production by the quench through a QCP.

Coarse-Grained Molecular Dynamics Simulation of the Aggregation Properties of Multiheaded Cationic Surfactants in Water

The aggregation property of multiheaded surfactants has been investigated by constant pressure molecular dynamics (MD) simulation in aqueous medium. The model multiheaded surfactants contain more than one headgroup (x = 2, 3, and 4) for a single tail group. This increases the hydrophilic charge progressively over the hydrophobic tail which has dramatic consequences in the aggregation behavior. In particular, we have looked at the change in the aggregation property such as critical micellar concentration (cmc), aggregation number, and size of the micelles for the multiheaded surfactants in water. We find with increasing number of headgroups of the Multiheaded surfactants that the cmc values increase and the aggregation numbers as well as the size of the micelles decrease. These trends are in agreement with the experimental findings as reported earlier with x = 1, 2, and 3. We also predict the aggregation properties of multiheaded surfactant With four headgroups (x = 4) for which no experimental studies exist yet.

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.

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.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.