A molecular dynamics investigation of the structure of PbO---PbF2 glasses*1

The structure of PbO---PbF2 glasses has been studied using molecular dynamics (MD). The existence of [OPb4] structural units is observed over the entire glass-forming range, in conformity with a model proposed earlier based on various structural investigations of this system. Various other features of the structural model are also supported by the MD calculations.

Molecular basis of the kidney filtration barrier : Role of the nephrin protein complex

The kidney filtration barrier consists of fenestrated endothelial cell layer, glomerular basement membrane and slit diaphragm (SD), the specialized junction between glomerular viscelar epithelial cells (podocytes). Podocyte injury is associated with the development of proteinuria, and if not reversed the injury will lead to permanent deterioration of the glomerular filter. The early events are characterized by disruption of the integrity of the SD, but the molecular pathways involved are not fully understood. Congenital nephrotic syndrome of the Finnish type (CNF) is caused by mutations in NPHS1, the gene encoding the SD protein nephrin. Lack of nephrin results in loss of the SD and massive proteinuria beginning before birth. Furthermore, nephrin expression is decreased in acquired human kidney diseases including diabetic nephropathy. This highlights the importance of nephrin and consequently SD in regulating the kidney filtration function. However, the precise molecular mechanism of how nephrin is involved in the formation of the SD is unknown. This thesis work aimed at clarifying the role of nephrin and its interaction partners in the formation of the SD. The purpose was to identify novel proteins that associate with nephrin in order to define the essential molecular complex required for the establishment of the SD. The aim was also to decipher the role of novel nephrin interacting proteins in podocytes. Nephrin binds to nephrin-like proteins Neph1 and Neph2, and to adherens junction protein P-cadherin. These interactions have been suggested to play a role in the formation of the SD. In this thesis work, we identified densin as a novel interaction partner for nephrin. Densin was localized to the SD and it was shown to bind to adherens junction protein beta-catenin. Furthermore, densin was shown to behave in a similar fashion as adherens junction proteins in cell-cell contacts. These results indicate that densin may play a role in cell adhesion and, therefore, may contribute to the formation of the SD together with nephrin and adherens junction proteins. Nephrin was also shown to bind to Neph3, which has been previously localized to the SD. Neph3 and Neph1 were shown to induce cell adhesion alone, whereas nephrin needed to trans-interact with Neph1 or Neph3 from the opposite cell surface in order to make cell-cell contacts. This was associated with the decreased tyrosine phosphorylation of nephrin. These data extend the current knowledge of the molecular composition of the nephrin protein complex at the SD and also provide novel insights of how the SD may be formed. This thesis work also showed that densin was up-regulated in the podocytes of CNF patients. Neph3 was up-regulated in nephrin deficient mouse kidneys, which share similar podocyte alterations and lack of the SD as observed in CNF patients podocytes. These data suggest that densin and Neph3 may have a role in the formation of morphological alterations in podocytes detected in CNF patients. Furthermore, this thesis work showed that deletion of beta-catenin specifically from adult mouse podocytes protected the mice from the development of adriamycin-induced podocyte injury and proteinuria compared to wild-type mice. These results show that beta-catenin play a role in the adriamycin induced podocyte injury. Podocyte injury is a hallmark in many kidney diseases and the changes observed in the podocytes of CNF patient share characteristics with injured podocytes observed in chronic kidney diseases. Therefore, the results obtained in this thesis work suggest that densin, Neph3 and beta-catenin participate in the molecular pathways which result in morphological alterations commonly detected in injured podocytes in kidney diseases.

Animal model and molecular interactions of Cln5

Neuronal ceroid lipofuscinoses (NCLs) are a family of inherited pediatric neurodegenerative disorders, leading to retinal degeneration, death of selective neuronal populations and accumulation of autofluorscent ceroid-lipopigments. The clinical manifestations are generally similar in all forms. The Finnish variant late infantile neuronal ceroid lipofuscinosis (vLINCLFin) is a form of NCL, especially enriched in the Finnish population. The aim of this thesis was to analyse the brain pathology of vLINCLFin utilising the novel Cln5-/- mouse model. Gene expression profiling of the brains of already symptomatic Cln5-/- mice revealed that inflammation, neurodegeneration and defects in myelinization are the major characteristics of the later stages of the disease. Histological characterization of the brain pathology confirmed that the thalamocortical system is affected in Cln5-/- mice, similarly to the other NCL mouse models. However, whereas the brain pathology in all other analyzed NCL mice initiate in the thalamus and spread only months later to the cortex, we observed that the sequence of events is uniquely reversed in Cln5-/- mice; beginning in the cortex and spreading to the thalamus only months later. We could also show that even though neurodegeneration is inititated in the cortex, reactive gliosis and loss of myelin are evident in specific nuclei of the thalamus already in the 1 month old brain. To obtain a deeper insight into the disturbed metabolic pathways, we performed gene expression profiling of presymptomatic mouse brains. We validated these findings with immunohistological analyses, and could show that cytoskeleton and myelin were affected in Cln5-/- mice. Comparison of gene expression profiling results of Cln5-/- and Cln1-/- mice, further highlighted that these two NCL models share a common defective pathway, leading to disturbances in the neuronal growth cone and cytoskeleton. Encouraged by the evidence of this defected pathway, we analyzed the molecular interactions of NCL-proteins and observed that Cln5 and Cln1/Ppt1 proteins interact with each other. Furthermore, we demonstrated that Cln5 and Cln1/Ppt1 share an interaction partner, the F1-ATP synthase, potentially linking both vLINCLFIN and INCL diseases to disturbed lipid metabolism. In addition, Cln5 was shown to interact with other NCL proteins; Cln2, Cln3, Cln6 and Cln8, implicating a central role for Cln5 in the NCL pathophysiology. This study is the first to describe the brain pathology and gene expression changes in the Cln5-/- mouse. Together the findings presented in this thesis represent novel information of the disease processes and the molecular mechanisms behind vLINCLFin and have highlighted that vLINCLFin forms a very important model to analyze the pathophysiology of NCL diseases.

High-throughput prediction of eucalypt lignin syringyl/guaiacyl content using multivariate analysis: a comparison between mid-infrared, near-infrared, and Raman spectroscopies for model development

BACKGROUND: In order to rapidly and efficiently screen potential biofuel feedstock candidates for quintessential traits, robust high-throughput analytical techniques must be developed and honed. The traditional methods of measuring lignin syringyl/guaiacyl (S/G) ratio can be laborious, involve hazardous reagents, and/or be destructive. Vibrational spectroscopy can furnish high-throughput instrumentation without the limitations of the traditional techniques. Spectral data from mid-infrared, near-infrared, and Raman spectroscopies was combined with S/G ratios, obtained using pyrolysis molecular beam mass spectrometry, from 245 different eucalypt and Acacia trees across 17 species. Iterations of spectral processing allowed the assembly of robust predictive models using partial least squares (PLS). RESULTS: The PLS models were rigorously evaluated using three different randomly generated calibration and validation sets for each spectral processing approach. Root mean standard errors of prediction for validation sets were lowest for models comprised of Raman (0.13 to 0.16) and mid-infrared (0.13 to 0.15) spectral data, while near-infrared spectroscopy led to more erroneous predictions (0.18 to 0.21). Correlation coefficients (r) for the validation sets followed a similar pattern: Raman (0.89 to 0.91), mid-infrared (0.87 to 0.91), and near-infrared (0.79 to 0.82). These statistics signify that Raman and mid-infrared spectroscopy led to the most accurate predictions of S/G ratio in a diverse consortium of feedstocks. CONCLUSION: Eucalypts present an attractive option for biofuel and biochemical production. Given the assortment of over 900 different species of Eucalyptus and Corymbia, in addition to various species of Acacia, it is necessary to isolate those possessing ideal biofuel traits. This research has demonstrated the validity of vibrational spectroscopy to efficiently partition different potential biofuel feedstocks according to lignin S/G ratio, significantly reducing experiment and analysis time and expense while providing non-destructive, accurate, global, predictive models encompassing a diverse array of feedstocks.

Discovering hidden structures in molecular data using a Bayesian partition model approach

Advancements in the analysis techniques have led to a rapid accumulation of biological data in databases. Such data often are in the form of sequences of observations, examples including DNA sequences and amino acid sequences of proteins. The scale and quality of the data give promises of answering various biologically relevant questions in more detail than what has been possible before. For example, one may wish to identify areas in an amino acid sequence, which are important for the function of the corresponding protein, or investigate how characteristics on the level of DNA sequence affect the adaptation of a bacterial species to its environment. Many of the interesting questions are intimately associated with the understanding of the evolutionary relationships among the items under consideration. The aim of this work is to develop novel statistical models and computational techniques to meet with the challenge of deriving meaning from the increasing amounts of data. Our main concern is on modeling the evolutionary relationships based on the observed molecular data. We operate within a Bayesian statistical framework, which allows a probabilistic quantification of the uncertainties related to a particular solution. As the basis of our modeling approach we utilize a partition model, which is used to describe the structure of data by appropriately dividing the data items into clusters of related items. Generalizations and modifications of the partition model are developed and applied to various problems. Large-scale data sets provide also a computational challenge. The models used to describe the data must be realistic enough to capture the essential features of the current modeling task but, at the same time, simple enough to make it possible to carry out the inference in practice. The partition model fulfills these two requirements. The problem-specific features can be taken into account by modifying the prior probability distributions of the model parameters. The computational efficiency stems from the ability to integrate out the parameters of the partition model analytically, which enables the use of efficient stochastic search algorithms.

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.

Diaryl dihydropyrazole-3-carboxamides with significant in vivo antiobesity activity related to CB1 receptor antagonism: Synthesis, biological evaluation, and molecular modeling in the homology model

A number of analogues of diaryl dihydropyrazole-3-carboxamides have been synthesized. Their activities were evaluated for appetite suppression and body weight reduction in animal models. Depending on the chemical modification of the selected dihydropyrazole scaffold, the lead compoundsthe bisulfate salt of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 26 and the bisulfate salt of (−)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 30showed significant body weight reduction in vivo, which is attributed to their CB1 antagonistic activity and exhibited a favorable pharmacokinetic profile. The molecular modeling studies also showed interactions of two isomers of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 9 with CB1 receptor in the homology model similar to those of N-piperidino-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4-methyl-3-pyrazole-carboxamide (rimonabant) 1 and 4S-(−)-3-(4-chlorophenyl)-N-methyl-N‘-[(4-chlorophenyl)-sulfonyl]-4-phenyl-4,5-dihydro-1H-pyrazole-1-carboxamidine (SLV-319) 2.

Inflammation-driven bone formation in a mouse model of ankylosing spondylitis: Sequential not parallel processes

Background Ankylosing spondylitis (AS) is an immune-mediated arthritis particularly targeting the spine and pelvis and is characterised by inflammation, osteoproliferation and frequently ankylosis. Current treatments that predominately target inflammatory pathways have disappointing efficacy in slowing disease progression. Thus, a better understanding of the causal association and pathological progression from inflammation to bone formation, particularly whether inflammation directly initiates osteoproliferation, is required. Methods The proteoglycan-induced spondylitis (PGISp) mouse model of AS was used to histopathologically map the progressive axial disease events, assess molecular changes during disease progression and define disease progression using unbiased clustering of semi-quantitative histology. PGISp mice were followed over a 24-week time course. Spinal disease was assessed using a novel semi-quantitative histological scoring system that independently evaluated the breadth of pathological features associated with PGISp axial disease, including inflammation, joint destruction and excessive tissue formation (osteoproliferation). Matrix components were identified using immunohistochemistry. Results Disease initiated with inflammation at the periphery of the intervertebral disc (IVD) adjacent to the longitudinal ligament, reminiscent of enthesitis, and was associated with upregulated tumor necrosis factor and metalloproteinases. After a lag phase, established inflammation was temporospatially associated with destruction of IVDs, cartilage and bone. At later time points, advanced disease was characterised by substantially reduced inflammation, excessive tissue formation and ectopic chondrocyte expansion. These distinct features differentiated affected mice into early, intermediate and advanced disease stages. Excessive tissue formation was observed in vertebral joints only if the IVD was destroyed as a consequence of the early inflammation. Ectopic excessive tissue was predominantly chondroidal with chondrocyte-like cells embedded within collagen type II- and X-rich matrix. This corresponded with upregulation of mRNA for cartilage markers Col2a1, sox9 and Comp. Osteophytes, though infrequent, were more prevalent in later disease. Conclusions The inflammation-driven IVD destruction was shown to be a prerequisite for axial disease progression to osteoproliferation in the PGISp mouse. Osteoproliferation led to vertebral body deformity and fusion but was never seen concurrent with persistent inflammation, suggesting a sequential process. The findings support that early intervention with anti-inflammatory therapies will be needed to limit destructive processes and consequently prevent progression of AS.

Anomalous reaction-diffusion as a model of nonexponential DNA escape kinetics

We show that data from recent experiments carried out on the kinetics of DNA escape from alpha-hemolysin nanopores [M. Wiggin, C. Tropini, C. T. Cossa, N. N. Jetha, and A. Marziali, Biophys. J. 95, 5317 (2008)] may be rationalized by a model of chain dynamics based on the anomalous diffusion of a particle moving in a harmonic well in the presence of a delta function sink. The experiments of Wiggin found, among other things, that the occasional occurrence of unusually long escape times in the distribution of chain trapping events led to nonexponential decays in the survival probability, S(t), of the DNA molecules within the nanopore. Wiggin ascribed this nonexponentiality to the existence of a distribution of trapping potentials, which they suggested was theresult of stochastic interactions between the bases of the DNA and the amino acids located on the surface of the nanopore. Based on this idea, they showed that the experimentally determined S(t) could be well fit in both the short and long time regimes by a function of the form (1+t/tau)(-alpha) (the so called Becquerel function). In our model, S(t) is found to be given by a Mittag-Leffler function at short times and by a generalized Mittag-Leffler function at long times. By suitable choice of certain parameter values, these functions are found to fit the experimental S(t) even better than the Becquerel function. Anomalous diffusion of DNA within the trap prior to escape over a barrier of fixed height may therefore provide a second, plausible explanation of the data, and may offer fresh perspectives on similar trapping and escape problems.

An analytical model for performance evaluation of multimedia applications over EDCA in an IEEE 802.11e WLAN

We extend the modeling heuristic of (Harsha et al. 2006. In IEEE IWQoS 06, pp 178 - 187) to evaluate the performance of an IEEE 802.11e infrastructure network carrying packet telephone calls, streaming video sessions and TCP controlled file downloads, using Enhanced Distributed Channel Access (EDCA). We identify the time boundaries of activities on the channel (called channel slot boundaries) and derive a Markov Renewal Process of the contending nodes on these epochs. This is achieved by the use of attempt probabilities of the contending nodes as those obtained from the saturation fixed point analysis of (Ramaiyan et al. 2005. In Proceedings ACM Sigmetrics, `05. Journal version accepted for publication in IEEE TON). Regenerative analysis on this MRP yields the desired steady state performance measures. We then use the MRP model to develop an effective bandwidth approach for obtaining a bound on the size of the buffer required at the video queue of the AP, such that the streaming video packet loss probability is kept to less than 1%. The results obtained match well with simulations using the network simulator, ns-2. We find that, with the default IEEE 802.11e EDCA parameters for access categories AC 1, AC 2 and AC 3, the voice call capacity decreases if even one streaming video session and one TCP file download are initiated by some wireless station. Subsequently, reducing the voice calls increases the video downlink stream throughput by 0.38 Mbps and file download capacity by 0.14 Mbps, for every voice call (for the 11 Mbps PHY). We find that a buffer size of 75KB is sufficient to ensure that the video packet loss probability at the QAP is within 1%.

Position preference and diffusion path of an oxygen ion in apatite-type lanthanum silicate La9.33Si6O26: A density functional study

Using density functional theory, we investigated the position preference and diffusion mechanisms of interstitial oxygen ions in lanthanum silicate La9.33Si6O26, which is an apatite-structured oxide and a promising candidate electrolyte material for solid oxide fuel cells. The reported lanthanum vacancies were explicitly taken into account by theoretically determining their arrangement with a supercell model. The most stable structures and the formation energies of oxygen interstitials were determined for each charged state. It was found that the double-negatively charged state is stable over a wide range of the Fermi level, and that the excess oxygen ions form split interstitials with the original oxygen ions, while the neutral and the single-negatively charged states preferably form molecular oxygen. These species were found near the lanthanum vacancy site. The theoretically determined migration pathway along the c-axis essentially follows an interstitialcy mechanism. The obtained migration barrier is sensitive to the charge state, and is also affected by the lanthanum vacancy. The barrier height of the double-negatively charged state was calculated to be 0.58 eV for the model structure, which is consistent with the measured activation energy.

Limiting the Extent of the RDN1 Heterochromatin Domain by a Silencing Barrier and Sir2 Protein Levels in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, transcriptional silencing occurs at the cryptic mating-type loci (HML and HMR), telomeres, and ribosomal DNA ( rDNA; RDN1). Silencing in the rDNA is unusual in that polymerase II (Pol II) promoters within RDN1 are repressed by Sir2 but not Sir3 or Sir4. rDNA silencing unidirectionally spreads leftward, but the mechanism of limiting its spreading is unclear. We searched for silencing barriers flanking the left end of RDN1 by using an established assay for detecting barriers to HMR silencing. Unexpectedly, the unique sequence immediately adjacent to RDN1, which overlaps a prominent cohesin binding site (CARL2), did not have appreciable barrier activity. Instead, a fragment located 2.4 kb to the left, containing a tRNA(Gln) gene and the Ty1 long terminal repeat, had robust barrier activity. The barrier activity was dependent on Pol III transcription of tRNA(Gln), the cohesin protein Smc1, and the SAS1 and Gcn5 histone acetyltransferases. The location of the barrier correlates with the detectable limit of rDNA silencing when SIR2 is overexpressed, where it blocks the spreading of rDNA heterochromatin. We propose a model in which normal Sir2 activity results in termination of silencing near the physical rDNA boundary, while tRNA(Gln) blocks silencing from spreading too far when nucleolar Sir2 pools become elevated.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Dimethyl sulfide and other biogenic volatile organic compound emissions from branching coral and reef seawater: Potential sources of secondary aerosol over the Great Barrier Reef

Volatile organic compounds (VOCs) in the headspace of bubble chambers containing branches of live coral in filtered reef seawater were analysed using gas chromatography with mass spectrometry (GC-MS). When the coral released mucus it was a source of dimethyl sulfide (DMS) and isoprene; however, these VOCs were not emitted to the chamber headspace from mucus-free coral. This finding, which suggests that coral is an intermittent source of DMS and isoprene, was supported by the observation of occasional large pulses of atmospheric DMS (DMSa) over Heron Island reef on the southern Great Barrier Reef (GBR), Australia, in the austral winter. The highest DMSa pulse (320 ppt) was three orders of magnitude less than the DMS mixing ratio (460 ppb) measured in the headspace of a dynamically purged bubble chamber containing a mucus-coated branch of Acropora aspera indicating that coral reefs can be strong point sources of DMSa. Static headspace GC-MS analysis of coral fragments identified mainly DMS and seven other minor reduced sulfur compounds including dimethyl disulfide, methyl mercaptan, and carbon disulfide, while coral reef seawater was an indicated source of methylene chloride, acetone, and methyl ethyl ketone. The VOCs emitted by coral and reef seawater are capable of producing new atmospheric particles < 15 nm diameter as observed at Heron Island reef. DMS and isoprene are known to play a role in low-level cloud formation, so aerosol precursors such as these could influence regional climate through a sea surface temperature regulation mechanism hypothesized to operate over the GBR.