Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Small dense particles in the retina observable by spectral-domain optical coherence tomography in age-related macular degeneration

To observe detailed changes in neurosensory retinal structure after anti-VEGF upload in age-related macular degeneration (AMD), by using spectral domain optical coherence tomography (SD-OCT).

Optical properties of a nanomatch-like plasmonic structure

The optical properties of a match-like plasmonic nanostructure are numerically investigated using full-wave finite-difference time-domain analysis in conjunction with dispersive material models. This work is mainly motivated by the developed technique enabling reproducible fabrication of nanomatch structures as well as the growing applications that utilize the localized field enhancement in plasmonic nanostructures. Our research revealed that due to the pronounced field enhancement and larger resonance tunabilities, some nanomatch topologies show potentials for various applications in the field of, e.g., sensing as well as a novel scheme for highly reproducible tips in scanning near field optical microscopy, among others. Despite the additional degrees of freedom that are offered by the composite nature of the proposed nanomatch topology, the paper also reflects on a fundamental complication intrinsic to the material interfaces especially in the nanoscale: stoichiometric mixing. We conclude that the specificity in material modeling will become a significant issue in future research on functionalized composite nanostructures.

Effects on choroidal neovascularization after anti-VEGF Upload using intravitreal ranibizumab, as determined by spectral domain-optical coherence tomography

It is unclear whether anti-VEGF monotherapy in age-related macular degeneration (AMD) achieves morphologic CNV regression or only stops further CNV growth. In this study, spectral domain-optical coherence tomography (SD-OCT) was used to image CNV structure before and after anti-VEGF treatment.

Structure and function of the glucose PTS transporter from Escherichia coli

The glucose transporter IICB of the Escherichia coli phosphotransferase system (PTS) consists of a polytopic membrane domain (IIC) responsible for substrate transport and a hydrophilic C-terminal domain (IIB) responsible for substrate phosphorylation. We have overexpressed and purified a triple mutant of IIC (mut-IIC), which had recently been shown to be suitable for crystallization purposes. Mut-IIC was homodimeric as determined by blue native-PAGE and gel-filtration, and had an eyeglasses-like structure as shown by negative-stain transmission electron microscopy (TEM) and single particle analysis. Glucose binding and transport by mut-IIC, mut-IICB and wildtype-IICB were compared with scintillation proximity and in vivo transport assays. Binding was reduced and transport was impaired by the triple mutation. The scintillation proximity assay allowed determination of substrate binding, affinity and specificity of wildtype-IICB by a direct method. 2D crystallization of mut-IIC yielded highly-ordered tubular crystals and made possible the calculation of a projection structure at 12Å resolution by negative-stain TEM. Immunogold labeling TEM revealed the sidedness of the tubular crystals, and high-resolution atomic force microscopy the surface structure of mut-IIC. This work presents the structure of a glucose PTS transporter at the highest resolution achieved so far and sets the basis for future structural studies.

Fic domain-catalyzed adenylylation: insight provided by the structural analysis of the type IV secretion system effector BepA

Numerous bacterial pathogens subvert cellular functions of eukaryotic host cells by the injection of effector proteins via dedicated secretion systems. The type IV secretion system (T4SS) effector protein BepA from Bartonella henselae is composed of an N-terminal Fic domain and a C-terminal Bartonella intracellular delivery domain, the latter being responsible for T4SS-mediated translocation into host cells. A proteolysis resistant fragment (residues 10-302) that includes the Fic domain shows autoadenylylation activity and adenylyl transfer onto Hela cell extract proteins as demonstrated by autoradiography on incubation with α-[(32)P]-ATP. Its crystal structure, determined to 2.9-Å resolution by the SeMet-SAD method, exhibits the canonical Fic fold including the HPFxxGNGRxxR signature motif with several elaborations in loop regions and an additional β-rich domain at the C-terminus. On crystal soaking with ATP/Mg(2+), additional electron density indicated the presence of a PP(i) /Mg(2+) moiety, the side product of the adenylylation reaction, in the anion binding nest of the signature motif. On the basis of this information and that of the recent structure of IbpA(Fic2) in complex with the eukaryotic target protein Cdc42, we present a detailed model for the ternary complex of Fic with the two substrates, ATP/Mg(2+) and target tyrosine. The model is consistent with an in-line nucleophilic attack of the deprotonated side-chain hydroxyl group onto the α-phosphorus of the nucleotide to accomplish AMP transfer. Furthermore, a general, sequence-independent mechanism of target positioning through antiparallel β-strand interactions between enzyme and target is suggested.

Crystal structure of the Leishmania major phosphodiesterase LmjPDEB1 and insight into the design of the parasite-selective inhibitors

Human leishmaniasis is a major public health problem in many countries, but chemotherapy is in an unsatisfactory state. Leishmania major phosphodiesterases (LmjPDEs) have been shown to play important roles in cell proliferation and apoptosis of the parasite. Thus LmjPDE inhibitors may potentially represent a novel class of drugs for the treatment of leishmaniasis. Reported here are the kinetic characterization of the LmjPDEB1 catalytic domain and its crystal structure as a complex with 3-isobutyl-1-methylxanthine (IBMX) at 1.55 A resolution. The structure of LmjPDEB1 is similar to that of human PDEs. IBMX stacks against the conserved phenylalanine and forms a hydrogen bond with the invariant glutamine, in a pattern common to most inhibitors bound to human PDEs. However, an extensive structural comparison reveals subtle, but significant differences between the active sites of LmjPDEB1 and human PDEs. In addition, a pocket next to the inhibitor binding site is found to be unique to LmjPDEB1. This pocket is isolated by two gating residues in human PDE families, but constitutes a natural expansion of the inhibitor binding pocket in LmjPDEB1. The structure particularity might be useful for the development of parasite-selective inhibitors for the treatment of leishmaniasis.

Two dimensional Lattice Boltzmann simulation of fluid flow through an idealized micro-structure of disordered media

This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

The crystal structure of the platelet activator aggretin reveals a novel (alphabeta)2 dimeric structure

Aggretin is a C-type lectin purified from Calloselasma rhodostoma snake venom. It is a potent activator of platelets, resulting in a collagen-like response by binding and clustering platelet receptor CLEC-2. We present here the crystal structure of aggretin at 1.7 A which reveals a unique tetrameric quaternary structure. The two alphabeta heterodimers are arranged through 2-fold rotational symmetry, resulting in an antiparallel side-by-side arrangement. Aggretin thus presents two ligand binding sites on one surface and can therefore cluster ligands in a manner reminiscent of convulxin and flavocetin. To examine the molecular basis of the interaction with CLEC-2, we used a molecular modeling approach of docking the aggretin alphabeta structure with the CLEC-2 N-terminal domain (CLEC-2N). This model positions the CLEC-2N structure face down in the "saddle"-shaped binding site which lies between the aggretin alpha and beta lectin-like domains. A 2-fold rotation of this complex to generate the aggretin tetramer reveals dimer contacts for CLEC-2N which bring the N- and C-termini into the proximity of each other, and a series of contacts involving two interlocking beta-strands close to the N-terminus are described. A comparison with homologous lectin-like domains from the immunoreceptor family reveals a similar but not identical dimerization mode, suggesting this structure may represent the clustered form of CLEC-2 capable of signaling across the platelet membrane.

Crystal structure of the yeast eIF4A-eIF4G complex: an RNA-helicase controlled by protein-protein interactions

Translation initiation factors eIF4A and eIF4G form, together with the cap-binding factor eIF4E, the eIF4F complex, which is crucial for recruiting the small ribosomal subunit to the mRNA 5' end and for subsequent scanning and searching for the start codon. eIF4A is an ATP-dependent RNA helicase whose activity is stimulated by binding to eIF4G. We report here the structure of the complex formed by yeast eIF4G's middle domain and full-length eIF4A at 2.6-A resolution. eIF4A shows an extended conformation where eIF4G holds its crucial DEAD-box sequence motifs in a productive conformation, thus explaining the stimulation of eIF4A's activity. A hitherto undescribed interaction involves the amino acid Trp-579 of eIF4G. Mutation to alanine results in decreased binding to eIF4A and a temperature-sensitive phenotype of yeast cells that carry a Trp579Ala mutation as its sole source for eIF4G. Conformational changes between eIF4A's closed and open state provide a model for its RNA-helicase activity.

Geology and Geological Structure of the Alberta and Saskatchewan Plains

Alberta, responsible for ninety per cent of Canada's output, had, by 1947 entered into her fifth year of production decline. Only ten per cent of Canada's oil requirements were secured from home fields. Ninety per cent had to be imported, mainly from the United States. How long could imports be maintained on present levels? During the year, the United States had started rationing; in one sector of its domain. Would this become general? If so, what was the answer for Canada?

A novel FIP1L1-PDGFRA mutant destabilizing the inactive conformation of the kinase domain in chronic eosinophilic leukemia/hypereosinophilic syndrome

BACKGROUND: The Fip1-like-1-platelet-derived growth factor receptor alpha (FIP1L1-PDGFRA) gene fusion is a common cause of chronic eosinophilic leukemia (CEL)/hypereosinophilic syndrome (HES), and patients suffering from this particular subgroup of CEL/HES respond to low-dose imatinib therapy. However, some patients may develop imatinib resistance because of an acquired T674I mutation, which is believed to prevent drug binding through steric hindrance. METHODS: In an imatinib resistant FIP1L1-PDGFRA positive patient, we analyzed the molecular structure of the fusion gene and analyzed the effect of several kinase inhibitors on FIP1L1-PDGFRA-mediated proliferative responses in vitro. RESULTS: Sequencing of the FIP1L1-PDGFRA fusion gene revealed the occurrence of a S601P mutation, which is located within the nucleotide binding loop. In agreement with the clinical observations, imatinib did not inhibit the proliferation of S601P mutant FIP1L1-PDGFRA-transduced Ba/F3 cells. Moreover, sorafenib, which has been described to inhibit T674I mutant FIP1L1-PDGFRA, failed to block S601P mutant FIP1L1-PDGFRA. Structural modeling revealed that the newly identified S601P mutated form of PDGFRA destabilizes the inactive conformation of the kinase domain that is necessary to bind imatinib as well as sorafenib. CONCLUSIONS: We identified a novel mutation in FIP1L1-PDGFRA resulting in both imatinib and sorafenib resistance. The identification of novel drug-resistant FIP1L1-PDGFRA variants may help to develop the next generation of target-directed compounds for CEL/HES and other leukemias.

Reanalysis-driven climate simulation over CORDEX North America domain using the Canadian Regional Climate Model, version 5: model performance evaluation

The performance of reanalysis-driven Canadian Regional Climate Model, version 5 (CRCM5) in reproducing the present climate over the North American COordinated Regional climate Downscaling EXperiment domain for the 1989–2008 period has been assessed in comparison with several observation-based datasets. The model reproduces satisfactorily the near-surface temperature and precipitation characteristics over most part of North America. Coastal and mountainous zones remain problematic: a cold bias (2–6 °C) prevails over Rocky Mountains in summertime and all year-round over Mexico; winter precipitation in mountainous coastal regions is overestimated. The precipitation patterns related to the North American Monsoon are well reproduced, except on its northern limit. The spatial and temporal structure of the Great Plains Low-Level Jet is well reproduced by the model; however, the night-time precipitation maximum in the jet area is underestimated. The performance of CRCM5 was assessed against earlier CRCM versions and other RCMs. CRCM5 is shown to have been substantially improved compared to CRCM3 and CRCM4 in terms of seasonal mean statistics, and to be comparable to other modern RCMs.