Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.

Understanding Molecular Orbitals; Sigma Orbitals

Various plots of sigma molecular orbitals in diatomic molecules are discussed.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

Molecular Genetic Basis of Opposite Gene by Environment Interactions in Reading Disability and Attention Deficit/Hyperactivity Disorder

The goal of this study is to better understand the genetic basis of Reading Disability (RD) and Attention Deficit Hyperactivity Disorder (ADHD) by examining molecular G x E interactions with parental education for each disorder. Research indicates that despite sharing genetic risk factors, RD and ADHD are influenced by different types of G x E interactions with parental education - a diathesis stress interaction in the case of ADHD and a bioecological interaction in RD. In order to resolve this apparent paradox, we conducted a preliminary study using behavioral genetic methods to test for G x E interactions in RD and the inattentive subtype of ADHD (ADHD-I) in the same sample of monozygotic and dizygotic Colorado Learning Disabilities Research Center same-sex twin pairs (DeFries et al., 1997), and our findings were consistent with the literature. We posited a genetic hypothesis for this opposite pattern of interactions, which suggests that only genes specific to each disorder enter into these opposite interactions, not the shared genes underlying their comorbidity. This study sought to further investigate this paradox using molecular genetics methods. We examined multiple candidate genes identified for RD or related language phenotypes and those identified for ADHD for G x E interactions with parental education. The specific aims of this study were as follows: 1) partition known risk alleles for RD and/or related language phenotypes and ADHD-I into those which are pleiotropic and non-pleiotropic by testing each risk allele for association with both RD and ADHD-I, 2) explore the main effects of parental education on both RD and ADHD-I, 3) address G-E correlations, and 4) conduct exploratory G x E interaction analyses in order to test the genetic hypothesis. Analyses suggested a number of pleiotropic genes that influence both RD and ADHD; however, results did not remain after correcting for multiple comparisons. Although exploratory G x E interaction findings were not significant after multiple comparison correction, results suggested a G x E interaction in the bioecological direction with KIAA0319, parental education, and ADHD-I. Given the limited power in the current study, replication of these findings with larger samples is necessary.

Dinâmica molecular de sistemas iônicos poliatômicos: modelos polarizável e não-polarizável

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito.

A novel and rational approach towards designing aluminium chelators

Aluminium (Al) is known to be neurotoxic and has been associated with the aetiology of Alzheimer's Disease. To date, only desferrioxamine (DFO), a trihydroxamic acid siderophore has been used in the clinical environment for the removal of Al from the body. However, this drug is expensive, orally inactive and is associated with many side effects. These studies employed a theoretical approach, with the use of quantum mechanics (QM) via semi-empirical molecular orbital (MO) calculations, and a practical approach using U87-MG glioblastoma cells as a model for evaluating the influence of potential chelators on the passage of aluminium into cells. Preliminary studies involving the Cambridge Structural Database (CSD) identified that Al prefers binding to bidentate ligands in a 3:1 manner, whereby oxygen was the exclusive donating atom. Statistically significant differences in M-O bond lengths when compared to other trivalent metal ions such as Fe3+ were established and used as an acceptance criterion for subsequent MO calculations. Of the semi-empirical methods parameterised for Al, the PM3 Hamiltonian was found to give the most reliable final optimised geometries of simple 3:1 Al complexes. Consequently the PM3 Hamiltonian was used for evaluating the Hf of 3:1 complexes with more complicated ligands. No correlation exists between published stability constants and individual parameters calculated via PM3 optimisations, although investigation of the dicarboxylates reveals a correlation of 0.961 showing promise for affinity prediction of closely related ligands. A simple and inexpensive morin spectrofluorescence assay has been developed and optimised producing results comparable to atomic absorption spectroscopy methods for the quantitative analysis of Al. This assay was used in subsequent in vitro models, initially on E. coli, which indicated that Al inhibits the antimicrobial action of ciprofloxacin, a potent quinolone antibiotic. Ensuing studies using the second model, U87-MG cells, investigated the influence of chelators on the transmembrane transport of Al, identifying 1,2-diethylhydroxypyridin-4-one as a ligand showing greatest potential for chelating Al in the clinical situation. In conclusion, these studies have explored semi-empirical MO Hamiltonians and an in-vitro U87-MG cell line, both as possible methods for predicting effective chelators of Al.

Encapsulated membrane proteins:a simplified system for molecular simulation

Over the past 50 years there has been considerable progress in our understanding of biomolecular interactions at an atomic level. This in turn has allowed molecular simulation methods employing full atomistic modeling at ever larger scales to develop. However, some challenging areas still remain where there is either a lack of atomic resolution structures or where the simulation system is inherently complex. An area where both challenges are present is that of membranes containing membrane proteins. In this review we analyse a new practical approach to membrane protein study that offers a potential new route to high resolution structures and the possibility to simplify simulations. These new approaches collectively recognise that preservation of the interaction between the membrane protein and the lipid bilayer is often essential to maintain structure and function. The new methods preserve these interactions by producing nano-scale disc shaped particles that include bilayer and the chosen protein. Currently two approaches lead in this area: the MSP system that relies on peptides to stabilise the discs, and SMALPs where an amphipathic styrene maleic acid copolymer is used. Both methods greatly enable protein production and hence have the potential to accelerate atomic resolution structure determination as well as providing a simplified format for simulations of membrane protein dynamics.

Development, validation and uncertainty analysis of quantitative structure and activity relationship models for Log P of disinfection by -products

Hydrophobicity as measured by Log P is an important molecular property related to toxicity and carcinogenicity. With increasing public health concerns for the effects of Disinfection By-Products (DBPs), there are considerable benefits in developing Quantitative Structure and Activity Relationship (QSAR) models capable of accurately predicting Log P. In this research, Log P values of 173 DBP compounds in 6 functional classes were used to develop QSAR models, by applying 3 molecular descriptors, namely, Energy of the Lowest Unoccupied Molecular Orbital (ELUMO), Number of Chlorine (NCl) and Number of Carbon (NC) by Multiple Linear Regression (MLR) analysis. The QSAR models developed were validated based on the Organization for Economic Co-operation and Development (OECD) principles. The model Applicability Domain (AD) and mechanistic interpretation were explored. Considering the very complex nature of DBPs, the established QSAR models performed very well with respect to goodness-of-fit, robustness and predictability. The predicted values of Log P of DBPs by the QSAR models were found to be significant with a correlation coefficient R2 from 81% to 98%. The Leverage Approach by Williams Plot was applied to detect and remove outliers, consequently increasing R 2 by approximately 2% to 13% for different DBP classes. The developed QSAR models were statistically validated for their predictive power by the Leave-One-Out (LOO) and Leave-Many-Out (LMO) cross validation methods. Finally, Monte Carlo simulation was used to assess the variations and inherent uncertainties in the QSAR models of Log P and determine the most influential parameters in connection with Log P prediction. The developed QSAR models in this dissertation will have a broad applicability domain because the research data set covered six out of eight common DBP classes, including halogenated alkane, halogenated alkene, halogenated aromatic, halogenated aldehyde, halogenated ketone, and halogenated carboxylic acid, which have been brought to the attention of regulatory agencies in recent years. Furthermore, the QSAR models are suitable to be used for prediction of similar DBP compounds within the same applicability domain. The selection and integration of various methodologies developed in this research may also benefit future research in similar fields.

Enzymatic differences between the endophyte Guignardia mangiferae (Botryosphaeriaceae) and the citrus pathogen G. citricarpa

The endophyte Guignardia mangiferae is closely related to G. citricarpa, the causal agent of citrus black spot; for many years these species had been confused with each other. The development of molecular analytical methods has allowed differentiation of the pathogen G. citricarpa from the endophyte G. mangiferae, but the physiological traits associated with pathogenicity were not described. We examined genetic and enzymatic characteristics of Guignardia spp strains; G. citricarpa produces significantly greater amounts of amylases, endoglucanases and pectinases, compared to G. mangiferae, suggesting that these enzymes could be key in the development of citrus black spot. Principal component analysis revealed pectinase production as the main enzymatic characteristic that distinguishes these Guignardia species. We quantified the activities of pectin lyase, pectin methylesterase and endopolygalacturonase; G. citricarpa and G. mangiferae were found to have significantly different pectin lyase and endopolygalacturonase activities. The pathogen G. citricarpa is more effective in pectin degradation. We concluded that there are significant physiological differences between the species G. citricarpa and G. mangiferae that could be associated with differences in pathogenicity for citrus plants.

Disparity between Multilocus Enzyme Electrophoresis, Microsatellite Markers and Pulsed-Field Gel Electrophoresis in epidemiological tracking of Candida albicans

Various molecular systems are available for epidemiological, genetic, evolutionary, taxonomic and systematic studies of innumerable fungal infections, especially those caused by the opportunistic pathogen C. albicans. A total of 75 independent oral isolates were selected in order to compare Multilocus Enzyme Electrophoresis (MLEE), Electrophoretic Karyotyping (EK) and Microsatellite Markers (Simple Sequence Repeats - SSRs), in their abilities to differentiate and group C. albicans isolates (discriminatory power), and also, to evaluate the concordance and similarity of the groups of strains determined by cluster analysis for each fingerprinting method. Isoenzyme typing was performed using eleven enzyme systems: Adh, Sdh, M1p, Mdh, Idh, Gdh, G6pdh, Asd, Cat, Po, and Lap (data previously published). The EK method consisted of chromosomal DNA separation by pulsed-field gel electrophoresis using a CHEF system. The microsatellite markers were investigated by PCR using three polymorphic loci: EF3, CDC3, and HIS3. Dendrograms were generated by the SAHN method and UPGMA algorithm based on similarity matrices (S(SM)). The discriminatory power of the three methods was over 95%, however a paired analysis among them showed a parity of 19.7-22.4% in the identification of strains. Weak correlation was also observed among the genetic similarity matrices (S(SM)(MLEE) x S(SM)(EK) x S(SM)(SSRs)). Clustering analyses showed a mean of 9 +/- 12.4 isolates per cluster (3.8 +/- 8 isolates/taxon) for MLEE, 6.2 +/- 4.9 isolates per cluster (4 +/- 4.5 isolates/taxon) for SSRs, and 4.1 +/- 2.3 isolates per cluster (2.6 +/- 2.3 isolates/taxon) for EK. A total of 45 (13%), 39(11.2%), 5 (1.4%) and 3 (0.9%) clusters pairs from 347 showed similarity (Si) of 0.1-10%, 10.1-20%, 20.1-30% and 30.1-40%, respectively. Clinical and molecular epidemiological correlation involving the opportunistic pathogen C. albicans may be attributed dependently of each method of genotyping (i.e., MLEE, EK, and SSRs) supplemented with similarity and grouping analysis. Therefore, the use of genotyping systems that give results which offer minimum disparity, or the combination of the results of these systems, can provide greater security and consistency in the determination of strains and their genetic relationships. (C) 2010 Elsevier B.V. All rights reserved.

A rapid single-step multiplex method for discriminating between Trichogramma (Hymenoptera: Trichogrammatidae) species in Australia

Inaccurate species identification confounds insect ecological studies. Examining aspects of Trichogramma ecology pertinent to the novel insect resistance management strategy for future transgenic cotton, Gossypium hirsutum L., production in the Ord River Irrigation Area (ORIA) of Western Australia required accurate differentiation between morphologically similar Trichogramma species. Established molecular diagnostic methods for Trichogramma identification use species-specific sequence difference in the internal transcribed spacer (ITS)-2 chromosomal region; yet, difficulties arise discerning polymerase chain reaction (PCR) fragments of similar base pair length by gel electrophoresis. This necessitates the restriction enzyme digestion of PCR-amplified ITS-2 fragments to readily differentiate Trichogramma australicum Girault and Trichogramma pretiosum Riley. To overcome the time and expense associated with a two-step diagnostic procedure, we developed a “one-step” multiplex PCR technique using species-specific primers designed to the ITS-2 region. This approach allowed for a high-throughput analysis of samples as part of ongoing ecological studies examining Trichogramma biological control potential in the ORIA where these two species occur in sympatry.

A first-principles density-functional calculation of the electronic and vibrational structure of the key melanin monomers

We report first-principles density-functional calculations for hydroquinone (HQ), indolequinone (IQ), and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of biomacromolecules with important biological functions (including photoprotection) and with the potential for certain bioengineering applications. We have used the difference of self-consistent fields method to study the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, HL. We show that HL is similar in IQ and SQ, but approximately twice as large in HQ. This may have important implications for our understanding of the observed broadband optical absorption of the eumelanins. The possibility of using this difference in HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to nondestructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behavior of the eumelanins.

Anti-inflammatory activity and possible mechanism of extract from Mikania laevigata in carrageenan-induced peritonitis

Objectives The aim was to test the potential use of an extract of Mikania laevigata (popularly known in Brazil as guaco), made from leaves harvested in different months of the year, oil neutrophil migration after all inflammatory Stimulus and investigate the underlying molecular mechanisms. Methods We examined the effect of guaco on vascular permeability and leucocyte function in carrageenan-induced peritonitis in mice. Key findings Our results demonstrated that guaco extract administered subcutaneously (3 mg/kg) decreased the vascular permeability and also leucocyte rolling and adhesion to the inflamed tissues by a mechanism dependent on nitric oxide. Specifically, inhibitors of nitric oxide synthase remarkably abrogated the guaco extract-mediated suppression of neutrophil migration to the inflammatory site. In addition, guaco extract-mediated suppression of neutrophil migration appeared to be dependent on the production of the cytokines interleukin-1 beta and tumour necrosis factor-alpha. One of the major constituents of the guaco extract, coumarin, was able to inhibit the neutrophil migration towards the inflammatory focus. Conclusions In conclusion the anti-inflammatory effect induced by guaco extract may be by inhibition of pro-inflammatory cytokine production at the inflammatory site.

Analysis of marine bivalve shellfish from the fish market in Santos city, Sao Paulo state, Brazil, for Toxoplasma gondii

The aim of this study was to determine if Toxoplasma gondii are present in oysters (Crassostrea rhizophorae) and mussels (Mytella guyanensis) under natural conditions using a bioassay in mice and molecular detection methods. We first compared two standard protocols for DNA extraction, phenol-chloroform (PC) and guanidine-thiocyanate (GT), for both molluscs. A total of 300 oysters and 300 mussels were then acquired from the fish market in Santos city, Sao Paulo state, Brazil, between March and August of 2008 and divided into 60 groups of 5 oysters and 20 groups of 15 mussels. To isolate the parasite, five mice were orally inoculated with sieved tissue homogenates from each group of oysters or mussels. For molecular detection of T. gondii, DNA from mussels was extracted using the PC method and DNA from oysters was extracted using the GT method. A nested-PCR (Polymerase Chain Reaction) based on the amplification of a 155 bp fragment from the B1 gene of T. gondii was then performed. Eleven PCR-RFLP (Restriction Fragment Length Polymorphism) markers, SAG1, SAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1, CS3 and Apico, were used to genotype positive samples. There was no isolation of the parasite by bioassay in mice. T. gondii was not detected in any of the groups of mussels by nested-PCR. DNA of T. gondii was apparently detected by nested-PCR in 2 groups of oysters (3.3%). Genotyping of these two positive samples was not successful. The results suggest that oysters of the species C. rhizophorae, the most common species from the coast of Sao Paulo, can filter and retain T. gondii oocysts from the marine environment. Ingestion of raw oysters as a potential transmission source of T. gondii to humans and marine mammals should be further investigated. (C) 2010 Elsevier B.V. All rights reserved.