|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 44, 32471-32479, November 2, 2007
Functional Analysis of Transmembrane Domain 2 of the M1 Muscarinic Acetylcholine Receptor*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The crystal structure of bovine rhodopsin in the ground state (6) has provided a template for modeling the inactive state of M1 mAChRs (7). However, a model does not in itself yield the secret of how ligand binding promotes the transition from an inactive to an active conformation that is transduced into a signal. Alanine-scanning mutagenesis with phenotypic classification of the mutants is one means of fleshing out the structural skeleton with functional data (8).
Transmembrane domain 2 (TM2) is one of the core domains of the rhodopsin-like GPCRs. It contains the highly conserved sequence motif (N/S)LAX(A/S)D (9). (SLACAD in the M1 mAChR) within which Asp2.50 (in standard nomenclature (10)) is one of the most highly conserved residues in the GPCR superfamily. In rhodopsin, it is directly hydrogen bonded to Asn1.50 (corresponding to Asn43 in the M1 mAChR), and indirectly, via a water-mediated network, to Asn7.49 (Asn414,M1) (6, 11, 12). In the M1 mAChR, Asn substitution first established the critical role of Asp71 in signal transduction (13), a finding subsequently confirmed in the M2 (14) and M3 (15) subtypes as well as a host of other receptor types (16). Reciprocal mutagenesis experiments on a number of receptors (17-22) have suggested that pairing between positions 2.50 and 7.49 is important for receptor activation.
Functional analysis of other positions in TM2 remains patchy and incomplete. The most systematic studies have been a substituted cysteine scan from residues 2.47 to 2.68 in the D2 receptor (23), and random mutagenesis from residues 2.43 to 2.63 of the C5a receptor (24). In the former case, the accessibility of the mutant side chains to sulfhydryl reagents was the issue, whereas in the latter the criterion was preservation of receptor signaling, which defined important positions but did not directly address the function of the side chains. The recent use of a yeast genetic screen to isolate inactivating mutations in the M3 mAChR (25) identified a number of non-conservative mutations in TM2. Of a total of 20 mutants, none were Ala substitutions, and only 9 were amenable to characterization at the level of binding or function. This study therefore confirmed the functional importance of the highly conserved amino acids of TM2, but again did not probe the role of the side chains within the structure.
Several other studies have been guided by naturally occurring allelic variants, some of which show constitutive activity, e.g. M2.53V in the thyroid-stimulating hormone receptor (26); interestingly, the adjacent position 2.52 in the angiotensin 2 receptor has been reported to exhibit an altered orientation in a constitutively active mutant (27). Thus, sequences surrounding position 2.50 are important for receptor activation.
The aim of this study was to make a systematic analysis of the functional roles of all of the amino acid side chains from positions 2.40 (Asn61) to 2.60 (Leu81) of the M1 mAChR by alanine substitution mutagenesis (Ala itself was mutated to Gly). Additionally, we have substituted Asn43 (1.50). The replacement of an amino acid by alanine deletes the side chain beyond the
-carbon atom, ablating any interactions that it makes, leaving a small hole in the three-dimensional structure without introducing new contacts to confuse the functional effect (28). By measuring the changes in expression of receptors that are appropriately folded to bind ligands, and the effects on antagonist and agonist affinity and signal transduction we have been able to assign the interactions made by the amino acid side chains to different functional categories. We have used this classification to group them with residues from other parts of the receptor structure. In this way, we have been able to suggest functional roles for most of the important residues in TM2.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Materials—(-)-N-[3H]Methylscopolamine (84 Ci/mmol) and myo-D-[3H]inositol (80 Ci/mmol) were purchased from Amersham Biosciences. Unlabeled ligands were from Sigma. The QuikChangeTM kit was from Stratagene.
Mutagenesis and Expression—Briefly, residues of the rat M1 mAChR were mutated to Ala using the QuikChange method. Ala itself was mutated to Gly. Mutant receptors, cloned into the pCD expression vector, and validated by dideoxy sequencing, were transiently expressed in COS-7 cells by electroporation in 1-cm cuvettes using a Bio-Rad Gene Pulser (260 volts, 960 microfarads). Total cell membrane preparations (100,000 x g, 30 min) were made after 72 h as described previously (29). Where necessary, expression levels of poorly expressed mutant receptors were rescued by treatment of the cultured cells with atropine (10-6 M) for 48 h before washing and harvesting for membrane preparations, or PI turnover assays, as described previously (30).
Equilibrium Binding and Functional Assays—Membrane binding and PI assays were performed as described previously with small variations (29, 31). Briefly, binding of (-)-N-[3H] methylscopolamine ([3H]NMS) to membrane preparations (5-20 µg of membrane protein/ml) was measured at 30 °C in a buffer containing 20 mM Na-Hepes, 100 mM NaCl, and 1 mM MgCl2, pH 7.5, using radioligand concentrations ranging from 1 x 10-11 to 3 x 10-9 M with 2 concentrations per decade in a final volume of 1 ml using an incubation time of 2 h. Nonspecific binding was defined with 10-6 M atropine. Assays were performed in triplicate. The binding reaction was terminated by rapid filtration on a Brandel or TomTec cell harvester. The mean expression level of wild-type receptor binding sites for [3H]NMS was 1.3 ± 0.4 pmol/mg protein (mean ± S.D.). The binding of ACh was measured by inhibition of the binding of [3H]NMS (usually 3 x 10-10 M) using ACh concentrations ranging from 10-8 to 10-2 M with 2 concentrations per decade.
ACh-stimulated total PI breakdown was assayed for 30 min at 37 °C in Krebs-bicarbonate solution supplemented with 10 mM LiCl to inhibit inositol monophosphatases following pre-labeling of the transfected COS-7 cells in 12-well plates with [3H]inositol (1 µCi/ml) for 48 h (31). ACh concentrations were 10-10-10-3 M.
Whole cell binding capacity of [3H]NMS was measured by incubating transfected cells in 12-well plates with 0.5 ml of 10-9 M [3H]NMS in phosphate-buffered saline at 4 °C for 4 h. Nonspecific binding was defined with 10-6 M atropine. The cells were washed three times in ice-cold phosphate-buffered saline, and harvested in 0.5 ml of 1% Triton-X-100. The lysate was counted in 4.5 ml of scintillation fluid. The mean expression level of the wild-type receptor binding sites for [3H]NMS was 1.5 ± 0.3 pmol/mg of protein (mean ± S.D.). This corresponds to
105 [3H]NMS binding sites per cell based on high levels (>75%) of transfection assessed by immunocytochemistry with an antibody directed against the C-terminal 13 amino acids of the M1 mAChR that were similar to those reported previously (29).
Receptor Thermal Stability Assays—Membranes were resuspended in binding buffer and aliquots preincubated at a pre-determined 53 °C for time periods between 2 and 240 min. Control samples and membranes not being incubated at the time were kept on ice. After preincubation, the samples were diluted into binding buffer to give a concentration of 5-20 µg of protein/ml and re-homogenized. A single concentration of 10-9 M (final concentration) of [3H]NMS was mixed with membrane preparations, vehicle (binding buffer), or vehicle containing 10-6 M atropine to determine total and nonspecific binding as described above.
Data Analysis—Saturation binding curves for [3H]NMS were fitted to a one-site model of binding using SigmaPlot 8.0 to yield a total concentration of binding sites and an affinity constant. The results are expressed as pKd. Expression levels (RT) were normalized to simultaneously transfected wild-type controls. Inhibition curves for ACh were fitted to the Hill equation to yield a pIC50 value and a slope factor (nH). The pIC50 values were corrected for the Cheng-Prusoff shift, as necessary. Thermal inactivation curves were fitted to single (or where appropriate double) exponential functions yielding a rate constant for thermal decay.
PI dose-response curves were fitted to a four-parameter logistic function, yielding a pEC50 (-log M) value and basal and maximum responses (Basal, Emax). Slope factors were close to 1.0. For the wild-type receptor and most mutants the ACh-independent basal PI signal was in the range 300-600 disintegrations/min, and the mean ACh-stimulated maximum activity was 3-5 times the basal activity. Values of Basal and Emax for the mutants are expressed relative to wild-type values obtained from contemporaneous controls.
We have shown previously, by blockade with an irreversible antagonist, that the apparent ratio of functional wild-type M1 mAChR to G protein is about 20 after expression in COS-7 cells (29). The levels of expression of functional receptors are affected by mutations, and this can affect the basal activity, ACh potency, and maximum response obtained in PI assays. These variations were taken into account in calculating values of the signaling efficacy of the ACh-receptor complex (eA) using equations described previously (30, 33, 34) that are based on the free association of receptors and G-proteins in the cell membrane (35) as shown,
![]() |
when the Emax value of the mutant (as a fraction of that of the wild-type receptor) is greater than 0.9 or
![]() |
when the Emax of the mutant is less than 0.9. The two equations yield equivalent results, but when the potency ratio IC50/EC50 is large, Emax is close to 1; conversely, when Emax is substantially less than 0.9, IC50
EC50.
|
Statistical Analysis—Experiments were repeated at least 3 times. Values are tabulated as mean ± S.E. Statistical comparisons of affinity and rate constants for mutants and wild-type controls were carried out by one-way analysis of variance followed by Dunnett's post-hoc test. Where values, such as expression levels, had been normalized to the wild-type control, two-tailed t-tests were used to ascertain the level of significance of differences from the wild-type.
Molecular Modeling—Results were analyzed in the context of a rhodopsin-based homology model of the M1 mAChR, as described previously (7).
| RESULTS |
|---|
|
|
|---|
We have found that certain mutants that yield very low levels of [3H]NMS binding sites also give reduced or undetectable PI signaling, but that signaling can be rescued by atropine treatment of the transfected cells in culture, in parallel with the restoration of [3H]NMS binding (7, 30). This suggests that atropine treatment promotes the formation of fully folded receptors, and that only these can bind [3H]NMS and signal. Following the procedure described previously (30), culture of the cells transfected with the D71A, G75A, and M79A mutants with 10-6 M atropine prior to harvesting rescued their expression levels, raising them to 65 ± 6, 111 ± 10, and 88 ± 13% of wild-type, respectively (supplementary Table S1), thus enabling their binding and functional properties to be characterized; expression of the wild-type receptor increased by only 43% as a result of this procedure.
|
Acetylcholine Binding—ACh binding affinities were measured by competition with [3H]NMS and analyzed using the Hill equation. For the wild-type receptor, the pIC50 was 5.00 ± 0.03 and the slope factor (nH) was 0.88 ± 0.02. Under the assay conditions used here, there is no reproducible GTP effect on ACh affinity (36), and a more complex form of data analysis did not seem to be justified. The effects of the mutations on ACh affinity are summarized in Fig. 2a. Full details are given in supplementary Table S2.
Reductions in ACh affinity of about 5-fold were seen for A70G and I74A, with 2-fold reductions for S66A and L67A. In contrast, two of the mutations, L64A and N61A, yielded 5-7-fold increases in ACh affinity, whereas N80A gave a 3-fold increase. A 10-fold increase was also seen for N43A (pIC50 = 6.01 ± 0.09). The mutations did not significantly affect the slope factor for ACh binding.
Phosphoinositide Functional Response—The effects of the mutations on receptor function were assessed by measurement of the ACh-induced total PI response following labeling of the cells with [3H]inositol. Representative PI dose-response curves for a selection of mutants are shown in supplementary Fig. S1. Values for basal (ACh-independent) and maximal signaling were computed as a fraction of the values measured for contemporaneously transfected wild-type controls. The wild-type M1 mAChR typically gave a maximum ACh-stimulated signal equivalent to 3-5 times the basal activity (300-600 dpm) with a mean pEC50 value of 6.91 ± 0.03 (n = 34). In the case of D71A, G75A, and M79A, which gave low levels of expression of [3H]NMS binding sites, additional measurements were carried out after enhancement of the expression levels by atropine rescue. The atropine rescue protocol, applied to the wild-type receptor, decreased basal activity by 23% and ACh potency by about 3-fold with little effect on the Emax, possibly indicating the persistence of some residual atropine. Values for atropine-rescued mutants were related to the corresponding atropine-treated controls.
The mutant D71A showed no ACh-induced PI signal, either before or after atropine rescue of expression. In contrast, all of the other TM2 mutants, as well as N43A, gave Emax values between 76 and 150% of the wild-type value (see supplementary Table S3 for full details). None of the reductions were statistically significant, but A68G and C69A gave statistically significant increases in Emax of up to 50% (n = 6; p < 0.05, p < 0.01, respectively). Some of the mutants also showed statistically significant reductions in ACh-independent signaling relative to wild-type, namely N61A (83 ± 3%), S66A (75 ± 2%), A68G (83 ± 2%), L72A (76 ± 3%), and I74A (71 ± 8%); in the case of N61A, S66A, and L72A this might reflect reduced expression levels. C69A showed a slightly increased level of basal signaling (115 ± 4%; p < 0.05). Because the receptor-dependent component of the basal signal is about 25%, these effects were hard to quantitate more precisely.
Ala substitution of particular residues in the N-terminal (cytoplasmic) two thirds of TM2 caused reductions (defined as EC50,mutant/EC50,wild-type) in the signaling potency of ACh of up to 70-fold (Fig. 2b). The largest (>20-fold) effects were due to mutation of Ser66, Leu67, Ala70, and Ile74, whereas Tyr62 gave a 9-fold effect. We noted a periodic distribution of the mutational effect on PI signaling peaking at Asp71, whose mutation completely abolished signaling. The N43A mutant showed a 4-fold reduction in signaling potency (pEC50 = 6.27 ± 0.07).
|
|
This representation emphasizes the primary importance of the integrity of Asp71 for the mediation of PI signaling. Mutation of this residue led to zero signaling efficacy. Ten-fold reductions in efficacy also resulted from mutation of its neighboring residues Ala70, Leu67, and Ile74. Smaller reductions (3-5-fold) followed mutation of Tyr62, Leu64, Leu65, and Ser66. The increased ACh affinity observed for the L64A mutation (Fig. 2a) was not translated into increased PI potency (Fig. 2b), whereas the reduction in PI potency seen for S66A may be partly attributable to decreased expression (Fig. 1a) and decreased affinity for ACh (Fig. 2a). The N43A mutant showed an
10-fold reduction in its calculated signaling efficacy, also reflecting the failure of the increased binding affinity to be reflected in enhanced signaling potency.
Investigations of Mutant Receptor Stability—To further understand the origin of the reduced expression of functional binding sites found for the M79A and D71A mutants and the apparently slightly enhanced expression found for the L65A mutant, we performed studies of receptor thermal stability following a protocol modeled on that used by Rasmussen and co-workers (37).
Preliminary experiments showed that all of the receptors were stable at temperatures up to 40 °C, but that incubation at 60 °C led to 60% loss of binding of the wild-type receptor in 30 min. After further investigation, a temperature of 53 °C was chosen as suitable for simultaneous measurements on atropine-rescued D71A and M79A (with time points of 0, 2, 4, 8, 16, and 32 min) and wild-type and L65A (with time points of 0, 15, 30, 60, 120, and 240 min). Following the incubations, residual binding capacity was measured by incubation of the membranes with 10-9 M [3H]NMS for 60 min at 30 °C in the absence or presence of 10-6 M atropine to determine nonspecific binding. The results were fitted to single or double exponentials.
Most of the thermal inactivation curves were adequately fitted by single exponentials. An example is shown in Fig. 4. This shows that, at 53 °C, the D71A mutant decayed at a rate of 0.22 min-1, 55-fold faster than the wild-type receptor with a rate of 0.004 min-1, whereas the M79A mutant had an intermediate rate of 0.04 min-1. The inactivation rate of the L65A mutant was similar to that of the wild-type receptor. In the case of L65A, and wild-type receptor from atropine-pretreated cells, a small proportion (about 30%) of an initial faster inactivation phase was seen in some experiments (data not shown). The full set of results is presented in Table 1.
|
| DISCUSSION |
|---|
|
|
|---|
|
-carbon atom. The consequent effect on receptor function depends on interactions made by the deleted side chain. The extended ternary complex model of GPCR activation (41) provides a suitable predictive framework. Here the receptor exists in a conformational equilibrium between a ground state that does not productively bind the cognate G-protein, and an activated state that catalyzes GTP-GDP exchange. Ligands that bind to the receptor act by perturbing this equilibrium. Four canonical classes of mutant are predicted. First, a null phenotype when the target side chain does not make energetically significant interactions in either the ground or the activated state. Second, a stability phenotype when the side chain interactions are equally important in the ground state and the activated state. Their deletion reduces either the stability of the protein fold if they are solely intramolecular potentially leading to a reduction in functional expression, or the affinity of the ligand if they are intermolecular and involved in ligand anchoring. Third, an enhanced agonist affinity phenotype when the side chain makes intramolecular interactions in the ground state that are broken when the agonist-activated state is formed. A group of related effects ensues, not only a reduction in the stability of the ground state, but also an enhancement of agonist affinity, because binding energy no longer has to be used to rupture the intramolecular contacts. Ligand-independent basal signaling may also be promoted. Fourth, a reduced signaling phenotype when the side chain forms interactions in the activated state that are absent in the ground state. Their deletion increases the free energy difference between the ground and activated states thus reducing the signaling efficacy of the agonist-receptor complex. No major effects on ground state stability are expected, but there may be reductions in agonist affinity and basal signaling, depending on whether the interactions are intramolecular or intermolecular, with the ligand or G-protein. In more complex cases, a superposition of these fundamental phenotypes may occur.
In TM2 of the M1 mAChR, F63A, L65A, A68G, C69A, I73A, T76A, F77A, S78A, N80A, and L81A showed little functional effect, and can be classified as null positions. L72A (2.51), G75A (2.54), and M79A (2.58) showed pure functional expression phenotypes. N61A (2.40) and L64A (2.43) showed enhanced ACh affinity phenotypes; this characteristic was shared by N43A (1.50). L67A (2.46), A70G (2.49), and I74A (2.53) showed pure reduced signaling efficacy phenotypes. Finally D71A (2.50) showed major effects on signaling efficacy and expression level, consistent with a dual role for this key residue. The mutations caused a maximum 3-fold reduction in NMS and ACh affinity, probably ruling out a primary role in anchoring these ligands for any of the target residues.
To interpret the functional effects of the Ala substitutions, we have used a rhodopsin-based homology model of the M1 mAChR (7). Like the crystal structure of rhodopsin (42), this emphasizes the close spatial relationships that exist between TM1, TM2, and TM7.
Interestingly, Ala substitution of Pro415 and Tyr418 in the highly conserved NPXXY sequence of TM7 evoked an increased ACh affinity phenotype very similar to that caused by the mutations of Asn43, Asn61, and Leu64 (7). As summarized in Fig. 5a, this occurred without major changes in NMS affinity, and was accompanied by variable reductions in expression of functional binding sites and signaling efficacy. As shown in Fig. 6a, the clustering of these residues by phenotype is reflected in the context of the model by the existence of a hydrogen bond between the amide moiety of Asn61 and the hydroxyl group of Tyr418 (also proposed in the M3 mAChR (25)), reinforced by a van der Waals contact between Leu64 and the aromatic ring of Tyr418; homologous contacts exist in the high-resolution crystal structure of rhodopsin (6, 12, 42). Again guided by the rhodopsin crystal structure, Asn43 is modeled as forming a H-bond with the backbone carbonyl of Ser411 (TM7), and may be in van der Waals contact with the pyrrolidine ring of Pro415, as well as H-bonded to Asp71 (Fig. 6b). These proposed ground state contacts may explain the important destabilizing effect of mutating Asn43 on receptor expression. Some of them (e.g. with Asp71) may be maintained in the activated state to account for the accompanying reduction in signaling efficacy, whereas others (e.g. with TM7) may be destabilized. These observations suggest that there is a cluster of intra-molecular contacts between the inner sections of TM1, TM2, and TM7 that is broken or modified when ACh binds. Conformational linkage between this inner membrane domain and the ACh binding site may account for the increased ACh affinity observed after Ala substitution of the participating residues. This interpretation is consistent with a Cys substitution and cross-linking study suggesting that a rotation and translation of the inner part of TM7 occurs after agonist binding to the M3 mAChR (43).
|
|
|
In TM2, the dominant effects of Ala substitution on functional receptor expression level are in the extracellular portion of the sequence C-terminal to Asp71. The opposite is true in TM3 (30). Mutation of three residues reduced expression of receptor binding sites without further effects on their function. As shown in Fig. 8b, Leu72 projects toward TM1 at the level of Leu44. Gly75 also faces toward TM1 and may have a role in helix packing or kinking, which is a dominant feature in TM2 of rhodopsin. In contrast, Met79 faces toward TM7, where it may be in van der Waals contact with Tyr408 (7.43), a residue that contributes to anchoring ACh in the binding site, and whose mutation also reduces the functional receptor expression level (Fig. 5c) (7). The mutation of Met79 also reduced the thermal stability of the receptor. The absence of a mutational effect on signaling efficacy implies that the contacts between the extracellular parts of TM1, TM2, and TM7 are unlikely to be rearranged during the receptor activation cycle.
The group of residues that showed definite phenotypes is highly conserved throughout the cationic amine receptors. Position 2.50 (Asp71) is completely conserved. His is the only significant variant at position 2.40 (Asn61). The bulky aliphatic residues admit only conservative substitution by other bulky aliphatics or, occasionally, Phe; position 2.46 (Leu67) is almost completely conserved. The minimal side chains of positions 2.49 (Ala70) and 2.54 (Gly75) only tolerate other small substituents such as Gly, Ala, or Ser. These sequence positions are clearly subject to strong evolutionary constraints, suggesting that they fulfill a generic function. This view is supported by the outcome of a random mutagenesis study on the M3 mAChR in which a series of non-conservative substitutions isolated at positions 1.50, 2.40, 2.46, 2.49, 2.50, 2.51, 2.54, and 2.58 as well as several other positions in TM2 inactivated receptor expression or signaling (25). These considerations may reduce the functional information to be expected from further non-conservative mutagenesis of these residues (28).
In conclusion, the rhodopsin-based homology model of the M1 mAChR gives a good account of phenotypes from Ala substitution mutagenesis of TM2. We propose that ACh binding allosterically destabilizes contacts between the inner parts of TM1, TM2, and TM7. This is likely part of a cooperative rearrangement of inter-helical contacts that accompanies receptor activation. Formation of a direct H-bond between Asp71 (TM2) and Asn414 (TM7) may be a key to the formation of the activated state of the receptor. These residues are among the most highly conserved in the rhodopsin-like GPCRs. Bulky hydrophobic residues above and below Asp71 may help to direct and stabilize the formation of this bond. Contacts at the top of TM2 may help to stabilize the receptor fold. These appear not to be modified by receptor activation.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables S1-S3 and Fig. S1. ![]()
1 Supported by a Medical Research Council post-graduate studentship. ![]()
2 To whom correspondence should be addressed. Tel.: 44-208-816-2057; Fax: 44-208-906-4477; E-mail: ehulme{at}nimr.mrc.ac.uk.
3 The abbreviations used are: mAChR, muscarinic acetylcholine receptor; ACh, acetylcholine; TM, transmembrane domain; NMS, (-)-N-methylscopolamine; PI, phosphoinositide. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Lebon, C. J. Langmead, B. G. Tehan, and E. C. Hulme Mutagenic Mapping Suggests a Novel Binding Mode for Selective Agonists of M1 Muscarinic Acetylcholine Receptors Mol. Pharmacol., February 1, 2009; 75(2): 331 - 341. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |