Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Unraveling Dynamical Heterogeneity in the Ionic Liquid 1-Butyl-3-methylimidazolium Chloride

Heterogeneous dynamics within a time range of nanoseconds was investigated by molecular dynamics (MD) simulations of 1 -butyl-3-methylimidazolium chloride ([bmim]Cl). After identifying groups of fast and slow ions, it was shown that the separation between the location of the center of mass and the center of charge of cations, d(CMCC), is a signature of such difference in ionic mobility. The distance d(CMCC) can be used as a signature because it relaxes in the time window of the dynamical heterogeneity. The relationship between the parameter dcmcc and conformations of the side alkyl chain in [bmim] is discussed. Since the relatively slow relaxation of dcmcc is a consequence of coexisting polar and nonpolar domains in the bulk, the MD simulations reveal a subtle interplay between structural and dynamical heterogeneity in ionic liquids.

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.

Defeitos pontuais nos compostos intermetálicos ZrNi e Zr/sub 2/Ni estudados por dinâmica molecular

Empregamos a técnica de Dinâmica Molecular para estudar propriedades de defeitos pontuais nos compostos intermetálicos ZrNi e Zr2Ni. Descrevemos as configurações estáveis de defeitos e mecanismos de migração, assim como as energias envolvidas. Os potenciais interatômicos foram derivados do Embedded Atom Model. No intuito de levar em conta a variação de estequiometria causada pela presença de alguns tipos de defeitos em intermetálicos, apresentamos um método numérico que fornece a energia efetiva de formação de defeitos e aplicamos o método ao ZrNi e Zr2Ni. Os resultados mostraram que vacâncias são mais estáveis na sub—rede do Ni, com energia de formação de 0,83 e-0,61 eV em ZrNi e Zr2Ni, respectivamente. Vacâncias de Zr são instáveis em ambos compostos; elas decaem espontaneamente em pares anti—sítio e vacância de Ni. Configurações e energias de formação de intersticiais também foram calculadas e mostraram comportamentos similares. Em ZrNi, a migração de vacâncias ocorre preferencialmente nas direções [025] e [100], com as respectivas energias de migração 0,67 e 0,73 eV, e é um processo essencialmente bidimensional no plano (001). Em Zr,Ni, a migração de vacâncias é unidimensional, ocorrendo na direção [001], com energia de migração de 0,67 eV. Em ambos compostos a presença de defeitos de anti—sítio de Ni diminui a energia de migração da vacância de Ni em até 3 vezes e facilita a movimentação em três dimensões. Mecanismos de anel não são energeticamente eficientes em comparação com saltos diretos. As configurações estáveis de intersticiais em ambos compostos consistem em um átomo de Ni sobre o plano (001) entre dois vizinhos de Zr fora do plano. Intersticiais de Zr são instáveis e tendem a deslocar um átomo de Ni, ocupando seu sítio. Energias de deslocamento foram estudadas através de simulações de irradiação de ambos compostos. Durante o processo de colisão binária, um potencial universal ZBL foi usado para colisões a curta distância. Para distâncias intermediárias usamos um potencial de união arbitrário. Zr mostrou—se mais difícil de ser arrancado de seu sítio do que Ni. Encontramos valores de energia de deslocamento no intervalo de aproximadamente 29 eV até 546 eV. Alguns resultados experimentais são mostrados e apresentam boa concordância com os cálculos.

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

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Obtenção de sistemas microemulsionados e estudo de simulação por dinâmica molecular de sistemas micelares objetivando a veiculação de produtos naturais bioativos

Among the new drugs launched into the market since 1980, up to 30% of them belong to the class of natural products or they have semisynthetic origin. Between 40-70% of the new chemical entities (or lead compounds) possess poor water solubility, which may impair their commercial use. An alternative for administration of poorly water-soluble drugs is their vehiculation into drug delivery systems like micelles, microemulsions, nanoparticles, liposomes, and cyclodextrin systems. In this work, microemulsion-based drug delivery systems were obtained using pharmaceutically acceptable components: a mixture Tween 80 and Span 20 in ratio 3:1 as surfactant, isopropyl mirystate or oleic acid as oil, bidistilled water, and ethanol, in some formulations, as cosurfactants. Self-Microemulsifying Drug Delivery Systems (SMEDDS) were also obtained using propylene glycol or sorbitol as cosurfactant. All formulations were characterized for rheological behavior, droplet size and electrical conductivity. The bioactive natural product trans-dehydrocrotonin, as well some extracts and fractions from Croton cajucara Benth (Euphorbiaceae), Anacardium occidentale L. (Anacardiaceae) e Phyllanthus amarus Schum. & Thonn. (Euphorbiaceae) specimens, were satisfactorily solubilized into microemulsions formulations. Meanwhile, two other natural products from Croton cajucara, trans-crotonin and acetyl aleuritolic acid, showed poor solubility in these formulations. The evaluation of the antioxidant capacity, by DPPH method, of plant extracts loaded into microemulsions evidenced the antioxidant activity of Phyllanthus amarus and Anacardium occidentale extracts. For Phyllanthus amarus extract, the use of microemulsions duplicated its antioxidant efficiency. A hydroalcoholic extract from Croton cajucara incorporated into a SMEDDS formulation showed bacteriostatic activity against colonies of Bacillus cereus and Escherichia coli bacteria. Additionally, Molecular Dynamics simulations were performed using micellar systems, for drug delivery systems, containing sugar-based surfactants, N-dodecylamino-1-deoxylactitol and N-dodecyl-D-lactosylamine. The computational simulations indicated that micellization process for N-dodecylamino-1- deoxylactitol is more favorable than N-dodecyl-D-lactosylamine system.

Cupin: A candidate molecular structure for the Nep1-like protein family

Background: NEP1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis. Some NLPs induce a hypersensitive-like response in dicot plants though the basis for this response remains unclear. In addition, the spatial structure and the role of these highly conserved proteins are not known.Results: We predict a 3d-structure for the beta-rich section of the NLPs based on alignments, prediction tools and molecular dynamics. We calculated a consensus sequence from 42 NLPs proteins, predicted its secondary structure and obtained a high quality alignment of this structure and conserved residues with the two Cupin superfamily motifs. The conserved sequence GHRHDWE and several common residues, especially some conserved histidines, in NLPs match closely the two cupin motifs. Besides other common residues shared by dicot Auxin-Binding Proteins (ABPs) and NLPs, an additional conserved histidine found in all dicot ABPs was also found in all NLPs at the same position.Conclusion: We propose that the necrosis inducing protein class belongs to the Cupin superfamily. Based on the 3d-structure, we are proposing some possible functions for the NLPs.

Molecular cloning and biochemical characterization of a myotoxin inhibitor from Bothrops alternatus snake plasma

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

How C-terminal carboxyamidation alters the biological activity of peptides from the venom of the eumenine solitary wasp

Inflammatory peptides display different types of post-transcriptional modifications, such as C-terminal amidation, that alter their biological activity. Here we describe the structural and molecular dynamics features of the mast cell degranulating peptide, eumenine mastoparan-AF (EMP-AF-NH2), found in the venom of the solitary wasp, and of its carboxyl-free C-terminal form (EMP-AF-COO-) characterized by a reduced activity. Circular dichroism indicates that both peptides switch from a random coil conformation in water to a helical structure in TFE and SDS micelles. NMR data, in 30% TFE, reveal that the two peptides fold into an alpha-helix spanning most of their length, while they differ in terms of molecular rigidity. To understand the origins of the conformational flexibility observed in the case of EMP-AF-COO-, a 5 ns MD simulation was carried out for each peptide, in an explicit water/TFE environment. The results show that the two peptides differ in an H-bond between Leu14 NH2 and the backbone carbonyl of Ile11. The loss of that H-bond in EMP-AF-COO- leads to a significant modification of its structural dynamics. In fact, as evidenced by essential dynamics analysis, while EMP-AF-NH2 exists mainly as a rigid structure, EMP-AF-COO- presents two helical stretches that fluctuate in some sort of independent fashion. We conclude that the diverse biological activity of the two peptides is not simply due to the reduction of the net positive charge, as generally suggested, but also to a structural perturbation of the amphipathic alpha-helix that affects their ability to perturb the cell membrane.

SMase II, a new sphingomyelinase D from Loxosceles laeta venom gland: Molecular cloning, expression, function and structural analysis

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

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.

Preparation, characterization and spectroscopy of the europium diphenylphosphinate complex

The purpose of this work is to study the preparation and spectroscopic behavior of the europium diphenylphosphinate complex -Eu(DPP)3. Elemental and thermogravimetric analysis, powder X-ray diffractometry, and infrared spectroscopy were applied to characterize the formula of the final product and the sixfold coordination of the Eu3+ ion. Excitation and emission spectra have been recorded at liquid nitrogen and room temperatures. The 5D0→7F2 transition intensity decreases when T decreases in comparison to the 5D0→7F1 transition intensity. Molecular mechanic calculations were developed in order to obtain the spatial coordinates of the Eu3+ and ligand ions. The simple overlap model was used to calculate the total splitting of the 5D0→7F1 transition, 5D0→7F0/5D 0→7F2 ntensity ratio and the intensity parameters, Ωλ (λ=2 and 4). Good agreements between theoretical predictions and experimental results have been obtained with g=2/3 as the effective charge and α=0.8×10-24 cm3 as the isotropic polarizability of the oxygen. © 1998 Elsevier Science S.A.

Steric constraints as folding coadjuvant

Through the analyses of the Miyazawa-Jernigan matrix it has been shown that the hydrophobic effect generates the dominant driving force for protein folding. By using both lattice and off-lattice models, it is shown that hydrophobic-type potentials are indeed efficient in inducing the chain through nativelike configurations, but they fail to provide sufficient stability so as to keep the chain in the native state. However, through comparative Monte Carlo simulations, it is shown that hydrophobic potentials and steric constraints are two basic ingredients for the folding process. Specifically, it is shown that suitable pairwise steric constraints introduce strong changes on the configurational activity, whose main consequence is a huge increase in the overall stability condition of the native state; detailed analysis of the effects of steric constraints on the heat capacity and configurational activity are provided. The present results support the view that the folding problem of globular proteins can be approached as a process in which the mechanism to reach the native conformation and the requirements for the globule stability are uncoupled.