Paulo C. T. Souza

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 mechanism of peroxisome proliferator-activated receptor α activation by WY14643 : a new mode of ligand recognition and receptor stabilization

Peroxisome proliferator-activated receptors (PPARs) are members of a superfamily of nuclear transcription factors. They are involved in mediating numerous physiological effects in humans, including glucose and lipid metabolism. PPARα ligands effectively treat dyslipidemia and have significant antiinflammatory and anti-atherosclerotic activities. These effects and their ligand-dependent activity make nuclear receptors obvious targets for drug design. Here, we present the structure of the human PPARα in complex with WY14643, a member of fibrate class of drug, and a widely used PPAR activator. The crystal structure of this complex suggests that WY14643 induces activation of PPARα in an unusual bipartite mechanism involving conventional direct helix 12 stabilization and an alternative mode that involves a second ligand in the pocket. We present structural observations, molecular dynamics and activity assays that support the importance of the second site in WY14643 action. The unique binding mode of WY14643 reveals a new pattern of nuclear receptor ligand recognition and suggests a novel basis for ligand design, offering clues for improving the binding affinity and selectivity of ligand. We show that binding of WY14643 to PPARα was associated with antiinflammatory disease in a human corneal cell model, suggesting possible applications for PPARα ligands.

CHARMM Force Field Parameterization of Peroxisome Proliferator-Activated Receptor γ Ligands

The peroxisome proliferator-activated receptor γ (PPARγ) ligands are important therapeutic drugs for the treatment of type 2 diabetes, obesity and cardiovascular diseases. In particular, partial agonists and non-agonists are interesting targets to reduce glucose levels, presenting few side effects in comparison to full agonists. In this work, we present a set of CHARMM-based parameters of a molecular mechanics force field for two PPARγ ligands, GQ16 and SR1664. GQ16 belongs to the thiazolidinedione class of drugs and it is a PPARγ partial agonist that has been shown to promote the “browning” of white adipose tissue. SR1664 is the precursor of the PPARγ non-agonist class of ligands that activates PPARγ in a non-classical manner. Here, we use quantum chemical calculations consistent with the CHARMM protocol to obtain bonded and non-bonded parameters, including partial atomic charges and effective torsion potentials for both molecules. The newly parameterized models were evaluated by examining the behavior of GQ16 and SR1664 free in water and bound to the ligand binding pocket of PPARγ using molecular dynamics simulations. The potential parameters derived here are readily transferable to a variety of pharmaceutical compounds and similar PPARγ ligands.

Identification of a New Hormone-Binding Site on the Surface of Thyroid Hormone Receptor

Thyroid hormone receptors (TRs) are members of the nuclear receptor superfamily of ligand-activated transcription factors involved in cell differentiation, growth, and homeostasis. Although X-ray structures of many nuclear receptor ligand-binding domains (LBDs) reveal that the ligand binds within the hydrophobic core of the ligand-binding pocket, a few studies suggest the possibility of ligands binding to other sites. Here, we report a new x-ray crystallographic structure of TR-LBD that shows a second binding site for T3 and T4 located between H9, H10, and H11 of the TRα LBD surface. Statistical multiple sequence analysis, site-directed mutagenesis, and cell transactivation assays indicate that residues of the second binding site could be important for the TR function. We also conducted molecular dynamics simulations to investigate ligand mobility and ligand-protein interaction for T3 and T4 bound to this new TR surface-binding site. Extensive molecular dynamics simulations designed to compute ligand-protein dissociation constant indicate that the binding affinities to this surface site are of the order of the plasma and intracellular concentrations of the thyroid hormones, suggesting that ligands may bind to this new binding site under physiological conditions. Therefore, the second binding site could be useful as a new target site for drug design and could modulate selectively TR functions.

Molecular Recognition of PPARγ by Kinase Cdk5/p25: Insights from a Combination of Protein–Protein Docking and Adaptive Biasing Force Simulations

The peroxisome proliferator-activated receptor γ (PPARγ) is an important transcription factor that plays a major role in the regulation of glucose and lipid metabolisms and has, therefore, many implications in modern-life metabolic disorders such as diabetes, obesity, and cardiovascular diseases. Phosphorylation of PPARγ by the cyclin-dependent kinase 5 (Cdk5) has been recently proved to promote obesity and loss of insulin sensitivity. The inhibition of this reaction is currently being pursued to develop PPARγ ligands for type 2 diabetes treatments. The knowledge of the protein-protein interactions between Cdk5/p25 and PPARγ can be an important asset for better understanding of the molecular basis of the Cdk5-meditated phosphorylation of PPARγ and its inhibition. By means of a computational approach that combines protein-protein docking and adaptive biasing force molecular dynamics simulations, we obtained PPARγ-Cdk5/p25 structural models that are consistent with the mechanism of the enzymatic reaction and with overall structural features of the full length PPARγ-RXRα heterodimer bound to DNA. In addition to the active site, our model shows that the interacting regions between the two proteins should involve two distal docking sites, comprising the PPARγ Ω-loop and Cdk5 N-terminal lobe and the PPARγ β-sheet and Cdk5 C-terminal lobe. These sites are related to PPARγ transactivation and directly interact with PPARγ ligands. Our results suggest that β-sheets and Ω-loop stabilization promoted by PPARγ agonists could be important to inhibit Cdk5-mediated phosphorylation.

On the denaturation mechanisms of the ligand binding domain of thyroid hormone receptors

The ligand binding domain (LBD) of nuclear hormone receptors adopts a very compact, mostly alpha-helical structure that binds specific ligands with very high affinity. We use circular dichroism spectroscopy and high-temperature molecular dynamics simulations to investigate unfolding of the LBDs of thyroid hormone receptors (TRs). A molecular description of the denaturation mechanisms is obtained by molecular dynamics simulations of the TRalpha and TRbeta LBDs in the absence and in the presence of the natural ligand Triac. The simulations show that the thermal unfolding of the LBD starts with the loss of native contacts and secondary structure elements, while the structure remains essentially compact, resembling a molten globule state. This differs from most protein denaturation simulations reported to date and suggests that the folding mechanism may start with the hydrophobic collapse of the TR LBDs. Our results reveal that the stabilities of the LBDs of the TRalpha and TRbeta subtypes are affected to different degrees by the binding of the isoform selective ligand Triac and that ligand binding confers protection against thermal denaturation and unfolding in a subtype specific manner. Our simulations indicate two mechanisms by which the ligand stabilizes the LBD: (1) by enhancing the interactions between H8 and H11, and the interaction of the region between H1 and the Omega-loop with the core of the LBD, and (2) by shielding the hydrophobic H6 from hydration.

An alternative conformation of ERβ bound to estradiol reveals H12 in a stable antagonist position

The natural ligand 17β-estradiol (E2) is so far believed to induce a unique agonist-bound active conformation in the ligand binding domain (LBD) of the estrogen receptors (ERs). Both subtypes, ERα and ERβ, are transcriptionally activated in the presence of E2 with ERβ being somewhat less active than ERα under similar conditions. The molecular bases for this intriguing behavior are mainly attributed to subtype differences in the amino-terminal domain of these receptors. However, structural details that confer differences in the molecular response of ER LBDs to E2 still remain elusive. In this study, we present a new crystallographic structure of the ERβ LBD bound to E2 in which H12 assumes an alternative conformation that resembles antagonist ERs structures. Structural observations and molecular dynamics simulations jointly provide evidence that alternative ERβ H12 position could correspond to a stable conformation of the receptor under physiological pH conditions. Our findings shed light on the unexpected role of LBD in the lower functional response of ERβ subtype.

High performance collision cross section calculation—HPCCS

Since the commercial introduction of Ion Mobility coupled with Mass Spectrometry (IM‐MS) devices in 2003, a large number of research laboratories have embraced the technique. IM‐MS is a fairly rapid experiment used as a molecular separation tool and to obtain structural information. The interpretation of IM‐MS data is still challenging and relies heavily on theoretical calculations of the molecules collision cross section (CCS) against a buffer gas. Here, a new software (HPCCS) is presented, which performs CCS calculations using high perfomance computing techniques. Based on the trajectory method, HPCCS can accurately calculate CCS for a great variety of molecules, ranging from small organic molecules to large protein complexes, using helium or nitrogen as buffer gas with considerable gains in computer time compared to publicly available codes under the same level of theory. HPCCS is available as free software under the Academic Use License at https://github.com/cepid-cces/hpccs.