Munir S. Skaf

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

Thiazolidinediones (TZDs) act through peroxisome proliferator activated receptor (PPAR) γ 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γ 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γ 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γ LBD, stronger partial agonists with full length PPARγ and exhibit full blockade of PPARγ phosphorylation by cyclin-dependent kinase 5 (cdk5), linked to reversal of adipose tissue insulin resistance. MCFAs that bind PPARγ 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/β-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γ modulators with useful clinical profiles among natural products.

Molecular association between water and dimethyl sulfoxide in solution: A molecular dynamics simulation study

A molecular dynamics simulation study of the local structures and H-bond distribution for water–dimethyl sulfoxide (DMSO) mixtures over the entire composition range is presented. Analysis of several site-site pair distribution functions reveals that two well-defined kinds of aggregates characterize the molecular association between water and DMSO in solution. One of them, already identified through recent neutron diffraction experiments and computer simulations, consists of two water molecules H-bonded to the oxygen atom of a DMSO molecule, such that the angle between the two H-bonds is nearly tetrahedral. The other complex features a central water molecule and two DMSO, making H-bonds to water hydrogens. According to the simulation data, these molecular aggregates coexist with each other in the mixture, but their proportions change with composition. 1DMSO-2water complexes predominate over 2DMSO-1water aggregates for water-rich mixtures (water mole fractions >50%), whereas the opposite is true for DMSO-ri...

The pattern of xylan acetylation suggests xylan may interact with cellulose microfibrils as a twofold helical screw in the secondary plant cell wall of Arabidopsis thaliana.

The interaction between xylan and cellulose microfibrils is important for secondary cell wall properties in vascular plants; however, the molecular arrangement of xylan in the cell wall and the nature of the molecular bonding between the polysaccharides are unknown. In dicots, the xylan backbone of β-(1,4)-linked xylosyl residues is decorated by occasional glucuronic acid, and approximately one-half of the xylosyl residues are O-acetylated at C-2 or C-3. We recently proposed that the even, periodic spacing of GlcA residues in the major domain of dicot xylan might allow the xylan backbone to fold as a twofold helical screw to facilitate alignment along, and stable interaction with, cellulose fibrils; however, such an interaction might be adversely impacted by random acetylation of the xylan backbone. Here, we investigated the arrangement of acetyl residues in Arabidopsis xylan using mass spectrometry and NMR. Alternate xylosyl residues along the backbone are acetylated. Using molecular dynamics simulation, we found that a twofold helical screw conformation of xylan is stable in interactions with both hydrophilic and hydrophobic cellulose faces. Tight docking of xylan on the hydrophilic faces is feasible only for xylan decorated on alternate residues and folded as a twofold helical screw. The findings suggest an explanation for the importance of acetylation for xylan–cellulose interactions, and also have implications for our understanding of cell wall molecular architecture and properties, and biological degradation by pathogens and fungi. They will also impact strategies to improve lignocellulose processing for biorefining and bioenergy.

GQ-16, a novel peroxisome proliferator-activated receptor γ (PPARγ) ligand, promotes insulin sensitization without weight gain.

Background: PPARγ agonists improve insulin sensitivity but also evoke weight gain. Results: GQ-16 is a PPARγ partial agonist that blocks receptor phosphorylation by Cdk5 and improves insulin sensitivity in diabetic mice in the absence of weight gain. Conclusion: The unique binding mode of GQ-16 appears to be responsible for the compounds advantageous pharmacological profile. Significance: Similar compounds could have promise as anti-diabetic therapeutics. The recent discovery that peroxisome proliferator-activated receptor γ (PPARγ) 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γ activation. Further characterization of this compound shows GQ-16 to be an effective inhibitor of Cdk5-mediated phosphorylation of PPARγ. The structure of GQ-16 bound to PPARγ demonstrates that the compound utilizes a binding mode distinct from other reported PPARγ 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 β-sheet region of the receptor, presumably explaining the compounds efficacy in inhibiting Cdk5-mediated phosphorylation of Ser-273. Molecular dynamics simulations suggest that the partial agonist activity of GQ-16 results from the compounds weak ability to stabilize helix 12 in its active conformation. Our results suggest that the emerging model, whereby “ideal” PPARγ-based therapeutics stabilize the β-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γ modulators that retain antidiabetic actions while minimizing untoward effects.

Cellulose-Builder: A toolkit for building crystalline structures of cellulose†

Cellulose‐builder is a user‐friendly program that builds crystalline structures of cellulose of different sizes and geometries. The program generates Cartesian coordinates for all atoms of the specified structure in the Protein Data Bank format, suitable for using as starting configurations in molecular dynamics simulations and other calculations. Crystalline structures of cellulose polymorphs Iα, Iβ, II, and IIII of practically any size are readily constructed which includes parallelepipeds, plant cell wall cellulose elementary fibrils of any length, and monolayers. Periodic boundary conditions along the crystallographic directions are easily imposed. The program also generates atom connectivity file in PSF format, required by well‐known simulation packages such as NAMD, CHARMM, and others. Cellulose‐builder is based on the Bash programming language and should run on practically any Unix‐like platform, demands very modest hardware, and is freely available for download from ftp://ftp.iqm.unicamp.br/pub/cellulose‐builder.

Plant Biomass Recalcitrance: Effect of Hemicellulose Composition on Nanoscale Forces that Control Cell Wall Strength

Efficient conversion of lignocellulosic biomass to second-generation biofuels and valuable chemicals requires decomposition of resilient plant cell wall structure. Cell wall recalcitrance varies among plant species and even phenotypes, depending on the chemical composition of the noncellulosic matrix. Changing the amount and composition of branches attached to the hemicellulose backbone can significantly alter the cell wall strength and microstructure. We address the effect of hemicellulose composition on primary cell wall assembly forces by using the 3D-RISM-KH molecular theory of solvation, which provides statistical-mechanical sampling and molecular picture of hemicellulose arrangement around cellulose. We show that hemicellulose branches of arabinose, glucuronic acid, and especially glucuronate strengthen the primary cell wall by strongly coordinating to hydrogen bond donor sites on the cellulose surface. We reveal molecular forces maintaining the cell wall structure and provide directions for genetic modulation of plants and pretreatment design to render biomass more amenable to processing.

Gaining ligand selectivity in thyroid hormone receptors via entropy.

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor β (TRβ) vs. TRα reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TRβ-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331β) in the TRβ ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3′-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TRβ selectivity. TR x-ray structures reveal better fit of ligand with the TRα LBC. The TRβ LBC, however, expands relative to TRα in the presence of Triac (549 Å3 vs. 461 Å3), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TRβ and permits greater flexibility of the Triac carboxylate group in TRβ than in TRα. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TRβ, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Mode of peroxisome proliferator-activated receptor gamma activation by luteolin.

The peroxisome proliferator-activated receptor γ (PPARγ) is a target for treatment of type II diabetes and other conditions. PPARγ 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γ but, unlike TZDs, does not promote adipocyte differentiation. However, previous reports suggested variously that luteolin is a PPARγ agonist or an antagonist. We show that luteolin exhibits weak partial agonist/antagonist activity in transfections, inhibits several PPARγ target genes in 3T3-L1 cells (LPL, ORL1, and CEBPα) and PPARγ-dependent adipogenesis, but activates GLUT4 to a similar degree as rosiglitazone, implying gene-specific partial agonism. The crystal structure of the PPARγ ligand-binding domain (LBD) reveals that luteolin occupies a buried ligand-binding pocket (LBP) but binds an inactive PPARγ LBD conformer and occupies a space near the β-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γ, and molecular dynamics simulations show that luteolin and myristic acid cooperate to stabilize the Ω-loop among H2′, H3, and the β-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γ-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γ modulators.

Solvation response of polar liquid mixtures: Water-dimethylsulfoxide

The solvation dynamics following the instantaneous creation of a positive or negative electronic charge in a previously neutral solute immersed in different water-dimethyl sulfoxide (DMSO) mixtures, spanning the entire composition range, is analyzed by molecular dynamics simulations. The solvation responses are strongly dependent on the sign of the solute charge, being considerably faster in the presence of cations for all mixtures considered. In terms of the composition dependence, the mixtures’ solvation response to the creation of the anion departs substantially from the pure solvents’, whereas for the cation, the mixtures’ responses are close to those exhibited by pure DMSO. In the case of anions, the mixture overall solvation time, defined as the time integral of the nonequilibrium response, can be as large as ten times the solvation time in pure DMSO, the slowest of the two cosolvents. The DMSO contribution to the mixtures’ solvation response may present an intriguing negative branch in the rotation...

Molecular dynamics simulations of dielectric properties of dimethyl sulfoxide: Comparison between available potentials

We present a molecular dynamics simulation study of the orientational correlations of molecular dipoles and dielectric properties of a series of recently proposed models for liquid dimethyl sulfoxide. The dipole–dipole correlations, as measured by the dipolar symmetry projections h110(r) and h112(r), are very similar for the models with four partially charged sites, but differ somewhat for the potential with no charges on the methyl groups. The effects of the charges on the methyl sites upon the local ordering of the molecular dipoles are discussed by comparing the h110(r) projections for two models which differ only in their site-charge assignments. The Kirkwood g factors for the simulated force fields turn out all greater than unity, in disagreement with some experimental estimates. The dielectric constants, however, are in good agreement with reported experimental data. In addition, we discuss the wave vector (k) dependence of the longitudinal static dielectric responses and also the single-particle mo...