Modeling and simulation of fibrinogen and its adsorption behavior

In this thesis different approaches for the modeling and simulation of the blood protein fibrinogen are presented. The approaches are meant to systematically connect the multiple time and length scales involved in the dynamics of fibrinogen in solution and at inorganic surfaces. The first part of the thesis will cover simulations of fibrinogen on an all atom level. Simulations of the fibrinogen protomer and dimer are performed in explicit solvent to characterize the dynamics of fibrinogen in solution. These simulations reveal an unexpectedly large and fast bending motion that is facilitated by molecular hinges located in the coiled-coil region of fibrinogen. This behavior is characterized by a bending and a dihedral angle and the distribution of these angles is measured. As a consequence of the atomistic detail of the simulations it is possible to illuminate small scale behavior in the binding pockets of fibrinogen that hints at a previously unknown allosteric effect. In a second step atomistic simulations of the fibrinogen protomer are performed at graphite and mica surfaces to investigate initial adsorption stages. These simulations highlight the different adsorption mechanisms at the hydrophobic graphite surface and the charged, hydrophilic mica surface. It is found that the initial adsorption happens in a preferred orientation on mica. Many effects of practical interest involve aggregates of many fibrinogen molecules. To investigate such systems, time and length scales need to be simulated that are not attainable in atomistic simulations. It is therefore necessary to develop lower resolution models of fibrinogen. This is done in the second part of the thesis. First a systematically coarse grained model is derived and parametrized based on the atomistic simulations of the first part. In this model the fibrinogen molecule is represented by 45 beads instead of nearly 31,000 atoms. The intra-molecular interactions of the beads are modeled as a heterogeneous elastic network while inter-molecular interactions are assumed to be a combination of electrostatic and van der Waals interaction. A method is presented that determines the charges assigned to beads by matching the electrostatic potential in the atomistic simulation. Lastly a phenomenological model is developed that represents fibrinogen by five beads connected by rigid rods with two hinges. This model only captures the large scale dynamics in the atomistic simulations but can shed light on experimental observations of fibrinogen conformations at inorganic surfaces.

Test sperimentali e modellazione di una torre di lavaggio per la filtrazione di particolato sottile in impianti di combustione delle biomasse

This study investigates biomass and particulate matter also known as PM produced from the combustion of a domestic boiler powered by mais and how to separate PM from the stream of smoke output from the boiler using wet scrubber with structured packing. Sperimentations show the inefficiency of the separator used, so we provide an optimization of the structured packing changing geometric parameters as angle of the bend or thickness of the channels. In order to obtain a higher separation efficiency we remove the structured packinkg and introduce a packed bed composed of spheres of polyethylene with a diameter of 3 mm.

pH-dependent distribution of chlorin e6 derivatives across phospholipid bilayers probed by NMR spectroscopy

The pH-dependent membrane adsorption and distribution of three chlorin derivatives, chlorin e6 (CE), rhodin G7 (RG), and monoaspartyl-chlorin e6 (MACE), in the physiological pH range (pH 6-8) were probed by NMR spectroscopy. Unilamellar vesicles consisting of dioleoyl-phosphatidyl-choline (DOPC) were used as membrane models. The chlorin derivatives were characterized with respect to their aggregation behavior, the pK(a) values of individual carboxylate groups, the extent of membrane adsorption, and their flip-flop rates across the bilayer membrane for pH 6-8. External membrane adsorption was found to be lower for RG than for CE and MACE. Both electrostatic interactions and the extent of aggregation seemed to be the main determinants of membrane adsorption. Rate constants for chlorin transfer across the membrane were found to correlate strongly with the pH of the surrounding medium, in particular, for CE and RG. In acidic solution, CE and RG transfer across the membrane was strongly accelerated, and in basic solution, all compounds were retained, mostly in the outer monolayer. In contrast, MACE flip-flop across the membrane remained very low even at pH 6. The protonation of ionizable groups is suggested to be a major determinant of chlorin transfer rates across the bilayer. pK(a) values of CE and RG were found to be between 6 and 8, and two of the carboxylate groups in MACE had pK(a) values below 6. For CE and RG, the kinetic profiles at acidic pH indicated that the initial fast membrane distribution was followed by secondary steps that are discussed in this article.

The contribution of surface residues to membrane binding and ligand transfer by the -tocopherol transfer protein ( -TTP)

Previous work has shown that the -tocopherol transfer protein ( -TTP) can bind to vesicular or immobilized phospholipid membranes. Revealing the molecular mechanisms by which -TTP associates with membranes is thought to be critical to understanding its function and role in the secretion of tocopherol from hepatocytes into the circulation. Calculations presented in the Orientations of Proteins in Membranes database have provided a testable model for the spatial arrangement of -TTP and other CRAL-TRIO family proteins with respect to the lipid bilayer. These calculations predicted that a hydrophobic surface mediates the interaction of -TTP with lipid membranes. To test the validity of these predictions, we used site-directed mutagenesis and examined the substituted mutants with regard to intermembrane ligand transfer, association with lipid layers and biological activity in cultured hepatocytes. Substitution of residues in helices A8 (F165A and F169A) and A10 (I202A, V206A and M209A) decreased the rate of intermembrane ligand transfer as well as protein adsorption to phospholipid bilayers. The largest impairment was observed upon mutation of residues that are predicted to be fully immersed in the lipid bilayer in both apo (open) and holo (closed) conformations such as Phe165 and Phe169. Mutation F169A, and especially F169D, significantly impaired -TTP-assisted secretion of -tocopherol outside cultured hepatocytes. Mutation of selected basic residues (R192H, K211A, and K217A) had little effect on transfer rates, indicating no significant involvement of nonspecific electrostatic interactions with membranes.

A tetrathiafulvalene-functionalized naphthalene diimide: synthesis, electrochemical and photophysical properties

A tetrathiafulvalene donor has been attached to the naphthalene diimide core via a rigid bridge affording a new planar molecular dyad. Its electronic properties have been studied experimentally by the combination of electrochemistry and UV-vis-NIR spectroscopy. Various electronic excited charge-transfer states are generated in different oxidation states, leading to almost full absorption in the visible to near-IR region with high extinction coefficients. The observed electronic properties are explained on the basis of density-functional-theory. In particular, the oxidized radical species show a strong tendency to undergo aggregation, in which the long-distance attractive interactions overcome the electrostatic repulsions. (C) 2011 Elsevier Ltd. All rights reserved.

Regulatory mechanism of the light-activable allosteric switch LOV-TAP for the control of DNA binding: a computer simulation study

The spatio-temporal control of gene expression is fundamental to elucidate cell proliferation and deregulation phenomena in living systems. Novel approaches based on light-sensitive multiprotein complexes have recently been devised, showing promising perspectives for the noninvasive and reversible modulation of the DNA-transcriptional activity in vivo. This has lately been demonstrated in a striking way through the generation of the artificial protein construct light-oxygen-voltage (LOV)-tryptophan-activated protein (TAP), in which the LOV-2-Jα photoswitch of phototropin1 from Avena sativa (AsLOV2-Jα) has been ligated to the tryptophan-repressor (TrpR) protein from Escherichia coli. Although tremendous progress has been achieved on the generation of such protein constructs, a detailed understanding of their functioning as opto-genetical tools is still in its infancy. Here, we elucidate the early stages of the light-induced regulatory mechanism of LOV-TAP at the molecular level, using the noninvasive molecular dynamics simulation technique. More specifically, we find that Cys450-FMN-adduct formation in the AsLOV2-Jα-binding pocket after photoexcitation induces the cleavage of the peripheral Jα-helix from the LOV core, causing a change of its polarity and electrostatic attraction of the photoswitch onto the DNA surface. This goes along with the flexibilization through unfolding of a hairpin-like helix-loop-helix region interlinking the AsLOV2-Jα- and TrpR-domains, ultimately enabling the condensation of LOV-TAP onto the DNA surface. By contrast, in the dark state the AsLOV2-Jα photoswitch remains inactive and exerts a repulsive electrostatic force on the DNA surface. This leads to a distortion of the hairpin region, which finally relieves its tension by causing the disruption of LOV-TAP from the DNA.

Structure-activity relationships from in vitro efficacies of the thiazolide series against the intracellular apicomplexan protozoan Neospora caninum

Nitazoxanide (2-acetolyloxy-N-(5-nitro 2-thiazolyl) benzamide; NTZ) represents the parent compound of a novel class of broad-spectrum anti-parasitic compounds named thiazolides. NTZ is active against a wide variety of intestinal and tissue-dwelling helminths, protozoa, enteric bacteria and a number of viruses infecting animals and humans. While potent, this poses a problem in practice, since this obvious non-selectivity can lead to undesired side effects in both humans and animals. In this study, we used real time PCR to determine the in vitro activities of 29 different thiazolides (NTZ-derivatives), which carry distinct modifications on both the thiazole- and the benzene moieties, against the tachyzoite stage of the intracellular protozoan Neospora caninum. The goal was to identify a highly active compound lacking the undesirable nitro group, which would have a more specific applicability, such as in food animals. By applying self-organizing molecular field analysis (SOMFA), these data were used to develop a predictive model for future drug design. SOMFA performs self-alignment of the molecules, and takes into account the steric and electrostatic properties, in order to determine 3D-quantitative structure activity relationship models. The best model was obtained by overlay of the thiazole moieties. Plotting of predicted versus experimentally determined activity produced an r2 value of 0.8052 and cross-validation using the "leave one out" methodology resulted in a q2 value of 0.7987. A master grid map showed that large steric groups at the R2 position, the nitrogen of the amide bond and position Y could greatly reduce activity, and the presence of large steric groups placed at positions X, R4 and surrounding the oxygen atom of the amide bond, may increase the activity of thiazolides against Neospora caninum tachyzoites. The model obtained here will be an important predictive tool for future development of this important class of drugs.

Pseudodominant Inheritance of Goitrous Congenital Hypothyroidism Caused by TPO Mutations: Molecular and in Silico Studies

Context and Objective: Most cases of goitrous congenital hypothyroidism (CH) from thyroid dyshormonogenesis 1) follow a recessive mode of inheritance and 2) are due to mutations in the thyroid peroxidase gene (TPO). We report the genetic mechanism underlying the apparently dominant inheritance of goitrous CH in a nonconsanguineous family of French Canadian origin. Design, Setting, and Participants: Two brothers identified by newborn TSH screening had severe hypothyroidism and a goiter with increased (99m)Tc uptake. The mother was euthyroid, but the father and two paternal uncles had also been diagnosed with goitrous CH. After having excluded PAX8 gene mutations, we hypothesized that the underlying defect could be TPO mutations. Results: Both compound heterozygous siblings had inherited a mutant TPO allele carried by their mother (c.1496delC; p.Pro499Argfs2X), and from their father, one brother had inherited a missense mutation (c.1978C-->G; p.Gln660Glu) and the other an insertion (c.1955insT; p.Phe653Valfs15X). The thyroid gland of one uncle who is a compound heterozygote for TPO mutations (p.Phe653Valfs15X/p.Gln660Glu) was removed because of concurrent multiple endocrine neoplasia type 2A. Immunohistochemistry revealed normal TPO staining, implying that Gln660Glu TPO is expressed properly. Modeling of this mutant in silico suggests that its three-dimensional structure is conserved, whereas the electrostatic binding energy between the Gln660Glu TPO and its heme group becomes repulsive. Conclusion: We report a pedigree presenting with pseudodominant goitrous CH due to segregation of three different TPO mutations. Although goitrous CH generally follows a recessive mode of inheritance, the high frequency of TPO mutations carriers may lead to pseudodominant inheritance.

Electrospun chitosan nanofibers for virus removal

Electrospinning uses electrostatic forces to create nanofibers that are far smaller than conventional fiber spinning process. Nanofibers made with chitosan were created and techniques to control fibers diameter and were well developed. However, the adsorption of porcine parvovirus (PPV) was low. PPV is a small, nonenveloped virus that is difficult to remove due to its size, 18-26 nm in diameter, and its chemical stability. To improve virus adsorption, we functionalized the nanofibers with a quaternized amine, forming N-[(2-hydroxy-3-trimethylammonium) propyl] chitosan chloride (HTCC). This was blended with additives to increase the ability to form HTCC nanofibers. The additives changed the viscosity and conductivity of the electrospinning solution. We have successfully synthesized and functionalized HTCC nanofibers that absorb PPV. HTCC blend with graphene have the ability to remove a minimum of 99% of PPV present in solution.

Laser Thomson scattering measurements of electron temperature and density in a Hall-effect plasma

Hall-effect thrusters (HETs) are compact electric propulsion devices with high specific impulse used for a variety of space propulsion applications. HET technology is well developed but the electron properties in the discharge are not completely understood, mainly due to the difficulty involved in performing accurate measurements in the discharge. Measurements of electron temperature and density have been performed using electrostatic probes, but presence of the probes can significantly disrupt thruster operation, and thus alter the electron temperature and density. While fast-probe studies have expanded understanding of HET discharges, a non-invasive method of measuring the electron temperature and density in the plasma is highly desirable. An alternative to electrostatic probes is a non-perturbing laser diagnostic technique that measures Thomson scattering from the plasma. Thomson scattering is the process by which photons are elastically scattered from the free electrons in a plasma. Since the electrons have thermal energy their motion causes a Doppler shift in the scattered photons that is proportional to their velocity. Like electrostatic probes, laser Thomson scattering (LTS) can be used to determine the temperature and density of free electrons in the plasma. Since Thomson scattering measures the electron velocity distribution function directly no assumptions of the plasma conditions are required, allowing accurate measurements in anisotropic and non-Maxwellian plasmas. LTS requires a complicated measurement apparatus, but has the potential to provide accurate, non-perturbing measurements of electron temperature and density in HET discharges. In order to assess the feasibility of LTS diagnostics on HETs non-invasive measurements of electron temperature and density in the near-field plume of a Hall thruster were performed using a custom built laser Thomson scattering diagnostic. Laser measurements were processed using a maximum likelihood estimation method and results were compared to conventional electrostatic double probe measurements performed at the same thruster conditions. Electron temperature was found to range from approximately 1 – 40 eV and density ranged from approximately 1.0 x 1017 m-3 to 1.3 x 1018 m-3 over discharge voltages from 250 to 450 V and mass flow rates of 40 to 80 SCCM using xenon propellant.

Dynamics of settling charged particles in turbulence : theory and experiments

It has been proposed that inertial clustering may lead to an increased collision rate of water droplets in clouds. Atmospheric clouds and electrosprays contain electrically charged particles embedded in turbulent flows, often under the influence of an externally imposed, approximately uniform gravitational or electric force. In this thesis, we present the investigation of charged inertial particles embedded in turbulence. We have developed a theoretical description for the dynamics of such systems of charged, sedimenting particles in turbulence, allowing radial distribution functions to be predicted for both monodisperse and bidisperse particle size distributions. The governing parameters are the particle Stokes number (particle inertial time scale relative to turbulence dissipation time scale), the Coulomb-turbulence parameter (ratio of Coulomb ’terminalar speed to turbulence dissipation velocity scale), and the settling parameter (the ratio of the gravitational terminal speed to turbulence dissipation velocity scale). For the monodispersion particles, The peak in the radial distribution function is well predicted by the balance between the particle terminal velocity under Coulomb repulsion and a time-averaged ’drift’ velocity obtained from the nonuniform sampling of fluid strain and rotation due to finite particle inertia. The theory is compared to measured radial distribution functions for water particles in homogeneous, isotropic air turbulence. The radial distribution functions are obtained from particle positions measured in three dimensions using digital holography. The measurements support the general theoretical expression, consisting of a power law increase in particle clustering due to particle response to dissipative turbulent eddies, modulated by an exponential electrostatic interaction term. Both terms are modified as a result of the gravitational diffusion-like term, and the role of ’gravity’ is explored by imposing a macroscopic uniform electric field to create an enhanced, effective gravity. The relation between the radial distribution functions and inward mean radial relative velocity is established for charged particles.

MULTISCALE EXAMINATION AND MODELING OF ELECTRON TRANSPORT IN NANOSCALE MATERIALS AND DEVICES

For half a century the integrated circuits (ICs) that make up the heart of electronic devices have been steadily improving by shrinking at an exponential rate. However, as the current crop of ICs get smaller and the insulating layers involved become thinner, electrons leak through due to quantum mechanical tunneling. This is one of several issues which will bring an end to this incredible streak of exponential improvement of this type of transistor device, after which future improvements will have to come from employing fundamentally different transistor architecture rather than fine tuning and miniaturizing the metal-oxide-semiconductor field effect transistors (MOSFETs) in use today. Several new transistor designs, some designed and built here at Michigan Tech, involve electrons tunneling their way through arrays of nanoparticles. We use a multi-scale approach to model these devices and study their behavior. For investigating the tunneling characteristics of the individual junctions, we use a first-principles approach to model conduction between sub-nanometer gold particles. To estimate the change in energy due to the movement of individual electrons, we use the finite element method to calculate electrostatic capacitances. The kinetic Monte Carlo method allows us to use our knowledge of these details to simulate the dynamics of an entire device— sometimes consisting of hundreds of individual particles—and watch as a device ‘turns on’ and starts conducting an electric current. Scanning tunneling microscopy (STM) and the closely related scanning tunneling spectroscopy (STS) are a family of powerful experimental techniques that allow for the probing and imaging of surfaces and molecules at atomic resolution. However, interpretation of the results often requires comparison with theoretical and computational models. We have developed a new method for calculating STM topographs and STS spectra. This method combines an established method for approximating the geometric variation of the electronic density of states, with a modern method for calculating spin-dependent tunneling currents, offering a unique balance between accuracy and accessibility.

Virus purification, detection and removal

The biopharmaceutical industry has a growing demand and an increasing need to improve the current virus purification technologies, especially as more and more vaccines are produced from cell-culture derived virus particles. Downstream purification strategies can be expensive and account for 70% of the overall manufacturing costs. The economic pressure and purification processes can be particularly challenging when the virus to be purified is small, as in our model virus, porcine parvovirus (PPV). Our efforts are focused on designing an easy, economical, scalable and efficient system for virus purification, and we focused on aqueous two-phase systems. Industry acceptable standards for virus vaccine recovery can be as low as 30% due to demand of high final titer, virus transduction inhibitors and presence of empty or defective virus capsids as impurities. We have overcome these shortcomings by recovering a high 64% of infectious virus using an aqueous two-phase system. We used high molecular weight polymer and citrate salt to achieve a good yield and eliminated the major contaminant bovine serum albumin. Viruses are also studied for ensuring pure and safe drinking water. Low pressure microfiltrations are continuously being investigated for water filters as they allow high permeate flux and low fouling. Viruses such as PPV are small enough to pass through the microporous membranes. Control of viruses in water is crucial for public health and we have designed an affinity based membrane filter to capture virus. Nanofibers have a high surface to volume ratio providing a highly accessible surface area for virus adsorption. Chitosan an insoluble, biocompatible and biodegradable polymer was used for adsorbing trimer peptide WRW. About 0.2 μmoles of cysteine terminal WRW peptide was conjugated to amine terminal chitosan using maleimide conjugation chemistry. We achieved 90-99% virus removal from water adjusted to a neutral pH. The virus removal from affinity based chitosan was attributed to electrostatic and hydrophobic driven binding effect.

The Effects of Zinc Sulphate, Manganese Sulphate, and Iron Sulphate on the Lead Storage Cell.

The lead storage battery as it is used today is made up of the pasted type plates of lead dioxide, the anode, and sponge lead, the cathode, and wooden or hard rubber separators, which serve to insulate these from one another. In manufacturing these, it is desirable to keep them free from impurities.

A MEMS Light Modulator Based on Diffractive Nanohole Gratings

We present the design, fabrication, and testing of a microelectromechanical systems (MEMS) light modulator based on pixels patterned with periodic nanohole arrays. Flexure-suspended silicon pixels are patterned with a two dimensional array of 150 nm diameter nanoholes using nanoimprint lithography. A top glass plate assembled above the pixel array is used to provide a counter electrode for electrostatic actuation. The nanohole pattern is designed so that normally-incident light is coupled into an in-plane grating resonance, resulting in an optical stop-band at a desired wavelength. When the pixel is switched into contact with the top plate, the pixel becomes highly reflective. A 3:1 contrast ratio at the resonant wavelength is demonstrated for gratings patterned on bulk Si substrates. The switching time is 0.08 ms and the switching voltage is less than 15V.