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.

Pressure Perturbation Calorimetry and Guest-Host Chemistry of NDMG and N-MAP with p-Sulfonatocalix[6]arene

The aim of this project is to provide an explanation for recently obtained binding constants for two similar guest molecules, NDMG and N-MAP, with a p-sulfonatocalix[6]arene host in ammonium acetate buffer. This work was done primarily using pressure perturbation calorimetry, which is a technique that determines the coefficient of thermal expansion, α, which is in turn related to the solute molecule's effect on the order of the surrounding water molecules. A series of experiments were designed to test the effects of suspected confounding variables on the validity of PPC data. PPC was then used to study NDMG and N-MAP in ammonium acetate buffer. NDMG exhibited a minimum in α as function of temperature, while N-MAP did not. This difference was theorized to be due to the formation of an intramolecular hydrogen bond in monocationic NDMG that would lower the heat capacity of the molecule and better distribute the molecule's charge. Computational work and nuclear magnetic resonance spectroscopy confirmed that monocationic, ring-closed NDMG has less concentrated charge and more constrained motion than monocationic, ring-open NDMG. This evidence supports the theory that monocationic NDMG forms an intramolecular hydrogen bond and that this may be responsible for the minimum in α. This difference may explain the differences in binding constants between NDMG and N-MAP.

A tellurium-based cathepsin B inhibitor: Molecular structure, modelling, molecular docking and biological evaluation

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo teórico QTAIM e DFT dos compostos de coordenação: efeito quelato, titanocenos e ligação química

This work is a study of coordination compounds by quantum theory of atoms in molecules (QTAIM), based on the topological analysis of the electron density of molecular systems, both theoretically and experimentally obtained. The coordination chemistry topics which were studied are the chelate effect, bent titanocene and chemical bond in coordination complexes. The chelate effect was investigated according to topological and thermodynamic parameters. The exchange of monodentate ligands on polydentate ligands from same transition metal increases the stability of the complex both from entropy and enthalpy contributions. In some cases, the latter had a higher contribution to the stability of the complex in comparison with entropy. This enthalpic contribution is explained according to topological analysis of the M-ligand bonds where polidentate complex had higher values of electron density of bond critical point, Laplacian of electron density of bond critical point and delocalization index (number of shared electrons between two atoms). In the second chapter, was studied bent titanocenes with bulky cyclopentadienyl derivative π-ligand. The topological study showed the presence of secondary interactions between the atoms of π-ligands or between atoms of π-ligand and -ligand. It was found that, in the case of titanocenes with small difference in point group symmetry and with bulky ligands, there was an nearly linear relationship between stability and delocalization index involving the ring carbon atoms (Cp) and the titanium. However, the titanocene stability is not only related to the interaction between Ti and C atoms of Cp ring, but secondary interactions also play important role on the stability of voluminous titanocenes. The third chapter deals with the chemical bond in coordination compounds by means of QTAIM. The quantum theory of atoms in molecules so far classifies bonds and chemical interactions in two categories: closed shell interaction (ionic bond, hydrogen bond, van der Waals interaction, etc) and shared interaction (covalent bond). Based on topological parameters such as electron density, Laplacian of electron density, delocalization index, among others, was classified the chemical bond in coordination compounds as an intermediate between closed shell and shared interactions

Crystal structure of a myotoxic Asp49-phospholipase A(2) with low catalytic activity: Insights into Ca2+ -independent catalytic mechanism

A myotoxic Asp49-phospholipase A(2) (Asp49-PLA(2)) with low catalytic activity (BthTX-II from Bothrops jararacussu venom) was crystallized and the molecular-replacement solution has been obtained with a dimer in the asymmetric unit. The quaternary structure of BthTX-II resembles the myotoxic Asp49-PLA2 PrTX-III (piratoxin III from B. pirajai venom) and all non-catalytic and myotoxic dimeric Lys49-PLA(2)s. Despite of this, BthTX-II is different from the highly catalytic and non-myotoxic BthA-I (acidic PLA(2) from B. jararacussu) and other Asp49-PLA(2)s. BthTX-II structure showed a severe distortion of calcium-binding loop leading to displacement of the C-terminal region. Tyr28 side chain, present in this region, is in an opposite position in relation to the same residue in the catalytic activity Asp49-PLA(2)s, making a hydrogen bond with the atom 0 delta 2 of the catalytically active Asp49, which should coordinate the calcium. This high distortion may also be confirmed by the inability of BthTX-II to bind Na+ ions at the Ca2+-binding loop, despite of the crystallization to have occurred in the presence of this ion. In contrast, other Asp49-PLA(2)s which are able to bind Ca2+ ions are also able to bind Na+ ions at this loop. The comparison with other catalytic, non-catalytic and inhibited PLA(2)s indicates that the BthTX-II is not able to bind calcium ions; consequently, we suggest that its low catalytic function is based on an alternative way compared with other PLA(2)s. (c) 2008 Elsevier B.V All rights reserved.

Crystal structure of a phospholipase A(2) homolog complexed with p-bromophenacyl bromide reveals important structural changes associated with the inhibition of myotoxic activity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Interaction of shikimic acid with shikimate kinase (Retracted Article. See vol 334, pg 967, 2005)

The crystal structure of shikimate kinase from Mycobacterium tuberculosis (MtSK) complexed with MgADP and shikimic acid (shikimate) has been determined at 2.3 Angstrom resolution, clearly revealing the amino acid residues involved in shikimate binding. In MtSK, the Glu61 strictly conserved in SK forms a hydrogen bond and salt-bridge with Arg58 and assists in positioning the guanidinium group of Arg58 for shikimate binding. The carboxyl group of shikimate interacts with Arg58, Gly81, and Arg136, and hydroxyl groups with Asp34 and Gly80. The crystal structure of MtSK-MgADP-shikimate will provide crucial information for elucidation of the mechanism of SK-catalyzed reaction and for the development of a new generation of drugs against tuberculosis. (C) 2004 Elsevier B.V. All rights reserved.

Structure of shikimate kinase from Mycobacterium tuberculosis reveals the binding of shikimic acid

Tuberculosis made a resurgence in the mid-1980s and now kills approximately 3 million people a year. The re-emergence of tuberculosis as a public health threat, the high susceptibility of HIV-infected persons and the proliferation of multi-drug-resistant strains have created a need to develop new drugs. Shikimate kinase and other enzymes in the shikimate pathway are attractive targets for development of non-toxic antimicrobial agents, herbicides and anti-parasitic drugs, because the pathway is essential in these species whereas it is absent from mammals. The crystal structure of shikimate kinase from Mycobacterium tuberculosis (MtSK) complexed with MgADP and shikimic acid ( shikimate) has been determined at 2.3 Angstrom resolution, clearly revealing the amino-acid residues involved in shikimate binding. This is the first three-dimensional structure of shikimate kinase complexed with shikimate. In MtSK, the Glu61 residue that is strictly conserved in shikimate kinases forms a hydrogen bond and salt bridge with Arg58 and assists in positioning the guanidinium group of Arg58 for shikimate binding. The carboxyl group of shikimate interacts with Arg58, Gly81 and Arg136 and the hydroxyl groups interact with Asp34 and Gly80. The crystal structure of MtSK-MgADP-shikimate will provide crucial information for the elucidation of the mechanism of the shikimate kinase-catalyzed reaction and for the development of a new generation of drugs against tuberculosis.

Three-dimensional structure of ribonuclease T-1 complexed with an isosteric phosphonate substrate analogue of GpU: Alternate substrate binding modes and catalysis

The X-ray crystal structure of a complex between ribonuclease T-1 and guanylyl(3'-6')-6'-deoxyhomouridine (GpcU) has been determined at 2.0 Angstrom resolution. This Ligand is an isosteric analogue of the minimal RNA substrate, guanylyl(3'-5')uridine (GpU), where a methylene is substituted for the uridine 5'-oxygen atom. Two protein molecules are part of the asymmetric unit and both have a GpcU bound at the active site in the same manner. The protein-protein interface reveals an extended aromatic stack involving both guanines and three enzyme phenolic groups. A third GpcU has its guanine moiety stacked on His92 at the active site on enzyme molecule A and interacts with GpcU on molecule B in a neighboring unit via hydrogen bonding between uridine ribose 2'- and 3'-OH groups. None of the uridine moieties of the three GpcU molecules in the asymmetric unit interacts directly with the protein. GpcU-active-site interactions involve extensive hydrogen bonding of the guanine moiety at the primary recognition site and of the guanosine 2'-hydroxyl group with His40 and Glu58. on the other hand, the phosphonate group is weakly bound only by a single hydrogen bond with Tyr38, unlike ligand phosphate groups of other substrate analogues and 3'-GMP, which hydrogen-bonded with three additional active-site residues. Hydrogen bonding of the guanylyl 2'-OH group and the phosphonate moiety is essentially the same as that recently observed for a novel structure of a RNase T-1-3'-GMP complex obtained immediately after in situ hydrolysis of exo-(S-p)-guanosine 2',3'-cyclophosphorothioate [Zegers et al. (1998) Nature Struct. Biol. 5, 280-283]. It is likely that GpcU at the active site represents a nonproductive binding mode for GpU [:Steyaert, J., and Engleborghs (1995) fur. J. Biochem. 233, 140-144]. The results suggest that the active site of ribonuclease T-1 is adapted for optimal tight binding of both the guanylyl 2'-OH and phosphate groups (of GpU) only in the transition state for catalytic transesterification, which is stabilized by adjacent binding of the leaving nucleoside (U) group.

Crystal Structure Determination of Mebendazole Form A Using High-Resolution Synchrotron X-Ray Powder Diffraction Data

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structural aspects, thermal behavior, and stability of a self-assembled supramolecular polymer derived from flunixin-meglumine supramolecular adducts

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo ab-initio da a-alanina em meio aquoso

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Trichloro-bridged diruthenium (III,III) complexes: X-ray isomorphous structures of [Ph3X=O center dot center dot center dot H center dot center dot center dot O=XPh3][Ru2Cl7(XPh3)(2)]center dot 0.5(CH2Cl2)(H2O)(X = As or P)

The triply chloro-bridged binuclear complexes [Ph3X=O...H...O=XPh3][Ru2Cl7(XPh3)(2)].0.5(CH2Cl2) (H2O) (X = As or P) were obtained from [RuCl3(XPh3)(2)DMA].DMA (DMA = dimethylacetamide) CH2Cl2/Et2O solution. The structures were characterized by X-ray diffraction studies. The complexes are formed from two Ru atoms bridged by three chloride anions. The two ruthenium atoms are also coordinated to two non-bridging Cl atoms and an AsPh3 or PPh3 ligand respectively. As an interesting feature, the cations of these complexes are protons, trapped in a very short hydrogen bond between two triphenylarsine or triphenylphosphine oxide molecules.

Active and exo-site inhibition of human factor Xa: Structure of des-Gla factor Xa inhibited by NAP5, a potent nematode anticoagulant protein from Ancylostoma caninum

Hookworms are hematophagous nematodes capable of growth, development and subsistence in living host systems such as humans and other mammals. Approximately one billion, or one in six, people worldwide are infected by hookworms causing gastrointestinal blood loss and iron deficiency anemia. The hematophagous hookworm Ancylostoma caninum produces a family of small, disulfide-linked protein anticoagulants (75-84 amino acid residues). One of these nematode anticoagulant proteins, NAP5, inhibits the amidolytic activity of factor Xa (fXa) with K-i = 43 pM, and is the most potent natural fXa inhibitor identified thus far. The crystal structure of NAP5 bound at the active site of gamma-carboxyglutamic acid domainless factor Xa (des-fXa) has been determined at 3.1 angstrom resolution, which indicates that Asp189 (fXa, S1 subsite) binds to Arg40 (NAP5, P1 site) in a mode similar to that of the BPTI/trypsin interaction. However, the hydroxyl group of Ser39 of NAP5 additionally forms a hydrogen bond (2.5 angstrom) with His57 NE2 of the catalytic triad, replacing the hydrogen bond of Ser195 OG to the latter in the native structure, resulting in an interaction that has not been observed before. Furthermore, the C-terminal extension of NAP5 surprisingly interacts with the fXa exosite of a symmetry-equivalent molecule forming a short intermolecular beta-strand as observed in the structure of the NAPc2/fXa complex. This indicates that NAP5 can bind to fXa at the active site, or the exosite, and to fX at the exosite. However, unlike NAPc2, NAP5 does not inhibit fVIIa of the fVIIa/TF complex. (c) 2007 Elsevier Ltd. All rights reserved.

Application of a new reverse Monte Carlo algorithm to polyatomic molecular systems. I. Liquid water

Using a new reverse Monte Carlo algorithm, we present simulations that reproduce very well several structural and thermodynamic properties of liquid water. Both Monte Carlo, molecular dynamics simulations and experimental radial distribution functions used as input are accurately reproduced using a small number of molecules and no external constraints. Ad hoc energy and hydrogen bond analysis show the physical consistency and limitations of the generated RMC configurations. (C) 2001 American Institute of Physics.