Binding Mechanism between Platelet Glycoprotein and Cyclic Peptide Elucidated by McMD-Based Dynamic Docking.

The cyclic peptide OS1 (amino acid sequence: CTERMALHNLC), which has a disulfide bond between both termini cysteine residues, inhibits complex formation between the platelet glycoprotein Ibα (GPIbα) and the von Willebrand factor (vWF) by forming a complex with GPIbα. To study the binding mechanism between GPIbα and OS1 and, therefore, the inhibition mechanism of the protein-protein GPIbα-vWF complex, we have applied our multicanonical molecular dynamics (McMD)-based dynamic docking protocol starting from the unbound state of the peptide. Our simulations have reproduced the experimental complex structure, although the top-ranking structure was an intermediary one, where the peptide was bound in the same location as in the experimental structure; however, the ß-switch of GPIbα attained a different conformation. Our analysis showed that subsequent refolding of the ß-switch results in a more stable binding configuration, although the transition to the native configuration appears to take some time, during which OS1 could dissociate. Our results show that conformational changes in the ß-switch are crucial for successful binding of OS1. Furthermore, we identified several allosteric binding sites of GPIbα that might also interfere with vWF binding, and optimization of the peptide to target these allosteric sites might lead to a more effective inhibitor, as these are not dependent on the ß-switch conformation.

Elucidation of binding mechanism, affinity, and complex structure between mWT1 tumor-associated antigen peptide and HLA-A*24:02.

We have applied our advanced computational and experimental methodologies to investigate the complex structure and binding mechanism of a modified Wilms' Tumor 1 (mWT1) protein epitope to the understudied Asian-dominant allele HLA-A*24:02 (HLA-A24) in aqueous solution. We have applied our developed multicanonical molecular dynamics (McMD)-based dynamic docking method to analyze the binding pathway and mechanism, which we verified by comparing the highest probability structures from simulation with our experimentally solved x-ray crystal structure. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that first an encounter complex is formed between the N-terminal's positive charge of the 9-residue mWT1 fragment peptide and a cluster of negative residues on the surface of HLA-A24, with the major histocompatibility complex (MHC) molecule preferring a predominantly closed conformation. The peptide first binds to this closed MHC conformation, forming an encounter complex, after which the binding site opens due to increased entropy of the binding site, allowing the peptide to bind to form the native complex structure. Further sequence and structure analyses also suggest that although the peptide loading complex would help with stabilizing the MHC molecule, the binding depends in a large part on the intrinsic affinity between the MHC molecule and the antigen peptide. Finally, our computational tools and analyses can be of great benefit to study the binding mechanism of different MHC types to their antigens, where it could also be useful in the development of higher affinity variant peptides and for personalized medicine.

Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL.

The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.

N-Terminal-Driven Binding Mechanism of an Antigen Peptide to Human Leukocyte Antigen-A*2402 Elucidated by Multicanonical Molecular Dynamic-Based Dynamic Docking and Path Sampling Simulations.

We have applied our advanced multicanonical molecular dynamics (McMD)-based dynamic docking methodology to investigate the binding mechanism of an HIV-1 Nef protein epitope to the Asian-dominant allele human leukocyte antigen (HLA)-A*2402. Even though pMHC complex formation [between a Major histocompatibility complex (MHC) class I molecule, which is encoded by an HLA allele, and an antigen peptide] is one of the fundamental processes of the adaptive human immune response, its binding mechanism has not yet been well studied, partially due to the high allelic variation of HLAs in the population. We have used our developed McMD-based dynamic docking method and have successfully reproduced the native complex structure, which is located near the free energy global minimum. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that the peptide binding is initially driven by the highly positively charged N-terminus of the peptide that is attracted to the various negatively charged residues on the MHC molecule's surface. Upon nearing the pocket, the second tyrosine residue of the peptide anchors the peptide by strongly binding to the B-site of the MHC molecule via hydrophobic driven interactions, resulting in a very strong bound complex structure. Our methodology can be effectively used to predict the bound complex structures between MHC molecules and their antigens to study their binding mechanism in close detail, which would help with the development of new vaccines against cancers, as well as viral infections such as HIV and COVID-19.

Flexibility and Cell Permeability of Cyclic Ras-Inhibitor Peptides Revealed by the Coupled Nosé-Hoover Equation.

Quantifying the cell permeability of cyclic peptides is crucial for their rational drug design. However, the reasons remain unclear why a minor chemical modification, such as the difference between Ras inhibitors cyclorasin 9A5 and 9A54, can substantially change a peptide's permeability. To address this question, we performed enhanced sampling simulations of these two 11-mer peptides using the coupled Nosé-Hoover equation (cNH) we recently developed. The present cNH simulations realized temperature fluctuations over a wide range (240-600 K) in a dynamic manner, allowing structural samplings that were well validated by nuclear Overhauser effect measurements. The derived structural ensembles were comprehensively analyzed by all-atom structural clustering, mapping the derived clusters onto principal components (PCs) that characterize the cyclic structure, and calculating cluster-dependent geometric and chemical properties. The planar-open conformation was dominant in aqueous solvent, owing to inclusion of the Trp side chain in the main-chain ring, while the compact-closed conformation, which favors cell permeation due to its compactness and high polarity, was also accessible. Conformation-dependent cell permeability was observed in one of the derived PCs, demonstrating that decreased cell permeability in 9A54 is due to the high free energy barrier separating the two conformations. The origin of the change in free energy surface was determined to be loss of flexibility in the modified residues 2-3, resulting from the increased bulkiness of their side chains. The derived molecular mechanism of cell permeability highlights the significance of complete structural dynamics surveys for accelerating drug development with cyclic peptides.

Dynamic Docking Using Multicanonical Molecular Dynamics: Simulating Complex Formation at the Atomistic Level.

Multicanonical molecular dynamics (McMD)-based dynamic docking has been applied to predict the native binding configurations for several protein receptors and their ligands. Due to the enhanced sampling capabilities of McMD, it can exhaustively sample bound and unbound ligand configurations, as well as receptor conformations, and thus enables efficient sampling of the conformational and configurational space, not possible using canonical MD simulations. As McMD samples a wide configurational space, extensive analysis is required to study the diverse ensemble consisting of bound and unbound structures. By projecting the reweighted ensemble onto the first two principal axes obtained via principal component analysis of the multicanonical ensemble, the free energy landscape (FEL) can be obtained. Further analysis produces representative structures positioned at the local minima of the FEL, where these structures are then ranked by their free energy. In this chapter, we describe our dynamic docking methodology, which has successfully reproduced the native binding configuration for small compounds, medium-sized compounds, and peptide molecules.

Cryptic-site binding mechanism of medium-sized Bcl-xL inhibiting compounds elucidated by McMD-based dynamic docking simulations.

We have performed multicanonical molecular dynamics (McMD) based dynamic docking simulations to study and compare the binding mechanism between two medium-sized inhibitors (ABT-737 and WEHI-539) that bind to the cryptic site of Bcl-xL, by exhaustively sampling the conformational and configurational space. Cryptic sites are binding pockets that are transiently formed in the apo state or are induced upon ligand binding. Bcl-xL, a pro-survival protein involved in cancer progression, is known to have a cryptic site, whereby the shape of the pocket depends on which ligand is bound to it. Starting from the apo-structure, we have performed two independent McMD-based dynamic docking simulations for each ligand, and were able to obtain near-native complex structures in both cases. In addition, we have also studied their interactions along their respective binding pathways by using path sampling simulations, which showed that the ligands form stable binding configurations via predominantly hydrophobic interactions. Although the protein started from the apo state, both ligands modulated the pocket in different ways, shifting the conformational preference of the sub-pockets of Bcl-xL. We demonstrate that McMD-based dynamic docking is a powerful tool that can be effectively used to study binding mechanisms involving a cryptic site, where ligand binding requires a large conformational change in the protein to occur.

Structural basis of mutants of PET-degrading enzyme from Saccharomonospora viridis AHK190 with high activity and thermal stability.

The cutinase-like enzyme from the thermophile Saccharomonospora viridis AHK190, Cut190, is a good candidate to depolymerize polyethylene terephthalate (PET) efficiently. We previously developed a mutant of Cut190 (S226P/R228S), which we designated as Cut190* that has both increased activity and stability and solved its crystal structure. Recently, we showed that mutation of D250C/E296C on one of the Ca2+ -binding sites resulted in a higher thermal stability while retaining its polyesterase activity. In this study, we solved the crystal structures of Cut190* mutants, Q138A/D250C-E296C/Q123H/N202H, designated as Cut190*SS, and its inactive S176A mutant, Cut190*SS_S176A, at high resolution. The overall structures were similar to those of Cut190* and Cut190*S176A reported previously. As expected, Cys250 and Cys296 were closely located to form a disulfide bond, which would assuredly contribute to increase the stability. Isothermal titration calorimetry experiments and 3D Reference Interaction Site Model calculations showed that the metal-binding properties of the Cut190*SS series were different from those of the Cut190* series. However, our results show that binding of Ca2+ to the weak binding site, site 1, would be retained, enabling Cut190*SS to keep its ability to use Ca2+ to accelerate the conformational change from the closed (inactive) to the open (active) form. While increasing the thermal stability, Cut190*SS could still express its enzymatic function. Even after incubation at 70°C, which corresponds to the glass transition temperature of PET, the enzyme retained its activity well, implying a high applicability for industrial PET depolymerization using Cut190*SS.

Molecular Interaction Mechanism of a 14-3-3 Protein with a Phosphorylated Peptide Elucidated by Enhanced Conformational Sampling.

Enhanced conformational sampling, a genetic-algorithm-guided multidimensional virtual-system coupled molecular dynamics, can provide equilibrated conformational distributions of a receptor protein and a flexible ligand at room temperature. The distributions provide not only the most stable but also semistable complex structures and propose a ligand-receptor binding process. This method was applied to a system consisting of a receptor protein, 14-3-3ε, and a flexible peptide, phosphorylated myeloid leukemia factor 1 (pMLF1). The results present comprehensive binding pathways of pMLF1 to 14-3-3ε. We identified four thermodynamically stable clusters of MLF1 on the 14-3-3ε surface and free-energy barriers among some clusters. The most stable cluster includes two high-density spots connected by a narrow corridor. When pMLF1 passes the corridor, a salt-bridge relay (switching) related to the phosphorylated residue of pMLF1 occurs. Conformations in one high-density spot are similar to the experimentally determined complex structure. Three-dimensional distributions of residues in the intermolecular interface rationally explain the binding constant changes resulting from the alanine mutation experiment for the residues. We also performed a simulation of nonphosphorylated peptide and 14-3-3ε, which demonstrated that the complex structure was unstable, suggesting that phosphorylation of the peptide is crucially important for binding to 14-3-3ε.

Accurate Prediction of Complex Structure and Affinity for a Flexible Protein Receptor and Its Inhibitor.

In order to predict the accurate binding configuration as well as the binding affinity for a flexible protein receptor and its inhibitor drug, enhanced sampling with multicanonical molecular dynamics (McMD) simulation and thermodynamic integration (TI) were combined as a general drug docking method. CDK2, cyclin-dependent kinase 2, is involved in the cell cycle regulation. Malfunctions in CDK2 can cause tumorigenesis, and thus it is a potential drug target. Here, we performed a long McMD simulation for docking the inhibitor CS3 to CDK2 starting from the unbound structure. Subsequently, a potential binding/unbinding pathway was given from the multicanonical ensemble, and the binding free energy was readily computed by TI along the pathway. Using this combination, the correct binding configuration of CS3 to CDK2 was obtained, and its affinity coincided well with the experimental value.

The Effect of Conformational Flexibility on Binding Free Energy Estimation between Kinases and Their Inhibitors.

Accurate prediction of binding affinities of drug candidates to their targets remains challenging because of protein flexibility in solution. Conformational flexibility of the ATP-binding site in the CDK2 and ERK2 kinases was identified using molecular dynamics simulations. The binding free energy (ΔG) of twenty-four ATP-competitive inhibitors toward these kinases was assessed using an alchemical free energy perturbation method, MP-CAFEE. However, large calculation errors of 2-3 kcal/mol were observed using this method, where the free energy simulation starts from a single equilibrated conformation. Here, we developed a new ΔG computation method, where the starting structure was set to multiconformations to cover flexibility. The calculation accuracy was successfully improved, especially for larger molecular size compounds, leading to reliable prediction of a broader range of drug candidates. The present study demonstrates that conformational flexibility of interactions between a compound and the glycine-rich loop in the kinases is a key factor in ΔG estimation.

Elastic properties of dynein motor domain obtained from all-atom molecular dynamics simulations.

Dyneins are large microtubule motor proteins that convert ATP energy to mechanical power. High-resolution crystal structures of ADP-bound cytoplasmic dynein have revealed the organization of the motor domain, comprising the AAA(+) ring, the linker, the stalk/strut and the C sequence. Recently, the ADP.vanadate-bound structure, which is similar to the ATP hydrolysis transition state, revealed how the structure of dynein changes upon ATP binding. Although both the ADP- and ATP-bound state structures have been resolved, the dynamic properties at the atomic level remain unclear. In this work, we built two models named 'the ADP model' and 'the ATP model', where ADP and ATP are bound to AAA1 in the AAA(+) ring, respectively, to observe the initial procedure of the structural change from the unprimed to the primed state. We performed 200-ns molecular dynamics simulations for both models and compared their structures and dynamics. The motions of the stalk, consisting of a long coiled coil with a microtubule-binding domain, significantly differed between the two models. The elastic properties of the stalk were analyzed and compared with the experimental results.

Structure of the entire stalk region of the Dynein motor domain.

Dyneins are large microtubule-based motor complexes that power a range of cellular processes including the transport of organelles, as well as the beating of cilia and flagella. The motor domain is located within the dynein heavy chain and comprises an N-terminal mechanical linker element, a central ring of six AAA+ modules of which four bind or hydrolyze ATP, and a long stalk extending from the AAA+ring with a microtubule-binding domain (MTBD) at its tip. A crucial mechanism underlying the motile activity of cytoskeletal motor proteins is precise coupling between the ATPase and track-binding activities. In dynein, a stalk region consisting of a long (~15nm) antiparallel coiled coil separates these two activities, which must facilitate communication between them. This communication is mediated by a small degree of helix sliding in the coiled coil. However, no high-resolution structure is available of the entire stalk region including the MTBD. Here, we have reported the structure of the entire stalk region of mouse cytoplasmic dynein in a weak microtubule-binding state, which was determined using X-ray crystallography, and have compared it with the dynein motor domain from Dictyostelium discoideum in a strong microtubule-binding state and with a mouse MTBD with its distal portion of the coiled coil fused to seryl-tRNA synthetase from Thermus thermophilus. Our results strongly support the helix-sliding model based on the complete structure of the dynein stalk with a different form of coiled-coil packing. We also propose a plausible mechanism of helix sliding together with further analysis using molecular dynamics simulations. Our results present the importance of conserved proline residues for an elastic motion of stalk coiled coil and imply the manner of change between high-affinity state and low-affinity state of MTBD.

NMR characterization of the interaction of the endonuclease domain of MutL with divalent metal ions and ATP.

MutL is a multi-domain protein comprising an N-terminal ATPase domain (NTD) and C-terminal dimerization domain (CTD), connected with flexible linker regions, that plays a key role in DNA mismatch repair. To expand understanding of the regulation mechanism underlying MutL endonuclease activity, our NMR-based study investigated interactions between the CTD of MutL, derived from the hyperthermophilic bacterium Aquifex aeolicus (aqMutL-CTD), and putative binding molecules. Chemical shift perturbation analysis with the model structure of aqMutL-CTD and circular dichroism results revealed that tight Zn(2+) binding increased thermal stability without changing secondary structures to function at high temperatures. Peak intensity analysis exploiting the paramagnetic relaxation enhancement effect indicated the binding site for Mn(2+), which shared binding sites for Zn(2+). The coexistence of these two metal ions appears to be important for the function of MutL. Chemical shift perturbation analysis revealed a novel ATP binding site in aqMutL-CTD. A docking simulation incorporating the chemical shift perturbation data provided a putative scheme for the intermolecular interactions between aqMutL-CTD and ATP. We proposed a simple and understandable mechanical model for the regulation of MutL endonuclease activity in MMR based on the relative concentrations of ATP and CTD through ATP binding-regulated interdomain interactions between CTD and NTD.

Theory for trivial trajectory parallelization of multicanonical molecular dynamics and application to a polypeptide in water.

Trivial trajectory parallelization of multicanonical molecular dynamics (TTP-McMD) explores the conformational space of a biological system with multiple short runs of McMD starting from various initial structures. This method simply connects (i.e., trivially parallelizes) the short trajectories and generates a long trajectory. First, we theoretically prove that the simple trajectory connection satisfies a detailed balance automatically. Thus, the resultant long trajectory is regarded as a single multicanonical trajectory. Second, we applied TTP-McMD to an alanine decapeptide with an all-atom model in explicit water to compute a free-energy landscape. The theory imposes two requirements on the multiple trajectories. We have demonstrated that TTP-McMD naturally satisfies the requirements. The TTP-McMD produces the free-energy landscape considerably faster than a single-run McMD does. We quantitatively showed that the accuracy of the computed landscape increases with increasing the number of multiple runs. Generally, the free-energy landscape of a large biological system is unknown a priori. The current method is suitable for conformational sampling of such a large system to reduce the waiting time to obtain a canonical ensemble statistically reliable.

Structural insight into PPARgamma activation through covalent modification with endogenous fatty acids.

Peroxisome proliferator-activated receptor (PPAR) gamma is a nuclear receptor that regulates lipid homeostasis, and several fatty acid metabolites have been identified as PPARgamma ligands. Here, we present four crystal structures of the PPARgamma ligand binding domain (LBD) covalently bound to endogenous fatty acids via a unique cysteine, which is reportedly critical for receptor activation. The structure analyses of the LBD complexed with 15-deoxy-Delta(12,14)-prostaglandin J(2) (15d-PGJ(2)) revealed that the covalent binding of 15d-PGJ(2) induced conformational changes in the loop region following helix H2', and rearrangements of the side-chain network around the created covalent bond in the LBD. Point mutations of these repositioned residues on the loop and helix H3 almost completely abolished PPARgamma activation by 15d-PGJ(2), indicating that the observed structural alteration may be crucial for PPARgamma activation by the endogenous fatty acid. To address the issue of partial agonism of endogenous PPARgamma ligands, we took advantage of a series of oxidized eicosatetraenoic acids (oxoETEs) as covalently bound ligands to PPARgamma. Despite similar structural and chemical properties, these fatty acids exhibited distinct degrees of transcriptional activity. Crystallographic studies, using two of the oxoETE/PPARgamma LBD complexes, revealed that transcriptional strength of each oxoETE is associated with the difference in the loop conformation, rather than the interaction between each ligand and helix H12. These results suggest that the loop conformation may be responsible for the modulation of PPARgamma activity. Based on these results, we identified novel agonists covalently bound to PPARgamma by in silico screening and a cell-based assay. Our crystallographic study of LBD complexed with nitro-233 demonstrated that the expected covalent bond is indeed formed between this newly identified agonist and the cysteine. This study presents the structural basis for the activation and modulation mechanism of PPARgamma through covalent modification with endogenous fatty acids.

In silico prediction of the chemical block of human ether-a-go-go-related gene (hERG) K+ current.

A variety of compounds with different chemical properties directly interact with the cardiac repolarizing K(+) channel encoded by the human ether-a-go-go-related gene (hERG). This causes acquired forms of QT prolongation, which can result in lethal cardiac arrhythmias including torsades de pointes one of the most serious adverse effects of various therapeutic agents. Prediction of this phenomenon will improve the safety of pharmacological therapy and also facilitate the process of drug development. Here we propose a strategy for the development of an in silico system to predict the potency of chemical compounds to block hERG. The system consists of two sequential processes. The first process is a ligand-based prediction to estimate half-maximal concentrations for the block of compounds inhibiting hERG current using the relationship between chemical features and activities of compounds. The second process is a protein-based prediction that comprises homology modeling of hERG, docking simulation of chemical-channel interaction, analysis of the shape of the channel pore cavity, and Brownian dynamics simulation to estimate hERG currents in the presence and absence of chemical blockers. Since each process is a combination of various calculations, the criterion for assessment at each calculation and the strategy to integrate these steps are significant for the construction of the system to predict a chemical's block of hERG current and also to predict the risk of inducing cardiac arrhythmias from the chemical information. The principles and criteria of elemental computations along this strategy are described.

Folding of the 25 residue Abeta(12-36) peptide in TFE/water: temperature-dependent transition from a funneled free-energy landscape to a rugged one.

The free-energy landscape of the Alzheimer beta-amyloid peptide Abeta(12-36) in a 40% (v/v) 2,2,2-trifluoroethanol (TFE)/water solution was determined by using multicanonical molecular dynamics simulations. Simulations using this enhanced conformational sampling technique were initiated from a random unfolded polypeptide conformation. Our simulations reliably folded the peptide to the experimental NMR structure, which consists of two linked helices. The shape of the free energy landscape for folding was found to be strongly dependent on temperature: Above 325 K, the overall shape was funnel-like, with the bottom of the funnel coinciding exactly with the NMR structure. Below 325 K, on the other hand, the landscape became increasingly rugged, with the emergence of new conformational clusters connected by low free-energy pathways. Finally, our simulations reveal that water and TFE solvate the polypeptide in different ways: The hydrogen bond formation between TFE and Abeta was enhanced with decreasing temperature, while that between water and Abeta was depressed.

Mutational analysis of block and facilitation of HERG current by a class III anti-arrhythmic agent, nifekalant.

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.

Alpha,beta-unsaturated ketone is a core moiety of natural ligands for covalent binding to peroxisome proliferator-activated receptor gamma.

Peroxisome proliferator-activated receptor gamma (PPARgamma) functions in various biological processes, including macrophage and adipocyte differentiation. Several natural lipid metabolites have been shown to activate PPARgamma. Here, we report that some PPARgamma ligands, including 15-deoxy-Delta12,14-prostaglandin J2, covalently bind to a cysteine residue in the PPARgamma ligand binding pocket through a Michael addition reaction by an alpha,beta-unsaturated ketone. Using rhodamine-maleimide as well as mass spectroscopy, we showed that the binding of these ligands is covalent and irreversible. Consistently, mutation at the cysteine residue abolished abilities of these ligands to activate PPARgamma, but not of BRL49653, a non-covalent synthetic agonist, indicating that covalent binding of the alpha,beta-unsaturated ketone in the natural ligands was required for their transcriptional activities. Screening of lipid metabolites containing the alpha,beta-unsaturated ketone revealed that several other oxidized metabolites of hydroxyeicosatetraenoic acid, hydroxyeicosadecaenoic acid, and prostaglandins can also function as novel covalent ligands for PPARgamma. We propose that PPARgamma senses oxidation of fatty acids by recognizing such an alpha,beta-unsaturated ketone as a common moiety.