Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

Theoretical investigation of atomic and electronic structures of Ga(2)O(3)(ZnO)(6)

Transparent conducting oxides (TCO) are widely used in technological applications ranging from photovoltaics to thin-film transparent field-effect transistors. In this work we report a first-principles investigation, based on density-functional theory, of the atomic and electronic properties of Ga(2)O(3)(ZnO)(6) (GZO(6)), which is a promising candidate to be used as host oxide for wide band gap TCO applications. We identify a low-energy configuration for the coherent distribution of the Ga and Zn atoms in the cation positions within the experimentally reported orthorhombic GZO(6) structure. Four Ga atoms are located in four-fold sites, while the remaining 12 Ga atoms in the unit cell form four shared Ga agglomerates (a motif of four atoms). The Zn atoms are distributed in the remaining cation sites with effective coordination numbers from 3.90 to 4.50. Furthermore, we identify the natural formation of twin-boundaries in GZO(6), which can explain the zigzag modulations observed experimentally by high-resolution transmission electron microscopy in GZO(n) (n=9). Due to the intrinsic twin-boundary formation, polarity inversion in the ZnO tetrahedrons is present which is facilitated by the formation of the Ga agglomerates. Our analysis shows that the formation of fourfold Ga sites and Ga agglomerates are stabilized by the electronic octet rule, while the distribution of Ga atoms and the formation of the twin-boundary help alleviate excess strain. Finally we identify that the electronic properties of GZO(6) are essentially determined by the electronic properties of ZnO, i.e., there are slight changes in the band gap and optical absorption properties.

Atomistic origins of the phase transition mechanism in Ge(2)Sb(2)Te(5)

The fast and reversible phase transition mechanism between crystalline and amorphous phases of Ge(2)Sb(2)Te(5) has been in debate for several years. Through employing first-principles density functional theory calculations, we identify a direct structural link between the metastable crystalline and amorphous phases. The phase transition is driven by the displacement of Ge atoms along the rocksalt [111] direction from stable octahedron to high energy unstable tetrahedron sites close to the intrinsic vacancy regions, which generates a high energy intermediate phase between metastable and amorphous phases. Due to the instability of Ge at the tetrahedron sites, the Ge atoms naturally shift away from those sites, giving rise to the formation of local-ordered fourfold motifs and the long-range structural disorder. Intrinsic vacancies, which originate from Sb(2)Te(3), lower the energy barrier for Ge displacements, and hence, their distribution plays an important role in the phase transition. The high energy intermediate configuration can be obtained experimentally by applying an intense laser beam, which overcomes the thermodynamic barrier from the octahedron to tetrahedron sites. The high figure of merit of Ge(2)Sb(2)Te(5) is achieved from the optimal combination of intrinsic vacancies provided by Sb(2)Te(3) and the instability of the tetrahedron sites provided by GeTe.

Electron theoretical investigation of the stability of the B2-TiFe compound

The metastable phase diagram of the BCC-based ordering equilibria in the Ti-Fe system has been calculated using a truncated cluster expansion, through the combination of FP-LAPW and cluster variation method (CVM) in the irregular tetrahedron cluster approximation. The results are compared with phenomenological CVM assessments of the system and suggest that the value for the experimental formation enthalpy of the B2-TiFe compound should be significantly more negative than the currently assessed value. (C) 2008 Elsevier B.V. All rights reserved.

Thermodynamic self-consistency issues related to the Cluster Variation Method: The case of the BCC Cr-Fe (Chromium-Iron) system

The Cluster Variation Method (CVM), introduced over 50 years ago by Prof. Dr. Ryoichi Kikuchi, is applied to the thermodynamic modeling of the BCC Cr-Fe system in the irregular tetrahedron approximation, using experimental thermochemical data as initial input for accessing the model parameters. The results are checked against independent data on the low-temperature miscibility gap, using increasingly accurate thermodynamic models, first by the inclusion of the magnetic degrees of freedom of iron and then also by the inclusion of the magnetic degrees of freedom of chromium. It is shown that a reasonably accurate description of the phase diagram at the iron-rich side (i.e. the miscibility gap borders and the Curie line) is obtained, but only at expense of the agreement with the above mentioned thermochemical data. Reasons for these inconsistencies are discussed, especially with regard to the need of introducing vibrational degrees of freedom in the CVM model. (C) 2008 Elsevier Ltd. All rights reserved.

QSAR Modeling of a Set of Pyrazinoate Esters as Antituberculosis Prodrugs

Tuberculosis is an infection caused mainly by Mycobacterium tuberculosis. A first-line antimycobacterial drug is pyrazinamide (PZA), which acts partially as a prodrug activated by a pyrazinamidase releasing the active agent, pyrazinoic acid (POA). As pyrazinoic acid presents some difficulty to cross the mycobacterial cell wall, and also the pyrazinamide-resistant strains do not express the pyrazinamidase, a set of pyrazinoic acid esters have been evaluated as antimycobacterial agents. In this work, a QSAR approach was applied to a set of forty-three pyrazinoates against M. tuberculosis ATCC 27294, using genetic algorithm function and partial least squares regression (WOLF 5.5 program). The independent variables selected were the Balaban index (I), calculated n-octanol/water partition coefficient (ClogP), van-der-Waals surface area, dipole moment, and stretching-energy contribution. The final QSAR model (N = 32, r(2) = 0.68, q(2) = 0.59, LOF = 0.25, and LSE = 0.19) was fully validated employing leave-N-out cross-validation and y-scrambling techniques. The test set (N = 11) presented an external prediction power of 73%. In conclusion, the QSAR model generated can be used as a valuable tool to optimize the activity of future pyrazinoic acid esters in the designing of new antituberculosis agents.

Assignment of the absolute configuration at the sulfur atom of thioridazine metabolites by the analysis of their chiroptical properties: The case of thioridazine 2-sulfoxide

Thioridazine (THD) is a commonly prescribed phenotiazine neuroleptic drug, which is extensively biotransformed in the organism producing as main metabolites sulfoxides and a sulfone by sulfur oxidation Significant differences have been observed in the activity of the THD enantiomers as well as for its main metabolites, and enantioselectivity phenomena have been proved in the metabolic pathway. Here the assignment of the absolute configuration at the sulfur atom of enantiomeric THD-2-sulfoxide (THD-2-SO) has been carried out by circular dichroism (CD) spectroscopy The stereoisomers were separated by HPLC on Chiralpak AS column, recording the CD spectra for the two collected enantiomeric fractions The theoretical electronic CD spectrum has been obtained by the TDDFT/B3LYP/6-31G*. as Boltzmann averaging of the contributions calculated for the most stable conformations of the drug The comparison of the simulated and experimental spectra allowed the absolute configuration at the sulfur atom of the four THD-2-SO stereoisomers to be assigned The developed method should be useful for a reliable correlation between stereochemistry and activity and/or toxicity

Ca(2+) binding to c-state of adenine nucleotide translocase (ANT)-surrounding cardiolipins enhances (ANT)-Cys(56) relative mobility: A computational-based mitochondrial permeability transition study

The oxidation of critical cysteines/related thiols of adenine nucleotide translocase (ANT) is believed to be an important event of the Ca(2+)-induced mitochondrial permeability transition (MPT), a process mediated by a cyclosporine A/ADP-sensitive permeability transition pores (PTP) opening. We addressed the ANT-Cys(56) relative mobility status resulting from the interaction of ANT/surrounding cardiolipins with Ca(2+) and/or ADP by means of computational chemistry analysis (Molecular Interaction Fields and Molecular Dynamics studies), supported by classic mitochondrial swelling assays. The following events were predicted: (i) Ca(2+) interacts preferentially with the ANT surrounding cardiolipins bound to the H4 helix of translocase, (ii) weakens the cardiolipins/ANT interactions and (iii) destabilizes the initial ANT-Cys(56) residue increasing its relative mobility. The binding of ADP that stabilizes the conformation ""m"" of ANT and/or cardiolipin, respectively to H5 and H4 helices, could stabilize their contacts with the short helix h56 that includes Cys(56), accounting for reducing its relative mobility. The results suggest that Ca(2+) binding to adenine nucleotide translocase (ANT)-surrounding cardiolipins in c-state of the translocase enhances (ANT)-Cys(56) relative mobility and that this may constitute a potential critical step of Ca(2+)-induced PTP opening. (C) 2009 Elsevier B.V. All rights reserved.

Using Computer-aided Drug Design and Medicinal Chemistry Strategies in the Fight Against Diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics. ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Structural and functional properties of Bp-LAAO, a new L-amino acid oxidase isolated from Bothrops pauloensis snake venom

An L-amino acid oxidase (Bp-LAAO) from Bothrops pauloensis snake venom was highly purified using sequential chromatography steps on CM-Sepharose, Phenyl-Sepharose CL4B, Benzamidine Sepharose and C18 reverse-phase HPLC. Purified Bp-LAAO showed to be a homodimeric acidic glycoprotein with molecular weight around 65 kDa under reducing conditions in SDS-PAGE. The best substrates for Bp-LAAO were L-Met, L-Leu, L-Phe and L-Ile and the enzyme showed a strong reduction of its catalytic activity upon L-Met and L-Phe substrates at extreme temperatures. Bp-LAAO showed leishmanicidal, antitumoral and bactericidal activities dose dependently. Bp-LAAO induced platelet aggregation in platelet-rich plasma and this activity was inhibited by catalase. Bp-LAAC-cDNA of 1548 bp codified a mature protein with 516 amino acid residues corresponding to a theoretical isoelectric point and molecular weight of 6.3 and 58 kDa, respectively. Additionally, structural and phylogenetic studies identified residues under positive selection and their probable location in Elp-LAAO and other snake venom LAAOs (svLAAOs). Structural and functional investigations of these enzymes can contribute to the advancement of toxinology and to the elaboration of novel therapeutic agents. (C) 2009 Elsevier Masson SAS. All rights reserved.

Computer-aided drug design and ADMET predictions for identification and evaluation of novel potential farnesyltransferase inhibitors in cancer therapy

We have used various computational methodologies including molecular dynamics, density functional theory, virtual screening, ADMET predictions and molecular interaction field studies to design and analyze four novel potential inhibitors of farnesyltransferase (FTase). Evaluation of two proposals regarding their drug potential as well as lead compounds have indicated them as novel promising FTase inhibitors, with theoretically interesting pharmacotherapeutic profiles, when Compared to the very active and most cited FTase inhibitors that have activity data reported, which are launched drugs or compounds in clinical tests. One of our two proposals appears to be a more promising drug candidate and FTase inhibitor, but both derivative molecules indicate potentially very good pharmacotherapeutic profiles in comparison with Tipifarnib and Lonafarnib, two reference pharmaceuticals. Two other proposals have been selected with virtual screening approaches and investigated by LIS, which suggest novel and alternatives scaffolds to design future potential FTase inhibitors. Such compounds can be explored as promising molecules to initiate a research protocol in order to discover novel anticancer drug candidates targeting farnesyltransferase, in the fight against cancer. (C) 2009 Elsevier Inc. All rights reserved.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

N-mesityl-C-acylketenimines: 1,5-sigmatropic shifts and electrocyclization to quinolines

Flash vacuum thermolysis (FVT) of triazoles 6a-c generates alpha-oxoketenimines 10, the ester 10a being isolable. FVT of pyrroledione 8 generates the isomeric imidoylketene 9a. Ketenes 9 and ketenimines 10 undergo thermal interconversion by 1,3-shifts of methoxy and dimethylamino groups under mild FVT conditions (ca. 350-400 degrees C). Both 9 and 10 are directly observable by IR spectroscopy at either 77 K or on Ar matrix isolation at 12 K. On FVT at temperatures above ca. 400 degrees C, the ketenimines 10 undergo a 1,5-H shift to o-quinoid imines 12/13, followed by electrocyclization to dihydroquinolines 14 (unobserved) and 15 (observed by NMR). The latter are easily oxidized to alkylquinoline-3-carboxylates or quinoline-3-carboxamides 16 by atmospheric oxygen. Ab initio calculations on model compounds 18-23 predict an energy barrier of ca. 38 kcal mol(-1) (161 kJ mol(-1)) for the 1,5-H shift in N-(o-methylphenyl)ketenimines via the transition state TS19 followed by an electrocyclization barrier to dihydroquinoline 23a via TS22a of ca. 16 kcal mol(-1).

Determination of chemical shielding tensor of an indole carbon and application of tryptophan orientation of a membrane peptide

Using tryptophan C-13-enriched at the C-4 (C epsilon(3)) of the indole, the orientation of the C epsilon(3) chemical shift tensor relative to the C epsilon(3)-H dipolar axis was determined from the C-13 chemical shift/C-13-H-1 dipolar 2D NMR powder pattern. The principal values obtained were 208, 137 and 15 ppm with sigma(33) perpendicular to the indole plane, and sigma(11) (least shielded direction) 5 degrees off the C epsilon(3)-H bond toward C xi(3). The side off the C epsilon(3)-H bond was determined by comparing the reduced chemical shift anisotropies obtained by solid-state NMR and from molecular dynamics calculations of [4-C-13] tryptophans in gramicidin A aligned in phospholipid membranes. (C) 1999 Elsevier Science B.V. All rights reserved.