Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

The classification of esterases: an important gene family involved in insecticide resistance - A review

The use of chemical insecticides continues to play a major role in the control of disease vector populations, which is leading to the global dissemination of insecticide resistance. A greater capacity to detoxify insecticides, due to an increase in the expression or activity of three major enzyme families, also known as metabolic resistance, is one major resistance mechanisms. The esterase family of enzymes hydrolyse ester bonds, which are present in a wide range of insecticides; therefore, these enzymes may be involved in resistance to the main chemicals employed in control programs. Historically, insecticide resistance has driven research on insect esterases and schemes for their classification. Currently, several different nomenclatures are used to describe the esterases of distinct species and a universal standard classification does not exist. The esterase gene family appears to be rapidly evolving and each insect species has a unique complement of detoxification genes with only a few orthologues across species. The examples listed in this review cover different aspects of their biochemical nature. However, they do not appear to contribute to reliably distinguish among the different resistance mechanisms. Presently, the phylogenetic criterion appears to be the best one for esterase classification. Joint genomic, biochemical and microarray studies will help unravel the classification of this complex gene family.

Substrate specificity of human kallikrein 2 (hK2) as determined by phage display technology.

Human glandular kallikrein 2 (hK2) is a trypsin-like serine protease expressed predominantly in the prostate epithelium. Recently, hK2 has proven to be a useful marker that can be used in combination with prostate specific antigen for screening and diagnosis of prostate cancer. The cleavage by hK2 of certain substrates in the proteolytic cascade suggest that the kallikrein may be involved in prostate cancer development; however, there has been very little other progress toward its biochemical characterization or elucidation of its true physiological role. In the present work, we adapt phage substrate technology to study the substrate specificity of hK2. A phage-displayed random pentapeptide library with exhaustive diversity was generated and then screened with purified hK2. Phages displaying peptides susceptible to hK2 cleavage were amplified in eight rounds of selection and genes encoding substrates were transferred from the phage to a fluorescent system using cyan fluorescent protein (derived from green fluorescent protein) that enables rapid determination of specificity constants. This study shows that hK2 has a strict preference for Arg in the P1 position, which is further enhanced by a Ser in P'1 position. The scissile bonds identified by phage display substrate selection correspond to those of the natural biological substrates of hK2, which include protein C inhibitor, semenogelins, and fibronectin. Moreover, three new putative hK2 protein substrates, shown elsewhere to be involved in the biology of the cancer, have been identified thus reinforcing the importance of hK2 in prostate cancer development.

Evolución diacrónica del sistema verbal de movimiento en español: entrar y salir (ss. XII-XV)

Esta investigación es una tesis doctoral que tiene como objetivo el estudio de la evolución semántica de los verbos entrar y salir desde el punto de vista diacrónico en español medieval mediante la aplicación de las herramientas teóricas proporcionadas por la Lingüística cognitiva. A través de un corpus de documentaciones extraídas del Corpus del nuevo diccionario histórico del español de la Real Academia Española y que abarcan el intervalo desde el siglo XIII al XV, se analizan los valores semánticos de ambos verbos con la finalidad de determinar cómo fue desarrollándose su uso, qué significados generaron otros nuevos y cómo se relacionan entre sí los significados de un mismo verbo. Para llevar a cabo esta tarea se ha partido de la Teoría de los prototipos aplicada a la categorización léxica, así como de la Teoría de la metáfora y la metonimia. Asimismo, se ha analizado la estructura sintáctica de cada ejemplo, la selección léxica de los argumentos verbales y la tradición discursiva a la que pertenece cada documentación. De esta forma, es posible, no sólo describir la evolución semántica de ambos verbos, sino completar la información sobre las causas de su desarrollo diacrónico. Finalmente, se establecen las relaciones semánticas de carácter diacrónico existentes entre ambos verbos en tanto que miembros de un mismo grupo léxico, los verbos de movimiento, y se corrobora su vinculación más allá de una relación de antonimia de tipo primario.

Modelling world investment markets using threshold conditional correlation models

In this paper we propose a parsimonious regime-switching approach to model the correlations between assets, the threshold conditional correlation (TCC) model. This method allows the dynamics of the correlations to change from one state (or regime) to another as a function of observable transition variables. Our model is similar in spirit to Silvennoinen and Teräsvirta (2009) and Pelletier (2006) but with the appealing feature that it does not suffer from the course of dimensionality. In particular, estimation of the parameters of the TCC involves a simple grid search procedure. In addition, it is easy to guarantee a positive definite correlation matrix because the TCC estimator is given by the sample correlation matrix, which is positive definite by construction. The methodology is illustrated by evaluating the behaviour of international equities, govenrment bonds and major exchange rates, first separately and then jointly. We also test and allow for different parts in the correlation matrix to be governed by different transition variables. For this, we estimate a multi-threshold TCC specification. Further, we evaluate the economic performance of the TCC model against a constant conditional correlation (CCC) estimator using a Diebold-Mariano type test. We conclude that threshold correlation modelling gives rise to a significant reduction in portfolio´s variance.

Full Exploration of the Diels–Alder Cycloaddition on Metallofullerenes M3N@C80 (M=Sc, Lu, Gd): The D5h versus Ih Isomer and the Influence of the Metal Cluster

In this work a detailed investigation of the exohedral reactivity of the most important and abundant endohedral metallofullerene (EMF) is provided, that is, Sc3N@Ih-C80 and its D5h counterpart Sc3N@D5h-C80, and the (bio)chemically relevant lutetium- and gadolinium-based M3N@Ih/D5h-C80 EMFs (M=Sc, Lu, Gd). In particular, we analyze the thermodynamics and kinetics of the Diels–Alder cycloaddition of s-cis-1,3-butadiene on all the different bonds of the Ih-C80 and D5h-C80 cages and their endohedral derivatives. First, we discuss the thermodynamic and kinetic aspects of the cycloaddition reaction on the hollow fullerenes and the two isomers of Sc3N@C80. Afterwards, the effect of the nature of the metal nitride is analyzed in detail. In general, our BP86/TZP//BP86/DZP calculations indicate that [5,6] bonds are more reactive than [6,6] bonds for the two isomers. The [5,6] bond D5h-b, which is the most similar to the unique [5,6] bond type in the icosahedral cage, Ih-a, is the most reactive bond in M3N@D5h-C80 regardless of M. Sc3N@C80 and Lu3N@C80 give similar results; the regioselectivity is, however, significantly reduced for the larger and more electropositive M=Gd, as previously found in similar metallofullerenes. Calculations also show that the D5h isomer is more reactive from the kinetic point of view than the Ih one in all cases which is in good agreement with experiments

Calorimetry of dehydrogenation and dangling-bond recombination in several hydrogenated amorphous silicon materials

Differential scanning calorimetry (DSC) was used to study the dehydrogenation processes that take place in three hydrogenated amorphous silicon materials: nanoparticles, polymorphous silicon, and conventional device-quality amorphous silicon. Comparison of DSC thermograms with evolved gas analysis (EGA) has led to the identification of four dehydrogenation processes arising from polymeric chains (A), SiH groups at the surfaces of internal voids (A'), SiH groups at interfaces (B), and in the bulk (C). All of them are slightly exothermic with enthalpies below 50 meV/H atoms , indicating that, after dissociation of any SiH group, most dangling bonds recombine. The kinetics of the three low-temperature processes [with DSC peak temperatures at around 320 (A),360 (A'), and 430°C (B)] exhibit a kinetic-compensation effect characterized by a linea relationship between the activation entropy and enthalpy, which constitutes their signature. Their Si-H bond-dissociation energies have been determined to be E (Si-H)0=3.14 (A), 3.19 (A'), and 3.28 eV (B). In these cases it was possible to extract the formation energy E(DB) of the dangling bonds that recombine after Si-H bond breaking [0.97 (A), 1.05 (A'), and 1.12 (B)]. It is concluded that E(DB) increases with the degree of confinement and that E(DB)>1.10 eV for the isolated dangling bond in the bulk. After Si-H dissociation and for the low-temperature processes, hydrogen is transported in molecular form and a low relaxation of the silicon network is promoted. This is in contrast to the high-temperature process for which the diffusion of H in atomic form induces a substantial lattice relaxation that, for the conventional amorphous sample, releases energy of around 600 meV per H atom. It is argued that the density of sites in the Si network for H trapping diminishes during atomic diffusion

Bitumens in the late Variscan hydrothermal vein-type uranium deposit of Pribram, Czech Republic: Sources, radiation-induced alteration, and relation to mineralization

The late Variscan (275-278 Ma) Pribram uranium deposit is one of the largest known accumulations of uraniferous bitumens in hydrothermal veins. The deposit extends along the northwestern boundary of the Central Bohemian pluton (345-335 Ma) with low-grade metamorphosed Late Proterozoic and unmetamorphosed Cambrian rocks. From a net uranium production of 41,742 metric tons (t), more than 6,000 t were extracted from bitumen-uraninite ores during 43 years of exploration and mining. Three morphological varieties of solid bitumen are recognized: globular, asphaltlike, and cokelike. While the globular bitumen is uranium free, the other two types are uraniferous. The amount of bitumen in ore veins gradually decreases toward the contact with the plutonic body and increases with depth. Two types of bitumen microtextures are recognized using high-resolution transmission electron microscopy: amorphous and microporous, the former being less common in uraniferous samples. A lower Raman peak area ratio (1,360/1,575 cm(-1)) in mineralized bitumens (0.9) compared with uranium-free samples (2.0) indicates a lower degree of microtextural organization in the latter The H/C and O/C atomic ratios in uranium-free bitumens (0.9-1.1 and 0.09, respectively) are higher than those in mineralized samples (H/C = 0.3-0.8, O/C = 0.03-0.09). The chloroform extractable matter yield is Very low in uranium-free bitumens (0.30-0.35% of the total organic carbon,TOC) and decreases with uranium content increase. The extracted solid uraniferous bitumen infrared spectra show depletion in aliphatic CH2 and CH3 groups compared to uranium-free samples. The concentration of oxygen-bearing functional groups relative to aromatic bonds in the IR spectra of uranium-free and mineralized bitumen, however, do not differ significantly. C-13 NMR confirmed than the aromaticity of a uraniferous sample is higher (F-ar = 0.61) than in the uranium-free bitumen (F-ar = 0.51). Pyrolysates from uraniferous and nonuraniferous bitumens do not differ significantly, being predominantly cresol, alkylphenols, alkylbenzenes, and alkylnaphthalenes. The liquid pyrolysate yield decreases significantly with increasing uranium content. The delta(13)C Values of bulk uranium-free bitumens and low-grade uraniferous, asphaltlike bitumens range from -43.6 to 52.3 per mil. High-grade, cokelike, uraniferous bitumens are more C-13 depleted (54.5 to -58.4 parts per thousand). In contrast to the very light isotopic ratios of the high-grade uraniferous cokelike bitumen bulk carbon, the individual n-alkanes and isoprenoids (pristane and phytane) extracted from the same sample are significantly C-13 enriched. The isotopic composition of the C13-24 n-alkanes extracted from the high-grade uraniferous sample (delta(13)C = -28.0 to 32.6 parts per thousand) are heavier compared with the same compounds in a uranium-free sample (delta(13)C = 31.9 to 33.8 parts per thousand). It is proposed that the bitumen source was the isotopically light (delta(13)C = 35.8 to 30.2 parts per thousand) organic matter of the Upper Proterozoic host rocks that were pyrolyzed during intrusion of the Central Bohemian pluton. The C-13- depleted pyrolysates were mobilized from the innermost part of the contact-metamorphic aureole, accumulated in structural traps in less thermally influenced parts of the sedimentary complex and were later extracted by hydrothermal fluids. Bitumens at the Pribram deposit are younger than the main part of the uranium mineralization and were formed through water-washing and radiation-induced polymerization of both the gaseous and liquid pyrolysates. Direct evidence for pyrolysate reduction of uranium in the hydrothermal system is difficult to obtain as the chemical composition of the original organic fluid phase was modified during water-washing and radiolytic alteration. However, indirect evidence-e.g., higher O/C atomic ratios in uranium-free bitumens (0.1) relative to the Upper Proterozoic source rocks (0.02-0.05), isotopically very light carbon in associated whewellite (delta(13)C = 31.7 to -28.4 parts per thousand), and the striking absence of bitumens in the pre-uranium, hematite stage of the mineralization-indicates that oxidation of organic fluids may have contributed to lowering of aO(2) and uraninite precipitation.

Critical analysis and extension of the Hirshfeld atoms in molecules

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting "Hirshfeld-I charges" correlate well with electrostatic potential derived atomic charges

The relevance of the Laplacian of intracule and extracule density distributions for analyzing electron-electron interactions in molecules

A topological analysis of intracule and extracule densities and their Laplacians computed within the Hartree-Fock approximation is presented. The analysis of the density distributions reveals that among all possible electron-electron interactions in atoms and between atoms in molecules only very few are located rigorously as local maxima. In contrast, they are clearly identified as local minima in the topology of Laplacian maps. The conceptually different interpretation of intracule and extracule maps is also discussed in detail. An application example to the C2H2, C2H4, and C2H6 series of molecules is presented

Determination of vibrational polarizabilities and hyperpolarizabilities using field-induced coordinates

An analytical set of field-induced coordinates is defined and is used to show that the vibrational degrees of freedom required to completely describe nuclear relaxation polarizabilities and hyperpolarizabilities is reduced from 3N-6 to a relatively small number. As this number does not depend upon the size of the molecule, the process provides computational advantages. A method is provided to separate anharmonic contributions from harmonic contributions as well as effective mechanical from electrical anharmonicity. The procedures are illustrated by Hartree-Fock calculations, indicating that anharmonicity can be very important

Electron localization function at the correlated level

The electron localization function (ELF) has been proven so far a valuable tool to determine the location of electron pairs. Because of that, the ELF has been widely used to understand the nature of the chemical bonding and to discuss the mechanism of chemical reactions. Up to now, most applications of the ELF have been performed with monodeterminantal methods and only few attempts to calculate this function for correlated wave functions have been carried out. Here, a formulation of ELF valid for mono- and multiconfigurational wave functions is given and compared with previous recently reported approaches. The method described does not require the use of the homogeneous electron gas to define the ELF, at variance with the ELF definition given by Becke. The effect of the electron correlation in the ELF, introduced by means of configuration interaction with singles and doubles calculations, is discussed in the light of the results derived from a set of atomic and molecular systems

The curvature of the conical intersection seam: an approximate second-order analysis

We present a method for analyzing the curvature (second derivatives) of the conical intersection hyperline at an optimized critical point. Our method uses the projected Hessians of the degenerate states after elimination of the two branching space coordinates, and is equivalent to a frequency calculation on a single Born-Oppenheimer potential-energy surface. Based on the projected Hessians, we develop an equation for the energy as a function of a set of curvilinear coordinates where the degeneracy is preserved to second order (i.e., the conical intersection hyperline). The curvature of the potential-energy surface in these coordinates is the curvature of the conical intersection hyperline itself, and thus determines whether one has a minimum or saddle point on the hyperline. The equation used to classify optimized conical intersection points depends in a simple way on the first- and second-order degeneracy splittings calculated at these points. As an example, for fulvene, we show that the two optimized conical intersection points of C2v symmetry are saddle points on the intersection hyperline. Accordingly, there are further intersection points of lower energy, and one of C2 symmetry - presented here for the first time - is found to be the global minimum in the intersection space

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed