| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Endocrinology and Internal Medicine (C.A.C., A.A.G., L.M.M., J.A.M., C.E.F.) and Faculty of Medicine and Center for Genomics and Bioinformatics (A.G., T.O.P.-A.), Pontificia Universidad Católica de Chile, 114-D Santiago, Chile; Department of Nephrology (E.T.L., M.d.P.H., M.P.R.), Hospital San Juan de Dios, Santiago, Chile; and Division of Endocrinology (D.G.R., C.E.G.-S.), University of Mississippi Medical Center and the G. V. Montgomery Veterans Affairs Medical Center, Jackson, Mississippi 39216
Address all correspondence and requests for reprints to: Carlos E. Fardella, M.D., Department of Endocrinology and Metabolism Faculty of Medicine, Pontificia Universidad Católica de Chile, Lira 85, 5th Floor Santiago, Chile. E-mail: cfardella{at}med.puc.cl.
| Abstract |
|---|
|
|
|---|
AAC) and a single nucleotide substitution C
T in intron 3. Using site-directed mutagenesis, we generated a mutant 11ßHSD2 cDNA containing the Asp223Asn mutation. Wild-type and mutant cDNA was transfected into Chinese hamster ovary cells and enzymatic activities were measured using radiolabeled cortisol and thin-layer chromatography. The mRNA and 11ßHSD2 protein were detected by RT-PCR and Western blot, respectively. Wild-type and mutant 11ßHSD2 protein was expressed in Chinese hamster ovary cells, but the mutant enzyme had only 6% of wild-type activity. In silico 3D modeling showed that Asp223Asn changed the enzymes surface electrostatic potential affecting the cofactor and substrate enzyme-binding capacity. The single substitution C
T in intron 3 (IVS3 + 14 C
T) have been previously reported that alters the normal splicing of pre-mRNA, given a nonfunctional protein. These findings may determine the full inactivation of this enzyme, explaining the biochemical profile and the early onset of hypertension seen in this patient. | Introduction |
|---|
|
|
|---|
The 11ßHSD2 is crucial in regulating mineralocorticoid action. Because the concentrations of cortisol in human blood are nearly 1000 times higher than that aldosterone, the discovery that glucocorticoids could bind to mineralocorticoid receptor posed the problem of how the aldosterone could exert specific and separate action from cortisol. It is now generally agree that 11ßHSD2 metabolizes active cortisol in cortisone and protects the mineralocorticoid receptor from glucocorticoid occupancy and thus defends the mineralocorticoid receptor from being overwhelmed by glucocorticoids (2).
The 11ßHSD2 enzyme is most abundantly expressed in tissues such as kidney and colon but is also expressed in pancreas, placenta, prostate, and gonads (3, 4). The enzyme catalyzes only the oxidase reaction (cortisol to cortisone), uses NAD+ as a cofactor, and has a low Michaelis-Menten constant (Km; 10100 nM) (5, 6). Thus, 11ßHSD2 found in the renal collecting ducts possesses all the features required to protect the mineralocorticoid receptor from occupancy by endogenous glucocorticoids: location in mineralocorticoid target cells, very high affinity for its substrate, and an ability to reduce steroids drastically to irreversible dehydrogenation. For these reasons, when there is a congenital absence of 11ßHSD2 (3) or the enzymes activity is inhibited by licorice (7) or carbenoxolone (8), the protective mechanism fails and cortisol gains inappropriate access to mineralocorticoid receptor, resulting in hypertension and hypokalemia.
The human HSD11B2 gene is in chromosome 16q22, consists of five exons spanning about 6.2 kb, and encodes a protein of 405 amino acids, with a predicted molecular mass of 42 kDa, which is only 14% identical to the type I enzyme (4, 9). The 11ßHSD2 belongs to short-chain dehydrogenase/reductase family, which had a 1530% sequence identity among them. However, the members of the family share a great similarity in their tertiary structures. The folding pattern of these proteins is a combination of
-helix and ß-sheets. Two important regions in the 11ßHSD2 are the cofactor-binding site (Gly-XXX-Gly-X-Gly motif) and steroid-binding site (Tyr-XXX-Lys motif). Initial modeling of the enzyme was performed with human 17 estradiol dehydrogenase and glucose dehydrogenase (Bacillus megaterium) whose tertiary structures correspond to patterns of the short-chain dehydrogenase/reductase family (4, 10, 11).
The syndrome of apparent mineralocorticoid excess (AME) is due to inactivating mutations in the HSD11B2 gene (5, 6). It is a rare disorder with an autosomal recessive inheritance and characterized by childhood onset of hypertension, hypokalemic alkalosis, and low plasma renin and aldosterone levels (12, 13, 14). This syndrome is also characterized by a high plasma cortisol to cortisone ratio or a high urinary tetrahydrocortisol plus allotetrahydrocortisol to tetrahydrocortisone (THF+aTHF/THE) ratio. Serum cortisol concentrations alone do not aid the diagnosis, although the plasma cortisol half-life is increased (13, 15, 16). In this study we described the clinical and molecular characterization of a patient with two homozygous mutations in the HSD11B2 gene, who associated the typical phenotype found in the syndrome of AME.
| Patients and Methods |
|---|
|
|
|---|
A 4-yr-old Chilean male was seen by one of us at the age of 3 yr. He was presented with increased thirst, polyuria, failure to thrive, height and weight less than the third percentile, hypertension (140/90 mm Hg), hypokalemia (serum potassium, 3.0 mEq/liter), and bilateral nephrocalcinosis. Plasma renin activity (PRA) and serum aldosterone (SA) were undetectable. Serum cortisol was normal (7.7 µg/dl), with a low cortisone value (0.53 µg/dl; normal value, 3.28 ± 1.09) in the presence of a high serum cortisol to cortisone (F/E) ratio of 14.5 (normal value, 2.8 ± 1.1). The patient was the only child of a 16-yr-old woman. The father and grandparents were unknown, and the mother refused to provide any information about the father. The patient and his mother came from a small and isolated village in northern Chile. The mother was normotensive and her PRA (3.1 ng/ml·h), SA (6.3 ng/dl), and potassium (4.0 mEq/liter) were normal, but the F/E ratio was slightly high (5.9). Human kidney proteins were obtained from the cortical region of renal tissue obtained from a biopsy of a healthy patient prepared for renal transplant. Informed consent was obtained from all participants in this study according to the guidelines of the Declaration of Helsinki, and the protocol approved by the Research Committee of the School of Medicine at Pontificia Universidad Católica de Chile.
Biochemical methods profile
SA was measured by RIA using a commercial kit (Diagnostic Products, Los Angeles, CA). Intra- and interassay variation for SA were 5.1% and 7.1%, respectively, and the normal range was 116 ng/dl. The PRA was determined as previously described (17), its intra- and interassay variations were 6.1% and 8.2%, respectively, and the normal range was 12.5 ng/ml·h (18). The lower limit of PRA determination was 0.1 ng/ml·h (19). Steroids were extracted with methane-dichloride. Cortisol was measured by RIA and cortisone by ELISA using highly specific antibody antihuman cortisone. Cortisol and cortisone samples were also measured by HPLC to confirm the values. Samples were analyzed by liquid chromatograph Waters 441 equipped with a reverse-phase column µBondapack C-18 (Waters, Milford, MA). Samples and steroids mixtures (containing known amounts of metabolites to be measured) were analyzed to 254 nm. HPLC mobile phase consists in isocratic gradient methanol:water (60:40) with 1% tetrahydrofuran (Fischer Chemicals, Pittsburgh, PA), flow rate 1.2 ml/min, and 1500
internal pressure. Retention times for cortisone, cortisol, and 6
-methylprednisolone (internal standard; Steraloids, Newport, RI) were 5.6, 6.7, and 9.2 min, respectively. Sensitivity of the chromatographic method for cortisone and cortisol was 0.5 µg/dl.
Analysis of genomic DNA
Genomic DNA of the patient and his mother was isolated from peripheral leukocytes. PCR amplification of five exons of HSD11B2 gene was performed in three segments using the oligonucleotides primers shown in Table 1
. Cycling conditions were 2 min at 94 C and 30 cycles of 1 min at 94 C, 1 min 30 sec at 60 C and 2 min at 72 C, and final extension 7 min at 72 C. Amplified gene fragments were purified in low-melting-point agarose and purified with Wizard PCR preps (Promega Corp., Madison, WI). Sequence analysis of the HSD11B2 gene was performed by fluorescent dideoxy chain terminator method in ABI Prism-377 DNA sequencer (Applied Biosystems, Foster City, CA). Sequences were matched with the published HSD11B2 gene sequence (GI 9989705) with the BLAST software (20). Restriction analysis with the BsmFI enzyme was performed to determinate the mutation haplotype of the patient and his mother. BsmFI recognizes and cuts the sequence GGGGACN14/N, which occurs only in the wild-type allele.
|
Full-length human HSD11B2 complementary DNA (GI 4504498) was cloned in a pCR3 expression vector (Invitrogen, Carlsbad, CA). The 11ßHSD2 Asp223Asn mutation in HSD11B2 was generated with the QuickChange site-directed mutagenesis kit (Stratagene, La Jolla, CA) using the following primers: GGGAGCCCAGCGGGGAACATGCCATATCCG (sense) and CGGATATGGCATGTTCCCCGCTGGGC-TCCC (antisense) (italicized base indicates the changed one). Mutant clones were screened by amplifying the mutation-flanking region by PCR with the following primers: (h-11HSD2-S-614 AGGCCACAATGAAGTAGTTGCT and h-11HSD2-AS-858 GTTCACAGCTGAATGTGTCCAT) and Advantaq Plus DNA polymerase (CLONTECH Laboratories, Inc., Palo Alto, CA), followed by enzymatic digestion with BslI. The fragment pattern was analyzed by electrophoresis in 4% metaphor agarose. The mutagenized cDNA was sequenced in its entirety to confirm the mutation and ensure that no other bases had been changed. Chinese hamster ovary K1 (CHO-K1) cells, grown on
MEM supplemented with 10% fetal bovine serum were split 1:16 and transferred to 6-well plates. Subsequently, 15 x 105 CHO-K1 cells were transfected with 3 µg plasmid DNA using Lipofectamine PLUS (Invitrogen) and grown in serum-free media for 3 h. The media were replaced with complete media and cells cultured for an additional 18 h before being subjected to the experimental protocols described below.
RT-PCR of HSD11B2 mRNA
Total RNA was extracted with TRIZOL LS reagent (Invitrogen). RNA (3 µg) was reverse transcribed using the SuperScript II kit (Invitrogen). Transfection efficiency was verified by RT-PCR within the exponential phase of the reaction, as described previously (21, 22), with the h-11HSD2-S614 and h-11HSD2-AS-858 primers.
11ßHSD2 immunodetection
Western blot analysis was performed with rabbit immunopurified polyclonal antiserum (anti-h11ßHSD2) in 1:2000 dilution, donated by Drs. Jonathan Seckl and Roger Brown (Edinburgh University, Edinburgh, UK) Briefly, 40 µg total proteins from the transfected and nontransfected cells, and 10 µg of the human kidney homogenate proteins were denatured and electrophoresed in 12% SDS-PAGE 2 h at 90 V and transferred to polyvinylidene fluoride membranes for 45 min at 100 V. Membranes were blocked with nonfat milk and incubated with polyclonal antiserum raised against human 11ßHSD2 for 18 h at 4 C, washed, and finally developed with anti-IgG-rabbit-horseradish-peroxidase (Bio-Rad Laboratories, Inc., Hercules, CA) (1:3000) using a chemiluminescent detection kit (Perkin-Elmer, Norwalk, CT).
11ßHSD2 enzyme activity assays
11ßHSD2 activity was analyzed by incubating transfected CHO cells with variable concentrations of [3H]-cortisol for 30 min and 2 h. Steroids were extracted and analyzed by thin-layer chromatography (TLC) in Silicagel 60 F-254 (Merck, Darmstadt, Germany), with a solvent mixture of hexane:ethylacetate (50:50). TLC was exposed to autoradiography films at -80 C for 1 wk and then developed to localize the standards and sample migration spots. Cortisol or cortisone migration spots on TLC plates were scrapped and extracted with isopropyl alcohol, and then activities of radiolabeled steroids were determined by scintillation spectrometry. To determinate the apparent kinetics parameters, Km, and maximal velocity (Vmax), whole CHO cells were incubated 30 min at 37 C with 2, 5, 10, 50, and 100 nM [3H]-cortisol and analyzed by Lineweaver-Burk plot (Hyper.exe version 1.02a; JS Esterby, Liverpool, UK). Cells were incubated for times (30 min) that resulted in less than 20% substrate conversion to ensure first-order kinetics as described previously (21).
Sequence alignment and protein modeling
11ßHSD2 sequence alignment was performed with the Vector NTI V.7.1 software (23, 24, 25) using human 17 estradiol dehydrogenase as query protein. Modeling of the core region (residues 80327 San Diego, CA) was performed using Modeler from Insight II software (Accelrys Inc., San Diego, CA) in a Silicon Graphics Origin 200 biprocessed server (Silicon Graphics, Mountain View, CA). The tertiary structure of 11ßHSD2 was manually corrected according to second structure consensus patterns, transferring coordinates for conflict regions from a previously generated 11ßHSD2 model, built from glucose dehydrogenase (Bacillus megaterium) structure. Model for transmembrane (TM) region (residues 180) gave a 3
-helix motif, which was corrected using the N-terminal of sintaxine 1A template. Core and TM region were independently modeled and manually assembled. The obtained structure was used as model for new coordinate assignment iteration until reaching the final structure (residues 1327). The C-terminal region (residues 328405) was not modeled. Mutant design methodology was substitution of Asp223 to Asn223 in the proposed model. Docking assays were performed with only core models, using natural cofactor NAD+ and substrate cortisol. We defined a 5 Å radius around the substrates atoms placed at binding sites vicinity in which it all were mobiles. Reaction mechanism used was sequentially ordered, first cofactor and then substrate. Energy minimization for docking assays of enzyme cofactor and subsequently enzyme-cofactor-substrate was performed during 15,000 iterations using the conjugate gradient method and consistent valence force field.
| Results |
|---|
|
|
|---|
The patient displayed typical features of AME such as low renin hypertension, polyuria, polydypsia, hypokalemia, failure to thrive, and nephrocalcinosis. PRA and SA were significantly decreased, indicating a suppression of the renin-aldosterone axis. Evidence of an impaired metabolism of serum cortisol to cortisone was given by the diminished amounts of serum cortisone and the increased F/E ratio. Treatment of the patient with spironolactone (25 mg/d) normalized the blood pressure (110/70 mm Hg) and electrolytes and improved growth progression. The heterozygous mother was normotensive and had no features of AME, except by the slightly high F/E ratio.
Molecular studies
Sequencing analysis of the HSD11B2 gene of the affected child revealed a homozygous single substitution C
T in intron 3 (IVS3 + 14 C
T; nt 7279) and a homozygous nucleotide change in the second codon of exon 4 (GAC
AAC, nt 7375) as shown in Fig. 1
. This nucleotide change corresponds to a novel mutation that results in an amino acid change from Asp223 to Asn223. Sequencing analysis of the mother showed a heterozygous haplotype for both mutations. To confirm the Asp223Asn mutation, we performed a restriction analysis using the BsmFI enzyme in the patient, mother, and 20 normotensive subjects. As shown in Fig. 2
, the patient had only the noncleaved 1100-bp fragment, indicating that the mutations are in both alleles. The mother showed the original 1100-bp fragment and two other bands of 810 and 290 bp, confirming that she was heterozygous for the exonic mutation. Enzymatic restriction in normotensives showed only the expected 810- and 290-bp fragments, indicating only the normal sequence (only two representative subjects, NT1 and NT2, are shown). We also studied several individuals with essential hypertension, and we were unable to detect the Asp223Asn mutation or the substitution in intron 3 (data not shown).
|
|
To determine the functional consequences of the Asp223Asn mutation, we constructed expression vectors to produce the normal and mutant 11ßHSD2 proteins in suitable mammalian cells. These vectors were used to transfect CHO cells, which do not normally express 11ßHSD2. RT-PCR analysis indicated comparable expression of HSD11B2 RNA from both constructs (data no shown). Western immunoblotting analysis of transfected CHO cells with either wild-type or mutant Asp223Asn vectors showed the presence of a protein approximately of 42 kDa (Fig. 3
), but less immunoreactivity was observed in lysates from cells transfected with mutant construct. By contrast, CHO cells transfected with the pCR3 vector alone and untransfected cells had no immunologically identifiable 11ßHSD2 protein.
|
|
We obtained a structure for wild-type protein (core and TM region, residues 1327), which computes a potential energy of 1642.17 Kcal/mol per Angstrom (Fig. 5
). Model validation was performed using a 3D structures database (Profiles-3D, Insight II; Accelrys Inc.). A high score obtained from the11ßHSD2 folding pattern (residues 1327) validated the proposed tertiary structure. In docking assays, NAD+ binds wild-type protein interacting with the cofactor-binding motif, in this case Gly-Cys-Asp-Ser-Gly-Phe-Gly, having a favorable reaction energy (-212.38 Kcal/mol per Angstrom). Distances for Gly (89)-P1NAD, Gly (93)-C17NAD, and Gly (95)-C18NAD were 4.49, 4.88, and 5.07 Å, respectively. Hydrogen bonding against Asp (91), Ser (92), and Phe (95) of the same motif suggests a right-binding site orientation. Cortisol docking assays converged in wild-type 11ßHSD2, in which cortisol interacts with the key atoms in the catalytic domain, Tyr (232)-H11-F and F-C9-C2-NAD, at distances of 2.49 and 4.57 Å, respectively. We assumed that energy contributions of the N-terminal region (3
-helix motif) and C-terminal region were not relevant for processes occurring at the active site. In the mutant 11ßHSD2, the docking assays for NAD+ and cortisol were not successful in spite of the folding model pattern similar to those found in the wild-type enzyme (root mean square deviation = 0.3506 Å/res) (Fig. 6A
). In the mutant, the binding site volumes for NAD+ and cortisol were smaller than those found in the wild-type enzyme (416 vs. 512 Å3, and 321 vs. 414 Å3, respectively). The calculated surface volume for NAD+ and cortisol were 420 and 297 Å3, respectively. NAD+ would be incapable of entering into the cofactor-binding site of mutant, owing to its size and the compression of the binding pocket. Surface electrostatic potential energy was higher in the mutant than the wild-type enzyme (816 vs. 718 Boltzmann constant of temperature per electrical charge, respectively), which is reflected in the new surface charge distribution (Fig. 6B
).
|
|
| Discussion |
|---|
|
|
|---|
The patients molecular analysis demonstrated the presence of two homozygous substitutions, the C
T change in intron 3 and the GAC
AAC mutation in exon 4 that results in the nonconservative amino acid change Asp223Asn. This is a novel mutation that would be translated to an enzyme with almost no dehydrogenase activity, impairing the conversion of cortisol to cortisone (approximately 6% of normal activity). Apparent kinetics parameters observed in wild-type enzyme were consistent with reports using similar systems (16, 21). Kinetics parameters, Km and Vmax, could not be calculated for the mutant enzyme owing to its very low activity.
Hydropathicity plots and tertiary structure shows that 11ßHSD2 has three successive hydrophobic segments (3
-helix) of approximately 20 amino acids each in the N-terminal region before the cofactor-binding domain. These could function as TM segments anchoring the enzyme to the membrane of the endoplasmic reticulum (13). In the 11ßHSD2 structure modeling, all measured cofactor and substrate distances are at minimal energy position for wild-type docking, which is in agreement with the proposed catalytic mechanism (3, 13), suggesting an adequate spatial disposition for atoms involved in the catalysis of the cortisol hydroxyl group (Fig. 6A
, left). Despite the mutant enzyme folding pattern similarity to the wild-type model, its modeling revealed significant changes, including energetic (changes in the enzymes surface electrostatic potential) and geometric (NAD+ and cortisol-binding site compression) ones. Therefore, we hypothesized that both kinds of changes, energetic and geometric, impair the effective binding of NAD+ and subsequently of cortisol to the catalytic domain, gating the enzyme inactivity.
The C
T change in intron 3 had been previously described by Mune et al. (5). This single base change, 14 bp distal to the exon-intron junction, causes skipping of exon 4, in which the sequence that codifies for the catalytic domain are localized, resulting in an inactive 11ßHSD2 enzyme. This mutation was postulated not to directly alter the splice donor or acceptor sites but induce significant change in the secondary structure of the pre-mRNA that may affect the splicing process (5). Thus, if some mutated RNA is translated, it will carry on the exonic mutation Asp223Asn that, as we just demonstrated, almost abolishes the enzymatic activity. Thus, in this patient both mutations cause the severe AME phenotype observed.
To date, more than 20 mutations have been identified in the 11ß-HSD2 gene, resulting in a functional impairment of the enzyme activity (13, 14, 28, 29, 30, 31, 32, 33, 34, 35). Almost all mutations described have been localized between exons 3 and 5 of the HSD11B2 gene, with the exception of the L114,
6nt mutation (26). In the South American population, only the homozygous Arg213Cys mutation had been reported. This mutation was found in two Chilean siblings with nonconsanguineous heterozygous parents (5, 34). Our patient is the second case described in Chile, but in this patient the possibility of consanguinity cannot be discarded because the patient comes from a small and isolated village in northern Chile and the mother refused to provide any information about the father and other family members.
The impaired activity of 11ßHSD2 was demonstrated in vivo by a high serum F/E ratio, which was several times higher than those found in a normotensive group with normal renin and serum aldosterone (Carvajal, C., L. Mosso, and C. Fardella, unpublished data). Measurement of cortisol and cortisone was performed by RIA and ELISA, respectively, and confirmed by HPLC; the correlation between both measurements was 0.92 and 0.77 (P < 0.05), respectively. There are only a few reports that use the serum cortisol and cortisone values as screening of 11ßHSD2 deficiency (22, 36, 37). In the study by Morineau et al. (22), the serum F/E ratio in affected patients was 35 times higher than that found in normotensive volunteers. This value is similar to those found in our patient, which was close to 5 times higher than that found in a similar control population. There are no reports comparing the classical urinary THF+aTHF/THE ratio with the serum F/E ratio. However, unpublished data (Ferrari, P.) suggest that THF+aTHF/THE ratio is quite similar to serum F/E ratio.
In summary, we report a case of severe hypertension with an AME phenotype that associates a new homozygous Asp223Asn mutation that abolished the enzymatic activity of 11ßHSD2. Likewise, new research tools, such as in silico assays demonstrated that the Asp223Asn mutation produces a disruption in the surface electrostatic potential and geometrical changes that could explain the impairment in cofactor binding in 11ßHSD2 enzyme and then the loss of enzyme activity. The finding of an intronic nucleotide transition in our patient corresponds to the already known change of C
T in intron 3 that cause skipping of exon 4, and impairing of Asp223Asn mutation can be expressed in the mutant enzyme. Both homozygous mutations induce important changes in the catalytic domain and cofactor binding, resulting in an inactive 11ßHSD2 enzyme, which explains the biochemical profile and early onset of hypertension seen in this patient.
| Acknowledgments |
|---|
| Footnotes |
|---|
Abbreviations: 11ßHSD2, 11ß-Hydroxysteroid dehydrogenase type 2; AME, apparent mineralocorticoid excess; CHO-K1, Chinese hamster ovary K1; F/E, cortisol to cortisone; Km, Michaelis-Menten constant; NAD+, nicotinamide adenine dinucleotide (oxidized form); PRA, plasma renin activity; SA, serum aldosterone; THF+aTHF/THE, tetrahydrocortisol plus allotetrahydrocortisol to tetrahydrocortisone; TLC, thin-layer chromatography; TM, transmembrane; Vmax, maximal velocity.
Received December 4, 2002.
Accepted February 16, 2003.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Morineau, V. Sulmont, R. Salomon, B. Fiquet-Kempf, X. Jeunemaitre, J. Nicod, and P. Ferrari Apparent Mineralocorticoid Excess: Report of Six New Cases and Extensive Personal Experience J. Am. Soc. Nephrol., November 1, 2006; 17(11): 3176 - 3184. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. O'Brien, D. Jia, L. I. Plotkin, T. Bellido, C. C. Powers, S. A. Stewart, S. C. Manolagas, and R. S. Weinstein Glucocorticoids Act Directly on Osteoblasts and Osteocytes to Induce Their Apoptosis and Reduce Bone Formation and Strength Endocrinology, April 1, 2004; 145(4): 1835 - 1841. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |