A comparison of the folding of two knotted proteins: YbeA and YibK

The extraordinary topology of proteins belonging to the alpha/beta-knot superfamily of proteins is unexpected, due to the apparent complexities involved in the formation of a deep trefoil knot in a polypeptide backbone. Despite this, an increasing number of knotted structures are being identified; how such proteins fold remains a mystery. Studies on the dimeric protein YibK from Haemophilus influenzae have led to the characterisation of its folding pathway in some detail. To complement research into the folding of YibK, and to address whether folding pathways are conserved for members of the alpha/beta-knot superfamily, the structurally similar knotted protein YbeA from Escherichia coli has been studied. A comprehensive thermodynamic and kinetic analysis of the folding of YbeA is presented here, and compared to that of YibK. Both fold via an intermediate state populated under equilibrium conditions that is monomeric and considerably structured. The unfolding/refolding kinetics of YbeA are simpler than those found for YibK and involve two phases attributed to the formation of a monomeric intermediate state and a dimerisation step. In contrast to YibK, a change in the rate-determining step on the unfolding pathway for YbeA is observed with a changing concentration of urea. Despite this difference, both proteins fold by a mechanism involving at least one sequential monomeric intermediate that has properties similar to that observed during the equilibrium unfolding. The rate of dimerisation observed for YbeA and YibK is very similar, as is the rate constant for formation of the kinetic monomeric intermediate that precedes dimerisation. The findings suggest that relatively slow folding and dimerisation may be common attributes of knotted proteins.     

Probing nature's knots: the folding pathway of a knotted homodimeric protein

The homodimeric protein YibK from Haemophilus influenzae belongs to a recently discovered superfamily of knotted proteins that has brought about a new protein-folding conundrum. Members of the alpha/beta-knot clan form deep trefoil knots in their native backbone structure, a topological feature that is currently unexplained in the protein-folding field. To help solve the puzzle of how a polypeptide chain can efficiently knot itself, the folding kinetics of YibK have been studied extensively and the results are reported here. Folding was monitored using probes for changes in both secondary and tertiary structure, and the monomer-dimer equilibrium was perturbed with a variety of solution conditions to allow characterisation of otherwise inaccessible states. Multiphasic kinetics were observed in the unfolding and refolding reactions of YibK, and under conditions where the dimer is favoured, dissociation and association were rate-limiting, respectively. A folding model consistent with all kinetic data is proposed: YibK appears to fold via two parallel pathways, partitioned by proline isomerisation events, to two distinct monomeric intermediates. These form a common third intermediate that is able to fold to native dimer. Kinetic simulations suggest that all intermediates are on-pathway. These results provide the valuable groundwork required to further understand how Nature codes for knot formation.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

Folding studies on a knotted protein

YibK is a 160 residue homodimeric protein belonging to the SPOUT class of methyltransferases. Proteins in this group all display a unique topological feature; the backbone polypeptide chain folds to form a deep trefoil knot. Such knotted structures were completely unpredicted, it being thought impossible for a protein to fold efficiently in this way. However, they are becoming more common and there are now a growing number of examples in the Protein Data Bank. These intriguing knotted structures represent a new and significant challenge in the field of protein folding. Here, we present an initial characterisation of the folding of YibK, one of the smallest knotted proteins to be identified. This is the first detailed folding study on a knotted protein to be reported. We have established conditions under which the protein can be denatured reversibly in vitro using urea, thereby showing that molecular chaperones are not required for the efficient folding of this protein. A series of equilibrium unfolding experiments were performed over a 400-fold range of protein concentration. Both secondary and tertiary structural probes show a single, protein concentration-dependent unfolding transition, and data are most consistent with a three-state equilibrium denaturation model involving a monomeric intermediate. Thermodynamic parameters obtained from the fit of the data to this model indicate that the intermediate is a stable species with appreciable secondary and tertiary structure; whether the topological knot remains in the intermediate state is still to be shown. Together, these results demonstrate that, despite its complex knotted structure, YibK is able to fold efficiently and behaves remarkably similarly to other dimeric proteins under equilibrium conditions.     

The folding mechanics of a knotted protein

An increasing number of proteins are being discovered with a remarkable and somewhat surprising feature, a knot in their native structures. How the polypeptide chain is able to "knot" itself during the folding process to form these highly intricate protein topologies is not known. Here we perform a computational study on the 160-amino-acid homodimeric protein YibK, which, like other proteins in the SpoU family of MTases, contains a deep trefoil knot in its C-terminal region. In this study, we use a coarse-grained C(alpha)-chain representation and Langevin dynamics to study folding kinetics. We find that specific, attractive nonnative interactions are critical for knot formation. In the absence of these interactions, i.e., in an energetics driven entirely by native interactions, knot formation is exceedingly unlikely. Further, we find, in concert with recent experimental data on YibK, two parallel folding pathways that we attribute to an early and a late formation of the trefoil knot, respectively. For both pathways, knot formation occurs before dimerization. A bioinformatics analysis of the SpoU family of proteins reveals further that the critical nonnative interactions may originate from evolutionary conserved hydrophobic segments around the knotted region.     

Crystal structure of class I acetohydroxy acid isomeroreductase from Pseudomonas aeruginosa

Acetohydroxy acid isomeroreductase (AHIR) is a key enzyme in the biosynthesis of branched-chain amino acids. We have determined the first crystal structure of a class I AHIR from Pseudomonas aeruginosa at 2.0 A resolution. Its dodecameric architecture of 23 point group symmetry is assembled of six dimeric units and dimerization is essential for the formation of the active site. The dimeric unit of P.aeruginosa AHIR partially superimposes with a three-domain monomer of spinach AHIR, a class II enzyme. This demonstrates that the so-called plant-specific insert in the middle of spinach AHIR is structurally and functionally equivalent to the C-terminal alpha-helical domain of P.aeruginosa AHIR, and the C-terminal alpha-helical domain was duplicated during evolution from the shorter, class I AHIRs to the longer, class II AHIRs. The dimeric unit of P.aeruginosa AHIR possesses a deep figure-of-eight knot, essentially identical with that in the spinach AHIR monomer. Thus, our work lowers the likelihood of the previous proposal that "domain duplication followed by exchange of a secondary structure element can be a source of such a knot in the protein structure" being correct.     

The monomeric dUTPase from Epstein-Barr virus mimics trimeric dUTPases

Deoxyuridine 5'-triphosphate pyrophosphatases (dUTPases) are ubiquitous enzymes cleaving dUTP into dUMP and pyrophosphate. They occur as monomeric, dimeric, or trimeric molecules. The trimeric and monomeric enzymes both contain the same five characteristic sequence motifs but in a different order, whereas the dimeric enzymes are not homologous. Monomeric dUTPases only occur in herpesviruses, such as Epstein-Barr virus (EBV). Here, we describe the crystal structures of EBV dUTPase in complex with the product dUMP and a substrate analog alpha,beta-imino-dUTP. The molecule consists of three domains forming one active site that has a structure extremely similar to one of the three active sites of trimeric dUTPases. The three domains functionally correspond to the subunits of the trimeric form. Domains I and II have the dUTPase fold, but they differ considerably in the regions that are not involved in the formation of the unique active site, whereas domain III has only little secondary structure.     

Cooperative mechanism of phosphorylation of the monomeric and dimeric forms of inorganic pyrophosphatase from baker's yeast

A comparative study of phosphorylation of native dimeric and artificial monomeric forms of inorganic pyrophosphatase and its fluoride-stabilized complex with PPi has been carried out. The maximal incorporation of Pi for the dimeric and monomeric proteins is 0.5 and 1 mole per mole of subunit, respectively. The saturation kinetic curves are suggestive of strong positive cooperative interactions. The value of the Hill coefficient (5.5) for the free dimeric enzyme drastically changes upon the active center blockage and/or transition to the monomeric enzyme. Acceleration of dephosphorylation induced by Pi in the presence of Mg2+ is observed only in the case of the dimeric protein. The data obtained indicate that phosphorylation of native dimeric pyrophosphatase occurs according to a "flip-flop" mechanism; the Pi binding in the active center exerts a strong influence on individual steps of the reaction.     

Expression of human topoisomerase I with a partial deletion of the linker region yields monomeric and dimeric enzymes that respond differently to camptothecin

Human topoisomerase I is a 765-residue protein composed of four major domains as follows: the unconserved and highly charged NH(2)-terminal domain, a conserved core domain, the positively charged linker region, and the highly conserved COOH-terminal domain containing the active site tyrosine. Previous studies of the domain structure revealed that near full topoisomerase I activity can be reconstituted in vitro by fragment complementation between recombinant polypeptides approximating the core and COOH-terminal domains. Here we demonstrate that deletion of linker residues Asp(660) to Lys(688) yields an active enzyme (topo70DeltaL) that purifies as both a monomer and a dimer. The dimer is shown to result from domain swapping involving the COOH-terminal and core domains of the two subunits. The monomeric form is insensitive to the anti-tumor agent camptothecin and distributive during in vitro plasmid relaxation assays, whereas the dimeric form is camptothecin-sensitive and processive. However, the addition of camptothecin to enzyme/DNA mixtures causes enhancement of SDS-induced breakage by both monomeric and dimeric forms of the mutant enzyme. The similarity of the dimeric form to the wild type enzyme suggests that some structural feature of the dimer is providing a surrogate linker. Yeast cells expressing topo70DeltaL were found to be insensitive to camptothecin.     

Deep knot structure for construction of active site and cofactor binding site of tRNA modification enzyme

The tRNA(Gm18) methyltransferase (TrmH) catalyzes the 2'-O methylation of guanosine 18 (Gua18) of tRNA. We solved the crystal structure of Thermus thermophilus TrmH complexed with S-adenosyl-L-methionine at 1.85 A resolution. The catalytic domain contains a deep trefoil knot, which mutational analyses revealed to be crucial for the formation of the catalytic site and the cofactor binding pocket. The tRNA dihydrouridine(D)-arm can be docked onto the dimeric TrmH, so that the tRNA D-stem is clamped by the N- and C-terminal helices from one subunit while the Gua18 is modified by the other subunit. Arg41 from the other subunit enters the catalytic site and forms a hydrogen bond with a bound sulfate ion, an RNA main chain phosphate analog, thus activating its nucleophilic state. Based on Gua18 modeling onto the active site, we propose that once Gua18 binds, the phosphate group activates Arg41, which then deprotonates the 2'-OH group for methylation.     

Molecular modeling of the dimeric structure of human lipoprotein lipase and functional studies of the carboxyl-terminal domain.

Lipoprotein lipase (LPL) plays a key role in lipid metabolism. Molecular modeling of dimeric LPL was carried out using insight ii based upon the crystal structures of human, porcine, and horse pancreatic lipase. The dimeric model reveals a saddle-shaped structure and the key heparin-binding residues in the amino-terminal domain located on the top of this saddle. The models of two dimeric conformations - a closed, inactive form and an open, active form - differ with respect to how surface-loop positions affect substrate access to the catalytic site. In the closed form, the surface loop covers the catalytic site, which becomes inaccessible to solvent. Large conformational changes in the open form, especially in the loop and carboxyl-terminal domain, allow substrate access to the active site. To dissect the structure-function relationships of the LPL carboxyl-terminal domain, several residues predicted by the model structure to be essential for the functions of heparin binding and substrate recognition were mutagenized. Arg405 plays an important role in heparin binding in the active dimer. Lys413/Lys414 or Lys414 regulates heparin affinity in both monomeric and dimeric forms. To evaluate the prediction that LPL forms a homodimer in a 'head-to-tail' orientation, two inactive LPL mutants - a catalytic site mutant (S132T) and a substrate-recognition mutant (W390A/W393A/W394A) - were cotransfected into COS7 cells. Lipase activity could be recovered only when heterodimerization occurred in a head-to-tail orientation. After cotransfection, 50% of the wild-type lipase activity was recovered, indicating that lipase activity is determined by the interaction between the catalytic site on one subunit and the substrate-recognition site on the other.     

Crystal structure of the dimeric phosphoenolpyruvate carboxykinase (PEPCK) from Trypanosoma cruzi at 2 A resolution.

ATP-dependent phosphoenolpyruvate carboxykinase (PEPCK) (ATP: oxaloacetate carboxylyase (transphosphorylating), EC 4.1.1.49) is a key enzyme involved in the catabolism of glucose and amino acids in the parasite Trypanosoma cruzi, the causative agent of Chagas' disease. Due to the significant differences in the amino acid sequence and substrate specificity of the human enzyme (PEPCK (GTP-dependent), EC 4.1.1.32), the parasite enzyme has been considered a good target for the development of new anti-chagasic drugs. We have solved the crystal structure of the recombinant PEPCK of T. cruzi up to 2.0 A resolution, characterised the dimeric organisation of the enzyme by solution small angle X-ray scattering (SAXS) and compared the enzyme structure with the known crystal structure of the monomeric PEPCK from Escherichia coli. The dimeric structure possesses 2-fold symmetry, with each monomer sharing a high degree of structural similarity with the monomeric structure of the E. coli PEPCK. Each monomer folds into two complex mixed alpha/beta domains, with the active site located in a deep cleft between the domains. The two active sites in the dimer are far apart from each other, in an arrangement that seems to permit an independent access of the substrates to the two active sites. All residues of the E. coli PEPCK structure that had been found to interact with substrates and metal cofactors have been found conserved and in a substantially equivalent spatial disposition in the T. cruzi PEPCK structure. No substrate or metal ion was present in the crystal structure. A sulphate ion from the crystallisation medium has been found bound to the active site. Solution SAXS data suggest that, in solutions with lower sulphate concentration than that used for the crystallisation experiments, the actual enzyme conformation may be slightly different from its conformation in the crystal structure. This could be due to a conformational transition upon sulphate binding, similar to the ATP-induced transition observed in the E. coli PEPCK, or to crystal packing effects. The present structure of the T. cruzi PEPCK will provide a good basis for the modelling of new anti-chagasic drug leads.     

The domain-swapped dimer of cyanovirin-N is in a metastable folded state: reconciliation of X-ray and NMR structures

The structure of the potent HIV-inactivating protein cyanovirin-N was previously found by NMR to be a monomer in solution and a domain-swapped dimer by X-ray crystallography. Here we demonstrate that, in solution, CV-N can exist both in monomeric and in domain-swapped dimeric form. The dimer is a metastable, kinetically trapped structure at neutral pH and room temperature. Based on orientational NMR constraints, we show that the domain-swapped solution dimer is similar to structures in two different crystal forms, exhibiting solely a small reorientation around the hinge region. Mutation of the single proline residue in the hinge to glycine significantly stabilizes the protein in both its monomeric and dimeric forms. By contrast, mutation of the neighboring serine to proline results in an exclusively dimeric protein, caused by a drastic destabilization of the monomer.     

Dissociation and reconstitution of bovine seminal RNAase: Construction of a hyperactive hybrid dimer

The quaternary structure of bovine seminal ribonuclease, the only dimeric protein in the superfamily of ribonucleases, is maintained both by noncovalent forces and by two intersubunit disulfides. The available monomeric derivatives of the enzyme may not be reassembled into dimers. They are catalytically active, but do not retain certain properties of the dimeric enzyme, such as: (i) the ability to respond cooperatively to increasing substrate concentrations in the rate-limiting reaction step; and (ii) the antitumor and immunosuppressive actions. In this report we describe the preparation of stable monomers of seminal ribonuclease which can be reassociated into covalent dimers indistinguishable from the native protein. With this procedure a hybrid dimer was constructed, made up of a native subunit associated to a subunit catalytically inactivated by selective alkylation of the active site His-119. This dimer was found to have enzymic properties typical of monomeric ribonucleases, such as a hyperbolic saturation curve in the hydrolytic rate-limiting step of the reaction. However, the hybrid dimer was one order-of-magnitude more active than the dimeric enzyme.     

A molecular dynamics investigation of mono and dimeric states of the outer membrane enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca(2+)) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca(2+)-bound dimer, and of the Ca(2+)-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.     

How the oxygen tolerance of a [NiFe]-hydrogenase depends on quaternary structure

‘Oxygen-tolerant’ [NiFe]-hydrogenases can catalyze H₂ oxidation under aerobic conditions, avoiding oxygenation and destruction of the active site. In one mechanism accounting for this special property, membrane-bound [NiFe]-hydrogenases accommodate a pool of electrons that allows an O₂ molecule attacking the active site to be converted rapidly to harmless water. An important advantage may stem from having a dimeric or higher-order quaternary structure in which the electron-transfer relay chain of one partner is electronically coupled to that in the other. Hydrogenase-1 from E. coli has a dimeric structure in which the distal [4Fe-4S] clusters in each monomer are located approximately 12 Å apart, a distance conducive to fast electron tunneling. Such an arrangement can ensure that electrons from H₂ oxidation released at the active site of one partner are immediately transferred to its counterpart when an O₂ molecule attacks. This paper addresses the role of long-range, inter-domain electron transfer in the mechanism of O₂-tolerance by comparing the properties of monomeric and dimeric forms of Hydrogenase-1. The results reveal a further interesting advantage that quaternary structure affords to proteins.     

Primary structure of cold-adapted alkaline phosphatase from a Vibrio sp. as deduced from the nucleotide gene sequence

Alkaline phosphatases (AP) are widely distributed in nature, and generally have a dimeric structure. However, there are indications that either monomeric or multimeric bacterial forms may exist. This paper describes the gene sequence of a psychrophilic marine Vibrio AP, previously shown to be particularly heat labile. The kinetic properties were also indicative of cold adaptation. The amino acid sequence of the Vibrio G15-21 AP reveals that the residues involved in the catalytic mechanism, including those ligating the metal ions, have precedence in other characterized APs. Compared with Escherichia coli AP, the two zinc binding sites are identical, whereas the metal binding site, normally occupied by magnesium, is not. Asp-153 and Lys-328 of E. coli AP are His-153 and Trp-328 in Vibrio AP. Two additional stretches of amino acids not present in E. coli AP are found inserted close to the active site of the Vibrio AP. The smaller insert could be accommodated within a dimeric structure, assuming a tertiary structure similar to E. coli AP. In contrast the longer insert would most likely protrude into the interface area, thus preventing dimer formation. This is the first primary structure of a putative monomeric AP, with indications as to the basis for a monomeric existence. Proximity of the large insert loop to the active site may indicate a surrogate role for the second monomer, and may also shape the catalytic as well as stability characteristics of this enzyme.     

Formation of Cystine Slipknots in Dimeric Proteins

We consider mechanical stability of dimeric and monomeric proteins with the cystine knot motif. A structure based dynamical model is used to demonstrate that all dimeric and some monomeric proteins of this kind should have considerable resistance to stretching that is significantly larger than that of titin. The mechanisms of the large mechanostability are elucidated. In most cases, it originates from the induced formation of one or two cystine slipknots. Since there are four termini in a dimer, there are several ways of selecting two of them to pull by. We show that in the cystine knot systems, there is strong anisotropy in mechanostability and force patterns related to the selection. We show that the thermodynamic stability of the dimers is enhanced compared to the constituting monomers whereas machanostability is either lower or higher.     

Urea-induced inactivation, dissociation, and unfolding of the allosteric phosphofructokinase from Escherichia coli

The influence of urea on the allosteric phosphofructokinase from Escherichia coli has been studied by measuring the changes in enzymatic activity, protein fluorescence, circular dichroism, and retention in size-exclusion chromatography. Tetrameric, dimeric, and monomeric forms of the protein can be discriminated by their elution from a high-performance liquid chromatography gel filtration column. Three successive steps can be detected during the urea-induced denaturation of phosphofructokinase: (i) the dissociation of the native tetramer into dimers which abolishes the activity; (ii) the dissociation of dimers into monomers which exposes the unique tryptophan, Trp-311, to the aqueous solvent; (iii) the unfolding of the monomers which disrupts most of the secondary structure. This pathway involves the ordered dissociation of the interfaces between subunits and supports a previous hypothesis (Deville-Bonne et al., 1989). Phosphofructokinase can be quantitatively renatured from urea solutions, provided that precautions are taken to avoid the aggregation of one insoluble monomeric state. The renaturation of phosphofructokinase from urea implies three steps: an initial folding reaction within the monomeric state is followed by two successive association steps. The faster association step restores the native fluorescence, and the slower regenerates the active enzyme. The renaturation and denaturation of phosphofructokinase correspond to the complex pathway: tetramer in equilibrium dimer in equilibrium folded monomer in equilibrium unfolded monomer. It is found that the subunit interface which forms the regulatory site is more stable and associates 40 times more rapidly than the subunit interface which forms the active site.     

Reversible sequestration of active site cysteines in a 2Fe-2S-bridged dimer provides a mechanism for glutaredoxin 2 regulation in human mitochondria

Human mitochondrial glutaredoxin 2 (GLRX2), which controls intracellular redox balance and apoptosis, exists in a dynamic equilibrium of enzymatically active monomers and quiescent dimers. Crystal structures of both monomeric and dimeric forms of human GLRX2 reveal a distinct glutathione binding mode and show a 2Fe-2S-bridged dimer. The iron-sulfur cluster is coordinated through the N-terminal active site cysteine, Cys-37, and reduced glutathione. The structures indicate that the enzyme can be inhibited by a high GSH/GSSG ratio either by forming a 2Fe-2S-bridged dimer that locks away the N-terminal active site cysteine or by binding non-covalently and blocking the active site as seen in the monomer. The properties that permit GLRX2, and not other glutaredoxins, to form an iron-sulfur-containing dimer are likely due to the proline-to-serine substitution in the active site motif, allowing the main chain more flexibility in this area and providing polar interaction with the stabilizing glutathione. This appears to be a novel use of an iron-sulfur cluster in which binding of the cluster inactivates the protein by sequestering active site residues and where loss of the cluster through changes in subcellular redox status creates a catalytically active protein. Under oxidizing conditions, the dimers would readily separate into iron-free active monomers, providing a structural explanation for glutaredoxin activation under oxidative stress.