Structural Domain Organization of Gastric H+,K+-ATPase and Its Rearrangement during the Catalytic Cycle*

(Received for publication, November 29, 1995, and in revised form, October 21, 1996)

María Gasset Dagger , José Laynez Dagger , Margarita Menéndez Dagger , Vincent Raussens § and Erik Goormaghtigh §par

From the Dagger  Instituto de Química-Física Rocasolano, CSIC, Serrano 119, 28006 Madrid, Spain and § Chimie-Physique des Macromolécules aux Interfaces, CP206/2, Universite Libre de Bruxelles, B-1050 Bruxelles, Belgium

ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
FOOTNOTES
REFERENCES


ABSTRACT

Differential scanning calorimetry has been used to characterized the thermal denaturation of gastric (H+,K+)-ATPase. The excess heat capacity function of (H+,K+)-ATPase in highly oriented gastric vesicles displays two peaks at 53.9 °C (Tm1) and 61.8 °C (Tm2). Its thermal denaturation is an irreversible process that does not exhibit kinetic control and can be resolved in two independent two-state processes. They can be assigned to two cooperative domains located in the cytoplasmic loops of the alpha -subunit, according to the disappearance of the endothermic signal upon removal of these regions by proteinase K digestion. Analysis of the thermal-induced unfolding of the enzyme trapped in different catalytic cycle intermediates has allowed us to get insight into the E1-E2 conformational change. In the E1 forms both transitions are always observed. As Tm1 is shifted to Tm2 by vanadate and ATP interaction, the unfolding mechanism changes from two independent to two sequential two-state transitions, revealing interdomain interactions. Stabilization of the E2 forms results in the disappearance of the second transition at saturation by K+, Mg2+-ATP, and Mg2+-vanadate as well as in significant changes in Tm2 and Delta H1. The catalytic domain melts following a process in which intermolecular interactions either in the native or in the unfolded state might be involved. Interestingly, the E2-vanadate-K+ form displays intermediate properties between the E1 and E2 conformational families.


INTRODUCTION

The (H+,K+)-ATPase (EC 3.6.1.36) is the electroneutral ion pump responsible for acid secretion in the gastric mucosa. In this process cytosolic H+ is exchanged 1:1 for extracellular K+. From a structural point of view this enzyme is a protomer composed of a catalytic beta -subunit (95 kDa) and a glycosylated beta -subunit (52 kDa). The beta -subunit, a membrane multispanning polypeptide, is responsible for the coupling of ATP hydrolysis with ion transport across the membrane, whereas the beta -subunit, a single membrane spanning polypeptide, is required for the proper assembly and targeting of the entire protein. According to the current models of (H+,K+)-ATPase (1), the beta -subunit has 795 amino acid residues exposed on the cytoplasmic side (i.e. extravesicular side), 172 amino acids embedded, and 66 amino acid residues located in the vesicle interior (intravesicular side). The beta -subunit is mainly intravesicular, with 225 amino acid residues inside the vesicle, 40 amino acid residues on the cytoplasmic side, and a single transmembrane segment of 27 residues.

The gastric (H+,K+)-ATPase is a member of the P-ATPase family and shares many features, including sequence homology and enzymatic mechanism, with other members such as the Ca2+-ATPase (2) and especially the (Na+,K+)-ATPase (3). Two major conformations, called E11 and E2, have been identified. In the E1 conformation the proton binding site is cytosolic, and the phosphorylated enzyme reacts with ADP to form ATP (4, 5). The E1-phosphorylated enzyme converts spontaneously to the E2-phosphorylated form of the enzyme, which is no longer sensitive to ADP but is rapidly dephosphorylated in the presence of K+ (6). While definite conformational changes take place during the catalytic cycle as demonstrated from fluorescence of fluorescein-labeled beta -subunit (7, 8) and from limited trypsin digestion pattern (5), little is known on the nature of these conformational changes. Amide hydrogen/deuterium exchange kinetics for Neurospora plasma membrane H+-ATPase have shown that at least 175 amino acid residues are shielded from the solvent in an E2 conformation (9). Membrane disposition of the H5-H6 hairpin of the homologous (Na+,K+)-ATPase beta -subunit has shown to be ligand-dependent (10).

Differential scanning calorimetry can provide information on the structural organization in cooperative domains and on domain-domain interactions in proteins (11-14). However, as discussed by others (11, 12, 14), the number of cooperative units does not always exhibit a straightforward relation with the number of structural domains identified by x-ray crystallography. While this approach has been widely used for water-soluble proteins, there are still only a few studies on membrane proteins (15). We have used this approach in the present study to investigate the structural organization into cooperative domains of the (H+,K+)-ATPase and its modulation by ligand binding. Proteinase K digestion of highly oriented gastric vesicles was used to locate the melting cooperative domains in the different topological regions of the protein. Analysis of the variations induced by ligand binding on the thermal denaturation of the enzyme has allowed us a tentative assignment of these cooperative domains in terms of their specific binding activities and has revealed the existence of structural rearrangements during the catalytic cycle.


EXPERIMENTAL PROCEDURES

Materials

The materials used throughout the whole experimentation were of the highest purity grade available. Proteinase K, vanadate, ATP, phenylmethanesulfonyl fluoride, nigericin, ouabain, and oligomycin were obtained from Sigma. SDS-PAGE reagents were of electrophoresis grade from Bio-Rad.

Gastric Vesicles Isolation and Purification

Gastric vesicles were isolated from hog gastric fundus by differential centrifugation and discontinuous sucrose density gradient ultracentrifugation as described previously (16, 17). The material collected at the 8-30% sucrose interface is referred to as gastric vesicles in this paper. SDS-PAGE reveals essentially a major band of 95 kDa, corresponding to the beta -subunit, and a smear of 60-90 kDa, due to the beta -subunit. This smear is only resolved as a single band after deglycosylation of the beta -subunit (see below).

ATPase Activity

ATPase activity was measured in 40 mM Hepes-Tris, pH 7.2, containing 2 mM ATP and 2 mM MgCl2 in presence or in absence of 20 mM KCl. Sucrose (8%) was present or absent as indicated to create isotonic or osmotic shock conditions, respectively. Reaction was performed at 37 °C for 15 min and stopped by addition of SDS (final concentration 1.75%). Inorganic phosphate released from ATP was quantified according to Stanton (18), except that coloration was developed with ascorbate.

Proteinase K Digestion

Vesicles containing the (H+,K+)-ATPase (2.5 mg ml-1 protein concentration) in 50 mM Hepes-Tris, pH 7.2, containing 8% sucrose (conservation buffer) were incubated at 37 °C with proteinase K in a protease/protein ratio of 1:4 (w/w). Proteolysis was halted by the addition of phenylmethanesulfonyl fluoride (6 mM final concentration). The digested vesicles were adjusted to 0.5 M NaCl and shaken on ice for a few minutes. Then the NaCl concentration was set at 0.25 M by adding 1 mM Hepes buffer, and the samples were centrifuged at 45,000 × g for 45 min. The resulting pellets were resuspended in the conservation buffer, pelleted again in the same conditions, and stored at -20 °C. Aliquots were removed after 5 min, 30 min, 1, 4, 7, 24, 48, and 72 h from the digestion mixture. Since ATPase digestion reaches a plateau after 1 h, as measured from the lipid/protein ratio in the reisolated vesicles, 1-h incubation periods were routinely used for the experiments described in this paper.

Differential Scanning Calorimetry

Calorimetric measurements were performed in a Microcal MC-2 differential scanning calorimeter (Microcal Inc., Northampton, MA) at a heating rate of 0.5 K min-1, unless otherwise stated and under an extra constant pressure of 2 atm. The standard DynaCp, CpCalc, DA-2, and Microcal Origin softwares were used for data acquisition and analysis. The excess heat capacity functions were obtained after base-line subtraction and correction for the instrument time response. In all cases thermal denaturation was found to be irreversible. Since buffer-buffer base lines had the same shape as rescans, the latter could be used as instrument base lines because the offset in Cp due to aggregation or other causes was unimportant (19). Gastric vesicle samples were prepared at a 0.5-2.0 mg/ml protein concentration in 50 mM Hepes-Tris, pH 7.2, containing 0.25 M sucrose, 1 mM EGTA, and the required ligands. The presence of 0.25 M sucrose in the buffer is essential to maintain the integrity of gastric vesicles. Samples requiring 2 mM MgCl2 or CaCl2 were prepared by centrifugation at 75,000 × g for 15 min after a 4-fold dilution in 50 mM Hepes, pH 7.2, and further resuspension of vesicles in an EGTA-depleted buffer containing 2 mM divalent cation chloride. Variations in the buffer pH with temperature were negligible in the 20-80 °C range (7.20 ± 0.15) and followed by a sharp drop to pH 6.65 at 95 °C. The addition of ionizable ligands at the concentration required in the present experiments did not cause additional modifications of pH. The (H+,K+)-ATPase Mr considered for calculations was 147,000 assuming a 90% of total protein abundance.

Analytical Methods

Purity of gastric vesicle (H+,K+)-ATPase as well as proteinase K digestion products was analyzed by SDS-PAGE on 7.5% acrylamide gels (20). Protein concentration was determined using the BCA kit (Pierce) or a modified Lowry assay (21), using bovine serum albumin as standard in both cases. Lipids were assayed with the MPR2 enzymatic test specific for the choline derivatives (Boehringer Mannheim), assuming that lipids carrying a choline group represent 50% of the total of the lipids (22). The beta -subunit of the (H+,K+)-ATPase was deglycosylated using N-glycosidase F (Boehringer Mannheim) under denaturing conditions. Gastric vesicles (20 ml of 1 mg/ml protein concentration solutions) were boiled for 3 min in SDS 1% (w/v) and beta -mercaptoethanol 1% (v/v). Then, 25 ml of n-octyl-beta -D-glucopyranoside 10% (w/v), 150 ml of phosphate-buffered saline/EDTA buffer (137 mM NaCl, 2.7 mM KCl, 1.47 mM KH2PO4, 8.1 mM Na2HPO4, and 20 mM EDTA), and 30 ml of N-glycosidase F (250 units/ml) were added before incubation at 37 °C during 24 h.


RESULTS

Thermal Denaturation of Gastric Vesicles

Gastric vesicle thermal denaturation has been studied by differential scanning calorimetry. The result of a typical scan of gastric vesicles in 50 mM Hepes-Tris, pH 7.05, containing 0.25 M sucrose and 1 mM EGTA at 0.5 K min-1 scan rate is shown in Fig. 1 (trace b). Such a result is representative of those obtained with at least three different enzyme preparations and varying the protein concentration in the 0.5-2.0 mg ml-1 range. The standard deviation of the mean in five independent experiments amounted to ± 25 kcal mol-1 for the overall Delta H (281 kcal mol-1) and ±0.2 °C in peak positions. The thermogram, in the 15-95 °C temperature range, is characterized by the presence of two overlapping peaks, at 53.9 Tm1 and 61.8 °C (Tm2), in agreement with previous studies (23). Since about 90% of gastric vesicle total protein is (H+,K+)-ATPase and the thermal-induced transitions occur within the temperature range described for most proteins, including (Ca2+,Mg2+)-ATPase from sarcoplasmic reticulum (24, 25), it seems likely that both peaks correspond to the (H+,K+)-ATPase denaturation. This assignment is further supported by the sensitivity of both peaks to (H+,K+)-ATPase-specific ligands (see below).


Fig. 1. Excess heat capacity function versus temperature of intact (a-c) and proteinase K-digested gastric vesicles after removal of soluble peptides (d). DSC traces were obtained at the following scanning rates: a, 0.33; b and d, 0.5; c, 1 K min-1. Dotted lines under the experimental traces a-c represent the result of the best fit of deconvolution analysis to two independent two-state transitions.
[View Larger Version of this Image (19K GIF file)]


Thermal denaturation of (H+,K+)-ATPase is an irreversible process under all conditions employed, as indicated by the fact that no endotherm is seen on rescanning the samples under any of the experimental conditions employed. Partial heating experiments show that both transitions are irreversible and proceed independently (data not shown). The effect of scanning rate on the heat capacity profile and thermodynamic parameters was studied (Fig. 1, traces A-C) in order to inspect a possible kinetic control in the denaturation process. Tm1 values are scanning rate independent, within the experimental uncertainty. Tm2 is initially shifted toward lower values as the heating rate decreases and then becomes constant at the lowest scanning rate tested (0.5 and 0.33 K min-1). This type of dynamic effect is consistent with the existence of a slow relaxation process between native and reversibly unfolded forms and excludes a kinetic control by the irreversible steps. Thus, application of equilibrium thermodynamics to the experimental traces would not lead to significant errors (26). In fact, irreversibility appears to be a general feature of membrane protein thermal denaturation arising from the aggregation of the unfolded state in the membrane plane (15, 27).

(H+,K+)-ATPase thermal-induced denaturation can be resolved into two independent two-state processes (Delta HvH/Delta H = 1) with enthalpy changes of 165 and 116 kcal mol-1 for the first and second transitions, respectively (Fig. 1, Table I). Identical values for the calorimetric and the van't Hoff enthalpies are found.

Table I.

Thermodynamic parameters of gastric (H+,K+)-ATPase thermal denaturation: effect of E1-promoting ligands

Tm is the temperature at which the maximum excess heat capacity occurs; Delta Hi is the calorimetric transition enthalpy; Delta HvHi is the van't Hoff enthalpy calculated as 4RT 2Cpmax/Delta Hi, where R is the gas constant. Protein concentration was 1 mg ml-1, except for Ca2+ experiments where 0.5 mg ml-1 was used. Thermodynamic parameters are given as the mean of at least three independent determinations. Maximum S.D. in the derived parameters amounted to: Tm1 ± 0.1 °C; Delta H 1 ± 10 kcal mol-1; Tm2 ± 0.1 °C; Delta H2 ± 15 kcal mol-1. Uncertainty derived from the deconvolution process, always lower than the S.D. obtained from different independent measurements, was used as a criteria to distinguish between random and sequential models.
Ligand Fitting model Tm1  Delta H1  Delta HvH1 Tm2  Delta H2  Delta HvH2

°C kcal/mol kcal/mol °C kcal/mol kcal/mol
None Random 53.9 165 165 61.8 116 116
van 34 µM Random 54.9 233 150 62.4 130 130
van 140 µM Random 56.6 220 181 62.1 123 123
van 1 mM Sequential 58.2 200 200 62.0 116 116
van 5 mM Sequential 60.4 196 196 62.0 120 120
Ca2+ 2 mM Sequential 55.2 196 197 60.1 120 120
Ca2+ 2 mM + ATP 0.27 mM Sequential 56.3 203 203 61.2 118 118
Ca2+ 2 mM + ATP 1 mM Sequential 57.2 202 202 62.2 108 108
EDTA 1 mM + ATP 1 mM Random 55.3 185 185 62.5 110 110

Origin of Thermal Transitions, ATPase Digestion with Proteinase K

To get insight into the protein regions involved in thermal transitions, (H+,K+)-ATPase was digested with proteinase K to differentiate the contributions of different topological regions (cytoplasmic, membrane-embedded, or vesicle internal parts). Accounting for asymmetric distribution of the alpha - and beta -subunits of (H+,K+)-ATPase across the membrane, treatment with proteases will also allow separation of the alpha -subunit contribution (alpha -subunit is the main cytoplasmic protruding fraction) from the beta -subunit, provided that gastric vesicles remain sealed during the isolation process. Activity measurements were carried out to check the integrity and orientation of the (H+,K+)-ATPase in the isolated vesicles. Enzymatic activities are 30 (expressed as µmol Pi/(mg·h)) in the absence of KCl (basal unspecific activity) and 122 in the presence of KCl after vesicle opening by osmotic shock by dilution in a sucrose-free assay medium. Therefore, the full access of ATP and K+ to both side of the membrane results in a supplementary activity of 92 µmol Pi/(mg·h). In the absence of osmotic shock, the supplementary activity was 6 µmol of Pi/(mg·h), accounting for the open vesicles present in the preparation (6.5%) and increases to 81 in the presence of the H+/K+ exchanger nigericin. These results indicate that at least 88% of the vesicles are sealed and exhibit the ATP binding site on the outside of the vesicle, which agrees with the data reported elsewhere (22, 27). Since the ATP binding site is located on the largest cytoplasmic loop of the alpha -subunit (loop C3 limited by H4 and H5 hydrophobic regions), the above results demonstrate that at least 88% of the isolated gastric vesicles exhibit that region directed outside the vesicle.

Before DSC experiments, gastric vesicles were first treated with proteinase K in iso-osmotic conditions. Separation from digested peptide fragments and protease followed, as described under "Experimental Procedures." According to the (H+,K+)-ATPase current model, a large fraction of the alpha -subunit should be removed by this treatment, whereas most of the beta -subunit should remain intact. SDS-PAGE analysis of the deglycosylated beta -subunit give molecular masses of 31 kDa for the untreated beta -subunit and 30.5 kDa (starting at Gly-31 (28)), without intensity loss, after proteinase K digestion. The beta -subunit band completely disappears when the digestion is performed under osmotic shock conditions (data not shown). Therefore, the protease protection observed in iso-osmotic conditions is due to the intravesicular disposition of the largest portion of the beta -subunit, which agrees with the above described model. Regarding the alpha -subunit, SDS-PAGE analysis reveals that its 95-kDa polypeptide chain practically disappears (more than 95%) upon proteinase K digestion, even in the absence of osmotic shock (data not shown). The largely cytoplasmic localization of this subunit is confirmed because no fragments with a mass greater than 20 kDa were detected.

No thermal transition was observed in the thermograms of proteinase K-digested gastric vesicles (Fig. 1, trace D), even at protein concentrations of 2.0 mg ml-1. This finding excludes an independent contribution of the largest part of the beta -subunit to the heat capacity profile of (H+,K+)-ATPase obtained with intact gastric vesicles. The contribution of the transmembrane stalk to the heat capacity function is also excluded, since the structure of this region is preserved in the protease treatment (28-30). Therefore, the thermal transitions observed in intact gastric vesicles should be assigned to the melting of two cooperative domains located in the cytoplasmic protruding segments of the alpha -subunit. This assignment is also supported by the finding of two different transitions in the sarcoplasmic reticulum (Ca2++Mg2+)-ATPase (24), homologous protein that does not contain beta -subunit.

Conformational Families of (H+,K+)-ATPase Associated to the Catalytic Cycle

As demonstrated for the (Na+,K+)-ATPase and other P-ATPases, the catalytic cycle of the (H+,K+)-ATPase is characterized by the presence of two families of conformations, E1 and E2, both in phosphorylated and dephosphorylated forms, with a conformational transition linked to cation transport. The E1 conformation presents a cytoplasmic-facing H+ site and exhibits a high affinity for ATP. After phosphorylation of the alpha -subunit by ATP, an E2 conformation appears, characterized by its high affinity for K+ ions and the low affinity for ATP and H+ (for a more complete description see Ref. 31). Among the different ligands potentially affecting the (H+,K+)-ATPase conformational equilibrium, K+ is known to promote the E2 conformation. Fluorescence of fluorescein-labeled alpha -subunit (7) and limited trypsin digestion (5) revealed qualitative differences in the tertiary structure of the protein when bound by different ligands. Secondary structure probed by Fourier transform infrared spectroscopy (32) did not show significant modification. However, it must be noted that some variation in the alpha -helical content had been observed before by circular dichroism (23). Since very little is known about the nature of these conformational changes, we have investigated here the effect of the E1-E2 ligand-induced transitions on the structural organization of the enzyme. In addition, the analysis of the effect of ligand binding on (H+,K+)-ATPase thermal denaturation can provide further insight into the functional role and, thus, the location of the unfolding cooperative units.

Thermal denaturation of ligand-bound (H+,K+)-ATPase is irreversible under all conditions tested and exhibits the same scanning rate dependence as the unligated protein. It should be mentioned that instrumental artifacts, resulting in exothermic peaks at high temperatures, may be observed at fast scanning rates in irreversible denaturations (33).

Thermal Denaturation of E1 Enzyme Forms

Vanadate binding drives the enzyme into the E1-van form. DSC scans of (H+,K+)-ATPase in the presence of increasing vanadate concentration are reported in Fig. 2, and the thermodynamic parameters are listed in Table I. Vanadate interaction produces a progressive increase in Tm1 that results in merging of both transitions above 0.14 mM concentration (Fig. 2). Therefore, deconvolution analysis was performed assuming either a random model for two independent elementary transitions (domains independently undergo a transition between the folded and denatured forms) or a sequential model for two sequential two-state transitions (domain stability is dependent on whether other domains are folded or unfolded) (13, 34, 35). At saturating vanadate concentrations, deconvolution assuming the sequential model gives best fit of the experimental data (Table I) revealing interdomain communication. However, when Tm1 and Tm2 are separated and the overlap is minimal (vanadate concentrations below 1 mM) consideration of one non-two-state transition is required for best fit. Taking into account the value of the Kd of vanadate (30 µM) (36), these deviations observed under non-saturating conditions probably arise from the reequilibration of the binding equilibrium coupled to protein denaturation, as the temperature is increased (37). Reequilibration of the binding and the unfolding processes with temperature causes a deviation of the first peak from the ideal two-state transition, and such distortion is also manifested in the apparent thermodynamic parameters derived, under these conditions, for the second peak. The observed effects upon vanadate interaction with the (H+,K+)-ATPase (increment of Tm1 and Delta H1, and invariability of Tm2 and Delta H2) is consistent with a direct interaction of this ligand with the cooperative domain involved in the lower temperature transition. Since vanadate is thought to bind at the phosphate site in the C3 loop of the alpha -subunit, near Asp-385 (23), it is reasonable to assign the first transition to a cooperative domain containing this region.


Fig. 2. Effect of vanadate on the thermal denaturation of (H+,K+)-ATPase. Vanadate total concentrations are as follows: 34 µM (a), 140 µM (b), 1 mM (c), and 5 mM (d). Solid line corresponds to the experimental trace, and the dotted line represents the best fit obtained from deconvolution analysis assuming two sequential two-state transitions for traces c and d.
[View Larger Version of this Image (20K GIF file)]


Binding of ATP in the presence of Ca2+ results in ATP hydrolysis and the appearance of the phosphorylated E1-P form of the enzyme, insensitive to K+ ions (7). This form exists in equilibrium with the E2-P conformation and can react with free ADP to dephosphorylate the initial E1 state. Addition of 2 mM Ca2+ modifies the thermal denaturation pattern of (H+,K+)-ATPase (Fig. 3, Table I). Again, the sequential model fits slightly better the experimental curve. The first transition is upward shifted, and its enthalpy change is increased up to 196 kcal mol-1. The second transition is downward shifted to 60.5 °C, with minor modifications in its enthalpy change. A more accurate analysis of Ca2+ effect was prevented by the appearance of a drop in the Cp function, immediately after the end of the second transition, upon increasing either cation or protein concentrations. With these limitations, the experimental results indicate a direct cation-enzyme interaction that affects both transitions. ATP, in the presence of 2 mM Ca2+, further stabilizes the first transition and counteracts the apparent Ca2+-induced destabilization of the second transition (Fig. 3, Table I). As in the previous cases, the experimental curves are better described in terms of a sequential model rather than a random model. Changes produced by Ca2+ and Ca2+-ATP on the first transition resemble those observed with vanadate and, therefore, could reflect the interaction of both ligands as cosubstrates with the catalytic domain contained in the alpha -subunit C3 loop. On the other hand, the antagonist action of both cosubstrates on the second transition could arise from a decrease in the free Ca2+ concentration due to cation chelation by ATP. A direct effect of ATP on the second transition cannot be discarded either because the main ATP binding is located on the second domain or because of the existence of multiple ATP binding sites (38).


Fig. 3. Excess heat capacity function versus temperature of (H+,K+)-ATPase in the presence of 2 mM CaCl2 (a), 2 mM CaCl2+ 0.27 mM ATP (b), 2 mM CaCl2 + 1 mM ATP (c), and 1 mM EDTA + 1 mM ATP (d). Curves displayed in dotted lines represent the result of the best fit obtained from deconvolution analysis. Experimental traces a-c obey a sequential model, and trace d was best fitted by a random model.
[View Larger Version of this Image (22K GIF file)]


Inclusion of EDTA in the medium permits uncoupling of ATP binding from enzyme phosphorylation and linked events. Addition of EDTA prevents ATP hydrolysis, by cosubstrate chelation, protecting the enzyme from total proteolysis (5). EDTA (1 mM) itself does not modify the thermal denaturation profile of the enzyme (data not shown). Addition of 1 mM ATP in the presence of 1 mM EDTA results in the appearance of an exotherm at low temperatures followed by the double peak endotherm characteristic of (H+,K+)-ATPase thermal unfolding (Fig. 3, Table I). The area of the exotherm, centered near 25 °C, is about -170 kcal mol-1, and its molecular basis is under investigation. ATP binding modifies the first transition in a manner similar to vanadate but to a lesser extent. The second transition is slightly stabilized.

Therefore, E1 forms of (H+,K+)-ATPase are characterized by the organization of the alpha -subunit cytoplasmic loops into at least two cooperative unfolding domains. Interdomain communication is evidenced by the analysis of their calorimetric transitions.

Thermal Denaturation of E2 Enzyme Forms

Binding of K+ drives the enzyme into the E2-K+ conformation. Fig. 4 and Table II show the effects of increasing concentrations of KCl on the heat capacity function of the protein and in the derived thermodynamic parameters. It should be noted that the observed effects are K+-specific since equal concentrations of NaCl do not produce significant modifications (Fig. 4, Table II). The effect of K+ on the first transition is saturable and shifts Tm1 to 57.4 ± 0.2 °C without modifications in the calorimetric enthalpy change, Delta H1 = 164 kcal mol-1. However, its Delta HvH increases, probably due to aggregation of the denatured form of this cooperative domain. The second peak undergoes major changes. Delta H2 continuously decreases with increasing K+ down to a negligible value. Tm2 remains nearly constant. The van't Hoff enthalpy change increases to 170 kcal mol-1, probably due to the aggregation of the denatured form. The effect of 10 mM KCl is reversed upon addition of CaCl2.


Fig. 4. Effect of K+ on the excess heat capacity function versus temperature of (H+,K+)-ATPase. The displayed traces were obtained in the presence of 100 mM NaCl (a), as control for ionic strength, 1 mM KCl (b), 10 mM KCl (c), 100 mM KCl (d), and 10 mM KCl + 2 mM CaCl2 (e). Curves displayed in dotted lines represent the result of the best fit obtained from deconvolution analysis.
[View Larger Version of this Image (21K GIF file)]


Table II.

Thermodynamic parameters of (H+,K+)-ATPase thermal denaturation: effect of E2 promoting ligands

Tm is the temperature at which the maximum excess heat capacity occurs; Delta Hi is the calorimetric transition enthalpy; Delta HvHi is the van't Hoff enthalpy calculated as 4RT 2Cpmax/Delta Hi, where R is the gas constant. Protein concentration used was 0.5 mg ml-1. Thermodynamic parameters are given as the mean of at least three independent determinations. Maximum S.D. in the derived parameters amounted to Tm1 ± 0.1 °C; Delta H1 ± 10 kcal mol-1; Tm2 ± 0.1 °C; Delta H2 ± 15 kcal mol-1. Uncertainty derived from the deconvolution process was always lower than the S.D. obtained from different independent measurements.
Ligand Tm1  Delta H1  Delta HvH1 Tm2  Delta H2  Delta HvH2

°C kcal/mol kcal/mol °C kcal/mol kcal/mol
NaCl 100 mM 54.2 165 160 61.8 115 115
KCl 1 mM 56.0 160 186 62.2 63 169
KCl 10 mM 57.4 162 222 62.2 44 170
KCl 100 mM 57.4 164 224
KCl 10 mM + CaCl2 2 mM 56.3 163 163 61.9 105 105
MgCl2 2 mM 54.0 170 170 58.0 100 100
MgCl2 2 mM + ATP 0.27 mM 55.7 195 195 57.6 60 91
MgCl2 2 mM + ATP 1 mM 56.8 197 197 56.9 40 87
MgCl2 2 mM + vana 1 mM 60.7 176 300
MgCl2 2 mM + van 1 mM + KCl 10 mM 63.0 189 302 56.4 107 107

a  van, vanadate.

Addition of ATP in the presence of Mg2+ blocks the E2-P form of the enzyme. As demonstrated above for Ca2+, Mg2+ alone binds to the enzyme, inducing significant protein conformational changes (5, 39). Fig. 5 shows the effect of 2 mM Mg2+ in the absence and in the presence of ATP. It must be mentioned that, in the presence of Mg2+, the experimental traces obtained at a protein concentration of 1 mg ml-1, in the absence or presence of additional ligands (see below), exhibited a constant abrupt drop in the post-transitional base line that is minimized at lower protein concentration. Therefore, protein concentrations of 0.5 mg ml-1 were employed in these cases. Mg2+ binding results in a notable destabilization of the second transition, whereas the first transition remains unaltered (Fig. 5 and Table II). Addition of Mg2+-ATP progressively increases Tm1 and Delta H1. The second transition is further downward shifted and its enthalpy change decreases without varying its van't Hoff enthalpy change.


Fig. 5. Excess heat capacity function versus temperature of (H+,K+)-ATPase in the presence of 2 mM MgCl2 (a), 2 mM MgCl2 + 0.27 mM ATP (b), 2 mM MgCl2 + 1 mM ATP (c), 2 mM MgCl2 + 1 mM vanadate (d), and 2 mM MgCl2 + 1 mM vanadate + 10 mM KCl (e). A, raw data after base-line subtraction. B, deconvolution analysis of the thermograms. Protein concentration used was 0.5 mg/ml. Curves displayed under the experimental trace in dotted lines represent the result of the best fit obtained from deconvolution analysis.
[View Larger Version of this Image (24K GIF file)]


In the presence of Mg2+ and vanadate, the conformational equilibrium is shifted to the E2-van conformation. The variation with temperature of the excess heat capacity curve is characterized by a single peak centered at 60.7 °C with a total enthalpy change of 176 kcal mol-1 (Fig. 5, Table II). The Delta HvH/Delta H ratio increases to 1.7 due to the high increase of the van't Hoff enthalpy, suggesting either intermolecular interactions between native catalytic domains or their association upon thermal denaturation. Differentiation between both possibilities is prevented by the appearance of an abrupt decay of the base line around 75 °C at high protein concentrations. Tetrameric arrangements of the native enzyme in the presence of Mg2+-vanadate have been reported (40).

Thermal denaturation of the E2-van-K+ form, stabilized by (2 mM) Mg2+ + (1 mM) vanadate + (10 mM) K+, is depicted in Fig. 5 (trace e). In this case the denaturation process is characterized by two peaks with Tm values of 63.0 and 56.4 °C (Table II). For the peak at 63 °C, the Delta HvH/Delta Hcal ratio is 1.9, suggesting again that the cooperative unit of the catalytic domain could be a dimer, under the above conditions. The E2-van-K+ form of the enzyme can be described as an intermediate conformation in which hallmarks of E2 (increased cooperativity of the catalytic domain melting) and of E1 (two calorimetric transitions) are seen.

In conclusion, the analysis of the thermal denaturation of the enzyme in E2 forms reveals major changes in both transitions, whereas ligands maintaining the ATPase in the E1 conformation mainly affect the first transition.


DISCUSSION

The excess heat capacity curves of gastric vesicles reveal the presence of at least two different peaks. The abundance of the (H+,K+)-ATPase in these membranes (more than 90% of protein content) and the sensitivity of these transitions to enzyme-specific ligands indicate that both endotherms arise from denaturation of the (H+,K+)-ATPase molecule. The total value of 281 kcal mol-1 measured for the enthalpy change, i.e. 1.9 cal/g (Mr 147,000), is much lower than the enthalpy change associated with the unfolding of water-soluble proteins (the average value 7.8-7 cal/g has been compiled in 41 from 42). It is even lower than the values of 2.9, 2.7, and 2.4 cal/g reported for the cytochrome oxidase from Paracoccus denitrificans, beef heart, and yeast, respectively (41, 43, 44). In the case of the oxidase, a protein fraction corresponding to the membrane-embedded segments (45% of the amino acid residues) was not denatured (41). The present results suggest that an even larger proportion of the (H+,K+)-ATPase could remain structured below 80 °C. In fact, Fourier transform infrared spectroscopy data obtained in our laboratory (28) indicate that up to 55% of the protein could be located in the membrane. While thermal denaturation of membrane proteins affects primarily the extramembrane regions (15), resulting in a reduced specific overall enthalpy change, the contribution of a partial unfolding of the regions responsible for the observed transitions cannot be discarded.

In the absence of ligands, (H+,K+)-ATPase thermal denaturation reveals the presence of two independent elementary transitions at 53.9 and 61.8 °C with enthalpy changes of 165 and 116 kcal mol-1, respectively. The origin of both transitions in terms of structural domains is difficult to assess a priori considering the complexity of the system that involves two different polypeptide chains and a variety of topological regions. The use of proteases to specifically cleave the protein fraction protruding from the membrane is legitimate only if the initial system is homogeneously oriented with respect to the membrane and if the structure of the so-isolated membrane region is maintained. The fact that rupture of the vesicles by hypotonic conditions reveals about the same ATPase activity as collapsing the electrochemical gradient by nigericin indicates that the ATP binding site of the ATPase is highly oriented toward the outside of the vesicles. Other works on membrane-embedded proteins support the idea that interactions between membrane-embedded helices stabilize their own assembly and that external connecting loops are not essential. Retinal binding and loop connections in bacteriorhodopsin were found to make a small contribution to stability by a variety of techniques (45). Rhodopsin can be proteolytically cleaved into three fragments which remain associated (46). Functional Escherichia coli lactose permease is obtained from genes which yield complementary fragments of the protein (47, 48). Dimerization of glycophorin A only depends on the presence of some specific amino acids (49) in its single alpha -helix transmembrane sequence (50). The fact that no thermal transition is observed in the proteolytically isolated membrane region of the protein is consistent with a high stability of the tertiary and secondary structures embedded in the membrane. This is due in part to the absence of competitor hydrogen bond donor or acceptor groups as already discussed elsewhere (41, 51-55). The total disappearance of calorimetric signal in (H+,K+)-ATPase after proteinase K treatment could then be attributed to the loss of those segments responsible for the excess heat capacity function observed in intact vesicles. Therefore, the two transitions can be assigned to two different cooperative domains located in the cytoplasmic loops of the alpha -subunit.

The experimental data also bring evidence for two major families of conformations based on their thermal denaturation mechanism. Under conditions in which the E1 form of the enzyme is promoted, as in the case of unligated E1, E1-van and E1-P, the thermal denaturation process can be described in terms of two two-state transitions. The first transition is sensitive to vanadate and ATP and is therefore suggested to contain a protein domain that encompasses the ATP and vanadate binding sites on the alpha -subunit (Table I). The strong stabilization of the first transition upon Mg2+-ATP and Mg2+-vanadate interaction could reflect the substrate or analog binding to this domain. The second transition is hardly affected. Yet some degree of interdomain communication can be concluded. As described below, there are also clues that ATP in the presence of Mg2+ could have a direct effect on the second domain.

Addition of those ligands promoting the E2 conformational family such as K+, Mg2+, Mg2+-ATP, Mg2+-vanadate, and Mg2+-vanadate-K+ dramatically modifies the previous picture. Importantly, the absence of effect of Na+ on the thermogram underlines the specificity for K+ and rules out any unspecific ionic strength effect on the protein conformation. Both transitions are affected by K+ binding. The second transition displays a progressive disappearance upon K+ addition until it completely vanishes from the scan at 100 mM KCl. It must be mentioned here that KCl was shown to protect the (H+,K+)-ATPase against trypsin digestion, whereas the enzyme is completely degraded into small fragments in its absence (5). The K+ binding site has been suggested to be on the H6 segment (31), a part of the H5-H6 hairpin, accessible to trypsin cleavage in the absence of ligands but not in the presence of K+ (30). This region can be found either membrane-embedded, when occluding cations are present, or in the aqueous phase in their absence in the (Na+,K+)-ATPase (56). Therefore, it is a possibility, although not demonstrated by the present experiments, that the second cooperative domain contains the H5-H6 hairpin. Whether or not the disappearance of the second calorimetric transition is due to membrane insertion of this domain remains to be demonstrated. Indeed, base-line problems in this region might distort the DSC scan profile. The saturable effect of K+ on the first transition parameters (Table I) might be the result of either a domain-domain interaction (if K+ binds the other domain) or an enzyme conformational change on going from E1 to E2. Alternatively, it could be the consequence of a direct interaction of K+ with this domain. In good agreement with fluorescence emission of the fluorescein-labeled enzyme (7) is the partial reversal of the E2 to E1 conformation observed here upon Ca2+ addition which results in the recovery of the calorimetric signal corresponding to thermal denaturation of the second transition. DSC data thereby confirm the direct competition of Ca2+ for Rb+ at the entrance of the cation occlusion site and its failure to replace Rb+ in inducing a conformation with a 19-kDa peptide protected against trypsin hydrolysis in the (Na+,K+)-ATPase (57). The reversibility of the first transition to the E1 state upon Ca2+ addition is not completely achieved. On the other hand, the thermal signal corresponding to the second transitions is 95% recovered (Delta H2 = 105 kcal mol-1) (Table II).

Induction of the E2 conformation by Mg2+, Mg2+-ATP, and Mg2+-vanadate also results in a destabilization of the second transition observed in free (H+,K+)-ATPase, which is accompanied by a progressive disappearance of the thermal transition and characterized by a reduction in the calorimetric enthalpy (Delta H2<Delta HvH2). Mg2+-ATP and Mg2+-vanadate produce essentially the same effect as K+ on the second domain but also modify the first domain, as expected from the observation that free ATP or vanadate mainly affect the first transition (Table I). The specific effect of ATP on the second transition in the presence of Mg2+ (but not in its absence) raises the possibility that ATP could interact with the second domain in the presence of Mg2+. Alternatively, Mg2+ could interact with bound ATP (58) and induce a protein structure rearrangement. In the E2-van-K+ the first domain is further stabilized probably as a result of the increased affinity of the enzyme for vanadate in the presence of K+ (36).

In conclusion, the study reported herein offers new points of view on the conformational multiplicity of P-ATPases, establishing the differences among forms in terms of the structural organization into cooperative domains of the alpha -subunit.


FOOTNOTES

*   This work was supported by Dirección General de Investigación Científica y Técnica (DGICYT, grant PB93-0114), Spain. The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
   F.R.I.A. fellow (Belgium).
par    Senior Research Associate of the National Fund for Scientific Research (Belgium). To whom the correspondence should be addressed: Tel.: 32-2-6505386; Fax: 32-2-6505113; E-mail: egoor{at}ulb.ac.be.
1    The abbreviations used are: E1, unligated form of the enzyme; E1-van, enzyme form obtained in the presence of vanadate; E1-P, phosphoenzyme form obtained in the presence of Ca2+; E2-K+, enzyme conformer blocked by K+; E2-P, phosphoenzyme generated in the presence of Mg2+; E2-van, conformer stabilized by vanadate in the presence of Mg2+; E2-van-K+, enzyme form obtained in the presence of Mg2+, vanadate and K+; PAGE, polyacrylamide gel electrophoresis; DSC, differential scanning calorimetry.

REFERENCES

  1. Sachs, G., Besancon, M., Shin, J. M., Mercier, F., Munson, K., and Hersey, S. (1992) J. Bioenerg. Biomembr. 24, 301-308 [Medline] [Order article via Infotrieve]
  2. McLennan, D. H., Brandl, C. J., Korczak, B., and Green, N. M. (1985) Nature 316, 696-700 [Medline] [Order article via Infotrieve]
  3. Shull, G. E., Schwartz, A., and Lingrel, J. B. (1985) Nature 316, 691-695 [Medline] [Order article via Infotrieve]
  4. Rabon, E., Sachs, G., Mardh, S., and Wallmark, B. (1982) Biochim. Biophys. Acta 688, 515-524 [Medline] [Order article via Infotrieve]
  5. Helmich-de Jong, M. L., van Emst-de Vries, S. E., and de Pont, J. J. H. H. M. (1987) Biochim. Biophys. Acta 905, 358-370 [Medline] [Order article via Infotrieve]
  6. Wallmark, B., Stewart, H. B., Rabon, E., Saccomani, G., and Sachs, G. (1980) J. Biol. Chem. 255, 5313-5319 [Free Full Text]
  7. Jackson, R. J., Mendlein, J., and Sachs, G. (1983) Biochim. Biophys. Acta 731, 9-15 [Medline] [Order article via Infotrieve]
  8. Faller, L. D., Diaz, R. A., Scheiner-Bobis, G., and Farley, R. A. (1991) Biochemistry 30, 3503-3510 [Medline] [Order article via Infotrieve]
  9. Goormaghtigh, E., Vigneron, L., Scarborough, G. A., and Ruysschaert, J.-M. (1994) J. Biol. Chem. 269, 27409-27413 [Abstract/Free Full Text]
  10. Lutsenko, S., Aderko, R., and Kaplan, J. H. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 7936-7940 [Abstract]
  11. Privalov, P. (1979) Adv. Protein Chem. 33, 167-241 [Medline] [Order article via Infotrieve]
  12. Privalov, P. (1982) Adv. Protein Chem. 35, 1-104 [Medline] [Order article via Infotrieve]
  13. Brandts, J. F., Hu, C. Q., Lin, N. L., and Mas, M. T. (1989) Biochemistry 28, 8588-8596 [Medline] [Order article via Infotrieve]
  14. Ramsay, G., and Freire, E. (1990) Biochemistry 29, 8677-8683 [Medline] [Order article via Infotrieve]
  15. Haltia, T., and Freire, E. (1995) Biochim. Biophys. Acta 1228, 1-27 [Medline] [Order article via Infotrieve]
  16. Soumarmon, A., Abastado, M., Bonfils, S., and Lewin, M. J. M. (1980) J. Biol. Chem. 255, 11682-11687 [Abstract/Free Full Text]
  17. Robert, J. C., Benkouka, F., Bayle, D., Hervatin, F., and Soumarmon, A. (1990) Biochim. Biophys. Acta 1024, 167-172 [Medline] [Order article via Infotrieve]
  18. Stanton, M. G. (1968) Anal. Biochem. 22, 27-34 [Medline] [Order article via Infotrieve]
  19. Conejero-Lara, F., Mateo, P. L., Avilés, F. X., and Sánchez-Ruiz, J. M. (1991) Biochemistry 30, 2067-2072 [Medline] [Order article via Infotrieve]
  20. Laemmli, U. K. (1970) Nature 227, 680-685 [Medline] [Order article via Infotrieve]
  21. Bensadou, A., and Weinstein, D. (1976) Anal. Biochem. 70, 241-250 [Medline] [Order article via Infotrieve]
  22. Swarts, H., Van Uem, T., Hoving, S., Fransen, J., and de Pont, J. (1991) Biochim. Biophys. Acta 1070, 283-292 [Medline] [Order article via Infotrieve]
  23. Sachs, G., Berglindh, T., Rabon, E., Stewart, H. B., Wallmark, B., and Saccomani, G. (1980) Ann. N. Y. Acad. Sci. 341, 312-334 [Abstract]
  24. Lepock, J. R., Rodahl, A. M., Zhang, C., Heynen, M. L., Waters, B., and Cheng, K. H. (1990) Biochemistry 29, 681-689 [Medline] [Order article via Infotrieve]
  25. Merino, J. M., Moller, J. V., and Gutiérrez-Merino, C. (1994) FEBS Lett. 343, 155-159 [CrossRef][Medline] [Order article via Infotrieve]
  26. Freire, E., van Osdol, W., Mayorga, O. L., and Sánchez-Ruiz, J. M. (1990) Annu. Rev. Biophys. Biophys. Chem. 19, 159-188 [CrossRef][Medline] [Order article via Infotrieve]
  27. Hall, K., Perez, G., Anderson, D., Gutierrez, C., Munson, K., Hersey, S., and Sachs, G. (1990) Biochemistry 29, 701-706 [Medline] [Order article via Infotrieve]
  28. 411467}Raussens, V., Ruysschaert, J. M., and Goormaghtigh, E. (1996) J. Biol. Chem. 271, in press {411467}
  29. Besancon, M., Shin, J. M., Mercier, F., Munson, K., Miller, M., Hersey, S., and Sachs, G. (1993) Biochemistry 32, 2345-2355 [Medline] [Order article via Infotrieve]
  30. Shin, J. M., and Sachs, G. (1994) J. Biol. Chem. 269, 8642-8646 [Abstract/Free Full Text]
  31. Rabbon, E. C., and Reuben, M. A. (1990) Annu. Rev. Physiol. 52, 321-344 [CrossRef][Medline] [Order article via Infotrieve]
  32. Mitchell, R. C., Haris, P. I., Fallowfield, C., Keeling, D. J., and Chapman, D. (1988) Biochim. Biophys. Acta 941, 31-38 [Medline] [Order article via Infotrieve]
  33. Burz, D. S., Allewell, N. M., Ghosaini, L., Hu, C. Q., and Sturtevant, J. M. (1990) Biophys. Chem. 37, 31-41 [CrossRef][Medline] [Order article via Infotrieve]
  34. Freire, E., and Biltonen, R. L. (1978) Biopolymers 17, 463-479
  35. Freire, E., and Biltonen, R. L. (1978) Biopolymers 17, 481-496
  36. Faller, L. D., Rabon, E., and Sachs, G. (1983) Biochemistry 22, 4676-4685 [Medline] [Order article via Infotrieve]
  37. Shrake, A., and Ross, P. D. (1990) J. Biol. Chem. 265, 5055-5059 [Abstract/Free Full Text]
  38. Faller, L. D. (1989) Biochemistry 28, 6771-6778 [Medline] [Order article via Infotrieve]
  39. Faller, L. D. (1990) Biochemistry 29, 3174-3186
  40. Herbert, H., Xian, Y., Hacksell, I., and Mardh, S. (1992) FEBS Lett. 299, 159-162 [CrossRef][Medline] [Order article via Infotrieve]
  41. Haltia, T., Semo, N., Arrondo, J. L. R., Goñi, F. M., and Freire, E. (1994) Biochemistry 33, 9731-9740 [Medline] [Order article via Infotrieve]
  42. Murphy, K. P., Bhakuni, V., Xie, D., and Freire, E. (1992) J. Mol. Biol. 227, 293-306 [Medline] [Order article via Infotrieve]
  43. Rigel, C. W., and Freire, E. (1987) Biochemistry 26, 4366-4371 [Medline] [Order article via Infotrieve]
  44. Morin, P., Diggs, D., Montgomery, D., and Freire, E. (1990) Biochemistry 29, 781-788 [Medline] [Order article via Infotrieve]
  45. Kahn, T. W., Sturtevant, J. M., and Engelman, D. M. (1992) Biochemistry 31, 8829-8839 [Medline] [Order article via Infotrieve]
  46. Litman, B. (1979) Photochem. Photobiol. 29, 671-677 [Medline] [Order article via Infotrieve]
  47. Bibi, E., and Kaback, R. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 4325-4329 [Abstract]
  48. Wrubel, W., Stochaj, U., Sonnewald, U., Theres, C., and Ehring, R. (1990) J. Bacteriol. 172, 5374-5381 [Medline] [Order article via Infotrieve]
  49. Lemmon, M. A., Flanagan, J. M., Hunt, J. F., Adair, B. D., Bormann, B.-J., Dempsey, C. E., and Engelman, D. M. (1992) J. Biol. Chem. 267, 7683-7689 [Abstract/Free Full Text]
  50. Challou, N., Goormaghtigh, E., Cabiaux, V., Conrath, K., and Ruysschaert, J. M. (1994) Biochemistry 33, 6902-6910 [Medline] [Order article via Infotrieve]
  51. Engelman, D. M., Steiz, T. A., and Goldman, A. (1986) Annu. Rev. Biophys. Biophys. Chem. 15, 321-353 [CrossRef][Medline] [Order article via Infotrieve]
  52. Zhang, Y.-P., Lewis, R. N. A. H., Hodges, R. S., and McElhaney, R. N. (1992) Biochemistry 31, 11572-11578 [Medline] [Order article via Infotrieve]
  53. Brouillette, C. G., Muccio, D. D., and Finney, T. K. (1987) Biochemistry 26, 7431-7438 [Medline] [Order article via Infotrieve]
  54. Shnyrov, V. L., and Mateo, P. L. (1993) FEBS Lett. 324, 237-240 [CrossRef][Medline] [Order article via Infotrieve]
  55. Oikawa, K., Lieberman, D. M., and Reithmeier, R. A. F. (1985) Biochemistry 24, 2843-2848 [Medline] [Order article via Infotrieve]
  56. Lutsenko, S., Anderko, R., and Kaplan, J. H. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 7936-7940 [Abstract]
  57. Shainskaya, A., and Karlish, S. J. D. (1994) J. Biol. Chem. 269, 10780-10789 [Abstract/Free Full Text]
  58. Grisham, C. M., Gupta, R. K., Barnett, R. E., and Mildvan, A. S. (1974) J. Biol. Chem. 249, 6738-6744 [Abstract/Free Full Text]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.