Thermodynamics of information transfer between subunits in oligomeric enzymes and kinetic cooperativity. 2. Thermodynamics of kinetic cooperativity

The principles of structural kinetics allow one to define the thermodynamic conditions that are sufficient to generate a certain type of kinetic behavior. If subunits are loosely coupled, that is if no quaternary constraint exists between them, the kinetic behavior of the polymeric enzyme is qualitatively defined by the behavior of an ideal dimer. The nature and the extent of the kinetic cooperativity are defined by the energy of interaction, delta G rho, between two subunits. This energy of interaction is that of an ideal dimer relative to that of the A2 and B2 states. This thermodynamic formulation of a given type of cooperativity holds whatever the degree of polymerization of the enzyme. Under these conditions of loose coupling between subunits, positive kinetic cooperativity cannot be associated with any sigmoidicity of the rate curve. The range of energy coupling where positive kinetic cooperativity must, of necessity, be observed becomes more and more narrow as the number of subunit interactions is increased. This range, however, is independent of the number of subunits. The same situation is not observed for negative cooperativity which appears to be independent of both the number of subunits and the number of subunit interactions. If the subunits are tightly coupled, that is if quaternary constraints exist between them, three thermodynamic parameters, delta G' rho, delta G lambda, delta G mu, are required to define the nature of kinetic cooperativity. delta G' rho is the free energy of an ideal strained dimer relative to that of strained A2 and B2 states. delta G lambda and delta G mu represent the difference of strain energies between conformations A and B and B and B relative to that existing between conformations A and A. One may determine in the parametric space (delta G' rho, delta G lambda, delta G mu) the boundaries between the sufficient conditions that generate a certain type of cooperativity and the lack of these conditions. The kinetic parameters of the rate equation are not all independent. A number of constraint conditions exist between them which depend upon the subunit design of the polymeric enzyme. The existence of these constraint conditions may be diagnostic of a certain type of subunit interactions.     

Evolution of regulatory enzymes towards functional simplicity

The potential kinetic complexity of polymeric regulatory enzymes does not seem to be often expressed in nature. Most of these enzymes exhibit in fact a rather simple kinetic behaviour. This functional simplicity is probably the consequence of constraints between rate constants or of blocking of some reaction steps. Functional simplicity is believed to have emerged in the course of neo-Darwinian evolution as a consequence of a trend towards an improved functional efficiency. Functional efficiency may be reached, in polymeric regulatory enzymes, when either of the two sets of conditions are met. The first set of conditions implies the occurrence of the unicity of enzyme conformation in any transition state, a loose coupling between subunits and an exact balance of the driving forces exerted by the enzyme in the forward and backward directions of the catalytic step. This situation results in constraints between rate constants which allow degenerescence of the steady state rate equation. The second set of conditions involves again the unicity of enzyme conformation in any of the transition states, associated with a tight coupling of subunits, and a driving force exerted by the enzyme much strongly in the forward than in the backward direction of the catalytic step. These conditions imply blocking of some reaction steps and again degenerescence of the corresponding rate equation. The most frequent types of quaternary structure and subunit interactions, namely loose coupling between subunits, and tight coupling associated with conservation of at least one symmetry axis, have probably emerged as molecular organizations, which precisely allow both functional efficiency and simplicity to occur. Indeed these situations probably represent the term of two different evolutionary trends. Therefore enzymes that have not reached this state usually exhibit more complex kinetic behaviour. Wavy curves, “bumps” and turning points may be considered as manifestations of the ancestral character of an enzyme.     

Subunit coupling and kinetic co-operativity of polymeric enzymes. Amplification, attenuation and inversion effects

The principles of structural kinetics, as applied to dimeric enzymes, allow us to understand how the strength of subunit coupling controls both substrate-binding co-operativity, under equilibrium conditions, and kinetic co-operativity, under steady state conditions. When subunits are loosely coupled, positive substrate-binding co-operativity may result in either an inhibition by excess substrate or a positive kinetic co-operativity. Alternatively, negative substrate-binding co-operativity is of necessity accompanied by negative kinetic co-operativity. Whereas the extent of negative kinetic co-operativity is attenuated with respect to the corresponding substrate-binding co-operativity, the positive kinetic co-operativity is amplified with respect to that of the substrate-binding co-operativity. Strong kinetic co-operativity cannot be generated by a loose coupling of subunits. If subunit is propagated to the other, the dimeric enzyme may display apparently surprising co-operativity effects. If the strain of the active sites generated by subunit coupling is relieved in the non-liganded and fully-liganded states, both substrate-binding co-operativity and kinetic co-operativity cannot be negative. If the strain of the active sites however, is not relieved in these states, negative substrate-binding co-operativity is accompanied by either a positive or a negative co-operativity. The possible occurrence of a reversal of kinetic co-operativity, with respect to substrate-binding co-operativity, is the direct consequence of quaternary constraints in the dimeric enzyme. Moreover, tight coupling between subunits may generate a positive kinetic co-operativity which is not associated with any substrate-binding co-operativity. In other words a dimeric enzyme may well bind the substrate in a non co-operative fashion and display a positive kinetic co-operativity generated by the strain of the active sites.     

Subunit interactions in enzyme transition states--antagonism between substrate binding and reaction rate

The principles of structural kinetics as applied to polymeric enzymes have been reinvestigated in order to take account of the probable existence of subunit interactions in the enzyme transition states. On the basis of simple and plausible postulates, structural rate equations have been derived for dimeric enzymes and compared to substrate binding isotherms. It then becomes possible to understand how subunit interactions affect substrate affinity and enzyme reaction rate. There exists an antagonism between substrate binding to the enzyme and the steady state rate of product appearance. If subunit interactions increase the rate of product appearance, they decrease the fractional saturation of the enzyme by the substrate. Alternatively, if they decrease the reaction velocity they increase the fractional saturation. This seemingly paradoxical effect is the direct consequence of subunit interactions occurring in both the ground and the transition states.     

Sequence determinants of quaternary structure in lumazine synthase

Riboflavin, an essential cofactor for all organisms, is biosynthesized in plants, fungi and microorganisms. The penultimate step in the pathway is catalyzed by the enzyme lumazine synthase. One of the most distinctive characteristics of this enzyme is that it is found in different species in two different quaternary structures, pentameric and icosahedral, built from practically the same structural monomeric unit. In fact, the icosahedral structure is best described as a capsid of twelve pentamers. Despite this noticeable difference, the active sites are virtually identical in all structurally studied members. Furthermore, the main regions involved in the catalysis are located at the interface between adjacent subunits in the pentamer. Thus, the two quaternary forms of the enzyme must meet similar structural requirements to achieve their function, but, at the same time, they should differ in the sequence traits responsible for the different quaternary structures observed. Here, we present a combined analysis that includes sequence-structure and evolutionary studies to find the sequence determinants of the different quaternary assemblies of this enzyme. A data set containing 86 sequences of the lumazine synthase family was recovered by sequence similarity searches. Seven of them had resolved three-dimensional structures. A subsequent phylogenetic reconstruction by maximum parsimony (MP) allowed division of the total set into two clusters in accord with their quaternary structure. The comparison between the patterns of three-dimensional contacts derived from the known three-dimensional structures and variation in sequence conservation revealed a significant shift in structural constraints of certain positions. Also, to explore the changes in functional constraints between the two groups, site-specific evolutionary rate shifts were analyzed. We found that the positions involved in icosahedral contacts suffer a larger increase in constraints than the rest. We found eight sequence sites that would be the most important icosahedral sequence determinants. We discuss our results and compare them with previous work. These findings should contribute to refinement of the current structural data, to the design of assays that explore the role of these positions, to the structural characterization of new sequences, and to initiation of a study of the underlying evolutionary mechanisms.     

Evolution of allosteric models for hemoglobin

We compare various allosteric models that have been proposed to explain cooperative oxygen binding to hemoglobin, including the two-state allosteric model of Monod, Wyman, and Changeux (MWC), the Cooperon model of Brunori, the model of Szabo and Karplus (SK) based on the stereochemical mechanism of Perutz, the generalization of the SK model by Lee and Karplus (SKL), and the Tertiary Two-State (TTS) model of Henry, Bettati, Hofrichter and Eaton. The preponderance of experimental evidence favors the TTS model which postulates an equilibrium between high (r)- and low (t)-affinity tertiary conformations that are present in both the T and R quaternary structures. Cooperative oxygenation in this model arises from the shift of T to R, as in MWC, but with a significant population of both r and t conformations in the liganded T and in the unliganded R quaternary structures. The TTS model may be considered a combination of the SK and SKL models, and these models provide a framework for a structural interpretation of the TTS parameters. The most compelling evidence in favor of the TTS model is the nanosecond - millisecond carbon monoxide (CO) rebinding kinetics in photodissociation experiments on hemoglobin encapsulated in silica gels. The polymeric network of the gel prevents any tertiary or quaternary conformational changes on the sub-second time scale, thereby permitting the subunit conformations prior to CO photodissociation to be determined from their ligand rebinding kinetics. These experiments show that a large fraction of liganded subunits in the T quaternary structure have the same functional conformation as liganded subunits in the R quaternary structure, an experimental finding inconsistent with the MWC, Cooperon, SK, and SKL models, but readily explained by the TTS model as rebinding to r subunits in T. We propose an additional experiment to test another key prediction of the TTS model, namely that a fraction of subunits in the unliganded R quaternary structure has the same functional conformation (t) as unliganded subunits in the T quaternary structure.     

The structure of 4-hydroxybenzoyl-CoA thioesterase from arthrobacter sp. strain SU

The 4-chlorobenzoyl-CoA dehalogenation pathway in certain Arthrobacter and Pseudomonas bacterial species contains three enzymes: a ligase, a dehalogenase, and a thioesterase. Here we describe the high resolution x-ray crystallographic structure of the 4-hydroxybenzoyl-CoA thioesterase from Arthrobacter sp. strain SU. The tetrameric enzyme is a dimer of dimers with each subunit adopting the so-called "hot dog fold" composed of six strands of anti-parallel beta-sheet flanked on one side by a rather long alpha-helix. The dimers come together to form the tetramer with their alpha-helices facing outwards. This quaternary structure is in sharp contrast to that previously observed for the 4-hydroxybenzoyl-CoA thioesterase from Pseudomonas species strain CBS-3, whereby the dimers forming the tetramer pack with their alpha-helices projecting toward the interfacial region. In the Arthrobacter thioesterase, each of the four active sites is formed by three of the subunits of the tetramer. On the basis of both structural and kinetic data, it appears that Glu73 is the active site base in the Arthrobacter thioesterase. Remarkably, this residue is located on the opposite side of the substrate-binding pocket compared with that observed for the Pseudomonas enzyme. Although these two bacterial thioesterases demonstrate equivalent catalytic efficiencies, substrate specificities, and metabolic functions, their quaternary structures, CoA-binding sites, and catalytic platforms are decidedly different.     

The transition model of RTK activation: A quantitative framework for understanding RTK signaling and RTK modulator activity

Here, we discuss the transition model of receptor tyrosine kinase (RTK) activation, which is derived from biophysical investigations of RTK interactions and signaling. The model postulates that (1) RTKs can interact laterally to form dimers even in the absence of ligand, (2) different unliganded RTK dimers have different stabilities, (3) ligand binding stabilizes the RTK dimers, and (4) ligand binding causes structural changes in the RTK dimer. The model is grounded in the principles of physical chemistry and provides a framework to understand RTK activity and to make predictions in quantitative terms. It can guide basic research aimed at uncovering the mechanism of RTK activation and, in the long run, can empower the search for modulators of RTK function.     

Thermodynamics of information transfer between subunits in oligomeric enzymes and kinetic cooperativity. 3. Information transfer between the subunits of chloroplast fructose bisphosphatase

The theory and the methods that have been described in the two preceding papers in this journal have been used to analyze the kinetic properties of chloroplast fructose bisphosphatase. The enzyme is a tetramer made up of apparently identical subunits and displays a sigmoidal kinetics with respect to its substrate, fructose bisphosphate. The free ionic species, magnesium and fructose bisphosphate bind to the enzyme and the chelate fructose-bisphosphate-magnesium does not affect the sigmoidicity of the rate curves. The Hill coefficient with respect to free fructose bisphosphate is equal to 2.3, which is indeed incompatible with the view that the enzyme behaves as a dimer of dimers. This conclusion is confirmed by direct analysis of the rate curve. On the basis of the sum of the residuals, their sum of squares, the standard error of the kinetic parameters of the equation, the kinetic scheme associated with a dimer of dimers may be ruled out. On the basis of the same criteria, the fit of an Adair equation to the rate data cannot be retained as satisfactory. This is a direct proof that neither the Monod nor the Koshland model can correctly fit these kinetic data. In fact the model that fits these data best is a structural kinetic scheme where information transfer occurs between each subunit and its three neighbors ('tetrahedral' mode of information transfer). The fit of these models to a large number of kinetic data allows one to compute the free energy profile during the successive binding processes of the four substrate molecules to the enzyme. Whereas the first two steps are associated with an increase of free energy, all the other subsequent steps are associated with a decrease of free energy.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

A range of catalytic efficiencies with avian retroviral protease subunits genetically linked to form single polypeptide chains.

Molecular modeling based on the crystal structure of the Rous sarcoma virus (RSV) protease dimer has been used to link the two identical subunits of this enzyme into a functional, single polypeptide chain resembling the nonviral aspartic proteases. Six different linkages were selected to test the importance of different interactions between the amino acids at the amino and carboxyl termini of the two subunits. These linkages were introduced into molecular clones of fused protease genes and the linked protease dimers were expressed in Escherichia coli and purified. Catalytically active proteins were obtained from the inclusion body fraction after renaturation. The linked protease dimers exhibited a 10-20-fold range in catalytic efficiencies (Vmax/Km) on peptide substrates. Both flexibility and ionic interactions in the linkage region affect catalytic efficiency. Some of the linked protease dimers were 2-3-fold more active than the nonlinked enzyme purified from bacteria, although substrate specificities were unchanged. Similar relative efficiencies were observed using a polyprotein precursor as substrate. Mutation of one catalytic Asp in the most active linked protease dimer inactivated the enzyme, demonstrating that these proteins function as single polypeptide chains rather than as multimers.     

Application of protein grey incidence degree measure to predict protein quaternary structural types

Many proteins are composed of two or more subunits, each associated with different polypeptide chains. The number and arrangement of subunits forming a protein are referred to as quaternary structure. It has been known for long that the functions of proteins are closely related to their quaternary structure. In this paper the grey incidence degree is introduced that can calculate the numerical relation between various components, expressed the similar or different degree between these components. We have demonstrated that introduction of the grey incidence degree can remarkably enhance the success rates in predicting the protein quaternary structural class. It is anticipated that the grey incidence degree can be also used to predict many other protein attributes, such as subcellular localization, membrane protein type, enzyme functional class, GPCR type, protease type, among many others.     

Multiple phosphate positions in the catalytic site of glycogen phosphorylase: structure of the pyridoxal-5''-pyrophosphate coenzyme-substrate analog.

The three-dimensional structure of an R-state conformer of glycogen phosphorylase containing the coenzyme-substrate analog pyridoxal-5'-diphosphate at the catalytic site (PLPP-GPb) has been refined by X-ray crystallography to a resolution of 2.87 A. The molecule comprises four subunits of phosphorylase related by approximate 222 symmetry. Whereas the quaternary structure of R-state PLPP-GPb is similar to that of phosphorylase crystallized in the presence of ammonium sulfate (Barford, D. & Johnson, L.N., 1989, Nature 340, 609-616), the tertiary structures differ in that the two domains of the PLPP-GPb subunits are rotated apart by 5 degrees relative to the T-state conformation. Global differences among the four subunits suggest that the major domains of the phosphorylase subunit are connected by a flexible hinge. The two different positions observed for the terminal phosphate of the PLPP are interpreted as distinct phosphate subsites that may be occupied at different points along the reaction pathway. The structural basis for the unique ability of R-state dimers to form tetramers results from the orientation of subunits with respect to the dyad axis of the dimer. Residues in opposing dimers are in proper registration to form tetramers only in the R-state.     

Purification and characterization of glutathione S-transferases P, S and N. Isolation from rat liver of Yb1 Yn protein, the existence of which was predicted by subunit hybridization in vitro.

The glutathione S-transferases are dimeric proteins and comprise subunits of Mr 25 500 (Ya), 26 500 (Yn), 27 000 (Yb1 and Yb2) and 28 500 (Yc). Enzymes containing Ya and/or Yc subunits have been isolated as have forms containing binary combinations of Yn, Yb1 and Yb2 subunits. To date only one enzyme, transferase S, has been described that is a YbYn heterodimer [Hayes & Chalmers (1983) Biochem. J. 215, 581-588]; the identity of the Yb monomer found in transferase S has not been reported previously. The identification and isolation of a YnYn dimer (transferase N) from rat testis is now described. This has enabled structural and functional comparisons to be made between Yb1, Yb2 and Yn monomers. Reversible dissociation experiments between the YnYn and Yb1Yb1 homodimers and between the YnYn and Yb2Yb2 homodimers demonstrated that Yn monomers can hybridize with both Yb1 and Yb2 monomers. Reversible dissociation of transferases N and C (Yb1Yb2) showed that both Yb1 and Yb2 monomers can hybridize with Yn monomers under competitive conditions. The hydridization data suggest that transferase S represents the Yb2Yn subunit combination. A knowledge of the elution position from chromatofocusing columns of the Yb1Yn hybrid that was formed in vitro enabled a purification scheme to be devised for an enzyme from rat liver (transferase P) believed to consist of Yb1Yn subunits. A comparison of the chromatographic behaviour of the YnYn, Yb1Yb1 and Yb2Yb2 dimers on chromatofocusing and hydroxyapatite columns with the behaviour of transferases P and S on the same matrices suggests these two enzymes may be identified as the Yb1Yn and Yb2Yn dimers respectively. The catalytic activities and the inhibitory effects of non-substrate ligands on transferases P and S are significantly different and again suggest they comprise Yb1 and Yn subunits and Yb2 and Yn subunits respectively; transferase P exhibits a 6-fold higher specific activity for 1,2-dichloro-4-nitrobenzene than does transferase S, whereas, conversely, transferase S possesses a 9-fold higher specific activity for trans-4-phenylbut-3-en-2-one than does transferase P. The quaternary structure of transferases P and S was verified by using peptide mapping and 'Western blotting' techniques.     

Cell surface-localized alkaline phosphatase of Escherichia coli as visualized by reaction product deposition and ferritin-labeled antibodies.

When cells of a wild-type Eschericia coli O8 strain bearing a complete lipopolysaccharide were incubated for alkaline phosphatase reaction product and examined by electron microscopy, the depostion of lead salts was to be observed primarily within the periplasmic space. A similar treatment of cells derived from this strain, which bears a highly abbreviated lipopolysaccharide, showed a mixed cell surface and periplasmic localization of reaction product, suggesting a surface association of a portion of the enzyme. To further explore this possibility, ferritin-antibody conjugates against the active enzyme and its irreversibly dissociated subunits were prepared and allowed to react with cells of both strains. The results obtained from these experiments revealed the presence of both the active enzyme and inactive subunits of the enzyme at the cell surface of the mutant strain. The evidence obtained offers further proof of the validity of the reaction product deposition technique and indicates that alkaline phosphatase may be associated with some component of the outer membrane in this organism. The observation of enzyme subunits at the cell surface further suggests that an association of these subunits with structural components of the cell envelope may provide a locus at which they may dimerize to form active enzyme.     

Brownian dynamics simulations of glycolytic enzyme subsets with F-actin

Previous Brownian dynamics (BD) simulations identified specific basic residues on fructose-1,6-bisphophate aldolase (aldolase) (I. V. Ouporov et al., Biophysical Journal, 1999, Vol. 76, pp. 17-27) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (I. V. Ouporov et al., Journal of Molecular Recognition, 2001, Vol. 14, pp. 29-41) involved in binding F-actin, and suggested that the quaternary structure of the enzymes may be important. Herein, BD simulations of F-actin binding by enzyme dimers or peptides matching particular sequences of the enzyme and the intact enzyme triose phosphate isomerase (TIM) are compared. BD confirms the experimental observation that TIM has little affinity for F-actin. For aldolase, the critical residues identified by BD are found in surface grooves, formed by subunits A/D and B/C, where they face like residues of the neighboring subunit enhancing their electrostatic potentials. BD simulations between F-actin and aldolase A/D dimers give results similar to the native tetramer. Aldolase A/B dimers form complexes involving residues that are buried in the native structure and are energetically weaker; these results support the importance of quaternary structure for aldolase. GAPDH, however, placed the critical residues on the corners of the tetramer so there is no enhancement of the electrostatic potential between the subunits. Simulations using GAPDH dimers composed of either S/H or G/H subunits show reduced binding energetics compared to the tetramer, but for both dimers, the sets of residues involved in binding are similar to those found for the native tetramer. BD simulations using either aldolase or GAPDH peptides that bind F-actin experimentally show complex formation. The GAPDH peptide bound to the same F-actin domain as did the intact tetramer; however, unlike the tetramer, the aldolase peptide lacked specificity for binding a single F-actin domain.     

Asymmetry of binding and physical assignments of CTP and ATP sites in aspartate transcarbamoylase.

The allosteric effectors of aspartate transcarbamoylase from Escherichia coli, CTP and ATP, associate with both the regulatory and the catalytic moieties of the enzyme. Studies with isolated, active subunits yield one binding site per regulatory dimer and one per catalytic trimer. Investigations of effector association with hybrid enzymes, containing either the three regulatory dimers or the two catalytic trimers in inactivated forms, indicate that the data obtained with isolated subunits can be used to analyze the binding patterns of these ligands to the native hexamer. Thus, the nonlinear Scatchard plots, characteristic of the binding of CTP and ATP to the native enzyme, can be interpreted in terms of three effector molecules associating with the regulatory subunits, and two binding to the catalytic moiety of the enzyme. Results with native protein in the presence of saturating concentrations of active site ligands support these assignments. The differences between the binding isotherms of CTP and ATP to the enzyme are due to their different affinities to the two types of subunits. The apparent half-of-the-site saturation of the regulatory moiety of aspartate transcarbamoylase supports the concept that this protein has a tendency to exist in an asymmetric state.     

Variability in quaternary association of proteins with the same tertiary fold: a case study and rationalization involving legume lectins

Legume lectins constitute a family of proteins in which small alterations arising from sequence variations in essentially the same tertiary structure lead to large changes in quaternary association. All of them are dimers or tetramers made up of dimers. Dimerization involves side-by-side or back-to-back association of the flat six-membered beta-sheets in the protomers. Variations within these modes of dimerization can be satisfactorily described in terms of angles defining the mutual disposition of the two subunits. In all tetrameric lectins, except peanut lectin, oligomerization involves the back-to-back association of side-by-side dimers. An attempt has been made to rationalize the observed modes of oligomerization in terms of hydrophobic surface area buried on association, interaction energy and shape complementarity, by constructing energy minimised models in each of which the subunit of one legume lectin is fitted in the quaternary structure of another. The results indicate that all the three indices favor and, thus, provide a rationale for the observed arrangements. However, the discrimination provided by buried hydrophobic surface area is marginal in a few instances. The same is true, to a lesser extent, about that provided by shape complementarity. The relative values of interaction energy turns out to be a still better discriminator than the other two indices. Variability in the quaternary association of homologous proteins is a widely observed phenomenon and the present study is relevant to the general problem of protein folding.     

The crystal structure of a plant 2C-methyl-D-erythritol 4-phosphate cytidylyltransferase exhibits a distinct quaternary structure compared to bacterial homologues and a possible role in feedback regulation for cytidine monophosphate

The homodimeric 2C-methyl-D-erythritol 4-phosphate cytidylyltransferase contributes to the nonmevalonate pathway of isoprenoid biosynthesis. The crystal structure of the catalytic domain of the recombinant enzyme derived from the plant Arabidopsis thaliana has been solved by molecular replacement and refined to 2.0 A resolution. The structure contains cytidine monophosphate bound in the active site, a ligand that has been acquired from the bacterial expression system, and this observation suggests a mechanism for feedback regulation of enzyme activity. Comparisons with bacterial enzyme structures, in particular the enzyme from Escherichia coli, indicate that whilst individual subunits overlay well, the arrangement of subunits in each functional dimer is different. That distinct quaternary structures are available, in conjunction with the observation that the protein structure contains localized areas of disorder, suggests that conformational flexibility may contribute to the function of this enzyme.