Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Dependence of Response Functions and Orbital Functionals on Occupation Numbers

Explicitly orbital-dependent approximations to the exchange-correlation energy functional of density functional theory typically not only depend on the single-particle Kohn-Sham orbitals but also on their occupation numbers in the ground-state Slater determinant. The variational calculation of the corresponding exchange-correlation potentials with the optimized effective potential (OEP) method therefore also requires a variation of the occupation numbers with respect to a variation in the effective single-particle potential, which is usually not taken into account. Here it is shown under which circumstances this procedure is justified.

Structural effects and hydrogen bonds on N,N `-di(methoxycarbonylsulfenyl)urea, [CH(3)OC(O)SNH](2)CO, studied by experimental and theoretical methods

Pure N,N`-di(methoxycarbonylsulfenyl)urea, [CH(3)OC(O)SNH](2)CO, is quantitatively prepared by the hydrolysis reaction of CH(3)OC(O)SNCO and characterized by (1)H NMR, GC-MS and FTIR spectroscopy techniques. Structural and conformational properties are analyzed using a combined approach with data obtained from X-ray diffraction, vibrational spectra and theoretical calculation methods. The IR and Raman spectra for normal and deuterated species are reported. The crystal structure of [CH(3)OC(O)SNH](2)CO was determined by X-ray diffraction methods. The substance crystallizes in the orthorhombic P2(1)2(1)2 space group with a = 9.524(2), b = 12.003(1), c = 4.481 (1) angstrom, and Z = 2 moieties in the unit cell. The molecule is sited on a twofold crystallographic axis (C(2)) parallel to c and shows the anti-anti conformation (S-N single bonds antiperiplanar with respect to the opposite C-N single bonds in sulfenyl-urea-sic group). Neighboring molecules are arranged in a chain motif that extends along the C(2)-axis and is held by bifurcated NH center dot center dot center dot O center dot center dot center dot HN intermolecular bonds. A local planar symmetry is observed in the crystal for the central -SN(H)C(O)N(H)S- skeleton. Experimental and calculated data allow to trace this structural feature to the occurrence of N-H center dot center dot center dot O=C hydrogen bonding interactions. Calculated vibrational and structural properties are in good agreement with the experimentally determined features. (C) 2008 Elsevier B.V. All rights reserved.

Sugarcane Phosphoribosyl Pyrophosphate Synthetase: Molecular Characterization of a Phosphate-independent PRS

Phosphoribosyl pyrophosphate synthetase (PRS-EC:2.7.6.1) is an important enzyme present in several metabolic pathways, thus forming a complex family of isoenzymes. However, plant PRS enzymes have not been extensively investigated. In this study, a sugarcane prs gene has been characterized from the Sugar Cane Expressed Sequence Tag Genome Project. This gene contains a 984-bp open reading frame encoding a 328-amino acid protein. The predicted amino acid sequence has 77% and 78% amino acid sequence identity to Arabidopsis thaliana and Spinacia oleracea PRS4, respectively. The assignment of sugarcane PRS as a phosphate-independent PRS isoenzyme (Class II PRS) is verified following enzyme assay and phylogenetic reconstruction of PRS homologues. To gain further insight into the structural framework of the phosphate independence of sugarcane PRS, a molecular model is described. This model reveals the formation of two conserved domains elucidating the structural features involved in sugarcane PRS phosphate independence. The recombinant PRS retains secondary structure elements and a quaternary arrangement consistent with known PRS homologues, based on circular dichroism measurements.

Orbital-polarization terms: From a phenomenological to a first-principles description of orbital magnetism in density-functional theory

Phenomenological orbital-polarizition (OP) terms have been repeatedly introduced in the single-particle equations of spin-density-functional theory, in order to improve the description of orbital magnetic moments in systems containing transition metal ions. Here we show that these ad hoc corrections can be interpreted as approximations to the exchange-correlation vector potential A(xc) of current-density functional theory (CDFT). This connection provides additional information on both approaches: phenomenological OP terms are connected to first-principles theory, leading to a rationale for their empirical success and a reassessment of their limitations and the approximations made in their derivation. Conversely, the connection of OP terms with CDFT leads to a set of simple approximations to the CDFT potential A(xc), with a number of desirable features that are absent from electron-gas-based functionals. (C) 2008 Wiley Periodicals, Inc.

Detailed molecular characterization of cord blood-derived endothelial progenitors

Objective. Given their involvement in pathological and physiological angiogenesis, there has been growing interest in understanding and manipulating endothellial progenitor cells (EPC) for therapeutic purposes. However, detailed molecular analysis of EPC before and during endothelial differentiation is lacking and is the subject of the present study. Materials and Methods. We report a detailed microarray gene-expression profile of freshly isolated (day 0) human cord blood (CB)-derived EPC (CD133(+)KDR(+) or CD34(+)KDR(+)), and at different time points during in vitro differentiation (early: day 13; late: day 27). Results. Data obtained reflect an EPC transcriptome enriched in genes related to stem/progenitor cells properties (chromatin remodeling, self-renewal, signaling, cytoskeleton organization and biogenesis, recruitment, and adhesion). Using a complementary DNA microarray enriched in intronic transcribed sequences, we observed, as well, that naturally transcribed intronic noncoding RNAs were specifically expressed at the EPC stage. Conclusion. Taken together, we have defined the global gene-expression profile of CB-derived EPC during the process of endothelial differentiation, which can be used to identify genes involved in different vascular pathologies. (C) 2008 ISEH - Society for Hematology and Stem Cells. Published by Elsevier Inc.

Host-guest system of 4-nerolidylcatechol in 2-hydroxypropyl-beta-cyclodextrin: preparation, characterization and molecular modeling

The interaction of 4-nerolidylcatechol (4-NRC), a potent antioxidant agent, and 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD) was investigated by the solubility method using Fourier transform infrared (FTIR) methods in addition to UV-Vis, (1)H-nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The inclusion complexes were prepared using grinding, kneading and freeze-drying methods. According to phase solubility studies in water a B(S)-type diagram was found, displaying a stoichiometry complexation of 2:1 (drug:host) and stability constant of 6494 +/- A 837 M(-1). Stoichiometry was established by the UV spectrophotometer using Job`s plot method and, also confirmed by molecular modeling. Data from (1)H-NMR, and FTIR, experiments also provided formation evidence of an inclusion complex between 4-NRC and HP-beta-CD. 4-NRC complexation indeed led to higher drug solubility and stability which could probably be useful to improve its biological properties and make it available to oral administration and topical formulations.

New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.

Expedient, accurate methods for the determination of the degree of substitution of cellulose carboxylic esters: Application of UV-vis spectroscopy (dye solvatochromism) and FTIR

Two techniques, namely UV-vis- and FTIR spectroscopy, have been employed in order to calculate the degree of substitution (DS) of cellulose carboxylic esters, including acetates, CAs, butyrates, CBs, and hexanoates, CHs. Regarding UV-vis spectroscopy, we have employed a novel approach, based on measuring the dependence of lambda(max) of the intra-molecular charge-transfer bands of polarity probes adsorbed on DS of the ester films (solvatochromism). Additionally, we have revisited the use of FTIR for DS determination. Several methods have been used in order to plot Beer`s law graph, namely: Absorption of KBr pellets, pre-coated with CA: reflectance (DRIFTS) of CAs-KBr solid-solid mixtures with, or without the use of 1.4-dicyanobenzene as an internal reference; reflectance of KBr powder pre-coated with CA. The methods indicated are simple, fast, and accurate, requiring much less ester than the titration method. The probe method is independent of the experimental variables examined. (c) 2010 Published by Elsevier Ltd.

Influence of neodymium ions on photocatalytic activity of TiO(2) synthesized by sol-gel and precipitation methods

Titanium dioxide with and without the addition of neodymium ions was prepared using sol-gel and precipitation methods. The resulting catalysts were characterized by thermal analysis, X-ray diffraction and BET specific surface area. Neodymium addition exerted a remarkable influence on the phase transition temperature and the surface properties of the TiO(2) matrix. TiO(2) samples synthesized by precipitation exhibit an exothermic event related from the amorphous to anatase phase transition at 510 degrees C, whereas in Nd-doped TiO(2) this transition occurred at 527 degrees C. A similar effect was observed in samples obtained using sol-gel method. The photocatalytic reactivity of the catalysts was evaluated by photodegradation of Remazol Black B (RB) under ultraviolet irradiation. Nd-doped TiO(2) showed enhanced photodegradation ability compared to undoped TiO(2) samples, independent of the method of synthesis. In samples obtained by sol-gel, RB decoloration was enhanced by 16% for TiO(2) doped with 0.5% neodymium ions. (C) 2010 Elsevier B.V. All rights reserved.

New Molecular Species of Potential Interest to Atmospheric Chemistry: Isomers on the [H, S(2), Br] Potential Energy Surface

This work reports a state-of-the-art theoretical characterization of four new sulfur-bromine species and five transition states on the [H, S(2), Br] potential energy surface. Our highest level theoretical approach employed the method coupled cluster singles and doubles with perturbative contributions of connected triples, CCSD(T), along with the series of correlation-consistent basis sets and with extrapolation to the complete basis set (CBS) limit in the optimization of the geometrical parameters and to quantify the energetic quantities. The structural and vibrational frequencies here reported are unique and represent the most accurate investigation to date of these species. The global minimum corresponds to a skewed structure HSSBr with a disulfide bond; this is followed by a pyramidal-like structure, SSHBr, 18.85 kcal/mol above the minimum. Much higher in energy, we found another skewed structure, HSBrS (50.29 kcal/mol), with one S-Br dative-type bond, and another pyramidal-like one, HBrSS (109.80 kcal/mol), with two S-Br dative-type bonds. The interconversion of HSSBr into SSHBr can occur via a transfer of either the hydrogen or the bromine atom but involves a very high barrier of about 43 kcal/mol. These molecules are potentially a new route of coupling the sulfur and bromine chemistry in the atmosphere, and conditions of high concentration of H(2)S like in volcanic eruptions might contribute to their formation. We note that HSSBr can act as a reservoir molecule for the reaction between the radicals HSS and Br. Also, an assessment of the methods DFT/B3LYP/CBS and MP2/CBS relative to CCSD(T)/CBS provides insights on the expected performance of these methods on the characterization of polysulfides and also of more complex systems containing disulfide bridges.

Insights into the Molecular Requirements for the Anti-obesity Activity of a Series of CB1 Ligands

Two-dimensional and 3D quantitative structure-activity relationships studies were performed on a series of diarylpyridines that acts as cannabinoid receptor ligands by means of hologram quantitative structure-activity relationships and comparative molecular field analysis methods. The quantitative structure-activity relationships models were built using a data set of 52 CB1 ligands that can be used as anti-obesity agents. Significant correlation coefficients (hologram quantitative structure-activity relationships: r 2 = 0.91, q 2 = 0.78; comparative molecular field analysis: r 2 = 0.98, q 2 = 0.77) were obtained, indicating the potential of these 2D and 3D models for untested compounds. The models were then used to predict the potency of an external test set, and the predicted (calculated) values are in good agreement with the experimental results. The final quantitative structure-activity relationships models, along with the information obtained from 2D contribution maps and 3D contour maps, obtained in this study are useful tools for the design of novel CB1 ligands with improved anti-obesity potency.

On the molecular mass of the extracellular hemoglobin of Glossoscolex paulistus: Analytical ultracentrifugation reexamination

The giant extracellular hemoglobin of Glossoscolex paulistus (HbGp) is constituted by Subunits containing heme groups with molecular masses (M) in the range of 15 to 19 kDa, monomers of 16 kDa (d), and trimers of 51 to 52 kDa (abc) linked by nonheme structures named linkers of 24 to 32 kDa (L). HbGp is homologous to Lumbricus terrestris hemoglobin (HbLt). Several reports propose M of HbLt in the range of 3.6 to 4.4 MDa. Based on subunits M determined by mass spectrometry and assuming HbGp stoichiometry of 12(abcd)(3)L(3) (Vinogradov model) plus 144 heme groups, a Value of M for HbGp oligomer of 3560 kDa can be predicted. This Value is nearly 500 kDa higher than the unique HbGp M Value reported in the literature. In the current work, sedimentation velocity analytical ultracentrifugation (AUC) experiments were performed to obtain M for HbGp in oxy and cyano-met forms. s(20,w)(0), values of 58.1 +/- 0.2 S and 59.6 +/- 0.2 S, respectively, for the two oxidation forms were obtained. The ratio between sedimentation and diffusion coefficients supplied values for M of approximately 3600 100 and 3700 100 kDa for oxy and cyano-met HbGp forms, respectively. An independent determination of the partial specific volume, V(bar), for HbGp was performed based on density measurements, providing a value of 0.764 +/- 0.008, in excellent agreement with the estimates from SEDFIT software. Our results show total consistency between M obtained by AUC and recent partial characterization by mass spectrometry. Therefore, HbGp possesses M very close to that of HbLt, suggesting an oligomeric assembly in agreement with the Vinogradov model. (c) 2008 Elsevier Inc. All rights reserved.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Development of a molecular test to determine the vitality status of Norway spruce (Picea abies) seedlings during frozen storage

In boreal forest regions, a great portion of forest tree seedlings are stored indoors in late autumn to prevent seedlings from outdoor winter damage. For seedlings to be able to survive in storage it is crucial that they store well and can cope with the dark and cold storage environment. The aim of this study was to search for genes that can determine the vitality status of Norway spruce (Picea abies (L.) Karst.) seedlings during frozen storage. Furthermore, the sensitivity of the ColdNSure (TM) test, a gene activity test that predicts storability was assessed. The storability of seedlings was tested biweekly by evaluating damage with the gene activity test and the electrolyte leakage test after freezing seedlings to -25 A degrees C (the SELdiff-25 method). In parallel, seedlings were frozen stored at -3 A degrees C. According to both methods, seedlings were considered storable from week 41. This also corresponded to the post storage results determined at the end of the storage period. In order to identify vitality indicators, Next Generation Sequencing (NGS) was performed on bud samples collected during storage. Comparing physiological post storage data to gene analysis data revealed numerous vitality related genes. To validate the results, a second trial was performed. In this trial, gene activity was better in predicting seedling storability than the conventional freezing test; this indicates a high sensitivity level of this molecular assay. For multiple indicators a clear switch between damaged and vital seedlings was observed. A collection of indicators will be used in the future development of a commercial vitality test.