Operator expansion: definition and molecular integral applications

Es defineix l'expansió general d'operadors com una combinació lineal de projectors i s'exposa la seva aplicació generalitzada al càlcul d'integrals moleculars. Com a exemple numèric, es fa l'aplicació al càlcul d'integrals de repulsió electrònica entre quatre funcions de tipus s centrades en punts diferents, i es mostren tant resultats del càlcul com la definició d'escalat respecte a un valor de referència, que facilitarà el procés d'optimització de l'expansió per uns paràmetres arbitraris. Es donen resultats ajustats al valor exacte

Els límits i les relacions entre els peixos acanthopterygii: filogènia molecular de mugilomorpha i atherinomorpha

Sovint, la sistemàtica, basada principalment en caràcters morfològics, no es correspon amb els processos evolutius relacionats amb l'aparició dels grups d'organismes. En l'actualitat, la utilització de les dades moleculars es fa indispensable per a una revisió i millora de la classificació biològica de diversos organismes, com els peixos Acanthopterygii. A la sèrie Mugilomorpha la incongruència entre la taxonomia i la filogènia sorgeix de l'elevada semblança morfològica trobada per part dels seus membres. Pel que fa referència a la sèrie Atherinomorpha, la problemàtica principal resideix en determinar la seva proximitat evolutiva respecte a la sèrie anterior i en establir les relacions filogenètiques dins de la mateixa. Per tant, s'hi ha volgut estimar tant la divergència genètica dins de cada sèrie com inferir les relacions filogenètiques entre ambdues mitjançant la seqüenciació directa del DNA de les regions mitocondrials corresponents al tRNA-Phe, 12S rRNA, COI, cytb, tRNA-Thr, tRNA-Pro i regió control.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Fast training of self organizing maps for the visual exploration of molecular compounds

Visual exploration of scientific data in life science area is a growing research field due to the large amount of available data. The Kohonen’s Self Organizing Map (SOM) is a widely used tool for visualization of multidimensional data. In this paper we present a fast learning algorithm for SOMs that uses a simulated annealing method to adapt the learning parameters. The algorithm has been adopted in a data analysis framework for the generation of similarity maps. Such maps provide an effective tool for the visual exploration of large and multi-dimensional input spaces. The approach has been applied to data generated during the High Throughput Screening of molecular compounds; the generated maps allow a visual exploration of molecules with similar topological properties. The experimental analysis on real world data from the National Cancer Institute shows the speed up of the proposed SOM training process in comparison to a traditional approach. The resulting visual landscape groups molecules with similar chemical properties in densely connected regions.

Notations and conventions in molecular spectroscopy: part 1. General spectroscopic notation

The field of Molecular Spectroscopy was surveyed in order to determine a set of conventions and symbols which are in common use in the spectroscopic literature. This document, which is Part I in a series, establishes the notations and conventions used for general spectroscopic notations and deals with quantum mechanics, quantum numbers (vibrational states, angular momentum and energy levels), spectroscopic transitions, and miscellaneous notations (e.g. spectroscopic terms). Further parts will follow, dealing inter alia with symmetry notation, permutation and permutation-inversion symmetry notation, vibration-rotation spectroscopy and electronic spectroscopy.

Molecular characterization of Iranian wheat stripe virus shows its taxonomic position as a distinct species in the genus Tenuivirus

The full lengths of three genome segments of Iranian wheat stripe virus (IWSV) were amplified by reverse transcription (RT) followed by polymerase chain reaction (PCR) using a primer complementary to tenuivirus conserved terminal sequences. The segments were sequenced and found to comprise 3469, 2337, and 1831 nt, respectively. The gene organization of these segments is similar to that of other known tenuiviruses, each displaying an ambisense coding strategy. IWSV segments, however, are different from those of other viruses with respect to the number of nucleotides and deduced amino acid sequence for each ORF. Depending on the segment, the first 16-22 nt at the 5' end and the first 16 nt at the 3' end are highly conserved among IWSV and rice hoja blanca virus (RHBV), rice stripe virus (RSV) and maize stripe virus ( MStV). In addition, the first 15-18 nt at the 5' end are complementary to the first 16-18 nt at the 3' end. Phylogenetic analyses showed close similarity and a common ancestor for IWSV, RHBV, and Echinochloa hoja blanca virus (EHBV). These findings confirm the position of IWSV as a distinct species in the genus Tenuivirus.

Molecular characterization of Iranian wheat stripe virus shows its taxonomic position as a distinct species in the genus Tenuivirus

The full lengths of three genome segments of Iranian wheat stripe virus (IWSV) were amplified by reverse transcription (RT) followed by polymerase chain reaction (PCR) using a primer complementary to tenuivirus conserved terminal sequences. The segments were sequenced and found to comprise 3469, 2337, and 1831 nt, respectively. The gene organization of these segments is similar to that of other known tenuiviruses, each displaying an ambisense coding strategy. IWSV segments, however, are different from those of other viruses with respect to the number of nucleotides and deduced amino acid sequence for each ORF. Depending on the segment, the first 16-22 nt at the 5' end and the first 16 nt at the 3' end are highly conserved among IWSV and rice hoja blanca virus (RHBV), rice stripe virus (RSV) and maize stripe virus ( MStV). In addition, the first 15-18 nt at the 5' end are complementary to the first 16-18 nt at the 3' end. Phylogenetic analyses showed close similarity and a common ancestor for IWSV, RHBV, and Echinochloa hoja blanca virus (EHBV). These findings confirm the position of IWSV as a distinct species in the genus Tenuivirus.

Molecular evidence from Ciona intestinalis for the evolutionary origin of vertebrate sensory placodes

Cranial sensory placodes are focused areas of the head ectoderm of vertebrates that contribute to the development of the cranial sense organs and their associated ganglia. Placodes have long been considered a key character of vertebrates, and their evolution is proposed to have been essential for the evolution of an active predatory lifestyle by early vertebrates. Despite their importance for understanding vertebrate origins, the evolutionary origin of placodes has remained obscure. Here, we use a panel of molecular markers from the Six, Eya, Pax, Dach, FoxI, COE and POUIV gene families to examine the tunicate Ciona intestinalis for evidence of structures homologous to vertebrate placodes. Our results identify two domains of Ciona ectoderm that are marked by the genetic cascade that regulates vertebrate placode formation. The first is just anterior to the brain, and we suggest this territory is equivalent to the olfactoty/adenohypophyseal placodes of vertebrates. The second is a bilateral domain adjacent to the posterior brain and includes cells fated to form the atrium and atrial siphon of adult Ciona. We show this bares most similarity to placodes fated to form the vertebrate acoustico-lateralis system. We interpret these data as support for the hypothesis that sensory placodes did not arise de novo in vertebrates, but evolved froth pre-existing specialised areas of ectoderm that contributed to sensory organs in the common ancestor of vertebrate and tunicates. Published by Elsevier Inc.

Molecular structures of benzoic acid and 2-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of benzoic acid (C6H5COOH) and 2-hydroxybenzoic acid (C6H4OHCOOH) have been determined in the gas phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d, p)) predict two conformers for benzoic acid, one which is 25.0 kJ mol(-1) (MP2) lower in energy than the other. In the low-energy form, the carboxyl group is coplanar with the phenyl ring and the O-H group eclipses the C=O bond. Theoretical calculations (HF and MP2/6-311+ G(d, p)) carried out for 2-hydroxybenzoic acid gave evidence for seven stable conformers but one low-energy form (11.7 kJ mol-1 lower in energy (MP2)) which again has the carboxyl group coplanar with the phenyl ring, the O-H of the carboxyl group eclipsing the C=O bond and the C=O of the carboxyl group oriented toward the O-H group of the phenyl ring. The effects of internal hydrogen bonding in 2-hydroxybenzoic acid can be clearly observed by comparison of pertinent structural parameters between the two compounds. These differences for 2-hydroxybenzoic acid include a shorter exocyclic C-C bond, a lengthening of the ring C-C bond between the substituents, and a shortening of the carboxylic single C-O bond.

Molecular structures of 3-hydroxybenzoic acid and 4-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of 3-hydroxybenzoic acid and 4-hydroxybenzoic acid have been determined by gas-phase electron diffraction using results from quantum chemical calculations to inform the choice of restraints applied to some of the structural parameters. The results from the study presented here demonstrate that resonance hybrids are not as helpful in rationalizing the structures of 2-, 3-, and 4-hydroxybenzoic acids as are models based upon electrostatic effects.

Quantum mechanics studies of the tautomers of dlhydrodiacetylformoin, an important Maillard intermediate and the synthesis of furaneol

Conformational analyses have been carried out on the acyclic and cyclic forms of dihydrodiacetylformoin, an important Maillard intermediate and precursor for furaneol. For the acyclic forms, the 2,5-dicarbonyl isomers have the lowest energy, while for the cyclic forms, the 3-carbonyl are favoured over the 4-carbonyl isomers. The likely path for cyclisation is investigated and it is shown that the favoured path is dependent upon the relative chiralities of the carbon atoms and in particular that the reaction proceeds more readily if C2 and C3 have different chiralities. After cyclisation, the reaction path to produce furaneol proceeds via the loss of a water molecule. This reaction has been studied with a model including two water molecules and a hydroxide anion and shows relatively low-energy barriers. (C) 2008 Published by Elsevier B.V.

A molecular and structural characterization of senescing Arabidopsis siliques and comparison of transcriptional profiles with senescing petals and leaves

Senescence of plant organs is a genetically controlled process that regulates cell death to facilitate nutrient recovery and recycling, and frequently precedes, or is concomitant with, ripening of reproductive structures. In Arabidopsis thaliana, the seeds are contained within a silique, which is itself a photosynthetic organ in the early stages of development and undergoes a programme of senescence prior to dehiscence. A transcriptional analysis of the silique wall was undertaken to identify changes in gene expression during senescence and to correlate these events with ultrastructural changes. The study revealed that the most highly up-regulated genes in senescing silique wall tissues encoded seed storage proteins, and the significance of this finding is discussed. Global transcription profiles of senescing siliques were compared with those from senescing Arabidopsis leaf or petal tissues using microarray datasets and metabolic pathway analysis software (MapMan). In all three tissues, members of NAC and WRKY transcription factor families were up-regulated, but components of the shikimate and cell-wall biosynthetic pathways were down-regulated during senescence. Expression of genes encoding ethylene biosynthesis and action showed more similarity between senescing siliques and petals than between senescing siliques and leaves. Genes involved in autophagy were highly expressed in the late stages of death of all plant tissues studied, but not always during the preceding remobilization phase of senescence. Analyses showed that, during senescence, silique wall tissues exhibited more transcriptional features in common with petals than with leaves. The shared and distinct regulatory events associated with senescence in the three organs are evaluated and discussed.

Quantum mechanics studies of the tautomers of diacetylformoin, an important maillard product and odorant

Diacetylformoin (3,4-dihydroxy-3-hexene-2,5-dione) has 16 tautomers, many with several possible conformations and all have been geometry optimised using quantum mechanics at the HF/6-31+G* level. Eleven structures have been identified with energies within 10 kcal mol(-1) of the minimum energy structure. Of these eight are acyclic and three cyclic. Calculations of NMR spectra have clarified the identity of the acyclic and cyclic structures found experimentally. The mechanism for cyclisation has been investigated and transition states obtained. The lowest energy reaction path requires the loss and gain of a proton during cyclisation. (c) 2006 Elsevier B.V. All rights reserved.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.