Estudos computacionais de esfingomielinases D: docking, dinâmica molecular e métodos híbridos QM/MM

Pós-graduação em Biofísica Molecular - IBILCE

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

Employing perichromism for probing the properties of carboxymethyl cellulose films: an expedient, accurate method for the determination of the degree of substitution of the biopolymer derivative

The properties of films of carboxymethyl cellulose, CMC, of different degree of substitution, DS, have been examined by the use of perichromic indicators (probes). The film properties that have been determined are: empirical polarity, E-T(33); "acidity", alpha; "basicity", beta; and dipolarity/polarizability, pi*. This has been achieved by employing the following perichromic probes: 4-nitroaniline, 4-nitroanisole, 4-nitro-N,N-dimethylaniline, and 2,6-dichloro-4-(2,4,6-triphenyl-pyridinium-1-yl)phenolate, WB. The correlations between both E-T(33)- or pi* and DS were found to be linear; that between beta and DS is a second order polynomial; no obvious correlation was found between alpha and DS. The polarities of CMC films are in the range of those of butyl alcohols. As models for CMC, we have employed cellulose plus CMC of high DS; oxidized cellulose with degree of oxidation = 0.5; sodium glucuronate. The former model behaved akin to CMC, but the plots of the perichromic properties versus DS showed different slopes/intercepts. FTIR data and molecular dynamics simulations on the solvation of WB have shown that this difference can be traced to more efficient hydrogen bonding between the film of the model and the probe. This affects the intra-molecular charge-transfer energy of the latter, leading to different responses to the variation of DS. Based on the excellent linear correlation between E-T(33) and DS, for CMC from different origins, we suggest that perichromism is a simple, accurate, and expedient alternative for the determination of DS of the biopolymer derivative.

Directed and elliptic flow of charged particles in Cu plus Cu collisions at root s(NN)=22.4 GeV

This paper reports results for directed flow v(1) and elliptic flow v(2) of charged particles in Cu + Cu collisions at root s(NN) = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4-GeV Cu + Cu collisions the prior observation that v1 is independent of the system size at 62.4 and 200 GeV and also extend the scaling of v(1) with eta/y(beam) to this system. The measured v(2)(p(T)) in Cu + Cu collisions is similar for root s(NN) throughout the range 22.4 to 200 GeV. We also report a comparison with results from transport model (ultrarelativistic quantum molecular dynamics and multiphase transport model) calculations. The model results do not agree quantitatively with the measured v(1)(eta), v(2)(p(T)), and v(2)(eta).

Hb S-Sao Paulo: A new sickling hemoglobin with stable polymers and decreased oxygen affinity

Hb S-Sao Paulo (SP) [HBB:c.20A > T p.Glu6Val: c.196A > G p.Lys65Glu] is a new double-mutant hemoglobin that was found in heterozygosis in an 18-month-old Brazilian male with moderate anemia. It behaves like Hb S in acid electrophoresis, isoelectric focusing and solubility testing but shows different behavior in alkaline electrophoresis, cation-exchange HPLC and RP-HPLC. The variant is slightly unstable, showed reduced oxygen affinity and also appeared to form polymers more stable than the Hb S. Molecular dynamics simulation suggests that the polymerization is favored by interfacial electrostatic interactions. This provides a plausible explanation for some of the reported experimental observations. (C) 2012 Elsevier Inc. All rights reserved.

Septin C-Terminal Domain Interactions: Implications for Filament Stability and Assembly

Septins form a conserved family of filament forming GTP binding proteins found in a wide range of eukaryotic cells. They share a common structural architecture consisting of an N-terminal domain, a central GTP binding domain and a C-terminal domain, which is often predicted to adopt a coiled-coil conformation, at least in part. The crystal structure of the human SEPT2/SEPT6/SEPT7 heterocomplex has revealed the importance of the GTP binding domain in filament formation, but surprisingly no electron density was observed for the C-terminal domains and their function remains obscure. The dearth of structural information concerning the C-terminal region has motivated the present study in which the putative C-terminal domains of human SEPT2, SEPT6 and SEPT7 were expressed in E. coli and purified to homogeneity. The thermal stability and secondary structure content of the domains were studied by circular dichroism spectroscopy, and homo- and hetero-interactions were investigated by size exclusion chromatography, chemical cross-linking, analytical ultracentrifugation and surface plasmon resonance. Our results show that SEPT6-C and SEPT7-C are able to form both homo- and heterodimers with a high alpha-helical content in solution. The heterodimer is elongated and considerably more stable than the homodimers, with a K (D) of 15.8 nM. On the other hand, the homodimer SEPT2-C has a much lower affinity, with a K (D) of 4 mu M, and a moderate alpha-helical content. Our findings present the first direct experimental evidence toward better understanding the biophysical properties and coiled-coil pairings of such domains and their potential role in filament assembly and stability.

Interactions between 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and gamma-Al2O3 (100) surface calculated by density functional theory

The main aim of this work is to investigate the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4C1Im]+[BF4]-) ionic liquid (IL) adsorption on the gamma-Al2O3 (100) by density functional theory calculations to try to rationalize the adsorption as an electrostatic phenomenon. Optimized geometries and interaction energies of IL one-monolayer on the gamma-Al2O3 were obtained on high surface coverage (one cationanion pair per 94.96 nm2). A study of dispersion force was made to estimate its contribution to the adsorption. Overall, the process is ruled by electrostatic interaction between ions and surface. Adsorption of the anion [BF4]- and cation [C4C1Im]+ was also studied by Bader charge analysis and charge density difference for supported and unsupported situations. It is suggested that the IL ions have their charges maintained with significant anion cloud polarization inward to the acid aluminum sites. (c) 2012 Wiley Periodicals, Inc.

The Carmaphycins: New Proteasome Inhibitors Exhibiting an alpha,beta-Epoxyketone Warhead from a Marine Cyanobacterium

Two new peptidic proteasome inhibitors were isolated as trace components from a Curacao collection of the marine cyanobacterium Symploca sp. Carmaphycin A (1) and carmaphycin B (2) feature a leucine-derived a,beta-epoxyketone warhead directly connected to either methionine sulfoxide or methionine sulfone. Their structures were elucidated on the basis of extensive NMR and MS analyses and confirmed by total synthesis, which in turn provided more material for further biological evaluations. Pure carmaphycins A and B were found to inhibit the beta 5 subunit (chymotrypsin-like activity) of the S. cerevisiae 20S proteasome in the low nanomolar range. Additionally, they exhibited strong cytotoxicity to lung and colon cancer cell lines, as well as exquisite antiproliferative effects in the NCI60 cell-line panel. These assay results as well as initial structural biology studies suggest a distinctive binding mode for these new inhibitors.

Medium Chain Fatty Acids Are Selective Peroxisome Proliferator Activated Receptor (PPAR) gamma Activators and Pan-PPAR Partial Agonists

Thiazolidinediones (TZDs) act through peroxisome proliferator activated receptor (PPAR) gamma to increase insulin sensitivity in type 2 diabetes (T2DM), but deleterious effects of these ligands mean that selective modulators with improved clinical profiles are needed. We obtained a crystal structure of PPAR gamma ligand binding domain (LBD) and found that the ligand binding pocket (LBP) is occupied by bacterial medium chain fatty acids (MCFAs). We verified that MCFAs (C8-C10) bind the PPAR gamma LBD in vitro and showed that they are low-potency partial agonists that display assay-specific actions relative to TZDs; they act as very weak partial agonists in transfections with PPAR gamma LBD, stronger partial agonists with full length PPAR gamma and exhibit full blockade of PPAR gamma phosphorylation by cyclin-dependent kinase 5 (cdk5), linked to reversal of adipose tissue insulin resistance. MCFAs that bind PPAR gamma also antagonize TZD-dependent adipogenesis in vitro. X-ray structure B-factor analysis and molecular dynamics (MD) simulations suggest that MCFAs weakly stabilize C-terminal activation helix (H) 12 relative to TZDs and this effect is highly dependent on chain length. By contrast, MCFAs preferentially stabilize the H2-H3/beta-sheet region and the helix (H) 11-H12 loop relative to TZDs and we propose that MCFA assay-specific actions are linked to their unique binding mode and suggest that it may be possible to identify selective PPAR gamma modulators with useful clinical profiles among natural products.

GQ-16, a Novel Peroxisome Proliferator-activated Receptor gamma (PPAR gamma) Ligand, Promotes Insulin Sensitization without Weight Gain

The recent discovery that peroxisome proliferator-activated receptor gamma (PPAR gamma) targeted anti-diabetic drugs function by inhibiting Cdk5-mediated phosphorylation of the receptor has provided a new viewpoint to evaluate and perhaps develop improved insulin-sensitizing agents. Herein we report the development of a novel thiazolidinedione that retains similar anti-diabetic efficacy as rosiglitazone in mice yet does not elicit weight gain or edema, common side effects associated with full PPAR gamma activation. Further characterization of this compound shows GQ-16 to be an effective inhibitor of Cdk5-mediated phosphorylation of PPAR gamma. The structure of GQ-16 bound to PPAR gamma demonstrates that the compound utilizes a binding mode distinct from other reported PPAR gamma ligands, although it does share some structural features with other partial agonists, such as MRL-24 and PA-082, that have similarly been reported to dissociate insulin sensitization from weight gain. Hydrogen/deuterium exchange studies reveal that GQ-16 strongly stabilizes the beta-sheet region of the receptor, presumably explaining the compound's efficacy in inhibiting Cdk5-mediated phosphorylation of Ser-273. Molecular dynamics simulations suggest that the partial agonist activity of GQ-16 results from the compound's weak ability to stabilize helix 12 in its active conformation. Our results suggest that the emerging model, whereby "ideal" PPAR gamma-based therapeutics stabilize the beta-sheet/Ser-273 region and inhibit Cdk5-mediated phosphorylation while minimally invoking adipogenesis and classical agonism, is indeed a valid framework to develop improved PPAR gamma modulators that retain antidiabetic actions while minimizing untoward effects.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

The Role of Molecular Conformation and Polarizable Embedding for One- and Two-Photon Absorption of Disperse Orange 3 in Solution

Solvent effects on the one- and two-photon absorption (IPA and 2PA) of disperse orange 3 (DO3) in dimethyl sulfoxide (DMSO) are studied using a discrete polarizable embedding (PE) response theory. The scheme comprises a quantum region containing the chromophore and an atomically granulated classical region for the solvent accounting for full interactions within and between the two regions. Either classical molecular dynamics (MD) or hybrid Car-Parrinello (CP) quantum/classical (QM/MM) molecular dynamics simulations are employed to describe the solvation of DO3 in DMSO, allowing for an analysis of the effect of the intermolecular short-range repulsion, long-range attraction, and electrostatic interactions on the conformational changes of the chromophore and also the effect of the solute-solvent polarization. PE linear response calculations are performed to verify the character, solvatochromic shift, and overlap of the two lowest energy transitions responsible for the linear absorption spectrum of DO3 in DMSO in the visible spectral region. Results of the PE linear and quadratic response calculations, performed using uncorrelated solute-solvent configurations sampled from either the classical or hybrid CP QM/MM MD simulations, are used to estimate the width of the line shape function of the two electronic lowest energy excited states, which allow a prediction of the 2PA cross-sections without the use of empirical parameters. Appropriate exchange-correlation functionals have been employed in order to describe the charge-transfer process following the electronic transitions of the chromophore in solution.

Correlations between Conjugation Length, Macromolecular Dynamics, and Photophysics of Phenylene-Vinylene/Aliphatic Multiblock Copolymers

This work reports a detailed spectroscopy study of a series of multiblock conjugated nonconjugated copolymers built by p-phenylene vinylene type units (PV) and octamethylene spacers, namely, poly(1,8-octanedioxy-2,6-dimethoxy-1,4-phenylene-1,2-ethenylene) (LaPPS18). The relative proportions of the PV and aliphatic segments were estimated on the basis of solid-state NMR and Raman spectroscopy. The overall structure was characterized by wide angle X-ray diffraction; H-1 wide-line dipolar chemical shift correlation (DIPSHIFT), and centerband-only detection of exchange (CODEX) NMR data, that together with glass transition temperatures allowed us to identify the groups involved in the molecular dynamics. These different structural properties were used to explain the photoluminescence properties in terms of peak position and spectral profile

Mode of Peroxisome Proliferator-Activated Receptor gamma Activation by Luteolin

The peroxisome proliferator-activated receptor gamma (PPAR gamma) is a target for treatment of type II diabetes and other conditions. PPAR gamma full agonists, such as thiazolidinediones (TZDs), are effective insulin sensitizers and anti-inflammatory agents, but their use is limited by adverse side effects. Luteolin is a flavonoid with anti-inflammatory actions that binds PPAR gamma but, unlike TZDs, does not promote adipocyte differentiation. However, previous reports suggested variously that luteolin is a PPAR gamma agonist or an antagonist. We show that luteolin exhibits weak partial agonist/antagonist activity in transfections, inhibits several PPAR gamma target genes in 3T3-L1 cells (LPL, ORL1, and CEBP alpha) and PPAR gamma-dependent adipogenesis, but activates GLUT4 to a similar degree as rosiglitazone, implying gene-specific partial agonism. The crystal structure of the PPAR gamma ligand-binding domain (LBD) reveals that luteolin occupies a buried ligand-binding pocket (LBP) but binds an inactive PPAR gamma LBD conformer and occupies a space near the beta-sheet region far from the activation helix (H12), consistent with partial agonist/antagonist actions. A single myristic acid molecule simultaneously binds the LBP, suggesting that luteolin may cooperate with other ligands to bind PPAR gamma, and molecular dynamics simulations show that luteolin and myristic acid cooperate to stabilize the Omega-loop among H2', H3, and the beta-sheet region. It is noteworthy that luteolin strongly suppresses hypertonicity-induced release of the pro-inflammatory interleukin-8 from human corneal epithelial cells and reverses reductions in transepithelial electrical resistance. This effect is PPAR gamma-dependent. We propose that activities of luteolin are related to its singular binding mode, that anti-inflammatory activity does not require H12 stabilization, and that our structure can be useful in developing safe selective PPAR gamma modulators.

Platinum-Based Nanoalloys PtnTM55-n (TM=Co, Rh, Au): A Density Functional Theory Investigation

Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.