Free fatty acid receptors:structural models and elucidation of ligand binding interactions

BACKGROUND: The free fatty acid receptors (FFAs), including FFA1 (orphan name: GPR40), FFA2 (GPR43) and FFA3 (GPR41) are G protein-coupled receptors (GPCRs) involved in energy and metabolic homeostasis. Understanding the structural basis of ligand binding at FFAs is an essential step toward designing potent and selective small molecule modulators.

RESULTS: We analyse earlier homology models of FFAs in light of the newly published FFA1 crystal structure co-crystallized with TAK-875, an ago-allosteric ligand, focusing on the architecture of the extracellular binding cavity and agonist-receptor interactions. The previous low-resolution homology models of FFAs were helpful in highlighting the location of the ligand binding site and the key residues for ligand anchoring. However, homology models were not accurate in establishing the nature of all ligand-receptor contacts and the precise ligand-binding mode. From analysis of structural models and mutagenesis, it appears that the position of helices 3, 4 and 5 is crucial in ligand docking. The FFA1-based homology models of FFA2 and FFA3 were constructed and used to compare the FFA subtypes. From docking studies we propose an alternative binding mode for orthosteric agonists at FFA1 and FFA2, involving the interhelical space between helices 4 and 5. This binding mode can explain mutagenesis results for residues at positions 4.56 and 5.42. The novel FFAs structural models highlight higher aromaticity of the FFA2 binding cavity and higher hydrophilicity of the FFA3 binding cavity. The role of the residues at the second extracellular loop used in mutagenesis is reanalysed. The third positively-charged residue in the binding cavity of FFAs, located in helix 2, is identified and predicted to coordinate allosteric modulators.

CONCLUSIONS: The novel structural models of FFAs provide information on specific modes of ligand binding at FFA subtypes and new suggestions for mutagenesis and ligand modification, guiding the development of novel orthosteric and allosteric chemical probes to validate the importance of FFAs in metabolic and inflammatory conditions. Using our FFA homology modelling experience, a strategy to model a GPCR, which is phylogenetically distant from GPCRs with the available crystal structures, is discussed.

A critical evaluation of Raman spectroscopy for the analysis of lipids: fatty acid methyl esters.

The work presented here is aimed at determining the potential and limitations of Raman spectroscopy for fat analysis by carrying out a systematic investigation of C-4-C-24 FAME. These provide a simple, well-characterized set of compounds in which the effect of making incremental changes can be studied over a wide range of chain lengths and degrees of unsaturation. The effect of temperature on the spectra was investigated over much larger ranges than would normally be encountered in real analytical measurements. It was found that for liquid FAME the best internal standard band was the carbonyl stretching vibration nu(C = O), whose position is affected by changes in sample chain length and physical state; in the samples studied here, it was found to lie between 1729 and 1748 cm(-1). Further, molar unsaturation could be correlated with the ratio of the nu(C = O) to either nu(C = C) or delta(H-C = ) with R-2 > 0.995. Chain length was correlated with the delta(CH2)(tw)/nu(C = O) ratio, (where "tw" indicates twisting) but separate plots for odd- and even-numbered carbon chains were necessary to obtain R-2 > 0.99 for liquid samples. Combining the odd- ani even-numbered carbon chain data in a single plot reduced the correlation to R-2 = 0.94-0.96, depending on the band ratios used. For molal unsaturation the band ratio that correlated linearly with unsaturation (R-2 > 0.99) was nu(C = C)/delta(CH2)(SC) (where "sc" indicates scissoring). Other band ratios show much more complex behavior with changes in chemical and physical structure. This complex behavior results from the fact that the bands do not arise from simple vibrations of small, discrete regions of the molecules but are due to complex motions of large sections of the FAME so that making incremental changes in structure does not necessarily lead to simple incremental changes in spectra.

Conformations, vibrational frequencies and Raman intensities of short-chain fatty acid methyl esters using DFT with 6-31G(d) and Sadlej pVTZ basis sets

Density functional calculations, using B3LPY/6-31G(d) methods, have been used to investigate the conformations and vibrational (Raman) spectra of three short-chain fatty acid methyl esters (FAMEs) with the formula CnH2nO2 (n = 3-5). In all three FAMEs, the lowest energy conformer has a simple 'all-trans' structure but there are other conformers, with different torsions about the backbone, which lie reasonably close in energy to the global minimum. One result of this is that the solid samples we studied do not appear to consist entirely of the lowest energy conformer. Indeed, to account for the 'extra' bands that were observed in the Raman data but were not predicted for the all-trans conformer, it was necessary to add-in contributions from other conformers before a complete set of vibrational assignments could be made. Provided this was done, the agreement between experimental Raman frequencies and 6-31G(d) values (after scaling) was excellent, RSD = 12.6 cm(-1). However, the agreement between predicted and observed intensities was much less satisfactory. To confirm the validity of the approach followed by the 6-3 1 G(d) basis set, we used a larger basis set, Sadlej pVTZ, and found that these calculations gave accurate Raman intensities and simulated spectra (summed from two different conformers) that were in quantitative agreement with experiment. In addition, the unscaled Sadlej pVTZ, and the scaled 6-3 1 G(d) calculations gave the same vibrational mode assignments for all bands in the experimental data. This work provides the foundation for calculations on longer-chain FAMEs (which are closer to those found as triglycerides in edible fats and oils) because it shows that scaled 6-3 1 G(d) calculations give equally accurate frequency predictions, and the same vibrational mode assignments, as the much more CPU-expensive Sadlej pVTZ basis set calculations.

Lipid content and fatty acid composition of the digital cushion of bulls offered different amounts of linseed

Previous research suggests that the digital cushion, a shock-absorbing structure in the claw, plays an important role in protecting cattle from lameness. This study aimed to assess the degree to which nutritional factors influence the composition of the digital cushion. This involved quantifying lipid content and fatty acid composition differences in digital cushion tissue from cattle offered diets with different amounts of linseed. Forty-six bulls were allocated to 1 of 4 treatments, which were applied for an average of 140 +/- 27 d during the finishing period. The treatments consisted of a linseed supplement offered once daily on top of the basal diet (grass silage:concentrate) at 0, 400, 800, or 1,200 g of supplement/animal per day. For each treatment, the concentrate offered was adjusted to ensure that total estimated ME intake was constant across treatments. Target BW at slaughter was 540 kg. Legs were collected in 3 batches after 120, 147 and 185 d on experiment. Six samples of the digital cushion were dissected from the right lateral hind claw of each animal. Lipids were extracted and expressed as a proportion of fresh tissue, and fatty acid composition of the digital cushion was determined by gas chromatography. Data were analyzed by ANOVA, with diet, location within the digital cushion, and their interactions as fixed effects and fat content (grams per 100 g of tissue) as a covariate. Linear or quadratic contrasts were examined. The lipid content of digital cushion tissue differed between sampling locations (P