Adrenocorticotropic Hormone (ACTH) Responses Require Actions of the Melanocortin-2 Receptor Accessory Protein on the Extracellular Surface of the Plasma Membrane

The melanocortin-2 (MC2) receptor is a G protein-coupled receptor that mediates responses to ACTH. The MC2 receptor acts in concert with the MC2 receptor accessory protein (MRAP) that is absolutely required for ACTH binding and signaling. MRAP has a single transmembrane domain and forms a highly unusual antiparallel homodimer that is stably associated with MC2 receptors at the plasma membrane. Despite the physiological importance of the interaction between the MC2 receptor and MRAP, there is little understanding of how the accessory protein works. The dual topology of MRAP has made it impossible to determine whether highly conserved and necessary regions of MRAP are required on the intracellular or extracellular face of the plasma membrane. The strategy used here was to fix the orientation of two antiparallel MRAP molecules and then introduce inactivating mutations on one side of the membrane or the other. This was achieved by engineering proteins containing tandem copies of MRAP fused to the amino terminus of the MC2 receptor. The data firmly establish that only the extracellular amino terminus (N_out) copy of MRAP, oriented with critical segments on the extracellular side of the membrane, is essential. The transmembrane domain of MRAP is also required in only the N_out orientation. Finally, activity of MRAP-MRAP-MC2-receptor fusion proteins with inactivating mutations in either MRAP or the receptor was rescued by co-expression of free wild-type MRAP or free wild-type receptor. These results show that the basic MRAP-MRAP-receptor signaling unit forms higher order complexes and that these multimers signal.     

Structural insight into the binding mechanism of ATP to EGFR and L858R,and T790M and L858R/T790 mutants

Abstract

The L858R mutation in EGFR is particularly responsive to small tyrosine kinase inhibitors (TKIs) such as gefitinib and erlotinib. This efficacy decreases due to drug resistance conferred by a second mutation, T790M, which subsequently produces a double mutant, L858R/T790M. Although this resistance was initially attributed to steric blocking by the T790M mutation, experimental studies have demonstrated that differences in the binding affinities of TKIs to T790M and L858R/T790M mutants are more a result of the increased sensitivity of these mutants to ATP than to a decrease in the affinity to TKIs. Regrettably, detailed information at the atomic level on the origins of the increased binding affinity of mutants for ATP is lacking. In this study, we have combined structural data and molecular dynamics simulations with the MMGBSA approach to determine how the L858R, T790M and L858R/T790 mutations impact the binding mechanism of ATP with respect to wild-type EGFR. Structural and energetic analyses provided novel information that helps to explain the increased affinity of ATP to T790M and L858R/T790 mutants with respect to L858R and wild-type systems. In addition, it was observed that dimerization of the wild-type and mutant systems exerts dissimilar effects on the ATP binding affinity characteristic of negative cooperativity.

Communicated by Ramaswamy H. Sarma     

Kinetics of pH-dependent interactions between PD-1 and PD-L1 immune checkpoint proteins from molecular dynamics

Immune checkpoint blockade of signaling pathways such as PD-1/PD-L1 has recently opened up a new avenue for highly efficient immunotherapeutic strategies to treat cancer. Since tumor microenvironments are characterized by lower pH (5.5-7.0), pH-dependent protein-ligand interactions can be exploited as efficient means to regulate drug affinity and specificity for a variety of malignancies. In this article, we investigate the mechanism and kinetics of pH-dependent binding and unbinding processes for the PD-1/PD-L1 checkpoint pair employing classical molecular dynamics simulations. Two representative pH levels corresponding to circumneutral physiological conditions of blood (pH 7.4) and acidic tumor microenvironment (pH 5.5) are considered. Our calculations demonstrate that pH plays a key role in protein-ligand interactions with small pH changes leading to several orders of magnitude increase in binding affinity. By identifying the binding pocket in the PD-1/PD-L1 complex, we show a pivotal role of the His68 protonation state of PD-1in the complex stabilization at low pH. The results on the reaction rate constants are in qualitative agreement with available experimental data. The obtained molecular details are important for further engineering of binding/unbinding kinetics to formulate more efficient immune checkpoint blockade strategies.     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Modeling and simulation of the human delta opioid receptor

A model for the human delta opioid receptor has been generated via sequence alignment, structure building using the crystal structure of bovine rhodopsin as a template, and refinement by molecular dynamics simulation. The model building suggested that, in addition to the previously postulated interaction between D128 and Y308, an internal salt bridge also exists between residues D128 and R192, both of which are conserved in all the opioid receptors. The model and salt bridge were then shown to be stable during a 20-nsec simulation in a lipid bilayer. It is therefore proposed that both of these interactions play a role in stabilizing the inactive state of the receptor. The model is also used in an effort to rationalize many of the mutational studies performed on delta opioid receptors, and to suggest a plausible explanation for the differences between known delta opioid agonists and antagonists.     

Structural activation pathways from dynamic olfactory receptor-odorant interactions

We have simulated an odor ligand's dynamic behavior in the binding region of an olfactory receptor (OR). Our short timescale computational studies (up to 200 ps) have helped identify unprecedented postdocking ligand behavior of ligands. From in vacuo molecular dynamics simulations of interactions between models of rat OR I7 and 10 aldehyde ligands, we have identified a dissociative pathway along which the ligand exits and enters the OR-binding pocket--a transit event. The ligand's transit through the receptor's binding region may mark the beginning of a signal transduction cascade leading to odor recognition. We have graphically traced the rotameric changes in key OR amino acid side chains during the transit. Our results have helped substantiate or refute previously held notions of amino acid contribution to ligand stability in the binding pocket. Our observations of ligand activity when compared to those of experimental (electroolfactogram response) OR-activation studies provide a view to predicting the stability of ligands in the binding pocket as a precursor to OR activation by the ligand.     

Computational and experimental approaches for assessing the interactions between the model calycin beta-lactoglobulin and two antibacterial fluoroquinolones

Norfloxacin and levofloxacin, two fluoroquinolones of different bulk, rigidity and hydrophobicity taken as model ligands, were docked to one apo and two holo crystallographic structures of bovine beta-lactoglobulin (BLG) using different computational approaches. BLG is a member of the lipocalin superfamily. Lipocalins show a typical b-barrel structure encompassing an internal cavity where small hydrophobic molecules are usually bound. Our studies allowed the identification of two putative binding sites in addition to the calyx. The rigid docking approximation resulted in strong repulsive forces when the ligands were docked into the calyx of the apo form. On the contrary, hindrance was not experienced in flexible docking protocols whether on the apo or on the holo BLG forms, due to allowance for side chain rearrangement. K(i) between 10(-7) and 10(-6) M were estimated for norfloxacin at pH 7.4, smaller than 10(-5) M for levofloxacin. Spectroscopic and electrophoretic techniques experimentally validated the occurrence of an interaction between norfloxacin and BLG. Changes in chemical shift and dynamic parameters were observed between the (19)F NMR spectra of the complex and of the ligand. A K(i) (ca 10(-7) M) comparable with the docking results was estimated through a NMR relaxation titration. Stabilization against unfolding was demonstrated by denaturant gradient gel electrophoresis on the complex versus apo BLG. NMR experimental evidence points to a very loose interaction for ofloxacin, the racemic mixture containing levofloxacin. Furthermore, we were able to calculate in silico K(i)'s comparable to the published experimental values for the complexes of palmitic and retinoic acid with BLG.     

Melanin content and MC1R function independently affect UVR‐induced DNA damage in cultured human melanocytes

Malignant transformation of melanocytes leads to melanoma, the most fatal form of skin cancer. Ultraviolet radiation (UVR)‐induced DNA photoproducts play an important role in melanomagenesis. Cutaneous melanin content represents a major photoprotective mechanism against UVR‐induced DNA damage, and generally correlates inversely with the risk of skin cancer, including melanoma. Melanoma risk is also determined by susceptibility genes, one of which is the melanocortin 1 receptor (MC1R) gene. Certain MC1R alleles are strongly associated with melanoma. We hereby present experimental evidence for the role of two melanoma risk factors, constitutive pigmentation, as assessed by total melanin, eumelanin and pheomelanin contents, and MC1R genotype and function, in determining the induction and repair of DNA photoproducts in cultured human melanocytes after irradiation with increasing doses of UVR. We found that total melanin and eumelanin contents (MC and EC) correlated inversely with the extent of UVR‐induced growth arrest, apoptosis and induction of cyclobutane pyrimidine dimers (CPD), but not with hydrogen peroxide release in melanocytes expressing functional MC1R. In comparison, melanocytes with loss‐of‐function MC1R, regardless of their MC or EC, sustained more UVR‐induced apoptosis and CPD, and exhibited reduced CPD repair. Therefore, MC, mainly EC, and MC1R function are independent determinants of UVR‐induced DNA damage in melanocytes.     

藏獒MC1R基因多态性与毛色表型的关联研究

目的:通过检测藏獒黑素皮质激素受体1(MC1R)基因的单链构象多态性(SSCP)在不同毛色群体中的分布,探讨MC1R基因多态性与毛色表型的相关性。方法:采用DNA测序技术,选择不同毛色藏獒的DNA为样本,根据GenBank发布的荷斯坦牛MC1R基因序列设计一对引物,采用PCR-SSCP技术分析MC1R基因在藏獒中的SSCP。结果:MC1R基因在藏獒中具有PCR-SSCP多态性,分别检测到3种基因型(AA、AB和BB);对MC1R基因多态性片段DNA克隆测序后发现,MC1R基因在编码区第313位存在单碱基突变(G→A),该突变导致第105位氨基酸发生由丙氨酸向苏氨酸的改变(T105A)。结论:MC1R基因的多态性与毛色性状不存在显著的相关性。     

Computational studies of LXR molecular interactions reveal an allosteric communication pathway

The liver X receptor, LXRα, is an important regulator of genes involved in metabolism and inflammation. The mechanism of communication between the cofactor peptide and the ligand in the ligand-binding pocket is a crucial and often discussed issue for the nuclear receptors (NRs), but such allosteric signaling pathways are difficult to detect and the transmission mechanism remains elusive. Here, we apply the anisotropic thermal diffusion method to the LXRα with bound coactivator and ligand. We detected a possible communication pathway between the coactivator peptide and the ligand. The signal is transmitted both through the receptor backbone and side chains. A key signaling residue is the first leucine in the cofactor peptide recognition motif LXXLL, which is conserved within the NR cofactors, suggesting a general mechanism for allosteric signaling. Furthermore, we studied the LXR receptor and cofactor molecular interactions in detail using molecular dynamics simulations. The protein-protein interaction patterns in the complexes of nine different cofactor peptides and holo-LXRα were characterized, revealing the importance of the receptor-cofactor charge clamp interaction. Specific, but infrequently occurring interactions were observed toward the cofactor peptide C-terminal residues. Thus, additional specificity between LXRα and its cofactors is likely to be found in molecular interactions outside the cofactor peptide or in other biological factors.     

The expression of MC4Rs in D1R neurons regulates food intake and locomotor sensitization to cocaine

While it is known that mice lacking melanocortin 4 receptor (MC4R) expression develop hyperphagia resulting in early‐onset obesity, the specific neural circuits that mediate this process remain unclear. Here, we report that selective restoration of MC4R expression within dopamine‐1 receptor‐expressing neurons [MC4R/dopamine 1 receptor (D1R) mice] partially blunts the severe obesity seen in MC4R‐null mice by decreasing meal size, but not meal frequency, in the dark cycle. We also report that both acute cocaine‐induced anorexia and the development of locomotor sensitization to repeated administration of cocaine are blunted in MC4R‐null mice and normalized in MC4R/D1R mice. Neuronal retrograde tracing identifies the lateral hypothalamic area as the primary target of MC4R‐expressing neurons in the nucleus accumbens. Biochemical studies in the ventral striatum show that phosphorylation of DARPP‐32^Thr‐34 and GluR1^Ser‐845 is diminished in MC4R‐null mice after chronic cocaine administration but rescued in MC4R/D1R mice. These findings highlight a physiological role of MC4R‐mediated signaling within D1R neurons in the long‐term regulation of energy balance and behavioral responses to cocaine.     

ACTH receptor: Ectopic expression, activity and signaling

Failure in obtaining expression of functional adrenocorticotropic hormone receptor (ACTHR, or melanocortin 2 receptor, MC2R) in non-adrenal cells has hindered molecular analysis of ACTH signaling pathways. Here, we ectopically expressed the mouse ACTHR in Balb/c mouse 3T3 fibroblasts to analyze ACTH signaling pathways involved in induction of fos and jun genes. Natural constitutive expression of the MC2R accessory protein (MRAP) in Balb3T3 and other mouse 3T3 fibroblasts (NIH, Swiss and 3T3-L1) renders these fibroblastic lines suitable for ectopic expression of ACTHR in its active form properly inserted into the plasma membrane at levels similar to those found in mouse Y1 adrenocortical tumor cells. The Y1 cell line is a cultured cell system well known for stably displaying normal adrenal specific metabolic pathways, ACTHR expression and ACTH functional responses. Thirty-nine sub-lines expressing ACTHR (3T3-AR transfectants) were selected for geneticin-resistance and clonally isolated after transfection of ACTHR-cDNA (in the pSVK3 mammalian plasmidial vector) into Balb3T3 fibroblasts. In addition, sixteen clonal sub-lines of Balb3T3 (3T3-0 transfectants) carrying the pSVK3 empty vector were likewise isolated. Fourteen 3T3-AR and four 3T3-0 clones were screened for response to ACTH₃₉ in comparison with Y1 adrenocortical cells. Eight 3T3-AR clones responded to ACTH₃₉ with activation of adenylate cyclase and induction of c-Fos protein, but the levels of, respectively, activation and induction were not strictly correlated. Other fos and jun genes were also induced by ACTH₃₉ in 3T3-AR transfectants, which express levels of ACTHR protein similar to parental Y1 cells. Signaling pathways relevant to c-Fos induction was extensively investigated in 3 clones: 3T3-AR01 and –07 and 3T3-04. In Y1 cells, specific inhibitors (H89/PKA; PD98059/MEK; Go6983/PKC and SP600125/JNK) show that signals initiated in the ACTH/ACTHR-system activate 4 pathways to induce the c-fos gene, namely: (a) cAMP/PKA/CREB; (b) MEK/ERK1/2; (c) PKC and d) JNK1/2. In 3T3-AR transfectants, both inhibitors PD98059 and Go6983 proved completely ineffective to inhibit c-Fos induction by ACTH₃₉, implying that MEK/ERK and PKC pathways are not involved in this process. On the other hand, SP600125 caused 85% inhibition of c-Fos induction by ACTH₃₉ and, in addition, ACTH₃₉ promotes JNK1/2 phosphorylation, suggesting that JNK is a major signaling pathway mediating c-Fos induction by ACTH₃₉ in these cells. In addiction, PKA inhibitor H89 also inhibits c-Fos induction in 3T3-AR7 cells by ACTH₃₉, implicating activation of the cAMP/PKA/CREB pathway in c-Fos induction by ACTH₃₉. However, the cAMP derivatives db-cAMP and 8Br-cAMP, do not promote CREB phosphorylation and c-Fos induction in parental Balb3T3 and 3T3-AR transfectants, confirming previous report by others. In conclusion, expression of active ACTHR in Balb3T3 fibroblasts renders these cells responsive to ACTH with activation of cAMP/PKA/CREB and JNK pathways and, also, induction of genes from the fos and jun families. These results show that Balb 3T3-AR sublines are useful cellular systems for genetic analysis of ACTH-signaling pathways. However, activation of cAMP/PKA/CREB and JNK pathways and induction of fos and jun genes are not yet sufficient to enable ACTH for interference in morphology, migration and proliferation of Balb3T3 fibroblasts as it does in Y1 adrenocortical cells.     

Relative stability of G‐quadruplex structures: Interactions between the human Bcl2 promoter region and derivatives of carbazole and diphenylamine

The bcl2 promoter region forms a G‐quadruplex structure, which is a crucial target for anticancer drug development. In this study, we provide theoretical predictions of the stability of different G‐quadruplex folds of the 23‐mer bcl2 promoter region and G‐quadruplex ligand. We take into account the whole G‐quadruplex structure, including bound‐cations and solvent effects, in order to compute the ligand binding free energy using molecular dynamics simulation. Two series of the carbazole and diphenylamine derivatives are used to screen for the most potent drug in terms of stabilization. The energy analysis identifies the predominant energy components affecting the stability of the various different G‐quadruplex folds. The energy associated with the stability of the G‐quadruplex‐K⁺ structures obtained displays good correlation with experimental T_m measurements. We found that loop orientation has an intrinsic influence on G‐quadruplex stability and that the basket structure is the most stable. Furthermore, parallel loops are the most effective drug binding site. Our studies also demonstrate that rigidity and planarity are the key structural elements of a drug that stabilizes the G‐quadruplex structure. BMVC‐4 is the most potential G‐quadruplex ligand. This approach demonstrates significant promise and should benefit drug design. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1038–1050, 2014.     

Theoretical evaluation and improvement on the potency of the rhodanine-based inhibitors for human serotonin N-acetyltransferase

Human serotonin N-acetyltransferase (hAANAT), included in the melatonin biosynthesis, plays a pivotal role in the regulation of the biological clock and the daily rhythm. In this research, a reliable model of hAANAT was first constructed by the homology modelling method. Then the inhibition mode of two representative rhodanine-based inhibitors was explored by molecular dynamics simulations and energy analyses. The results show that the inhibitor class could share a similar inhibition mechanism in which the carboxyl moiety is positioned in the Ac-CoA binding region while the other end spans the serotonin binding pocket. The interaction between the inhibitor's carboxyl and the enzyme seems to be more important according to the decomposition of binding free energy. Based on the proposed inhibition mode, the inhibitor's improvement was carried out to obtain a more potent compound. The newly designed inhibitor, with the larger binding free energy, exhibits the stronger interaction with the related residues of the enzyme by the added chemical groups. This work will shed light on the inhibition mechanism of the rhodanine-based inhibitors and promote the development of a new drug targeting hAANAT.     

How inaccuracies in protein structure models affect estimates of protein-ligand interactions: computational analysis of HIV-I protease inhibitor binding

The influence of possible inaccuracies that can arise during homology modeling of protein structures used for ligand binding studies were investigated with the molecular mechanics generalized Born surface area (MM-GBSA) method. For this, a family of well-characterized HIV-I protease-inhibitor complexes was used. Validation of MM-GBSA led to a correlation coefficient ranging from 0.72 to 0.93 between calculated and experimental binding free energies DeltaG. All calculated DeltaG values were based on molecular dynamics simulations with explicit solvent. Errors introduced into the protein structure through misplacement of side-chains during rotamer modeling led to a correlation coefficient between DeltaG(calc) and DeltaG(exp) of 0.75 compared with 0.90 for the correctly placed side chains. This is in contrast to homology models for members of the retroviral protease family with template structures ranging in sequence identity between 32% and 51%. For these protein models, the correlation coefficients vary between 0.84 and 0.87, which is considerably closer to the original protein (0.90). It is concluded that HIV-I low sequence identity with the template structure still allows creating sufficiently reliable homology models to be used for ligand-binding studies, although placement of the rotamers is a critical step during the modeling.     

Evaluation of the porcine melanocortin 4 receptor (MC4R) gene as a positional candidate for a fatness QTL in a cross between Landrace and Hampshire

Melanocortin 4 receptor (MC4R) is expressed in the appetite-regulating areas of the brain where it is central in the regulation of feed intake and energy balance. A mutation in MC4R causing an Asp298Asn substitution has been associated with fatness, high daily gain and feed intake in the pig. In a previously performed genome scan based on a Hampshire x Landrace cross, we detected one quantitative trait loci (QTL) affecting carcass fat/meat ratio and one QTL affecting the biceps femoris muscle, both close to the position of MC4R on porcine chromosome 1. In this study, the two lines were found to be close to fixation for alternative alleles of the Asp298Asn polymorphism. Additional QTL analyses supported our hypothesis of MC4R as a positional candidate gene but only for the fat/meat QTL. The Asp298Asn polymorphism was also evaluated as a selection target for daily gain in a Danish pig breeding population that included four breeds (Hampshire, Duroc, Landrace and Yorkshire). Over a 12-year period (1990-2002), a significant increase in the allele frequency of 298Asn was found in Landrace and Duroc, whereas a non-significant decrease in the 298Asn allele frequency was observed in Yorkshire. The Hampshire breed was fixed for the 298Asn allele in 1990. The high 298Asn allele frequencies in Hampshire, Landrace and Duroc are most likely due to selection for daily gain, whereas selection for daily gain in the Yorkshire breed apparently focuses on other loci.     

Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions

1. Download : Download high-res image (214KB)
2. Download : Download full-size image

     

敲减MC4R表达对牛胎儿成纤维细胞CMS系统关键因子的影响

为获得敲减黑素皮质素4受体（melanocortin 4 receptor,MC4R）基因的牛胎儿成纤维细胞,并探讨其在能量平衡神经调节系统中的作用,将构建成功并已鉴定为有效序列的短发夹状RNA (short hairpin RNA, shRNA)真核表达载体pGSH1 GFP MC4R,利用阳离子脂质体转染牛胎儿成纤维细胞并使用G418筛选稳定转染细胞株.利用实时荧光定量和Western印迹检测MC4R及中枢黑素皮质素系统(central melanocortin system, CMS)关键因子的表达水平变化.结果表明,在稳定转染的牛胎儿成纤维细胞系中, MC4R表达显著抑制,瘦蛋白(leptin)和阿黑色素原(POMC)表达下调,黑素皮质素拮抗物agouti相关蛋白(AGRP)和MC3R表达上调,而神经肽Y (NPY)表达无明显改变.综上所述,本研究成功获得了敲减MC4R基因表达的牛胎儿成纤维细胞.相关基因表达水平检测结果提示, MC4R的表达水平对CMS系统中的各关键基因的表达有不同的抑制或促进影响.     

Modeling the interaction between the N‐terminal domain of the tumor suppressor p53 and azurin

It is known that the half life of the tumor suppressor p53 can be increased by the interaction with the bacterial protein azurin, resulting in an enhanced anti‐tumoral activity. The understanding of the molecular mechanisms on the basis of this phenomenon can open the way to new anti‐cancer strategies. Some experimental works have given evidence of an interaction between p53 and azurin (AZ); however the binding regions of the proteins are still unknown. Recently, fluorescence studies have shown that p53 partakes in the binding with the bacterial protein by its N‐terminal (NT) domain. Here we have used a computational method to get insight into this interacting mode. The model that we propose for the best complex between AZ and p53 has been obtained from a rigid‐body docking, coupled with a molecular dynamics (MD) simulation, a free energy calculation, and validated by mutagenesis analysis. We have found a high degree of geometric fit between the two proteins that are kept together by several hydrophobic interactions and numerous hydrogen bonds. Interestingly, it has emerged that AZ binds essentially to the helices H_I and H_III of the p53 NT domain, which are also interacting regions for the foremost inhibitor of p53, MDM2. Copyright © 2009 John Wiley & Sons, Ltd.     

Modeling fatty acid delivery from intestinal fatty acid binding protein to a membrane

Intestinal fatty acid binding protein (IFABP) interacts with biological membranes and delivers fatty acid (FA) into them via a collisional mechanism. However, the membrane-bound structure of the protein and the pathway of FA transfer are not precisely known. We used molecular dynamics (MD) simulations with an implicit membrane model to determine the optimal orientation of apo- and holo-IFABP (bound with palmitate) on an anionic membrane. In this orientation, the helical portal region, delimited by the alphaII helix and the betaC-betaD and betaE-betaF turns, is oriented toward the membrane whereas the putative beta-strand portal, delimited by the betaB-betaC, betaF-betaG, betaH-betaI turns and the N terminus, is exposed to solvent. Starting from the MD structure of holo-IFABP in the optimal orientation relative to the membrane, we examined the release of palmitate via both pathways. Although the domains can widen enough to allow the passage of palmitate, fatty acid release through the helical portal region incurs smaller conformational changes and a lower energetic cost.