Protein fluorescence of nicotinamide nucleotide-dependent dehydrogenases

1. The decrease in the protein fluorescence (F) of Neurospora crassa glutamate dehydrogenase is linearly related to the increase in the fraction of the coenzyme sites occupied by NADPH (alpha) at pH6.35. Under these conditions NADPH causes this enzyme to dissociate to monomers. 2. There is a non-linear relationship of F to alpha for NADH binding to give the alcohol dehydrogenase-NADH-isobutyramide complex, the l-glycerol 3-phosphate dehydrogenase-NADH complex and the bovine glutamate dehydrogenase-NADH-glutamate complex. The non-linearity is accurately represented by F=[1-alpha(1-x)](n) where n is the number of NADH-binding sites per protein molecule. 3. The co-operative binding of GTP to bovine glutamate dehydrogenase in the presence of NADH gives a linear relationship between F and alpha. 4. The prediction from the equation F=[1-alpha(1-x)](n) that initial tangents to non-linear protein-fluorescence-quenching curves will intercept the fluorescence when alpha=1 at a value of total ligand concentration less than the sum of the concentration of binding sites in the solution plus the dissociation constant of ligand is quantitatively fulfilled. 5. Non-linear protein-fluorescence titrations may be used to obtain information about the distribution of ligand among the protein molecules in solution.     

DNA bis-intercalation: Theoretical calculation of binding curves

A statistical mechanical calculation of the binding properties of DNA bis-intercalators is presented, based on the sequence-generating function method of Lifson. The effects of binding by intercalation of one or both chromophores of a bifunctional intercalating agent are examined. The secular equation for a general model that includes the effects of neighbor (nearest and non-nearest) exclusion and/or cooperativity in the binding of both singly and doubly intercalated ligands is derived. Numerical results for binding curves are presented for a more restricted model in which each type of bound ligand rigorously excludes its nearest neighbor and the total number of sites covered by a doubly intercalated ligand is variable. At low values of free ligand concentration bis-intercalation dominates the binding process, while at high value of free ligand concentration, intercalation of only one chromophore per ligand becomes significant due to the unavailability of contiguous free sites required for bis-intercalation. Also, depending on the binding parameters, the free energy of the system can be lowered by a loss of doubly intercalated ligands in favor of singly intercalated ones. Corresponding to this transition in binding mode, the average number of sites occupied by a bound ligand decreases from that characteristic of bis-intercalation to that characteristic of mono-intercalation as free ligand concentration increases. An analysis of Scatchard plots describing bis-intercalation is presented.     

Ligand Depletion in vivo Modulates the Dynamic Range and Cooperativity of Signal Transduction

Biological signal transduction commonly involves cooperative interactions in the binding of ligands to their receptors. In many cases, ligand concentrations in vivo are close to the value of the dissociation constant of their receptors, resulting in the phenomenon of ligand depletion. Using examples based on rotational bias of bacterial flagellar motors and calcium binding to mammalian calmodulin, we show that ligand depletion diminishes cooperativity and broadens the dynamic range of sensitivity to the signaling ligand. As a result, the same signal transducer responds to different ranges of signal with various degrees of cooperativity according to its effective cellular concentration. Hence, results from in vitro dose-response analyses cannot be applied directly to understand signaling in vivo. Moreover, the receptor concentration is revealed to be a key element in controlling signal transduction and we propose that its modulation constitutes a new way of controlling sensitivity to signals. In addition, through an analysis of the allosteric enzyme aspartate transcarbamylase, we demonstrate that the classical Hill coefficient is not appropriate for characterizing the change in conformational state upon ligand binding to an oligomeric protein (equivalent to a dose-response curve), because it ignores the cooperativity of the conformational change for the corresponding equivalent monomers, which are generally characterized by a Hill coefficient . Therefore, we propose a new index of cooperativity based on the comparison of the properties of oligomers and their equivalent monomers.     

A model of cooperative rearrangements of biological membranes without changes in the affinity of receptor binding centers

A mathematical model is proposed for cooperative rearrangements induced by specific ligand in certain biomembrane domains. They are considered as the N-valent receptors undergoing rearrangement when n less than or equal to N ligand-binding receptor sites are occupied. The model predicts distinct sigmoidal dependence for change of some structural parameter on ligand concentration when the receptor site-ligand affinity remains constant as binding rises (positive cooperativity is absent).     

Characteristics of a pure endogenous activator of the gastric H+,K+-ATPase system. Evaluation of the role as a possible intracellular regulator

An endogenous activator capable of stimulating the gastric H+,K+-ATPase activity has been purified to homogeneity from dog and pig gastric cells and found to be a dimer of two identical 40-kDa subunits in the active state. Identical nature of the activator monomers was revealed by the detection of lysine as the sole N-terminal amino acid. The activator from one species can stimulate the H+,K+-ATPase from another species and vice versa. Such cross-activation is consistent with the striking similarities in the amino acid composition between the two species, suggesting considerable homology in the activator molecules from different species. The activator exhibited several unique features during modulation of the H+,K+-ATPase reaction. It appreciably enhances affinity of the H+,K+-ATPase for K+, known to increase turnover of the enzyme. To complement this K+ affinity, the activator also enhances ability of the H+,K+-ATPase to generate more transition state (E*.ATP) complex by increasing the entropy of activation (delta S++) of the system as revealed from an Arrhenius plot of the data on temperature activation. In addition, the activator shows both positive cooperativity and strong inhibition, depending on its concentration. Thus, up to the ratio of the H+,K+-ATPase and activator of about 1:2 (on the protein basis), the activator shows sigmoidal activation (Hill coefficient = 4.5), but beyond such concentration a strong inhibition was observed. Finally, Ca2+ at low (2-4 microM) concentration strongly inhibits the activator-stimulated H+,K+-ATPase. It is proposed that the activator may be acting as a link in the signal transducing cascade system between the intracellular second messenger (Ca2+) and the physiological response (gastric H+ transport).     

Anadara ovalis hemoglobins: distinct dissociation and ligand binding characteristics

The red cells of the arcid clam Anadara ovalis contain two electrophoretically distinct hemoglobins: Hb Major (Hb Ma) and Hb Minor (Hb Mi). The major component consists of two electrophoretically indistinguishable tetramers each composed of two heterodimers; the minor hemoblogin is a homodimer whose subunits are different from the tetramer. Functionally, Hb Ma has a higher P50, exhibits a concentration dependent oxygen affinity, has significant ligand cooperativity (n = 2.0), lacks a Bohr effect and is unaffected by ATP. HB Mi has a P50 which is lower and independent of hemoglobin concentration, shows appreciable cooperativity (n = 1.4) and exhibits no heterotropic effects. Both Hb Ma and Mi are resistant to dissociation in the presence of 1.0 M NaI, NaCl and guanidine-HCl but dissociate to monomers when converted to the aquamet but not the cyanmet derivative. The dissociation is completely inhibited by mercaptoethanol. The large number of reactive -SH groups (10-13 per tetramer) suggests that the monomerization is mediated by intra-subunit disulfide bridge formation.     

Decrease in stability of human albumin with increase in protein concentration

The stability (reflected in denaturation temperature, Td) of defatted human albumin monomer, monitored by differential scanning calorimetry, decreases with increasing protein concentration. This is shown to be compatible with a simple model in which reversible polymerization of denatured monomer promotes unfolding. This model also predicts an increase in transition cooperativity with decreasing protein concentration whereas experimentally cooperativity decreases because the rate of thermally induced polymerization of unfolded monomer is slow relative to the scan rate of the calorimeter. The denaturation of undefatted human albumin monomer, subsaturated with high affinity endogenous long-chain fatty acid (LCFA), was previously observed by differential scanning calorimetry to be a biphasic process. Td for the first endotherm, associated with the denaturation of LCFA-poor species, decreases with increasing protein concentration similar to that for defatted monomer whereas Td for the second endotherm, associated with denaturation of LCFA-rich species, is independent of concentration. The magnitude of the concentration dependence of Td relates directly to the extent of polymerization of denatured monomer, which decreases with increasing level of bound ligand. The bimodal thermogram observed for undefatted monomer persists upon simultaneous extrapolation of Td values to low concentration and low scan rate thereby demonstrating that this biphasic denaturation arising from ligand redistribution during denaturation is a true thermodynamic phenomenon and not an artifact of specific experimental conditions or the method used to induce denaturation.     

Phosphate ion partially relieves the cooperativity of effector binding in D-3-phosphoglycerate dehydrogenase without altering the cooperativity of inhibition

The binding of L-serine to phosphoglycerate dehydrogenase from E. coli displays elements of both positive and negative cooperativity. In addition, the inhibition of enzymatic activity by L-serine is also cooperative with Hill coefficients greater than 1. However, phosphate buffer significantly reduces the cooperative effects in serine binding without affecting the cooperativity of inhibition of activity. The maximal degree of inhibition and fluorescence quenching in Tris buffer occurs when an average of two serine binding sites out of four are occupied. This value increases to three out of the four sites at maximal levels of inhibition and quenching in phosphate buffer. The increase from two to three sites appears to be due to the ability of phosphate to reduce the site to site cooperative effects and render each ligand binding site less dependent on each other. The correlation between the level of inhibition and the fractional site occupancy indicates that in Tris buffer, one serine is bound to each interface at maximal effect. In the presence of phosphate, the order of binding appears to change so that both sites at one interface fill before the first site at the opposite interface is occupied. In each case, there is a good correlation between serine binding, conformational change at the regulatory site interfaces, and inhibition of enzyme activity. The observation that phosphate does not appear to have a similar effect on the cooperativity of inhibition of enzymatic activity suggests that there are two distinct cooperative pathways at work: one path between the four serine binding sites, and one path between the serine binding sites and the active sites.     

The relationship between effector binding and inhibition of activity in D-3-phosphoglycerate dehydrogenase

The binding of L-serine to phosphoglycerate dehydrogenase from Escherichia coli displays elements of both positive and negative cooperativity. At pH 7.5, approximately 2 mol of serine are bound per mole of tetrameric enzyme. A substantial degree of positive cooperativity is seen for the binding of the second ligand, but the binding of the third and fourth ligand display substantial negative cooperativity. The data indicate a state of approximately 50% inhibition when only one serine is bound and approximately 80-90% inhibition when two serines are bound. This is consistent with the tethered domain hypothesis that has been presented previously. Comparison of the data derived directly from binding stoichiometry to the binding constants determined from the best fit to the Adair equation, produce a close agreement, and reinforce the general validity of the derived binding constants. The data also support the conclusion that the positive cooperativity between the binding to the first and second site involves binding sites at opposite interfaces over 110 A apart. Thus, an order of binding can be envisioned where the binding of the first ligand initiates a conformational transition that allows the second ligand to bind with much higher affinity at the opposite interface. This is followed by the third ligand, which binds with lesser affinity to one of the two already occupied interfaces, and in so doing, completes a global conformational transition that produces maximum inhibition of activity and an even lower affinity for the fourth ligand, excluding it completely. Thus, maximal inhibition is accomplished with less than maximal occupancy of effector sites through a mechanism that displays strong elements of both positive and negative cooperativity.     

Screening for positive allosteric modulators: assessment of modulator concentration-response curves as a screening paradigm

Targeting allosteric binding sites represents a powerful mechanism for selectively modulating receptor function. The advent of functional assays as the screening method of choice is leading to an increase in the number of allosteric modulators identified. These include positive allosteric modulators that can increase the affinity of the orthosteric agonist and potentiate the evoked response. A common method for screening for positive allosteric modulators is to examine a concentration-response (C/R) curve to the putative modulator in the presence of a single, low concentration of agonist. The study reported here has used data simulations for positive allosteric modulators according to the allosteric ternary complex model to generate modulator C/R curves. The results are then compared to the mechanistic parameters used to simulate the data. It is clear from the simulations that the potency of a positive modulator C/R curve in a screening assay is the product of both its affinity and positive cooperativity. However, it is often difficult to tell which parameter dominates the response; not knowing the actual affinity or cooperativity of a ligand may have consequences for receptor selectivity. Further modeling demonstrates that the use and choice of single agonist concentration, as well as changes in the agonist curve Hill slope, can have significant effects on the modulator C/R curve. Finally, the quantitative relationship between modulator C/R curves and the allosteric ternary complex model is explored. These simulations emphasize the importance of careful interpretation of screening data and of conducting full mechanism of action studies for positive allosteric modulators.     

Cooperative ligand-lattice binding. Approximate Gaussian binding distribution

Nearest-neighbor cooperative binding of a ligand covering n sites and binding with equilibrium constant K and cooperativity factor omega to a large molecule with m binding sites (m much greater than n omega, n/omega) can be approximately described by a Gaussian distribution P(q-qmax), where q is the number of ligands bound and qmax the most probable value of q. The variance of the Gaussian is equal to the derivative dqmax/d ln(L), where L is the free ligand concentration. This variance, sigma 2, is a complicated function of qmax. However, in the limits of very large cooperativity, omega much greater than 1, very large anticooperativity, omega much less than 1, or noncooperativity, omega = 1, simpler expressions for sigma 2 can be given. For qmax = m/(n + 1), where the most probable number of bound ligands equals the number of free binding sites, sigma 2 has a particularly simple form: sigma 2 = 2m omega 1/2/(n + 1)3. The Gaussian and the infinite lattice approximations for the average number of ligands bound are good approximations only if sigma is much smaller than the number of binding sites. The variance may therefore provide an easy check on the validity of the infinite lattice approximation, which is commonly used to analyze experimental binding data.     

Yeast glyceraldehyde-3-phosphate dehydrogenase. Evidence that subunit cooperativity in catalysis can be controlled by the formation of a complex with phosphoglycerate kinase

Sepharose-bound tetrameric, dimeric and monomeric forms of yeast glyceraldehyde-3-phosphate dehydrogenase were prepared, as well as immobilized hybrid species containing (by selective oxidation of an active center cysteine residue with H2O2) one inactivated subunit per tetramer or dimer. The catalytic properties of these enzyme forms were compared in the forward reaction (glyceraldehyde-3-phosphate oxidation) and reverse reaction (1,3-bisphosphoglycerate reductive dephosphorylation) under steady-state conditions. In the reaction of glyceraldehyde-3-phosphate oxidation, immobilized monomeric and tetrameric forms exhibited similar specific activities. The hybrid-modified dimer contributed on half of the total activity of a native dimer. The tetramer containing one modified subunit possessed 75% of the activity of an unmodified tetramer. In the reaction of 1,3-bisphosphoglycerate reductive dephosphorylation, the specific activity of the monomeric enzyme species was nearly twice as high as that of the tetramer, suggesting that only one-half of the active centers of the oligomer were acting simultaneously. Subunit cooperativity in catalysis persisted in an isolated dimeric species. The specific activity of a monomer associated with a peroxide-inactivated monomer in a dimer was equal to that of an isolated monomeric species and twice as high as that of a native immobilized dimer. The specific activity of subunits associated with a peroxide-inactivated subunit in a tetramer did not differ from that of a native immobilized tetramer; this indicates that interdimeric interactions are involved in catalytic subunit cooperativity. A complex was formed between the immobilized glyceraldehyde-3-phosphate dehydrogenase and soluble phosphoglycerate kinase. Three monomers of phosphoglycerate kinase were bound per tetramer of the dehydrogenase and one per dimer. Evidence is presented that if the reductive dephosphorylation of 1,3-bisphosphoglycerate proceeds in the phosphoglycerate kinase - glyceraldehyde-3-phosphate dehydrogenase complex, all active sites of the latter enzyme act independently, i.e. subunit cooperativity is abolished.     

Unfolding pathway of the dimeric and tetrameric forms of phosphofructokinase-2 from Escherichia coli

Escherichia coli phosphofructokinase-2 (Pfk-2) is an oligomeric enzyme characterized by two kinds of interfaces: a monomer-monomer interface, critical for enzymatic activity, and a dimer-dimer interface formed upon tetramerization due to allosteric binding of MgATP. In this work, Pfk-2 was denatured by guanidine hydrochloride (GdnHCl) and the impact of ligand binding on the unfolding pathway of the dimeric and the tertrameric forms of the enzyme was examined. The unligated dimeric form unfolds and dissociates from 0.15 to 0.8 M GdnHCl without the accumulation of native monomers, as indicated by circular dichroism and size exclusion chromatography measurements. However, a monomeric intermediate with an expanded volume and residual secondary structure accumulates above 0.8 M GdnHCl. The dimeric fructose-6-P-enzyme complex shows a shift in the simultaneous dissociation and unfolding process to elevated GdnHCl concentrations (from 0.8 to 1.4 M) together with the expulsion of the ligand detected by intrinsic fluorescence measurements. The unfolding pathway of the tetrameric MgATP-enzyme complex shows the accumulation of a tetrameric intermediate with altered fluorescence properties at about 0.4 M GdnHCl. Above this concentration a sharp transition from tetramers to monomers, without the accumulation of either compact dimers or monomers, was detected by light scattering measurements. Indeed, the most populated species was a partially unfolded monomer about 0.7 M GdnHCl. On the basis of these results, we suggest that the subunit contacts are critical for the maintenance of the overall structure of Pfk-2 and for the binding of ligands, explaining the reported importance of the dimeric state for enzymatic activity.     

Spectroscopic and thermodynamic investigations on the binding of azure B to chondroitin sulfate and the structure of the metachromatic dye complex

The binding of azur B to chondroitin sulfate (CHS) was investigated using absorption spectroscopy. In aqueous solutions it is possible to distinguish three different dye species with absorption bands at 646, 597, and 555 nm. They are assigned to monomers, dimers, and higher aggregates of azure B, which become bound to CHS as the dye concentration (CD) increases. The short-wavelength band (555 nm) causes metachromasia in stained histological materials. When saturation occurs, the metachromatic azure B-CHS complex has a 1:1 composition, i.e., each anionic SO-4 and COO(-)-binding site of CHS binds one dye cation. The composition of the saturated metachromatic complex was determined by spectrophotometric and conductometric titration of CHS with azure B, while the SO-4 and COO- content of CHS was determined by conductometric titration of CHS-acid with NaOH. The binding isotherm of azure B to CHS was determined using gelpermeation chromatography. The isotherm can be described by the model of cooperative binding of ligands to linear biopolymers. We found good agreement between theoretical predictions and experimental findings in the range of 0 less than r less than 0.8 (r = the fraction of occupied binding sites). Using a Schwarz plot, we determined the binding constants of nucleation (Kn = 2.5 X 10(3) M-1) and aggregation (Kq = 1.2 X 10(5) M-1), as well as the cooperativity parameter (q = 50), T = 295 K. With increasing CD, the strong cooperativity of the dye binding favors the formation of metachromatic aggregates rather than monomers and dimers. From the temperature dependence of Kq we evaluated the standard binding enthalpy (delta Hoq = -20.0 kJ mol-1) and entropy (delta Soq = 29.7 JK-1 mol-1) of the cooperative dye binding. The binding was found to be strongly exothermic and accompanied by a thermodynamically favorable entropy increase, this being typical of hydrophobic interactions. Solid azure B-CHS complexes were prepared according to a special dialytic technique and were studied using a microspectrophotometer equipped with a polarizer and an analyzer. The metachromatic 1:1 complex has a broad, intense absorption band whose main peak occurs at 560 nm. This corresponds with the maximum of the metachromatic dye complex in aqueous solution, i.e. 555 nm. The CHS chains of the azure B-CHS complex can be mechanically aligned in a preferred direction (k). We were able to prepare excellently orientated and very fine dye-CHS films which were birefringent and dichroic - the more birefringent, the better the mechanical orientation.(ABSTRACT TRUNCATED AT 400 WORDS)     

The binding capacity is a probability density function.

The binding capacity of a system, or equivalently, the fluctuations of the number of ligands bound around the average value defined by the binding isotherm, can be regarded as a probability density function for the chemical potential of the ligand. The first moment of this density function is the mean ligand activity as defined by Wyman and gives the average free energy (in kT units) of binding per site. The second moment is directly related to the cooperativity of the system. These and higher moments can be obtained from numerical integration of experimental data in a direct way. An analytical expression for the moment generating function shows that the N independent coefficients of the partition function of a system containing N sites are uniquely defined by the first N moments of the binding capacity.     

Mathematical modeling suggests cooperative interactions between a disordered polyvalent ligand and a single receptor site

BACKGROUND: The CDK inhibitor Sic1 must be phosphorylated on at least six sites in order to allow its recognition by the SCF ubiquitin ligase subunit Cdc4. However, because Cdc4 appears to have only a single phospho-epitope binding site, the apparent cooperative dependence on the number of phosphorylation sites in Sic1 cannot be accounted for by traditional thermodynamic models of cooperativity. RESULTS: We develop a general kinetic model, which predicts an unexpected multiplicative increase in affinity as a function of ligand sites. This effect, termed allovalency, derives from a high local concentration of interaction sites moving independently of each other. Modeling of this interaction by a first exit time approach indicates that the probability of ligand rebinding increases exponentially with the number of sites. This type of interaction is relatively immune to loss of any one site and may be easily tuned to any given threshold by adjusting the properties of individual sites. CONCLUSIONS: The allovalency model suggests that a previously undescribed mechanism may underlie certain cooperative interactions. The widespread occurrence of flexible polyvalent ligands in biological systems suggests that this principle may be broadly applicable.     

Unraveling Subunit Cooperativity in Homotetrameric HCN2 Channels

In a multimeric receptor protein, the binding of a ligand can modulate the binding of a succeeding ligand. This phenomenon, called cooperativity, is caused by the interaction of the receptor subunits. By using a complex Markovian model and a set of parameters determined previously, we analyzed how the successive binding of four ligands leads to a complex cooperative interaction of the subunits in homotetrameric HCN2 pacemaker channels. The individual steps in the model were characterized by Gibbs free energies for the equilibria and activation energies, specifying the affinity of the binding sites and the transition rates, respectively. Moreover, cooperative free energies were calculated for each binding step in both the closed and the open channel. We show that the cooperativity sequence positive-negative-positive determined for the binding affinity is generated by the combined effect of very different cooperativity sequences determined for the binding and unbinding rates, which are negative-negative-positive and no-negative-no, respectively. It is concluded that in the ligand-induced activation of HCN2 channels, the sequence of cooperativity based on the binding affinity is caused by two even qualitatively different sequences of cooperativity that are based on the rates of ligand binding and unbinding.     

Ligand-dependent aggregation and cooperativity: a critique

Ligand-dependent site-site (or subunit-subunit) interactions provide the basis for explaining cooperativity in chemical reactions. Even in the simplest possible nonaggregating system, interpretation of the interactions in terms of structural details requires an explicit assumption (or model) for the binding of the ligand to the sites when there are no interactions. This paper develops in detail the processes by which aggregation will yield ligand-dependent cooperativity. Two conceptually distinct free energy differences may contribute to cooperativity in an aggregation reaction. One is the free energy difference in ligand binding between the monomer and the aggregate. The other is derived from ligand-dependent interactions between the sites of the aggregate. In this analysis an explicit distinction is made between the experimentally accessible constants and those derived from assumed models. Experimental measurements of an aggregation cycle in which all of the species in equilibrium are defined do not allow for an evaluation of the energies of interaction without some model (or assumption). In the analysis presented, an explicit assumption is employed relating the constant for binding of the ligand to the isolated monomer and the constant for the binding of the ligand to aggregate under conditions where there are no ligand-dependent interactions.     

Analysis of DNA-RecA protein interactions involving the protein self-association reaction

A method to analyze the DNA binding properties of RecA protein is developed to take into account the protein self-association reaction and is applied to reanalyze the interaction with chemically modified single-stranded DNA (epsilon-DNA). Protein oligomerization is investigated by static light-scattering measurements and analyzed with the accumulated strain model. A coupled equilibrium between DNA binding and self-association of RecA is resolved considering that the complex formed between an m polymer (where m represents the number of units of the polymer) and DNA is identical to the complex resulting from the cooperative binding of m monomers. The cooperativity parameter thus determined is about 10(4), which is more than 100 times higher than the apparent parameter estimated without consideration of the protein oligomerization. This extremely high figure is in good agreement with the formation of large clusters of complex observed by electron microscopy. The apparent DNA binding constant depends upon the ratio of the DNA binding affinity and the self-association constant. For this reason, the variation of the DNA binding constant with the salt concentration is amplified, and the number of ion pairs formed between DNA and RecA obtained from the apparent salt dependence (11 ion pairs/monomer) has been overestimated. Only 2 ion pairs may be formed.     

Immobilized active monomers of D-glyceraldehyde-3-phosphate dehydrogenase from rabbit skeletal muscles and their coenzyme-binding properties

Experimental conditions favouring the dissociation of tetrameric rabbit muscle D-glyceraldehyde-3-phosphate dehydrogenase into active monomers were elaborated. The urea-induced dissociation of the tetramer was shown to be a stepwise process (in 2 M urea only dimers are formed; an increase in urea concentration up to 3 M causes the splitting of the dimers into monomers). The specific activity of immobilized monomers in the glyceraldehyde-3-phosphate oxidation reaction does not differ from that of the parent immobilized tetrameric form. The tetrameric enzyme molecule binds the coenzyme with a negative cooperativity (the first two NAD+ molecules bind with KD below 0.1 microM; for the third and fourth molecules the dissociation constant was determined to be equal to 5.5 +/- 1.5 microM (50 mM medinal buffer, 10 mM sodium phosphate, pH 8.2). The cooperativity of NAD+ binding is preserved in the immobilized preparation of tetrameric dehydrogenase. The immobilized monomers bind NAD+ with KD of 1.6 +/- 1.0 microM. The experimental results are consistent with the hypothesis according to which the association of catalytically active subunits into a tetramer changes their coenzyme-binding properties in such a way that the first two NAD+ molecules bind more firmly to a tetramer than to a monomer, whereas the third and the fourth NAD+ molecules bind less firmly.