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.

Diagnostic Molecular Markers for Phosphine Resistance in U.S. Populations of Tribolium castaneum and Rhyzopertha dominica

Stored product beetles that are resistant to the fumigant pesticide phosphine (hydrogen phosphide) gas have been reported for more than 40 years in many places worldwide. Traditionally, determination of phosphine resistance in stored product beetles is based on a discriminating dose bioassay that can take up to two weeks to evaluate. We developed a diagnostic cleaved amplified polymorphic sequence method, CAPS, to detect individuals with alleles for strong resistance to phosphine in populations of the red flour beetle, Tribolium castaneum, and the lesser grain borer, Rhyzopertha dominica, according to a single nucleotide mutation in the dihydrolipoamide dehydrogenase (DLD) gene. We initially isolated and sequenced the DLD genes from susceptible and strongly resistant populations of both species. The corresponding amino acid sequences were then deduced. A single amino acid mutation in DLD in populations of T.castaneum and R.dominica with strong resistance was identified as P45S in T.castaneum and P49S in R.dominica, both collected from northern Oklahoma, USA. PCR products containing these mutations were digested by the restriction enzymes MboI and BstNI, which revealed presence or absence, respectively of the resistant (R) allele and allowed inference of genotypes with that allele. Seven populations of T.castaneum from Kansas were subjected to discriminating dose bioassays for the weak and strong resistance phenotypes. Application of CAPS to these seven populations confirmed the R allele was in high frequency in the strongly resistant populations, and was absent or at a lower frequency in populations with weak resistance, which suggests that these populations with a low frequency of the R allele have the potential for selection of the strong resistance phenotype. CAPS markers for strong phosphine resistance will help to detect and confirm resistant beetles and can facilitate resistance management actions against a given pest population.

Mathematical models on the impact of noise and dyadic molecular structures on the properties of a cardiac myocyte

In cardiac myocytes (heart muscle cells), coupling of electric signal known as the action potential to contraction of the heart depends crucially on calcium-induced calcium release (CICR) in a microdomain known as the dyad. During CICR, the peak number of free calcium ions (Ca) present in the dyad is small, typically estimated to be within range 1-100. Since the free Ca ions mediate CICR, noise in Ca signaling due to the small number of free calcium ions influences Excitation-Contraction (EC) coupling gain. Noise in Ca signaling is only one noise type influencing cardiac myocytes, e.g., ion channels playing a central role in action potential propagation are stochastic machines, each of which gates more or less randomly, which produces gating noise present in membrane currents. How various noise sources influence macroscopic properties of a myocyte, how noise is attenuated and taken advantage of are largely open questions. In this thesis, the impact of noise on CICR, EC coupling and, more generally, macroscopic properties of a cardiac myocyte is investigated at multiple levels of detail using mathematical models. Complementarily to the investigation of the impact of noise on CICR, computationally-efficient yet spatially-detailed models of CICR are developed. The results of this thesis show that (1) gating noise due to the high-activity mode of L-type calcium channels playing a major role in CICR may induce early after-depolarizations associated with polymorphic tachycardia, which is a frequent precursor to sudden cardiac death in heart failure patients; (2) an increased level of voltage noise typically increases action potential duration and it skews distribution of action potential durations toward long durations in cardiac myocytes; and that (3) while a small number of Ca ions mediate CICR, Excitation-Contraction coupling is robust against this noise source, partly due to the shape of ryanodine receptor protein structures present in the cardiac dyad.

Molecular evidence of Ebola Reston virus infection in Philippine bats

Background: In 2008-09, evidence of Reston ebolavirus (RESTV) infection was found in domestic pigs and pig workers in the Philippines. With species of bats having been shown to be the cryptic reservoir of filoviruses elsewhere, the Philippine government, in conjunction with the Food and Agriculture Organization of the United Nations, assembled a multi-disciplinary and multi-institutional team to investigate Philippine bats as the possible reservoir of RESTV. Methods: The team undertook surveillance of bat populations at multiple locations during 2010 using both serology and molecular assays. Results: A total of 464 bats from 21 species were sampled. We found both molecular and serologic evidence of RESTV infection in multiple bat species. RNA was detected with quantitative PCR (qPCR) in oropharyngeal swabs taken from Miniopterus schreibersii, with three samples yielding a product on conventional hemi-nested PCR whose sequences differed from a Philippine pig isolate by a single nucleotide. Uncorroborated qPCR detections may indicate RESTV nucleic acid in several additional bat species (M. australis, C. brachyotis and Ch. plicata). We also detected anti-RESTV antibodies in three bats (Acerodon jubatus) using both Western blot and ELISA. Conclusions: The findings suggest that ebolavirus infection is taxonomically widespread in Philippine bats, but the evident low prevalence and low viral load warrants expanded surveillance to elaborate the findings, and more broadly, to determine the taxonomic and geographic occurrence of ebolaviruses in bats in the region. © 2015 Jayme et al.

Neutron Scattering and Molecular Dynamics Evidence for Levitation Effect in Nanopores

Neutron Scattering and Molecular Dynamics Evidence for Levitation Effect in Nanopores ... Neutron scattering measurements and molecular dynamics simulations have been carried out on the three isomers of pentane (neopentane (neo), isopentane (iso), and n-pentane (n-)) adsorbed in zeolite NaY. ... In order to understand this surprising dependence, the dimensionless levitation parameter, γ, for atomic systems may be modified to suit molecular systems.

Molecular vibrations and conformation of primary dithiocarbamate ester. Infrared and NMR studies on S-methyl dithiocarbamate

The i.r. spectra of a primary dithiocarbamate ester namely, S-methyl dithiocarbamate (SMDTC) and its N-dideuterated compound have been measured between 4000 and 30 cm−1. Spectra in solution and at liquid nitrogen temperature have also been obtained. Assignment of all the fundamentals has been proposed and supported from a full normal coordinate analysis. The band assignments for SMDTC have been compared with those of related molecules and the characteristic bands of primary thioamides are derived. Conformational flexibility of SMDTC has been examined by i.r. and proton NMR spectroscopy. The hindered rotation around the C---N bond has been studied by a complete line shape analysis. The magnitude of ---NH2 and ---CH3 torsional barriers is also estimated from vibrational frequencies and force constants.

Molecular structures of cytidine-5-diphosphate and cytidine-5-diphospho-choline, and their role in intermediary metabolism

The nucleotide coenzyme cytidine-5-diphospho-choline is highly folded. The CMP-5 parts of the molecules in the crystal structure are strongly linked by metal ligation and hydrogen bonds leaving the phosphoryl-choline residues relatively free. Cytidine-5-diphosphoric acid exists as a zwitterion with N31 protonated. The P−O bond lengths from the anhydride bridging oxygen in the pyrophosphate are significantly different.

Molecular structure of nicotinamide adenine dinucleotide

NICOTINAMIDE adenine dinucleotide (NAD) has a fundamental role in metabolic processes as an electron transport molecule. Although its chemical structure was elucidated1 in 1934, its detailed conformation remains still to be established in spite of numerous physicochemical applications2. NAD analogues with a variety of substitutions on the bases are known to retain considerable activity of the natural coenzyme as long as the pyrophosphate diester group has been retained3,4. The geometry of this backbone moiety is therefore indispensable to our understanding of the conformation and function of the coenzyme. We have so far no experimental evidence on this in NAD or any other nucleotide coenzyme molecule. X-ray studies have been possible only on those analogues5,6 where the nicotinamide and adenine rings are linked by a trimethylene bridge. The results are conflicting and it is difficult to use them to provide a structural basis for the NAD molecule itself, particularly as the phosphate backbone is absent from these analogues.

An uncommon nucleotide conformation shown by molecular structure of deoxyuridine-5-phosphate and nucleic acid stereochemistry

CRYSTAL structure determinations of nucleic acid fragments have shown that several of the conformational features found in the monomeric building blocks are also manifested at the nucleic acid level. Stereochemical variations between thymine and uracil nucleotides are therefore of interest as they can provide a structural basis for some of the differences between the conformations of DNA and RNA. X-ray studies have so far not shown any major dissimilarities between these two nucleotide species although the sugar ring of deoxyribonucleotides is found to possess greater flexibility than that in ribonucleotides. We report here the molecular structure of deoxyuridine-5'-phosphate (dUMP-5') which is not a common monomer unit of DNAs as it is replaced by its 5-methyl analogue deoxythymidine-5'-phosphate (dTMP-5'). The investigation was undertaken to help determine whether or not this implied a fundamental difference between the geometries of these two molecules.

A molecular survey of Eimeria in chickens across Australia

Coccidiosis is a costly enteric disease of chickens caused by protozoan parasites of the genus Eimeria. Disease diagnosis and management is complicated since there are multiple Eimeria species infecting chickens and mixed species infections are common. Current control measures are only partially effective and this, combined with concerns over vaccine efficacy and increasing drug resistance, demonstrates a need for improved coccidiosis diagnosis and control. Before improvements can be made, it is important to understand the species commonly infecting poultry flocks in both backyard and commercial enterprises. The aim of this project was to conduct a survey and assessment of poultry Eimeria across Australia using genetic markers, and create a collection of isolates for each Eimeria species. A total of 260 samples (faecal or caecal) was obtained, and survey results showed that Eimeria taxa were present in 98% of commercial and 81% of backyard flocks. The distribution of each Eimeria species was widespread across Australia, with representatives of all species being found in every state and territory, and the Eimeria species predominating in commercial flocks differed from those in backyard flocks. Three operational taxonomic units also occurred frequently in commercial flocks highlighting the need to understand the impact of these uncharacterised species on poultry production. As Eimeria infections were also frequent in backyard flocks, there is a potential for backyard flocks to act as reservoirs for disease, especially as the industry moves towards free range production systems. This Eimeria collection will be an important genetic resource which is the crucial first step in the development of more sophisticated diagnostic tools and the development of new live vaccines which ultimately will provide savings to the industry in terms of more efficient coccidiosis management.

Molecular Breeding for Complex Adaptive Traits: How Integrating Crop Ecophysiology and Modelling Can Enhance Efficiency

Progress in crop improvement is limited by the ability to identify favourable combinations of genotypes (G) and management practices (M) in relevant target environments (E) given the resources available to search among the myriad of possible combinations. To underpin yield advance we require prediction of phenotype based on genotype. In plant breeding, traditional phenotypic selection methods have involved measuring phenotypic performance of large segregating populations in multi-environment trials and applying rigorous statistical procedures based on quantitative genetic theory to identify superior individuals. Recent developments in the ability to inexpensively and densely map/sequence genomes have facilitated a shift from the level of the individual (genotype) to the level of the genomic region. Molecular breeding strategies using genome wide prediction and genomic selection approaches have developed rapidly. However, their applicability to complex traits remains constrained by gene-gene and gene-environment interactions, which restrict the predictive power of associations of genomic regions with phenotypic responses. Here it is argued that crop ecophysiology and functional whole plant modelling can provide an effective link between molecular and organism scales and enhance molecular breeding by adding value to genetic prediction approaches. A physiological framework that facilitates dissection and modelling of complex traits can inform phenotyping methods for marker/gene detection and underpin prediction of likely phenotypic consequences of trait and genetic variation in target environments. This approach holds considerable promise for more effectively linking genotype to phenotype for complex adaptive traits. Specific examples focused on drought adaptation are presented to highlight the concepts.

Molecular orbital studies on some 4-substituted pyridine N-oxides and 4-nitroquinoline N-oxide

The electronic structures of a series of 4-substituted pyridine N-oxides and 4-nitroquinoline N-oxide are investigated using the simple Pariser-Parr-Pople (PPP), a modified PPP, IEH and MINDO/2 methods. The electronic absorption band maxima and dipole moments are calculated and compared with experimental values. The photoelectron spectra of these compounds are assigned. The nature of the N-oxide group is characterized using the orbital population distributions. The antifungal activity exhibited by some of these compounds is discussed in terms of the nucleophilic frontier electron densities, superdelocalizabilities and electron acceptor properties. The effect of the electron releasing as well as the electron withdrawing substituents on the physico-chemical properties is explained.

Crystal and molecular structure of allo-4-hydroxy-L-proline dihydrate

allo-4-Hydroxy-L-proline crystallizes from an aqueous solution as the dihydrate. The crystals are orthorhombic, space group P212121, with a=7.08 (2), b=22.13 (3), c= 5"20 (2) A,. The structure was solved by direct methods and refined by block-diagonal least squares. The final R for 733 observed reflexions is 0.054. The molecule exists as a zwitterion with hydroxyl and carboxyl groups cis to the pyrrolidine ring. The latter is puckered at the fl-carbon atom, which deviates by -0.54 A, from the best plane formed by the four remaining atoms. The molecules are held together by a network of hydrogen bonds, the water molecules playing a dominant role in the stability of the structure.

Sequence-dependent molecular conformation of polynucleotides: right and left-handed helices

Based upon a stereochemical guideline, two topologically distinct types of helicalduplexes have been deduced for a polynucleotide duplex with alternating purine pyrimidine sequence (PAPP): (a) right-handed uniform (RU) helix and (b) left-handed zig-zag (LZ) helix. Both structures have trinucleoside diphosphate as the basic unit wherein the purine pyrimidine fragment has a different conformation from the pyrimidine-purine fragment. Thus, RU and LZ helices represent two different classes of sequence-dependent molecular conformations for PAPP. The conformationalf eatures of an RU helix of PAPP in B-form and three LZ-helices for B-, D- and Z-forms are discussed.

Crystal and molecular structure of 15-bromolongibornane-8,9-dione

The crystal and molecular structure of the title compound (1) has been determined by the heavy-atom method from 1038 observed three-dimensional photographic data. Crystals are orthorhombic, with a = 20.07 ± 0.02, b= 10.05 ± 0.02, c= 7.31 ± 0.01 Å, space group P212121, with Z= 4. The structure was refined by block diagonal leastsquares to R 0.099. The conformation of the norbornane moiety is discussed. The seven-membered ring portion of the molecule adopts an approximate chair conformation. The packing of the molecules in the crystal is mainly a consequence of van der Waals interactions.