|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 281, Issue 32, 23066-23074, August 11, 2006
Heparin Modulates the 99-Loop of Factor IXa
EFFECTS ON REACTIVITY WITH ISOLATED KUNITZ-TYPE INHIBITOR DOMAINS*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The major inhibitor of fIXa in plasma is antithrombin, whose reactivity with fIXa essentially requires heparin (7-9). Heparin is known to bind to antithrombin and sterically alter its conformation to allow this serpin to react with its target (10-13). Heparin also binds to fIXa (14) allowing long chains of heparin to additionally catalyze the interaction of fIXa with antithrombin via the formation of bridged complexes where heparin acts as a "template." Recently, we have shown that low molecular weight heparin binding to fIXa enhances reactivity of fIXa with the Kunitz-type inhibitor BPTI (15), suggesting that oligosaccharide binding can also allosterically modulate the fIXa active site region. In this study we examine in greater detail the ability of heparin to modulate fIXa reactivity toward several isolated Kunitz-type inhibitor domains. We show that the modulatory effect of heparin can be completely abrogated by mutating a single amino acid residue in the 99-loop region of the extended fIXa active site cleft outside of the heparin binding exosite.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
, factor VIIa, factor XIa, and the factor X activator from Russell's Viper venom were purchased from Hematologic Technologies Inc. (Essex Junction, VT). Recombinant soluble tissue factor (the extracellular domain of tissue factor) was expressed and purified from bacteria as previously described (16). Factor Xa was prepared from plasma-derived factor X as previously described (17). Enoxaparin (Lovenox®) was purchased from Aventis Pharmaceuticals (Bridgewater, NJ). Purified heparin-derived oligosaccharides of 6, 10, 14, and 18 saccharide units (H6, H10, H14, and H18) were prepared and characterized essentially as described (18-20) and were a generous gift of Dr. Steven T. Olson, University of Illinois, Chicago, IL. Bovine serum albumin (Fraction V, fatty acid free) was from Calbiochem (La Jolla, CA), and ethylene glycol was from Fisher Scientific. The chromogenic substrate CBS 31.39 (CH3SO2-D-LGR-pNA) was purchased from Diagnostica Stago (Parsippany, NJ). All other reagents were of the highest quality available.
Construction and Expression of Recombinant InhibitorsAppropriate expression clones encoded for: BPTI (59 amino acids) (21, 22), PN2-KPI (61 amino acids corresponding to residues 285-344 of Protease Nexin-2) (23), TFPI-K1 (58 amino acids corresponding to residues 50-107 of TFPI) (24), and TFPI-K2 (59 amino acids corresponding to residues 121-178 of TFPI) (24). Each construct was directionally cloned into pET11a (Novagen) and verified by sequencing. Inhibitors were expressed as inclusion bodies in Escherichia coli strain BL21(DE3). Transformed bacterial cells were first grown to log phase at 37 °C in TB media containing 50 µg/ml carbenicillin. Protein expression was induced by addition of isopropyl 1-thio-
-D-galactopyranoside to 0.5 mM (0.1 mM for TFPI-K1) and the cells were allowed to continue growing for 4 h at 37°C.
Inclusion bodies were isolated essentially as described (25) and solubilized with 6 M guanidine HCl containing 20 mM dithiothreitol, 50 mM Tris-HCl, pH 8.0, and 1 mM EDTA to obtain a total protein concentration of roughly 20 mg/ml. The solution was then clarified by centrifugation (16,000 x g for 30 min) and oxidative refolding (26) of each protein preparation was performed by rapid dilution into 20 volumes of buffer containing 50 mM Tris-HCl, pH 8.0, 1 M guanidine HCl, 1 mM EDTA, 2.5 mM oxidized glutathione (Sigma), and 1 mM dithiothreitol. The diluted protein solution was incubated at room temperature for 6 h with slow stirring for completion of protein refolding followed by exhaustive dialysis into an appropriate buffer for ionexchange chromatography.
Construction and Expression of Wild-type and Mutant fIXThe coding sequence for wild-type fIX in pBR322 (27) was a generous gift of Dr. Earl Davie (University of Washington). The fIX coding sequence was removed into the mammalian expression vector pcDNA3 (Invitrogen) and sequenced to verify the correct orientation. This construct (pFN04) was used for expression of wild-type fIX as well as PCR-based mutagenesis (28) to generate fIXK98A essentially as previously described for constructing fVII mutants (29). Expression constructs were transfected into human 293 cells using Lipofectin ® (Invitrogen) and high expressing clones isolated by limiting dilution.
Protein PurificationsPurification of refolded BPTI was accomplished by ion-exchange chromatography using Mono-S HR 5/5 (Amersham Biosciences) in 20 mM Tris-HCl, pH 8.0. The column was developed with a 0-1 M NaCl gradient and BPTI eluted as a single peak at roughly 0.43 M NaCl. The specific activity of recombinant BPTI preparations was equivalent to or better than that of commercial preparations of aprotinin (not shown). Purification of refolded PN2-KPI and TFPI-K2 was accomplished by ion-exchange chromatography using Mono-Q HR 5/5 (Amersham Biosciences) in 20 mM MES, pH 6.0. In both cases the column was developed with a 0-0.5 M NaCl gradient. PN2-KPI eluted at roughly 90 mM NaCl and TFPI-K2 eluted at roughly 50 mM NaCl. Purification of refolded TFPI-K1 was accomplished by affinity chromatography over a trypsin-agarose column. Trypsin-agarose was prepared by coupling 20 mg of L-1-tosylamido-2-phenylethyl chloromethyl ketone-treated trypsin (Worthington Biochemicals) to 2.5 ml of Affi-Gel 10 in 100 mM MOPS, pH 7.4, 10 mM CaCl2, and 100 µg/ml leupeptin overnight at 4 °C. This was followed by blocking non-reacted sites with 1 M ethanolamine-HCl, pH 8.0. Refolded TFPI-K1 was loaded onto the prepared trypsin-agarose column in 50 mM Tris-HCl, pH 7.5, 100 mM NaCl. The column was extensively washed with the same buffer before eluting the inhibitor with 10 mM HCl, pH 2.1, 100 mM NaCl. The pH of the eluted fractions was immediately neutralized with 0.02 volumes of 2 M Tris buffer. All inhibitors were judged >95% pure by SDS-PAGE.
Recombinant human wild-type fIX and fIXK98A were isolated from 293 cell supernatants using a combination of ionexchange and heparin affinity chromatography. A 10-fold concentrate of cell supernatant was diluted 2-fold with deionized water to reduce the ionic strength before loading a 150-ml DEAE FF Sepharose (Amersham Biosciences) column equilibrated in 25 mM sodium citrate, pH 6.0, 33 mM NaCl, and 1 mM benzamidine. After loading, the column was extensively washed in the same buffer before elution of the fIX protein with a 0.033-0.4 M NaCl gradient over 10 column volumes. The fIX protein peak was identified by clotting activity, pooled, and dialyzed versus 50 mM Tris-HCl, pH 7.5, 100 mM NaCl before heparin affinity chromatography using either POROS ® HE2 (Applied Biosystems) or HiPrepTM Heparin FF 16/10 (Amersham Biosciences) and eluting with a NaCl gradient. Wild-type fIX and fIXK98A both eluted as single peaks at roughly 0.46 M NaCl.
Wild-type and mutant fIX proteins were activated with the purified factor X activator from the venom of Russell's viper, which also cleaves fIX after Arg180 to generate active enzyme (fIXa
). The activated enzyme was purified away from the venom protease by subsequent heparin affinity chromatography essentially as described above using HiTrapTM Heparin HP (Amersham Biosciences). Although the activation peptide remains attached to the light chain of fIXa
, this enzyme retains 100% amidolytic activity compared with fIXa
(activation peptide proteolytically removed) and is comparable with fIXa
in kinetics of inhibition by Kunitz-type inhibitors (see "Results"). Unless indicated otherwise, the fIXa
form was used in experiments.
Clotting AssaysCoagulant activities of wild-type and mutant fIXa proteins were assayed by a standard single-stage clotting assay using a Coag-a-mate XM (Organon Teknika) coagulometer, fIX-deficient plasma (George King Biomedical), and APTT Reagent (Sigma).
Reactive Site Titration of InhibitorsThe active concentration of inhibitor preparations was determined by reactive site titration essentially as described (23) using 10 nM active site-titrated trypsin (30) and S-2222 substrate (Chromogenix, Milano, Italy) to measure residual trypsin activity after a 15-min incubation period. This method assumes a 1:1 stoichiometry of inhibitor and trypsin. Amino acid analysis performed on an initial PN2-KPI preparation indicated an equivalent concentration as that determined by reactive site titration (not shown).
Active Site Titration of fIXa EnzymesActive concentrations of wild-type and mutant fIXa preparations were determined by active site titration using biotin-EGR-ck (Hematologic Technologies Inc.) essentially as described (31). Briefly, wild-type or mutant fIXa (roughly 5 µM) were incubated with 150 µM biotin-EGR-ck in 50 mM Tricine, pH 8.0, 200 mM NaCl, 10 mM CaCl2, and 30% ethylene glycol for 24 h at room temperature. Biotin-EGR-fIXa was then quantified by enzyme-linked immunosorbent assay using a goat anti-fIX polyclonal capture antibody and alkaline phosphatase-conjugated streptavidin. These were compared with a standard line made using native fIXa
to determine active concentrations for all fIXa enzyme preparations.
|
![]() |
Slow-binding Enzyme Inhibition AnalysisKinetic parameters for slow-binding inhibition were obtained using Equations 2-5 as derived by others (33, 34) and are described in Schemes Ia and Ib (Fig. 1). In both of these schemes the non-covalent enzyme-inhibitor complex, EI, isomerizes into EI*; either after significant formation of the EI complex (Scheme Ia) or without significant accumulation of the EI complex (Scheme Ib). In both cases, the overall inhibition constant describing generation of EI* is defined as Ki*, which is equivalent to Ki,eq obtained from equilibrium experiments (above). In cases following Scheme Ia, the parameter Ki* can be further broken down to obtain Ki, which describes the establishment of the initial "loose" EI complex, k5 and k6 (see below).
Final reaction conditions were the same as described above for equilibrium studies. In these assays, however, the fIXa was preincubated with or without enoxaparin for 15 min at 25 °C in the reaction mixture before the combined addition of inhibitor and substrate at time 0. The absorbance at 405 nm was then monitored for up to 30 min in a Thermo-Max microplate reader (Molecular Devices) set at 25 °C to monitor substrate hydrolysis using KINEMAX software (written and kindly provided by Dr. Jolyon Jesty, SUSB, Stony Brook, NY). Data for each generated curve were fitted with the following integrated rate equation describing slow-binding inhibition,
![]() |
where A is the absorbance at 405 nm at any time, t. Fits of progress curves with Equation 2 yield values for Ao (the initial absorbance at t = 0), vo (the initial rate of substrate hydrolysis), vs (the steady-state rate of substrate hydrolysis), and kobs (the apparent first-order rate constant for inhibition).
For analyses using Scheme Ia, values of k6 (the reverse rate constant for EI* isomerization) were determined from progress curves above using the following relationship.
![]() |
Values of k5 and initial Ki (defined as k4/k3) were then obtained from secondary plots of kobs versus I using the following hyperbolic equation.
![]() |
For analyses using Scheme Ib, values of k10 were obtained from progress curves also using Equation 3. However, in these cases vo does not vary with inhibitor concentration and a plot of kobs versus I yields a straight line, indicating conditions where Ki (1+ S/Km) >> I. Thus for Scheme Ib EI formation is insignificant and EI* can be considered formed directly from E + I. For these cases the following linear equation is applicable for obtaining an estimate of k9, the apparent second-order on-rate constant,
![]() |
where the y intercept reflects k10 and the slope of the line is equal to k9/(1 + S/Km). Alternatively, k9 can be obtained from Ki*, which is equivalent to Ki,eq in Equation 1, using the relationship k9 = k10/Ki*.
Although fits with Equations 4 or 5 yield estimates of k6 or k10, respectively, the values reported herein were obtained from Equation 3 using the more accurate fits of progress curves to Equation 2 and then verified in fits with Equations 4 or 5. Experimental values of S as well as experimentally determined values of Km (defined in the traditional manner as (k2 + k7)/k1 in Schemes Ia and Ib) were used as necessary in all fitting procedures. All regression procedures were performed using Slide-WritePlus 6.0 (Advanced Graphics Software), which uses the Levenberg-Marquardt algorithm.
|
| RESULTS |
|---|
|
|
|---|
|
are compared in Fig. 2. As expected, fIXa exhibited remarkable specificity toward these inhibitors despite their high homology. Of the inhibitors examined, PN2-KPI showed the highest level of reactivity (Ki,eq = 10 µM), followed by TFPI-K2 (Ki,eq = 336 µM), BPTI (Ki,eq > 500 µM), and TFPI-K1 (Ki,eq > 1mM). Consistent with our previous observations, enoxaparin was able to enhance the reactivity of fIXa with BPTI more than 10-fold (Ki,eq = 46 µM). Surprisingly, however, this same level of enhancement by enoxaparin was not observed with any of the other inhibitors examined; TFPI-K2 and PN2-KPI each showed only a small, but consistent, enhancement in reactivity with enoxaparin (1.7- and 1.4-fold, respectively; Ki,eq values of 203 and 7 µM) and TFPI-K1 showed no measurable enhancement in reactivity with enoxaparin. The highly basic nature of BPTI compared with the other isolated Kunitz domains along with its ability to bind to heparin (albeit weakly; Kd = 172 µM (15)) raised the potential that enoxaparin, although short (15 saccharide units; H15), may retain some capacity to facilitate the interaction of BPTI with fIXa via a bridging-type mechanism. Although unlikely based on previous equilibrium kinetic studies and the level of enoxaparin used in these experiments (10 µM; or 0.06 x Kd for BPTI binding versus 78 x Kd for fIXa binding), this issue was examined by using increasing concentrations of enoxaparin as well as progressively smaller heparin oligosaccharides; H18, H14, H10, and H6 (18-20). As shown in Fig. 3, the typical bell-shaped profile for bridging-type mechanisms was not observed at enoxaparin concentrations ranging from 1 nM to >100 µM.In addition, and of greater significance, is the observation that progressively smaller oligosaccharides do not lose the ability to enhance reactivity of fIXa. These results along with previous kinetic studies support the ability of heparin to modulate fIXa reactivity via a mechanism other than bridging, and are consistent with allosteric modulation of the fIXa protease domain.
|
Based on this hypothesis, we examined a mutant of fIXa in which Lys98 was mutated to Ala (fIXK98A). Wild-type and mutant forms of fIXa were expressed in human 293 cells and purified to >95% homogeneity as judged by SDS-PAGE (Fig. 5A). The fIXK98A mutant was found to retain 100% clotting activity compared with wild-type fIXa (Fig. 5B) and upon activation retained near normal amidolytic activity toward CBS 31.39 substrate (Table 1).
|
(not shown) and are well described by Equation 2. Values of k6 subsequently obtained using Equation 3 showed no significant differences in the absence or presence of enoxaparin; k6 = 2.7 (±1.5) x 10-3 s-1 without enoxaparin and 1.2 (±0.2) x 10-3 s-1 with enoxaparin.
|
From fits of the progress curves in Fig. 6, A and B, with Equation 2, secondary plots of kobs versus BPTI for WT fIXa show little or no curvature (Fig. 6C). This suggests a very high initial Ki value and an indeterminate value for k5. These observations as well as the lack of change in vo with increasing inhibitor concentrations (cf. Fig. 6, A and B) are diagnostic for Scheme Ib. The data of Fig. 6C were thus fitted with Equation 5 to obtain estimates for the apparent second-order rate constant; k9 = 5.4 M-1 s-1 without enoxaparin and 23 M-1 s-1 with enoxaparin. The value obtained above for k6 (k10 in Scheme Ib) is consistent with the values obtained for k10 using Equation 5; k10 = 1.3 x 10-3 s-1 without enoxaparin and 1.2 x 10-3 s-1 with enoxaparin.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
38 Å and that of a hexasaccharide is expected to be
23 Å at full extension. This is compared with the measured length of the identified heparin binding site on fIXa (
35 Å) based on mutational studies (14).
|
It is important to note that the reduction in the effect of heparin with the fIXaK98A mutant is not merely due to reduced heparin binding because fIXaK98A retained the ability to bind to heparin-Sepharose and eluted at the same salt concentration as WT fIXa during purification procedures. In support of this, preliminary experiments performed by titrating fIXaK98A with enoxaparin in the presence of 100 µM BPTI yielded results consistent with the high nanomolar affinity for enoxaparin previously observed with fIXa
(Kd = 128 nM (15)). No further increase in inhibition of fIXaK98A by BPTI was observed when the enoxaparin level was increased from 1 to 10 µM (not shown). Because all experiments were conducted using 10 µM enoxaparin, it seems reasonable to assume that fIXaK98A was saturated in these experiments.
Whereas the data for fIXa with BPTI in Fig. 6C are consistent with Scheme Ib, this scheme is essentially a simplified version of Scheme Ia with a very large value for Ki (33, 34). Examination of the Scheme Ia kinetic constants in Table 3 for fIXa with PN2-KPI as well as those for fIXaK98A with both inhibitors shows fairly consistent values for k5 and k6. This suggests that once formed, isomerization of the EI complex to EI* is essentially the same for any of the enzyme-inhibitor pairs examined. Based on this it seems reasonable to tentatively extend the value of k5 to the WT fIXa-BPTI pair and re-examine the data of Fig. 6C with respect to Scheme Ia. This results in fits to Equation 4 shown in Fig. 6C as dashed lines and yields initial Ki values of 2.1 x 10-3 M and 4.1 x 10-4 M in the absence and presence of enoxaparin, respectively. This 5-fold difference in Ki along with
2-fold difference in k6 (k10) for this enzyme-inhibitor pair (Table 2) would seem to account for the
10-fold effect of enoxaparin observed in Ki,eq (Table 3). Unfortunately, the rather large errors in the values of k6 (k10) preclude definitive conclusions concerning potential effects of enoxaparin on this rate constant.
Regardless, these results are consistent with a very weak initial interaction of fIXa with BPTI. Enoxaparin binding to the heparin binding exosite in fIXa at least in part acts to allosterically modulate the 99-loop of fIXa in a manner that facilitates this initial interaction. The results with BPTI are in contrast to the interaction of fIXa with TFPI-K1, TFPI-K2, or PN2-KPI. The former shows no effect of heparin binding, whereas the latter two show only small (but reproducible) responses to heparin binding. These results would suggest that these three Kunitz inhibitor domains are not as greatly hindered as BPTI by the 99-loop of fIXa. Indeed, mutation of Lys98 enhanced the inhibition by TFPI-K1 roughly 8-fold, TFPI-K2 roughly 9-fold, and by PN2-KPI only 6-7-fold. Interestingly, however, the effect of enoxaparin on fIXaK98A reactivity toward both TFPI-K1 and TFPI-K2 was unchanged, whereas that for PN2-KPI was completely abrogated. Thus, there is some disparity in the role of Lys98 and additional effects of heparin binding are likely.
Previous studies have demonstrated the ability of full-length TFPI to inhibit fIXa with an affinity of roughly 300 nM (40). It is thus interesting that neither the isolated TFPI-K1 nor isolated TFPI-K2 domains retain very high reactivity with fIXa. Preliminary studies indicated that our preparations of isolated inhibitor domains were of good quality and activity. This would seem to indicate that the folding and activity of these inhibitor domains are correct. Taken as a whole, the results here would seem to imply that either the connecting regions of TFPI are important for its activity toward fIXa or other structural/topological elements are involved. Nonetheless, it is likely that the TFPI-K2 domain is the inhibitor domain of TFPI most likely responsible for interaction with fIXa. Although the isolated third Kunitz repeat of TFPI was not examined in this study, this domain is likely not an active inhibitor domain per se (41).
The ability of heparin to allosterically modulate the active site of fIXa demonstrates the ability of fIXa to respond to binding at this exosite. This observation is made more interesting by the observation by Sheehan et al. (42, 43) that the heparin binding exosite in part represents a fVIIIa interactive site on fIXa. Thus, whereas heparin binding in itself may inhibit fIXa coagulant activity via steric hindrance of fVIIIa binding, fVIIIa binding to this exosite may inversely act to allosterically modulate fIXa in an as yet undefined manner. Although previous studies found no effect of fVIIIa on fIXa inhibition by the isolated PN2-KPI domain (44), this is consistent with our observations here using enoxaparin in place of fVIIIa and does not preclude potential effects of fVIIIa toward other inhibitors. It thus remains possible that occupation of this exosite by fVIIIa, like heparin, results in inhibitor-specific modulatory effects. Alternatively, the modulating effect of heparin on the 99-loop may be simply due to heparin-specific electrostatic forces that are introduced by heparin binding in close proximity to Lys98 (cf. Fig. 4). These forces may or may not be mimicked by fVIIIa binding. Further clarification of these issues must await future studies.
The allosteric modulation of fIXa by enoxaparin is somewhat reminiscent of the effect of thrombomodulin on the interaction of thrombin with BPTI, where binding of thrombomodulin alters the conformation of one of the specificity loops (60-loop) at the mouth of the active site of thrombin, resulting in enhanced reactivity with BPTI (45). Interestingly, in that study the chondroitin sulfate moiety of thrombomodulin (which binds to the heparin binding site on thrombin; anion binding exosite 2) further enhanced the inhibition of an E192Q mutant of thrombin but not wild-type thrombin. Comparison of the sequence of fIXa with thrombin reveals that fIXa contains a Gln at position 192, similar to the thrombin E192Q mutant and fXa. In addition, previous studies have revealed that mutation of the homologous residue in fVIIa (Lys192) to Gln enhances its reactivity with BPTI (29), and that mutation of Gln192 in fIXa altered its reactivity toward TFPI (40). Neither of these studies examined the potential effect of heparin on reactivity with these Kunitz inhibitors. These studies are intriguing and seem to implicate residue 192 in potential heparin responsiveness along with the 99-loop. The investigation of potential interplay between residue 192 and the 99-loop of fIXa in its response to heparin binding is beyond the scope of the present study, however, and must await future studies.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Table S1 and Figs. S1-S3. ![]()
1 To whom correspondence should be addressed: Biomedical Research Lab C8, The University of Texas Health Science Center at Tyler, 11937 US Hwy. 271, Tyler, TX 75708. Tel.: 903-877-7678; Fax: 903-877-7679; E-mail: Pierre.Neuenschwander{at}uthct.edu.
2 The abbreviations used are: FIXa, factor IXa; BPTI, basic pancreatic trypsin inhibitor; PN2-KPI, the isolated Kunitz-type inhibitor domain from protease nexin-2; TFPI, tissue factor pathway inhibitor; TFPI-K1, the isolated first Kunitz-type inhibitor domain of TFPI; TFPI-K2, the isolated second Kunitz-type inhibitor domain of TFPI; WT, wild-type; MES, 4-morpholineethanesulfonic acid; MOPS, 4-morpholinepropanesulfonic acid; Tricine, N-[2-hydroxy-1,1-bis(hydroxymethyl)ethyl]glycine. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. Wu, C. Eigenbrot, W.-C. Liang, S. Stawicki, S. Shia, B. Fan, R. Ganesan, M. T. Lipari, and D. Kirchhofer Structural insight into distinct mechanisms of protease inhibition by antibodies PNAS, December 11, 2007; 104(50): 19784 - 19789. [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 |