Parity violation in molecular magnetic resonance properties

The Standard Model of particle physics consists of the quantum electrodynamics (QED) and the weak and strong nuclear interactions. The QED is the basis for molecular properties, and thus it defines much of the world we see. The weak nuclear interaction is responsible for decays of nuclei, among other things, and in principle, it should also effects at the molecular scale. The strong nuclear interaction is hidden in interactions inside nuclei. From the high-energy and atomic experiments it is known that the weak interaction does not conserve parity. Consequently, the weak interaction and specifically the exchange of the Z^0 boson between a nucleon and an electron induces small energy shifts of different sign for mirror image molecules. This in turn will make the other enantiomer of a molecule energetically favorable than the other and also shifts the spectral lines of the mirror image pair of molecules into different directions creating a split. Parity violation (PV) in molecules, however, has not been observed. The topic of this thesis is how the weak interaction affects certain molecular magnetic properties, namely certain parameters of nuclear magnetic resonance (NMR) and electron spin resonance (ESR) spectroscopies. The thesis consists of numerical estimates of NMR and ESR spectral parameters and investigations of the effects of different aspects of quantum chemical computations to them. PV contributions to the NMR shielding and spin-spin coupling constants are investigated from the computational point of view. All the aspects of quantum chemical electronic structure computations are found to be very important, which makes accurate computations challenging. Effects of molecular geometry are also investigated using a model system of polysilyene chains. PV contribution to the NMR shielding constant is found to saturate after the chain reaches a certain length, but the effects of local geometry can be large. Rigorous vibrational averaging is also performed for a relatively small and rigid molecule. Vibrational corrections to the PV contribution are found to be only a couple of per cents. PV contributions to the ESR g-tensor are also evaluated using a series of molecules. Unfortunately, all the estimates are below the experimental limits, but PV in some of the heavier molecules comes close to the present day experimental resolution.

Johnalcornia gen. et. comb. nov., and nine new combinations in Curvularia based on molecular phylogenetic analysis

An examination of ex-type and authentic cultures of 34 species of Bipolaris and Curvularia by phylogenetic analysis of four loci (EF-1α, GAPDH, ITS and LSU) resulted in nine new combinations in Curvularia, as well as new synonymies for some species of Bipolaris and Curvularia. Lectotypes are designated for Bipolaris secalis and Curvularia richardiae, and an epitype is designated for Curvularia crustacea. A new monotypic genus, Johnalcornia, is introduced to accommodate Bipolaris aberrans, which clusters sister to the newly described Porocercospora. Johnalcornia differs morphologically from this taxon by producing distinctive conidia-like chlamydospores as well as comparatively thick-walled, geniculate conidiophores, with conidiogenous cells that have conspicuous scars. Johnalcornia further differs from related genera by forming the second conidial septum in the apical cell.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.

Molecular and crystal structure of deoxyguanosine 5'-phosphate

RECENT crystallographic studies of the dinucleosides ApU (ref. 1) and GpC (ref. 2) have given experimental proof for the base pairing arrangement proposed by Watson and Crick for the DNA double helix3. Another striking feature of this structure relates to the torsional angle about the C5'-C4' bond in the phosphate−sugar backbone chain. In the Crick and Watson model4, this conformation is gauche−trans (GT). Crystal structures of 5'-nucleotides, dinucleosides and dinucleotides so far studied, however, have shown only the gauche−gauche (GG) conformation about this bond. The GG conformer is also the only one found in the refined models of the proposed structure of the double helical nucleic acids and polynucleotides5−7. The only nucleotide with a GT conformation is 6-azauridine-5'-phosphate8 which is not a normal monomer unit of nucleic acids. It is also reported that 5'-dGMP assumes preferentially GT conformation in solution9.

Mechanism of adsorption effects on electrode kinetics: Part I: Selection rules in the case of low coverage by molecular surfactants

It is shown that the effect of adsorption of inert molecules on electrode reaction rates is completely accounted for, by introducing into the rate equation, adsorption-induced changes in both the effective electrode area as well as in the electrostatic potential at the reaction site with an additional term for the noncoulombic interaction between the reactant and the adsorbate. The electrostatic potential at the reaction site due to the adsorbed layer is calculated using a model of discretely-distributed molecules in parallel orientation when adsorbed on the electrode with an allowance for thermal agitation. The resulting expression, which is valid for the limiting case of low coverages, is used to predict the types of molecular surfactants that are most likely to be useful for acceleration and inhibition of electrode reactions.

Conformation And C-N Rotational Barrier Of Dithiocarbamate Esters Nmr And Molecular-Orbital Studies

The variable temperature proton and ambient temperature carbon-13 NMR spectra of S-methyl dithiocarbamate esters have been recorded. The results of the theoretical energy calculations (CNDO/2 and EHT types) together with the experimental data have been interpreted in terms of the molecular conformations. The barrier heights for the rotation about the thioamide C—N bond are calculated using the CNDO/2 method and the results are discussed in terms of the computed charge densities and bond orders.

Polymer supported metal complexes as catalysts for oxidation of thiosalts by molecular oxygen IV. Quaternised poly(4-vinyl pyridine) complexes with Cu2+ as template

Copper(II) complexes of quaternised poly(4-vinylpyridine) (PVP) of different degrees of quaternisation and copper content have been prepared by crosslinking the polymer with 1,2-dibromoethane in the presence of Cu2+ ion as template. The stability constant of the PVP---Cu(II) complexes is found to increase with the degree of crosslinking quaternisation of the resin, but the rate at which Cu2+ is adsorbed by the resin decreases. An optimum combination of both stability and rate can be achieved with a moderate degree (31%) of crosslinking. A kinetic study reveals that quaternisation increases significantly the catalytic activity of the complex for the oxidation of S2O2−3 by O2 compared with PVP----Cu(II) without quaternisation, but it deactivates the complex for the oxidation of both S3O2−6 and S4O2−6. The batch reactor oxidation kinetics at pH 2.16, where the rate is observed to be maximum, is well explained by the Langmuir—Hinshelwood model assuming the coordination of both O2 and thioanion to Cu(II) as a precursor to the oxidation reaction.

Molecular diversity of Chickpea chlorotic dwarf virus in Sudan: High rates of intra-species recombination – a driving force in the emergence of new strains

In Sudan Chickpea chlorotic dwarf virus (CpCDV, genus Mastrevirus, family Geminiviridae) is an important pathogen of pulses that are grown both for local consumption, and for export. Although a few studies have characterised CpCDV genomes from countries in the Middle East, Africa and the Indian subcontinent, little is known about CpCDV diversity in any of the major chickpea production areas in these regions. Here we analyse the diversity of 146 CpCDV isolates characterised from pulses collected across the chickpea growing regions of Sudan. Although we find that seven of the twelve known CpCDV strains are present within the country, strain CpCDV-H alone accounted for ∼73% of the infections analysed. Additionally we identified four new strains (CpCDV-M, -N, -O and -P) and show that recombination has played a significant role in the diversification of CpCDV, at least in this region. Accounting for observed recombination events, we use the large amounts of data generated here to compare patterns of natural selection within protein coding regions of CpCDV and other dicot-infecting mastrevirus species.

Proton NMR and infrared investigations of secondary dithiocarbamate ester: molecular conformation and vibrational assignment of S-methyl N-methyl dithiocarbamate

Rotational isomerism of S-methyl N-methyl dithiocarbamate (MMDTC) has been investigated by means of variable temperature proton NMR and i.r. spectroscopy. The i.r. spectra of MMDTC as neat, solution and at sub-ambient temperatures have been examined. Normal vibrational analysis of all the fundamentals of MMDTC has been carried out, the vibrational assignment has been compared with those of related secondary thioamides to note the consistency in the assignments and to obtain the pattern characteristic of the secondary thioamide vibrations.

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.