Modeling oligomers with Cn or Dn symmetry: application to CAPRI target 10

The abundance of oligomeric proteins makes them a frequent target for structure prediction. However, homologous proteins sometimes adopt different oligomerization states, rendering the prediction of structures of whole oligomers beyond the scope of comparative modeling. This obstacle can be overcome by combining comparative modeling of the single subunit of an oligomer with docking techniques, designed for predicting subunit-subunit interfaces. We present here algorithms for predicting the structures of homo-oligomers with C(n) or D(n) (n > 2) symmetry. The prediction procedure includes a symmetry-restricted docking step followed by a C(n) or D(n) oligomer-forming step, in which the dimers from the docking step are assembled to oligomers. The procedure is applied to each of the crystallographically independent subunits in 8 C(n) and 3 D(n) oligomers, producing very accurate predictions. It is further applied to a single monomer of the tick-borne encephalitis virus coat protein E (Target 10 of the CAPRI experiment). The predicted trimer ranked 30, obtained via rigid-body geometric-hydrophobic docking followed by C(n) oligomer formation, is very similar to the experimentally observed trimer formed by domain II of this protein. Furthermore, the predicted trimer formed from the separated domain I is also close to the experimental structure.     

Formation and carbon monoxide‐dependent dissociation of Allochromatium vinosum cytochrome c′ oligomers using domain‐swapped dimers

The number of artificial protein supramolecules has been increasing; however, control of protein oligomer formation remains challenging. Cytochrome c′ from Allochromatium vinosum (AVCP) is a homodimeric protein in its native form, where its protomer exhibits a four‐helix bundle structure containing a covalently bound five‐coordinate heme as a gas binding site. AVCP exhibits a unique reversible dimer–monomer transition according to the absence and presence of CO. Herein, domain‐swapped dimeric AVCP was constructed and utilized to form a tetramer and high‐order oligomers. The X‐ray crystal structure of oxidized tetrameric AVCP consisted of two monomer subunits and one domain‐swapped dimer subunit, which exchanged the region containing helices αA and αB between protomers. The active site structures of the domain‐swapped dimer subunit and monomer subunits in the tetramer were similar to those of the monomer subunits in the native dimer. The subunit–subunit interactions at the interfaces of the domain‐swapped dimer and monomer subunits in the tetramer were also similar to the subunit–subunit interaction in the native dimer. Reduced tetrameric AVCP dissociated to a domain‐swapped dimer and two monomers upon CO binding. Without monomers, the domain‐swapped dimers formed tetramers, hexamers, and higher‐order oligomers in the absence of CO, whereas the oligomers dissociated to domain‐swapped dimers in the presence of CO, demonstrating that the domain‐swapped dimer maintains the CO‐induced subunit dissociation behavior of native ACVP. These results suggest that protein oligomer formation may be controlled by utilizing domain swapping for a dimer–monomer transition protein.     

Quaternary structure of pig liver phosphofructokinase

The quaternary structure of phosphofructokinase from pig liver has been studied by electron microscopy. Particles ranging in size from tetramers to long flexible chains of tetramers were commonly observed. Phosphofructokinase tetramers are square planar and approximately 110 A on a side; individual subunits are roughly spherical, with a mean radius of 28 A. Chains are formed by end-to-end association of tetramers rather than by tetramer stacking. The geometry of association implies that phosphofructokinase tetramers possess D2 symmetry, with distinct isologous bonding domains for dimer, tetramer, and chain formation.     

Reaction pathway for the quaternary structure change in hemoglobin

We perform a computer simulation of the quaternary structure change during the allosteric transition of hemoglobin. The simulation is based on a docking procedure by which αβ dimers of human hemoglobin are associated into tetramers after being rotated in various orientations. The stability of tetramers thus reconstituted is estimated from the values of a simplified energy function describing nonbonded interactions and from the area of the surface buried in dimer–dimer contacts (their interface area), which we take to represent stabilizing interactions and solvent contribution. A systematic analysis of tetramers reconstituted with twofold symmetry reveals that when the dimers have the R tertiary structure, only tetramers having R-like quaternary structures are stable. When the dimers have the T tertiary structure, they may associate into T-like tetramers or a variety of quaternary structures ranging from T to near R, thus tracing a plausible reaction pathway for the allosteric transition. We subject intermediates of this pathway to energy refinement with rigid αβ dimers. The refinement demonstrates that symmetrical structures are more stable than non symmetrical ones. A detailed analysis of dimer–dimer contacts in intermediates shows how close packing is maintained over large interfaces throughout the quaternary structure change, especially in the “switch region” of contact between the C helix of α-chains and the FG corner of β-chains.     

Molecular characterization of peripherin-2 and rom-1 mutants responsible for digenic retinitis pigmentosa

Peripherin-2 and Rom-1 are homologous tetraspanning membrane proteins that assemble into noncovalent tetramers and higher order disulfide-linked oligomers implicated in photoreceptor disc morphogenesis. Individuals who coinherit a L185P peripherin-2 mutation and a null or G113E rom-1 mutation are afflicted with retinitis pigmentosa, whereas individuals who inherit only one defective gene are normal. We examined the expression, subunit assembly, and disulfide-mediated oligomerization of L185P and L185A peripherin-2 and L188P Rom-1 by velocity sedimentation, co-immunoprecipitation, and cross-linking. These mutants formed noncovalent dimers under disulfide-reducing conditions but failed to assemble into core tetramers. Under nonreducing conditions, L185P dimers formed disulfide-linked tetramers but not higher order oligomers. L185P coassembled with wild-type peripherin-2 and Rom-1 to form tetramers and higher order disulfide-linked oligomers characteristic of the wild-type proteins. The G113E Rom-1 mutant expressed 20-fold lower than wild-type Rom-1, indicating that it behaves mechanistically as a null allele. We conclude that Leu(185) of peripherin-2 (Leu(188) of Rom-1) is critical for tetramer but not dimer formation and that the core tetramer has 2-fold symmetry. Peripherin-2-containing tetramers are required for higher order disulfide-linked oligomer formation. The level of these oligomers is critical for stable photoreceptor disc formation and the digenic retinitis pigmentosa disease phenotype.     

Construction of a stable dimer of Bacillus stearothermophilus lactate dehydrogenase.

A molecular graphics analysis of the features which prevent cytosolic malate dehydrogenase dimers from forming tetramers was evaluated by its success in predicting the synthesis of a version of the LDH framework which is a stable dimer. Surface residues responsible for malate dehydrogenases being dimers were revealed by superimposing the structures of two dimers of pig cytosolic malate dehydrogenase on one homologous tetramer of L-lactate dehydrogenase from Bacillus stearothermophilus. Four regions were identified as composing the P-axis dimer-dimer interface. Two regions of the dimer were surface loops that collided when built as a tetramer: a large loop (residues 203-207, KNOBI) and a small loop (residues 264-269, KNOBII), and these were candidates to explain the dimeric character of malate dehydrogenase. The analysis was tested by constructing a synthetic B. stearothermophilus lactate dehydrogenase (KNOBI) containing the large malate dehydrogenase loop (residues 203-207 being AYIKLQAKE, and extra four amino acids). The new construct was thermotolerant (90 degrees C) and enzymically active with kcat and KM (pyruvate) values similar to those of the wild-type enzyme. However, whereas the allosteric activator fructose 1,6-bisphosphate decreased KM 100 times for wild type, it had no influence on KNOBI. The molecular volumes of 1-120 microM concentrations of the construct were measured by time-resolved decay of tryptophan fluorescence anisotropy and by gel filtration. Both methods showed the molecular weight of wild type increased from dimer to tetramer with Kd about 20 microM dimer. KNOBI remained a dimer under these conditions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Pyridoxal phosphate inhibits dynamic subunit interchange among serine hydroxymethyltransferase tetramers

Cytoplasmic serine hydroxymethyltransferase (cSHMT) is a tetrameric, pyridoxal phosphate (PLP)-dependent enzyme that catalyzes the reversible interconversion of serine and tetrahydrofolate to glycine and methylenetetrahydrofolate. The enzyme has four active sites and is best described as a dimer of obligate dimers. Each monomeric subunit within the obligate dimer contributes catalytically important amino acid residues to both active sites. To investigate the interchange of subunits among cSHMT tetramers, a dominant-negative human cSHMT enzyme (DNcSHMT) was engineered by making three amino acid substitutions: K257Q, Y82A, and Y83F. Purified recombinant DNcSHMT protein was catalytically inactive and did not bind 5-formyltetrahydrofolate. Coexpression of the cSHMT and DNcSHMT proteins in bacteria resulted in the formation of heterotetramers with a cSHMT/DNcSHMT subunit ratio of 1. Characterization of the cSHMT/DNcSHMT heterotetramers indicates that DNcSHMT and cSHMT monomers randomly associate to form tetramers and that cSHMT/DNcSHMT obligate dimers are catalytically inactive. Incubation of recombinant cSHMT protein with recombinant DNcSHMT protein did not result in the formation of hetero-oligomers, indicating that cSHMT subunits do not exchange once the tetramer is assembled. However, removal of the active site PLP cofactor does permit exchange of obligate dimers among preformed cSHMT and DNcSHMT tetramers, and the formation of heterotetramers from cSHMT and DNcSHMT homodimers does not affect the activity of the cSHMT homodimers. The results of these studies demonstrate that PLP inhibits dimer exchange among cSHMT tetramers and suggests that cellular PLP concentrations may influence the stability of cSHMT protein in vivo.     

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.     

Crystallization, crystal structure analysis and preliminary molecular model of the bilin binding protein from the insect Pieris brassicae

The bilin binding protein of the butterfly Pieris brassicae has been prepared, crystallized and its crystal structure determined at high resolution using film and FAST area detector intensity data. The crystallographic asymmetric unit contains a tetramer of identical subunits with a molecular weight of about 90,000. The crystal structure was determined by isomorphous replacement. Use was made of the molecular symmetry to improve phases. A molecular interpretation of the electron density distribution and partial tracing of the polypeptide chain was possible without amino acid sequence information, as the fold is very similar to retinol binding protein. It is characterized by a beta-barrel formed by two orthogonal beta-sheets and an alpha-helix. The bilin pigment seems to be bound within the beta-barrel analogously to retinol in retinol binding protein. The tetramer in the crystal has C2 symmetry and is a dimer of dimers of quasi-equivalent subunits.     

Predicting oligomeric assemblies: N-mers a primer

Multi-protein complexes play key roles in many biological processes. However, since the structures of these assemblies are hard to resolve experimentally, the detailed mechanism of how they work cooperatively in the cell has remained elusive. Similarly, recent advances on in silico prediction of protein-protein interactions have so far avoided this difficult problem. In this paper, we present a general algorithm to predict molecular assemblies of homo-oligomers. Given the number of N-mers and the 3D structure of one monomer, the method samples all the possible symmetries that N-mers can be assembled. Based on a scoring function that clusters the low free energy structures at each binding interface, the algorithm predicts the complex structure as well as the symmetry of the protein assembly. The method is quite general and does not involve any free parameters. The algorithm has been implemented as a public server and integrated to the protein-protein complex prediction server ClusPro. Using this application, we validated predictions for trimers, tetramers (discriminating between dimer of dimers and 4-fold symmetry structures), pentamers and hexamers (discriminating between trimer of dimers, dimer of trimers, and 6-fold symmetry structures), for a total of 107 assemblies. For 85% of the multimers, the server predicts the complex structure within an average rms deviation of 2A from the full crystal. For complexes that involve more than one binding interface, the cluster size at each surface provides a strong indication as to which interface forms first. With improving scoring functions and computer power, our multimer docking approach could be used as a framework to address the more general problem of multi-protein assemblies.     

Conformational flexibility of the acetylcholinesterase tetramer suggested by x-ray crystallography.

Acetylcholinesterase, a polymorphic enzyme, appears to form amphiphilic and nonamphiphilic tetramers from a single splice variant; this suggests discrete tetrameric arrangements where the amphipathic carboxyl-terminal sequences can be either buried or exposed. Two distinct, but related crystal structures of the soluble, trypsin-released tetramer of acetylcholinesterase from Electrophorus electricus were solved at 4.5 and 4.2 A resolution by molecular replacement. Resolution at these levels is sufficient to provide substantial information on the relative orientations of the subunits within the tetramer. The two structures, which show canonical homodimers of subunits assembled through four-helix bundles, reveal discrete geometries in the assembly of the dimers to form: (a) a loose, pseudo-square planar tetramer with antiparallel alignment of the two four-helix bundles and a large space in the center where the carboxyl-terminal sequences may be buried or (b) a compact, square nonplanar tetramer that may expose all four sequences on a single side. Comparison of these two structures points to significant conformational flexibility of the tetramer about the four-helix bundle axis and along the dimer-dimer interface. Hence, in solution, several conformational states of a flexible tetrameric arrangement of acetylcholinesterase catalytic subunits may exist to accommodate discrete carboxyl-terminal sequences of variable dimensions and amphipathicity.     

Low-resolution structures of thyroid hormone receptor dimers and tetramers in solution

High-resolution X-ray structures of thyroid hormone (TH) receptor (TR) DNA and ligand binding domains (DBD and LBD) have yielded significant insights into TR action. Nevertheless, the TR DBD and LBD act in concert to mediate TH effects upon gene expression, and TRs form multiple oligomers; however, structures of full-length TRs or DBD-LBD constructs that would clarify these influences are not available. Here, we report low-resolution X-ray structures of the TRbeta DBD-LBD construct in solution which define the shape of dimers and tetramers and likely positions of the DBDs and LBDs. The holo TRbeta DBD-LBD construct forms a homodimer with LBD-DBD pairs in close contact and DBDs protruding from the base in the same direction. The DBDs are connected to the LBDs by crossed extended D domains. The apo hTRbeta DBD-LBD construct forms tetramers that resemble bulged cylinders with pairs of LBD dimers in a head-to-head arrangement with DBD pairs packed tightly against the LBD core. Overall, there are similarities with our previous low-resolution structures of retinoid X receptors, but TRs exhibit two unique features. First, TR DBDs are closely juxtaposed in the dimer and tetramer forms. Second, TR DBDs are closely packed against LBDs in the tetramer, but not the dimer. These findings suggest that TRs may be able to engage in hitherto unknown interdomain interactions and that the D domain must rearrange in different oligomeric forms. Finally, the data corroborate our suggestion that apo TRs form tetramers in solution which dissociate into dimers upon hormone binding.     

The dihedral symmetry of the p53 tetramerization domain mandates a conformational switch upon DNA binding.

The p53 tumor suppressor forms stable tetramers, whose DNA binding activity is allosterically regulated. The tetramerization domain is contained within the C-terminus (residues 323-355) and its three-dimensional structure exhibits dihedral symmetry, such that a p53 tetramer can be considered a dimer of dimers. Under conditions where monomeric p53 fails to bind DNA, we studied the effects of p53 C-terminal mutations on DNA binding. Residues 322-355 were sufficient to drive DNA binding of p53 as a tetramer. Within this region residues predicted by the three-dimensional structure to stabilize tetramerization, such as Arg337 and Phe341, were critical for DNA binding. Furthermore, substitution of Leu344 caused p53 to dissociate into DNA binding-competent dimers, consistent with the location of this residue at the dimer-dimer interface. The p53 DNA site contains two inverted repeats juxtaposed to a second pair of inverted repeats. Thus, the four repeats exhibit cyclic-translation symmetry and cannot be recognized simultaneously by four dihedrally symmetric p53 DNA binding domains. The discrepancy may be resolved by flexible linkers between the p53 DNA binding and tetramerization domains. When these linkers were deleted p53 exhibited novel DNA binding properties consistent with an inability to recognize four contiguous DNA repeats. Allosteric regulation of p53 DNA binding may involve repositioning the DNA binding domains from a dihedrally symmetric state to a DNA-bound asymmetric state.     

Near-planar solution structures of mannose-binding lectin oligomers provide insight on activation of lectin pathway of complement

The complement system is a fundamental component of innate immunity that orchestrates complex immunological and inflammatory processes. Complement comprises over 30 proteins that eliminate invading microorganisms while maintaining host cell integrity. Protein-carbohydrate interactions play critical roles in both the activation and regulation of complement. Mannose-binding lectin (MBL) activates the lectin pathway of complement via the recognition of sugar arrays on pathogenic surfaces. To determine the solution structure of MBL, synchrotron x-ray scattering and analytical ultracentrifugation experiments showed that the carbohydrate-recognition domains in the MBL dimer, trimer, and tetramer are positioned close to each other in near-planar fan-like structures. These data were subjected to constrained modeling fits. A bent structure for the MBL monomer was identified starting from two crystal structures for its carbohydrate-recognition domain and its triple helical region. The MBL monomer structure was used to identify 10-12 near-planar solution structures for each of the MBL dimers, trimers, and tetramers starting from 900 to 6,859 randomized structures for each. These near-planar fan-like solution structures joined at an N-terminal hub clarified how the carbohydrate-recognition domain of MBL binds to pathogenic surfaces. They also provided insight on how MBL presents a structural template for the binding and auto-activation of the MBL-associated serine proteases to initiate the lectin pathway of complement activation.     

Regulatory mutations in Sin recombinase support a structure-based model of the synaptosome

The resolvase Sin regulates DNA strand exchange by assembling an elaborate interwound synaptosome containing catalytic and regulatory Sin tetramers, and an architectural DNA-bending protein. The crystal structure of the regulatory tetramer was recently solved, providing new insights into the structural basis for regulation. Here we describe the selection and characterization of two classes of Sin mutations that, respectively, bypass or disrupt the functions of the regulatory tetramer. Activating mutations, which allow the catalytic tetramer to assemble and function independently at site I (the crossover site), were found at ∼20% of residues in the N-terminal domain. The most strongly activating mutation (Q115R) stabilized a catalytically active synaptic tetramer in vitro . The positions of these mutations suggest that they act by destabilizing the conformation of the ground-state site I-bound dimers, or by stabilizing the altered conformation of the active catalytic tetramer. Mutations that block activation by the regulatory tetramer mapped to just two residues, F52 and R54, supporting a functional role for a previously reported crystallographic dimer–dimer interface. We suggest how F52/R54 contacts between regulatory and catalytic subunits might promote assembly of the active catalytic tetramer within the synaptosome.     

The crystal structure of the catalytic domain of the site-specific recombination enzyme gamma delta resolvase at 2.7 A resolution

The crystal structure of the catalytic domain of the site-specific recombination enzyme gamma delta resolvase has been determined at 2.7 A resolution. Its first 120 amino acids form a central five-stranded, beta-pleated sheet surrounded by five alpha helices. In one of the four dyad-related dimers, the two active site Ser-10 residues are 19 A apart, perhaps close enough to contact and become covalently linked to the DNA at the recombination site. This dimer also forms the only closely packed tetramer found in the crystal. The subunit interface at a second dyad-related dimer is more extensive and more highly conserved among the homologous recombinases; however, its active site Ser-10 residues are more than 30 A apart. Side chains, identified by mutations that eliminate catalysis but not DNA binding, are located on the subunit surface near the active site serine and at the interface between a third dyad-related pair of subunits of the tetramer.     

The synthesis of oligomers of oxetane-based dipeptide isosteres derived from L-rhamnose or D-xylose.

Routes to oligomers (dimers, tetramers, hexamers) of five oxetane-based dipeptide isosteres have been established. Methyl 2,4-anhydro-5-azido-5-deoxy-L-rhamnonate 'monomer' led, by coupling the corresponding carboxylic acid and amine, to a 'dimer'. Reverse-aldol ring-opening occurred on attempted saponification of the dimer, so all further oligomerization was performed using TBDMS C-3 hydroxyl protection. The silyl protected L-rhamnonate monomer led in turn to the dimer (via the monomer acid and amine), the tetramer (via the dimer acid and amine) and finally the hexamer (via the tetramer acid and dimer amine). In each case the acids were obtained through saponification of the respective methyl esters and the amines were obtained by hydrogenation of the azides; coupling was TBTU-mediated. Essentially the same strategy was employed on equivalent D-lyxonate, 6-deoxy-L-altronate, 6-deoxy-D-gulonate and D-fuconate dipeptide isosteres to give the respective dimers, tetramers and hexamers.     

Amphiphilic and Nonamphiphilic Forms of Bovine and Human Dopamine β-Hydroxylase

We show that human and bovine dopamine beta-hydroxylases (DBH) exist under three main molecular forms: a soluble nonamphiphilic form and two amphiphilic forms. Sedimentation in sucrose gradients and electrophoresis under nondenaturing conditions, by comparison with acetylcholinesterase (AChE), suggest that the three forms are tetramers of the DBH catalytic subunit and bind either no detergent, one detergent micelle, or two detergent micelles. By analogy with the Gna4 and Ga4 AChE forms, we propose to call the nonamphiphilic tetramer Dna4 and the amphiphilic tetramers Da4I and Da4II. In addition to the major tetrameric forms, DBH dimers occur as very minor species, both amphiphilic and nonamphiphilic. Reduction under nondenaturing conditions leads to a partial dissociation of tetramers into dimers, retaining their amphiphilic character. This suggests that the hydrophobic domain is not linked to the subunits through disulfide bonds. The two amphiphilic tetramers are insensitive to phosphatidylinositol phospholipase C, but may be converted into soluble DBH by proteolysis in a stepwise manner; Da4II----Da4I----Dna4. Incubation of soluble DBH with various phospholipids did not produce any amphiphilic form. Several bands corresponding to the catalytic subunits of bovine DBH were observed in sodium dodecyl sulfate-polyacrylamide gel electrophoresis, but this multiplicity was not simply correlated with the amphiphilic character of the enzyme. In the case of human DBH, we observed two bands of 78 and 84 kDa. As previously reported by others, the presence of the heavy subunit characterizes the amphiphilic forms of the enzyme. We discuss the nature of the hydrophobic domain, which could be an uncleaved signal peptide, and the organization of the different amphiphilic and nonamphiphilic DBH forms. We present two models in which dimers may possess either one hydrophobic domain or two domains belonging to each subunit; in both cases, a single detergent micelle would be bound per dimer.     

Crystal structure of sorbitol dehydrogenase

Sorbitol dehydrogenase (SDH) is a distant relative to the alcohol dehydrogenases (ADHs) with sequence identities around 20%. SDH is a tetramer with one zinc ion per subunit. We have crystallized rat SDH and determined the structure by molecular replacement using a tetrameric bacterial ADH as search object. The conformation of the bound coenzyme is extended and similar to NADH bound to mammalian ADH but the interactions with the NMN-part have several differences with those of ADH. The active site zinc coordination in SDH is significantly different than in mammalian ADH but similar to the one found in the bacterial tetrameric NADP(H)-dependent ADH of Clostridiim beijerinckii. The substrate cleft is significantly more polar than for mammalian ADH and a number of residues are ideally located to position the sorbitol molecule in the active site. The SDH molecule can be considered to be a dimer of dimers, with subunits A-B and C-D, where the dimer interactions are similar to those in mammalian ADH. The tetramers are composed of two of these dimers, which interact with their surfaces opposite the active site clefts, which are accessible on the opposite side. In contrast to the dimer interactions, the tetramer-forming interactions are small with only few hydrogen bonds between side-chains.