Tailoring heme -thiolate proteins into efficient biocatalysts with high specificity and selectivity

Cytochrome P450 monooxygenases, one of the most important classes of heme-thiolate proteins, have attracted considerable interest in the biochemical community because of its catalytic versatility, substrate diversity and great number in the superfamily. Although P450s are capable of catalyzing numerous difficult oxidation reactions, the relatively low stability, low turnover rates and the need of electron-donating cofactors have limited their practical biotechnological and pharmaceutical applications as isolated enzymes. The goal of this study is to tailor such heme-thiolate proteins into efficient biocatalysts with high specificity and selectivity by protein engineering and to better understand the structure-function relationship in cytochromes P450. In the effort to engineer P450cam, the prototype member of the P450 superfamily, into an efficient peroxygenase that utilizes hydrogen peroxide via the “peroxide-shunt” pathway, site-directed mutagenesis has been used to elucidate the critical roles of hydrophobic residues in the active site. Various biophysical, biochemical and spectroscopic techniques have been utilized to investigate the wild-type and mutant proteins. Three important P450cam variants were obtained showing distinct structural and functional features. In P450camV247H mutant, which exhibited almost identical spectral properties with the wild-type, it is demonstrated that a single amino acid switch turned the monooxygenase into an efficient preoxidase by increasing the peroxidase activity nearly one thousand folds. In order to tune the distal pocket of P450cam with polar residues, Leu 246 was replaced with a basic residue, lysine, resulting in a mutant with spectral features identical to P420, the inactive species of P450. But this inactive-species-like mutant showed catalytic activities without the facilitation of any cofactors. By substituting Gly 248 with a histidine, a novel Cys-Fe-His ligation set was obtained in P450cam which represented the very rare case of His ligation in heme-thiolate proteins. In addition to serving as a convenient model for hemoprotein structural studies, the G248H mutant also provided evidence about the nature of the axial ligand in cytochrome P420 and other engineered hemoproteins with thiolate ligations. Furthermore, attempts have been made to replace the proximal ligand in sperm whale myoglobin to construct a heme-thiolate protein model by mimicking the protein environment of cytochrome P450cam and chloroperoxidase.

Tailoring Heme-Thiolate Proteins into Efficient Biocatalysts with High Specificity and Selectivity

Cytochrome P450 monooxygenases, one of the most important classes of heme-thiolate proteins, have attracted considerable interest in the biochemical community because of its catalytic versatility, substrate diversity and great number in the superfamily. Although P450s are capable of catalyzing numerous difficult oxidation reactions, the relatively low stability, low turnover rates and the need of electron-donating cofactors have limited their practical biotechnological and pharmaceutical applications as isolated enzymes. The goal of this study is to tailor such heme-thiolate proteins into efficient biocatalysts with high specificity and selectivity by protein engineering and to better understand the structure-function relationship in cytochromes P450. In the effort to engineer P450cam, the prototype member of the P450 superfamily, into an efficient peroxygenase that utilizes hydrogen peroxide via the “peroxide-shunt” pathway, site-directed mutagenesis has been used to elucidate the critical roles of hydrophobic residues in the active site. Various biophysical, biochemical and spectroscopic techniques have been utilized to investigate the wild-type and mutant proteins. Three important P450cam variants were obtained showing distinct structural and functional features. In P450camV247H mutant, which exhibited almost identical spectral properties with the wild-type, it is demonstrated that a single amino acid switch turned the monooxygenase into an efficient preoxidase by increasing the peroxidase activity nearly one thousand folds. In order to tune the distal pocket of P450cam with polar residues, Leu 246 was replaced with a basic residue, lysine, resulting in a mutant with spectral features identical to P420, the inactive species of P450. But this inactive-species-like mutant showed catalytic activities without the facilitation of any cofactors. By substituting Gly 248 with a histidine, a novel Cys-Fe-His ligation set was obtained in P450cam which represented the very rare case of His ligation in heme-thiolate proteins. In addition to serving as a convenient model for hemoprotein structural studies, the G248H mutant also provided evidence about the nature of the axial ligand in cytochrome P420 and other engineered hemoproteins with thiolate ligations. Furthermore, attempts have been made to replace the proximal ligand in sperm whale myoglobin to construct a heme-thiolate protein model by mimicking the protein environment of cytochrome P450cam and chloroperoxidase.

Mastigatory efficiency in denture wears with bilateral balanced occlusion and canine guidance

The maintenance of masticatory function is especially important for patients wearing complete dentures due to their limitations. Thus, the bilateral balanced occlusal concept is used to achieve greater masticatory efficiency. However, a critical review of the literature reveals that there is not sufficient scientific evidence to support bilateral balanced occlusion as the most appropriate occlusal concept in complete dentures. Therefore, the aim of this study was to evaluate the masticatory efficiency in complete dentures wearers with bilateral balanced occlusion and canine guidance. A double-blinded controlled crossover clinical trial was conducted. The sample was composed by 24 edentulous patients who wore sets of complete dentures with both occlusal concepts during equal periods of 3 months. Objective data were collected through the masticatory efficiency test performed by the colorimetric method with the beads, in which capsules of a synthetic material enclosing fuchsine-containing granules were used. Subjective data were recorded by patient's ratings of their chewing function. No significant statistical difference was found for masticatory efficiency (p=0.095) between the two occlusal concepts studied. The results suggest that bilateral balanced occlusion does not improve the masticatory efficiency in complete denture wearers.

Mastigatory efficiency in denture wears with bilateral balanced occlusion and canine guidance

The maintenance of masticatory function is especially important for patients wearing complete dentures due to their limitations. Thus, the bilateral balanced occlusal concept is used to achieve greater masticatory efficiency. However, a critical review of the literature reveals that there is not sufficient scientific evidence to support bilateral balanced occlusion as the most appropriate occlusal concept in complete dentures. Therefore, the aim of this study was to evaluate the masticatory efficiency in complete dentures wearers with bilateral balanced occlusion and canine guidance. A double-blinded controlled crossover clinical trial was conducted. The sample was composed by 24 edentulous patients who wore sets of complete dentures with both occlusal concepts during equal periods of 3 months. Objective data were collected through the masticatory efficiency test performed by the colorimetric method with the beads, in which capsules of a synthetic material enclosing fuchsine-containing granules were used. Subjective data were recorded by patient's ratings of their chewing function. No significant statistical difference was found for masticatory efficiency (p=0.095) between the two occlusal concepts studied. The results suggest that bilateral balanced occlusion does not improve the masticatory efficiency in complete denture wearers.

Regional flood frequency analysis for Newfoundland and Labrador using the L-Moments index-flood method

The L-moments based index-flood procedure had been successfully applied for Regional Flood Frequency Analysis (RFFA) for the Island of Newfoundland in 2002 using data up to 1998. This thesis, however, considered both Labrador and the Island of Newfoundland using the L-Moments index-flood method with flood data up to 2013. For Labrador, the homogeneity test showed that Labrador can be treated as a single homogeneous region and the generalized extreme value (GEV) was found to be more robust than any other frequency distributions. The drainage area (DA) is the only significant variable for estimating the index-flood at ungauged sites in Labrador. In previous studies, the Island of Newfoundland has been considered as four homogeneous regions (A,B,C and D) as well as two Water Survey of Canada's Y and Z sub-regions. Homogeneous regions based on Y and Z was found to provide more accurate quantile estimates than those based on four homogeneous regions. Goodness-of-fit test results showed that the generalized extreme value (GEV) distribution is most suitable for the sub-regions; however, the three-parameter lognormal (LN3) gave a better performance in terms of robustness. The best fitting regional frequency distribution from 2002 has now been updated with the latest flood data, but quantile estimates with the new data were not very different from the previous study. Overall, in terms of quantile estimation, in both Labrador and the Island of Newfoundland, the index-flood procedure based on L-moments is highly recommended as it provided consistent and more accurate result than other techniques such as the regression on quantile technique that is currently used by the government.

Hyperspectral Image Classification with Nonlinear Methods

This thesis introduces two related lines of study on classification of hyperspectral images with nonlinear methods. First, it describes a quantitative and systematic evaluation, by the author, of each major component in a pipeline for classifying hyperspectral images (HSI) developed earlier in a joint collaboration [23]. The pipeline, with novel use of nonlinear classification methods, has reached beyond the state of the art in classification accuracy on commonly used benchmarking HSI data [6], [13]. More importantly, it provides a clutter map, with respect to a predetermined set of classes, toward the real application situations where the image pixels not necessarily fall into a predetermined set of classes to be identified, detected or classified with.

The particular components evaluated are a) band selection with band-wise entropy spread, b) feature transformation with spatial filters and spectral expansion with derivatives c) graph spectral transformation via locally linear embedding for dimension reduction, and d) statistical ensemble for clutter detection. The quantitative evaluation of the pipeline verifies that these components are indispensable to high-accuracy classification.

Secondly, the work extends the HSI classification pipeline with a single HSI data cube to multiple HSI data cubes. Each cube, with feature variation, is to be classified of multiple classes. The main challenge is deriving the cube-wise classification from pixel-wise classification. The thesis presents the initial attempt to circumvent it, and discuss the potential for further improvement.

Experimental Study of Storage Ring Free-Electron Laser with Novel Capabilities

The Duke Free-electron laser (FEL) system, driven by the Duke electron storage ring, has been at the forefront of developing new light source capabilities over the past two decades. In 1999, the Duke FEL demonstrated the first lasing of a storage ring FEL in the vacuum ultraviolet (VUV) region at $194$ nm using two planar OK-4 undulators. With two helical undulators added to the outboard sides of the planar undulators, in 2005 the highest FEL gain ($47.8\%$) of a storage ring FEL was achieved using the Duke FEL system with a four-undulator configuration. In addition, the Duke FEL has been used as the photon source to drive the High Intensity $\gamma$-ray Source (HIGS) via Compton scattering of the FEL beam and electron beam inside the FEL cavity. Taking advantage of FEL's wavelength tunability as well as the adjustability of the energy of the electron beam in the storage ring, the nearly monochromatic $\gamma$-ray beam has been produced in a wide energy range from $1$ to $100$ MeV at the HIGS. To further push the FEL short wavelength limit and enhance the FEL gain in the VUV regime for high energy $\gamma$-ray production, two additional helical undulators were installed in 2012 using an undulator switchyard system to allow switching between the two planar and two helical undulators in the middle section of the FEL system. Using different undulator configurations made possible by the switchyard, a number of novel capabilities of the storage ring FEL have been developed and exploited for a wide FEL wavelength range from infrared (IR) to VUV. These new capabilities will eventually be made available to the $\gamma$-ray operation, which will greatly enhance the $\gamma$-ray user research program, creating new opportunities for certain types of nuclear physics research.

With the wide wavelength tuning range, the FEL is an intrinsically well-suited device to produce lasing with multiple colors. Taking advantage of the availability of an undulator system with multiple undulators, we have demonstrated the first two-color lasing of a storage ring FEL. Using either a three- or four-undulator configuration with a pair of dual-band high reflectivity mirrors, we have achieved simultaneous lasing in the IR and UV spectral regions. With the low-gain feature of the storage ring FEL, the power generated at the two wavelengths can be equally built up and precisely balanced to reach FEL saturation. A systematic experimental program to characterize this two-color FEL has been carried out, including precise power control, a study of the power stability of two-color lasing, wavelength tuning, and the impact of the FEL mirror degradation. Using this two-color laser, we have started to develop a new two-color $\gamma$-ray beam for scientific research at the HIGS.

Using the undulator switchyard, four helical undulators installed in the beamline can be configured to not only enhance the FEL gain in the VUV regime, but also allow for the full polarization control of the FEL beams. For the accelerator operation, the use of helical undulators is essential to extend the FEL mirror lifetime by reducing radiation damage from harmonic undulator radiation. Using a pair of helical undulators with opposite helicities, we have realized (1) fast helicity switching between left- and right-circular polarizations, and (2) the generation of fully controllable linear polarization. In order to extend these new capabilities of polarization control to the $\gamma$-ray operation in a wide energy range at the HIGS, a set of FEL polarization diagnostic systems need to be developed to cover the entire FEL wavelength range. The preliminary development of the polarization diagnostics for the wavelength range from IR to UV has been carried out.

Fast 2D/3D object representation with growing neural gas

This work presents the design of a real-time system to model visual objects with the use of self-organising networks. The architecture of the system addresses multiple computer vision tasks such as image segmentation, optimal parameter estimation and object representation. We first develop a framework for building non-rigid shapes using the growth mechanism of the self-organising maps, and then we define an optimal number of nodes without overfitting or underfitting the network based on the knowledge obtained from information-theoretic considerations. We present experimental results for hands and faces, and we quantitatively evaluate the matching capabilities of the proposed method with the topographic product. The proposed method is easily extensible to 3D objects, as it offers similar features for efficient mesh reconstruction.

Computational solid mechanics using a vertex-based finite volume method

A number of research groups are now developing and using finite volume (FV) methods for computational solid mechanics (CSM). These methods are proving to be equivalent and in some cases superior to their finite element (FE) counterparts. In this paper we will describe a vertex-based FV method with arbitrarily structured meshes, for modelling the elasto-plastic deformation of solid materials undergoing small strains in complex geometries. Comparisons with rational FE methods will be given.

Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.

Issues and strategies in the parallelisation of unstructured multiphysics code

Unstructured mesh based codes for the modelling of continuum physics phenomena have evolved to provide the facility to model complex interacting systems. Such codes have the potential to provide a high performance on parallel platforms for a small investment in programming. The critical parameters for success are to minimise changes to the code to allow for maintenance while providing high parallel efficiency, scalability to large numbers of processors and portability to a wide range of platforms. The paradigm of domain decomposition with message passing has for some time been demonstrated to provide a high level of efficiency, scalability and portability across shared and distributed memory systems without the need to re-author the code into a new language. This paper addresses these issues in the parallelisation of a complex three dimensional unstructured mesh Finite Volume multiphysics code and discusses the implications of automating the parallelisation process.

Modelling continuum mechanics phenomena using three dimensional unstructured meshes on massively parallel processors

The difficulties encountered in implementing large scale CM codes on multiprocessor systems are now fairly well understood. Despite the claims of shared memory architecture manufacturers to provide effective parallelizing compilers, these have not proved to be adequate for large or complex programs. Significant programmer effort is usually required to achieve reasonable parallel efficiencies on significant numbers of processors. The paradigm of Single Program Multi Data (SPMD) domain decomposition with message passing, where each processor runs the same code on a subdomain of the problem, communicating through exchange of messages, has for some time been demonstrated to provide the required level of efficiency, scalability, and portability across both shared and distributed memory systems, without the need to re-author the code into a new language or even to support differing message passing implementations. Extension of the methods into three dimensions has been enabled through the engineering of PHYSICA, a framework for supporting 3D, unstructured mesh and continuum mechanics modeling. In PHYSICA, six inspectors are used. Part of the challenge for automation of parallelization is being able to prove the equivalence of inspectors so that they can be merged into as few as possible.

Partition alignment in three dimensional unstructured mesh multi-physics modelling

Unstructured mesh codes for modelling continuum physics phenomena have evolved to provide the facility to model complex interacting systems. Parallelisation of such codes using single Program Multi Data (SPMD) domain decomposition techniques implemented with message passing has been demonstrated to provide high parallel efficiency, scalability to large numbers of processors P and portability across a wide range of parallel platforms. High efficiency, especially for large P requires that load balance is achieved in each parallel loop. For a code in which loops span a variety of mesh entity types, for example, elements, faces and vertices, some compromise is required between load balance for each entity type and the quantity of inter-processor communication required to satisfy data dependence between processors.

Symmetric Interior Penalty DG Methods for the Compressible Navier-Stokes Equations I: Method Formulation

In this article we consider the development of discontinuous Galerkin finite element methods for the numerical approximation of the compressible Navier-Stokes equations. For the discretization of the leading order terms, we propose employing the generalization of the symmetric version of the interior penalty method, originally developed for the numerical approximation of linear self-adjoint second-order elliptic partial differential equations. In order to solve the resulting system of nonlinear equations, we exploit a (damped) Newton-GMRES algorithm. Numerical experiments demonstrating the practical performance of the proposed discontinuous Galerkin method with higher-order polynomials are presented.

Near-infraread observation of supernovae with the Nordic Optical Telescope (NOT)

Supernova (SN) is an explosion of a star at the end of its lifetime. SNe are classified to two types, namely type I and II through the optical spectra. They have been categorised based on their explosion mechanism, to core collapse supernovae (CCSNe) and thermonuclear supernovae. The CCSNe group which includes types IIP, IIn, IIL, IIb, Ib, and Ic are produced when a massive star with initial mass more than 8 M⊙ explodes due to a collapse of its iron core. On the other hand, thermonuclear SNe originate from white dwarfs (WDs) made of carbon and oxygen, in a binary system. Infrared astronomy covers observations of astronomical objects in infrared radiation. The infrared sky is not completely dark and it is variable. Observations of SNe in the infrared give different information than optical observations. Data reduction is required to correct raw data from for example unusable pixels and sky background. In this project, the NOTCam package in the IRAF was used for the data reduction. For measuring magnitudes of SNe, the aperture photometry method with the Gaia program was used. In this Master’s thesis, near-infrared (NIR) observations of three supernovae of type IIn (namely LSQ13zm, SN 2009ip and SN2011jb), one type IIb (SN2012ey), in addition to one type Ic (SN2012ej) and type IIP (SN 2013gd) are studied with emphasis on luminosity and colour evolution. All observations were done with the Nordic Optical Telescope (NOT). Here, we used the classification by Mattila & Meikle (2001) [76], where the SNe are differentiated by the infrared light curves into two groups, namely ’ordinary’ and ’slowly declining’. The light curves and colour evolution of these supernovae were obtained in J, H and Ks bands. In this study, our data, combined with other observations, provide evidence to categorize LSQ13zm, SN 2012ej and SN 2012ey as being part of the ordinary type. We found interesting NIR behaviour of SN 2011jb, which lead it to be classified as a slowly declining type.