Relations among several nuclear and electronic density functional reactivity indexes

A set of connections among several nuclear and electronic indexes of reactivity in the framework of the conceptual Density Functional Theory by using an expansion ofthe energy functional in terms of the total number of electrons and the normal coordinates within a canonical ensemble was derived. The relations obtained provided explicit links between important quantities related to the chemical reactivity of a system. This paper particularly demonstrates that the derivative of the electronic energy with respect to the external potential of a system in its equilibrium geometry was equal to the negative of the nuclear repulsion derivative with respect to the external potential

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Intramolecular basis set superposition error effects on the planarity of benzene and other aromatic molecules: A solution to the problem

Recently, the surprising result that ab initio calculations on benzene and other planar arenes at correlated MP2, MP3, configuration interaction with singles and doubles (CISD), and coupled cluster with singles and doubles levels of theory using standard Pople’s basis sets yield nonplanar minima has been reported. The planar optimized structures turn out to be transition states presenting one or more large imaginary frequencies, whereas single-determinant-based methods lead to the expected planar minima and no imaginary frequencies. It has been suggested that such anomalous behavior can be originated by two-electron basis set incompleteness error. In this work, we show that the reported pitfalls can be interpreted in terms of intramolecular basis set superposition error (BSSE) effects, mostly between the C–H moieties constituting the arenes. We have carried out counterpoise-corrected optimizations and frequency calculations at the Hartree–Fock, B3LYP, MP2, and CISD levels of theory with several basis sets for a number of arenes. In all cases, correcting for intramolecular BSSE fixes the anomalous behavior of the correlated methods, whereas no significant differences are observed in the single-determinant case. Consequently, all systems studied are planar at all levels of theory. The effect of different intramolecular fragment definitions and the particular case of charged species, namely, cyclopentadienyl and indenyl anions, respectively, are also discussed

EU cohesion aid to Spain: a data set. Part II: 1994-99 planning period

In this paper we construct a data set on EU cohesion aid to Spain during the planning period 1994-99. The data are disaggregated by region, year and function and attempt to approximate the timing of actual executed expenditure on assisted projects.

Drugs, schizotypy and cognition: Cognitive attenuations seem more consistently associated with substance use than schizotypal symptoms

Psychosis is a debilitating disease, causing harm to the individual and society. Since early detection of the disease is associated with a more benign course, factors are warranted that enable the early detection of psychosis. In the present thesis we will be focusing on two potential risk factors, namely schizotypy and drug use. The schizotypy concept, originally developed by Meehl (1962), states that schizophrenia symptoms exist on a spectrum, with symptoms ranging from the most severe in patients with schizophrenia to the least affected individual in the general population. Along the schizophrenia spectrum cognitive impairments are commonly found, for instance reduced hemispheric asymmetry or frontal lobe functions. The second risk factor (drug use), affects similar cognitive functions as those attenuated along the schizophrenia spectrum, and drug use is elevated in schizophrenia and people scoring high on schizotypy. Therefore, we set out to investigate whether cognitive attenuations formerly allocated to schizotypal symptoms could have been influenced by elevated substance use in this population. To test this idea, we assessed various drugs (nicotine, cannabis, mephedrone, general substance dependence) and schizotypy symptoms (O-LIFE), and measured either hemispheric asymmetry of function (left hemisphere dominance for language, and right hemisphere dominance for facial processing) or functions largely relying on the frontal lobes (such as cognitive flexibility, working memory, verbal short-term memory, verbal learning and verbal fluency). Results of all studies suggest that it is mostly drugs, and not schizotypy in general that predict cognitive functioning. Therefore, cognitive attenuations subscribed to schizotypy dimensions are likely to have been affected by enhanced drug use. Future studies should extend the list of potential risk factors (e.g. depression and IQ) to acquire a comprehensive overview of the most reliable predictors of disadvantageous cognitive profiles.

Esfingolipids i senyal•lització cel.lular: l'(E)-2-hexadecenal com a possible lípid bioactiu

L’esfingosina-1-fosfat (S1P) és un lípid bioactiu amb funcions crucials en la biologia cel•lular. Entre aquestes, la seva activitat mitogènica i citoprotectora són les més estudiades. L’S1P és catabolitzada intracel•lularment mitjançant l’esfingosina-1-fosfat liasa (SGPL1) per generar (E)-2-hexadecenal i fosforiletanolamina. L’objectiu d’aquest projecte és explorar si l’(E)-2-hexadecenal és realment un catabòlit innocu o bé si, pel seu caràcter acceptor de Michael, és capaç de reaccionar amb pèptids o proteïnes específics. Aquesta interacció podria traduïr-se en funcions biològiques determinades, algunes de les quals són possiblement atribuïdes a l’esfingosina-1-fosfat com a tal. Per poder explorar el potencials adductes proteïcs amb l’aldehid, s’han emprat, sobre cèl•lules HeLa que sobreexpressen SGPL1, sondes anàlegs a esfingosina i esfinganina (i els seus derivats fosforil•lats) que presenten una funció azida en la posició omega de la cadena esfingoide. Aquestes, mitjançant química click sense coure, s’han fet reaccionar amb una molècula que presenta un dibenzociclooctí unit a biotina DBCObiotina). Després d’aïllar les proteïnes així biotinilades amb una reïna d’estreptavidina, aquestes es van separar per electroforesi. Les bandes proteïques observades es van extreure del gel i es van digerir amb tripsina, per posteriorment analitzar els pèptids per MALDI-TOF, el que permetria l’identificació de proteïnes a partir de “peptide mass fingerprinting”. Lamentablement, a la fi d’aquest contracte, encara no s’ha pogut identificar cap proteïna que s’uneixi a l’aldehid alliberat per la reacció de l’esfingosina-1- fosfat liasa. No obstant, durant aquest temps s’ha millorat el mètode per detectar aquests adductes proteïcs. Per això, si la recerca continua en aquesta línia, properament es podria saber amb certesa si existeixen o no aquestes interaccions covalents entre determinades proteïnes i l’(E)-2-hexadecenal.

Review article: Intestinal epithelia and barrier functions.

The mucosal epithelia of the digestive tract acts as a selective barrier, permeable to ions, small molecules and macromolecules. These epithelial cells aid the digestion of food and absorption of nutrients. They contribute to the protection against pathogens and undergo continuous cell renewal which facilitates the elimination of damaged cells. Both innate and adaptive defence mechanisms protect the gastrointestinal-mucosal surfaces against pathogens. Interaction of microorganisms with epithelial cells triggers a host response by activating specific transcription factors which control the expression of chemokines and cytokines. This host response is characterized by the recruitment of macrophages and neutrophils at the site of infection. Disruption of epithelial signalling pathways that recruit migratory immune cells results in a chronic inflammatory response. The adaptive defence mechanism relies on the collaboration of epithelial cells (resident sampling system) with antigen-presenting and lymphoid cells (migratory sampling system); in order to obtain samples of foreign antigen, these samples must be transported across the barriers without affecting the integrity of the barrier. These sampling systems are regulated by both environmental and host factors. Fates of the antigen may differ depending on the way in which they cross the epithelial barrier, i.e. via interaction with motile dendritic cells or epithelial M cells in the follicle-associated epithelium.

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

Comparative gene finding in chicken indicates that we are closing in on the set of multi-exonic widely-expressed human genes

The recent availability of the chicken genome sequence poses the question of whether there are human protein-coding genes conserved in chicken that are currently not included in the human gene catalog. Here, we show, using comparative gene finding followed by experimental verification of exon pairs by RT–PCR, that the addition to the multi-exonic subset of this catalog could be as little as 0.2%, suggesting that we may be closing in on the human gene set. Our protocol, however, has two shortcomings: (i) the bioinformatic screening of the predicted genes, applied to filter out false positives, cannot handle intronless genes; and (ii) the experimental verification could fail to identify expression at a specific developmental time. This highlights the importance of developing methods that could provide a reliable estimate of the number of these two types of genes.

Splitting statistical potentials into meaningful scoring functions: testing the prediction of near-native structures from decoy conformations

Background: Recent advances on high-throughput technologies have produced a vast amount of protein sequences, while the number of high-resolution structures has seen a limited increase. This has impelled the production of many strategies to built protein structures from its sequence, generating a considerable amount of alternative models. The selection of the closest model to the native conformation has thus become crucial for structure prediction. Several methods have been developed to score protein models by energies, knowledge-based potentials and combination of both.Results: Here, we present and demonstrate a theory to split the knowledge-based potentials in scoring terms biologically meaningful and to combine them in new scores to predict near-native structures. Our strategy allows circumventing the problem of defining the reference state. In this approach we give the proof for a simple and linear application that can be further improved by optimizing the combination of Zscores. Using the simplest composite score () we obtained predictions similar to state-of-the-art methods. Besides, our approach has the advantage of identifying the most relevant terms involved in the stability of the protein structure. Finally, we also use the composite Zscores to assess the conformation of models and to detect local errors.Conclusion: We have introduced a method to split knowledge-based potentials and to solve the problem of defining a reference state. The new scores have detected near-native structures as accurately as state-of-art methods and have been successful to identify wrongly modeled regions of many near-native conformations.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.