Amino acid residues that determine functional specificity of NADP- and NAD-dependent isocitrate and isopropylmalate dehydrogenases

Isocitrate and isopropylmalalte dehydrogenases are homologous enzymes important for the cell metabolism. They oxidize their substrates using NAD or NADP as cofactors. Thus, they have two specificities, towards the substrate and the cofactor, appearing in three combinations. Although many three-dimensional (3D) structures are resolved, identification of amino acids determining these specificities remains a challenge. We present computational identification and analysis of specificity-determining positions (SDPs). Besides many experimentally proven SDPs, we predict new SDPs, for example, four substrate-specific positions (103Leu, 105Thr, 337Ala, and 341Thr in IDH from E. coli) that contact the cofactor and may play a role in the recognition process.     

Positive selection and subfunctionalization of duplicated CCT chaperonin subunits

To reach a functional and energetically stable conformation, many proteins need molecular helpers called chaperonins. Among the group II chaperonins, CCT proteins provide crucial machinery for the stabilization and proper folding of several proteins in the cytosol of eukaryotic cells through interactions that are subunit-specific and geometry-dependent. CCT proteins are made up of eight different subunits, all with similar sequences, positioned in a precise arrangement. Each subunit has been proposed to have a specialized function during the binding and folding of the CCT protein substrate. Here, we demonstrate that functional divergence occurred after several CCT duplication events due to the fixation of amino acid substitutions by positive selection. Sites critical for ATP binding and substrate binding were found to have undergone positive selection and functional divergence predominantly in subunits that bind tubulin but not actin. Furthermore, we show clear functional divergence between CCT subunits that bind the C-terminal domains of actin and tubulin and those that bind the N-terminal domains. Phylogenetic analyses could not resolve the deep relationships between most subunits, except for the groups alpha/beta/eta and delta/epsilon, suggesting several almost simultaneous ancient duplication events. Together, the results support the idea that, in contrast to homo-oligomeric chaperonins such as GroEL, the high divergence level between CCT subunits is the result of positive selection after each duplication event to provide a specialized role for each CCT subunit in the different steps of protein folding.     

Sequence determinants of protein architecture

Delineation of the relationship between sequence and structure in proteins has proven elusive. Most studies of this problem use alignment methods and other approaches based on the characteristics of individual residues. It is demonstrated herein that the sequence‐structure relationship is determined in significant part by global characteristics of sequence organization. Information encoded in complete sequences is required to distinguish proteins in different architectural groups. It is found that the statistically significant differences between sequences encoding different architectures are encoded in a surprisingly small set of low‐wave‐number sequence periodicities. It would therefore appear that unexpected simplicity in an appropriately defined Fourier space may be an inherent characteristic of the sequences of folded proteins. Proteins 2013; 81:1681–1685. © 2013 Wiley Periodicals, Inc.     

The crystal structure of yeast CCT reveals intrinsic asymmetry of eukaryotic cytosolic chaperonins

The cytosolic chaperonin CCT is a 1‐MDa protein‐folding machine essential for eukaryotic life. The CCT interactome shows involvement in folding and assembly of a small range of proteins linked to essential cellular processes such as cytoskeleton assembly and cell‐cycle regulation. CCT has a classic chaperonin architecture, with two heterogeneous 8‐membered rings stacked back‐to‐back, enclosing a folding cavity. However, the mechanism by which CCT assists folding is distinct from other chaperonins, with no hydrophobic wall lining a potential Anfinsen cage, and a sequential rather than concerted ATP hydrolysis mechanism. We have solved the crystal structure of yeast CCT in complex with actin at 3.8 Å resolution, revealing the subunit organisation and the location of discrete patches of co‐evolving ‘signature residues’ that mediate specific interactions between CCT and its substrates. The intrinsic asymmetry is revealed by the structural individuality of the CCT subunits, which display unique configurations, substrate binding properties, ATP‐binding heterogeneity and subunit–subunit interactions. The location of the evolutionarily conserved N‐terminus of Cct5 on the outside of the barrel, confirmed by mutational studies, is unique to eukaryotic cytosolic chaperonins.     

Computational method for predicting protein functional sites with the use of specificity determinants

The currently available body of decoded amino acid sequences of various proteins exceeds manifold the experimental capabilities of their functional annotation. Therefore, in silico annotation using bioinformatics methods becomes increasingly important. Such annotation is actually a prediction; however, this can be an important starting point for further laboratory research. This work describes a new method for predicting functionally important protein sites, SDPsite, on the basis of identification of specificity determinants. The algorithm proposed utilizes a protein family aglinment and a phylogenetic tree to predict the conserved positions and specificity determinants, map them onto the protein structure, and search for clusters of the predicted positions. Comparison of the resulting predictions with experimental data and published predictions of functional sites by other methods demonstrates that the results of SDPsite agree well with experimental data and exceed the results obtained with the majority of previous methods. SDPsite is publicly available at http://bioinf.fbb.msu.ru/SDPsite.     

Interaction with the CCT chaperonin complex limits APOBEC3A cytidine deaminase cytotoxicity

The APOBEC3 cytidine deaminases are implicated as the cause of a prevalent somatic mutation pattern found in cancer genomes. The APOBEC3 enzymes act as viral restriction factors by mutating viral genomes. Mutation of the cellular genome is presumed to be an off‐target activity of the enzymes, although the regulatory measures for APOBEC3 expression and activity remain undefined. It is therefore difficult to predict circumstances that enable APOBEC3 interaction with cellular DNA that leads to mutagenesis. The APOBEC3A (A3A) enzyme is the most potent deaminase of the family. Using proteomics, we evaluate protein interactors of A3A to identify potential regulators. We find that A3A interacts with the chaperonin‐containing TCP‐1 (CCT) complex, a cellular machine that assists in protein folding and function. Importantly, depletion of CCT results in A3A‐induced DNA damage and cytotoxicity. Evaluation of cancer genomes demonstrates an enrichment of A3A mutational signatures in cancers with silencing mutations in CCT subunit genes. Together, these data suggest that the CCT complex interacts with A3A, and that disruption of CCT function results in increased A3A mutational activity.     

Versatile substrate protein recognition mechanism of the eukaryotic chaperonin CCT

Group II chaperonins, found in eukaryotic and archaeal organisms, recognize substrate proteins through diverse mechanisms that involve either hydrophobic‐ or electrostatic‐dominated interactions. This action is distinct from the universal substrate recognition mechanism of group I chaperonins, which bind a wide spectrum of non‐native proteins primarily through hydrophobic interactions. We use computational approaches to pinpoint the substrate protein binding sites of the γ‐subunit of the eukaryotic chaperonin CCT and to identify its interactions with the stringent substrate β‐tubulin. Protein–protein docking methods reveal intrinsic binding sites of CCT comprising a helical (HL) region, homologous to the GroEL‐binding site, and the helical protrusion (HP) region. We performed molecular dynamics simulations of the solvated CCTγ apical domain, β‐tubulin peptide‐CCTγ complexes, and isolated β‐tubulin peptides. We find that tubulin binds to CCTγ through an extensive interface that spans both the HL region and the HP region. HL interactions involve both hydrophobic and electrostatic contacts, while binding to the HP region is stabilized almost exclusively by a salt bridge network. On the basis of additional simulations of a β‐tubulin‐CCTγ complex that involves a reduced interface, centered onto the HP region, we conclude that this salt bridge network is the minimal stabilizing interaction required. Strong conservation of the charged amino acids that participate in the salt bridge network, Arg306 and Glu271, indicates a general mechanism across the nonidentical CCT subunits and group II chaperonins. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The 13 angstroms structure of a chaperonin GroEL-protein substrate complex by cryo-electron microscopy

The 13 angstroms resolution structures of GroEL bound to a single monomer of the protein substrate glutamine synthetase (GS(m)), as well as that of unliganded GroEL have been determined from a heterogeneous image population using cryo-electron microscopy (cryo-EM) coupled with single-particle image classification and reconstruction techniques. We combined structural data from cryo-EM maps and dynamic modeling, taking advantage of the known X-ray crystallographic structure and normal mode flexible fitting (NMFF) analysis, to describe the changes that occur in GroEL structure induced by GS(m) binding. The NMFF analysis reveals that the molecular movements induced by GS(m) binding propagate throughout the GroEL structure. The modeled molecular motions show that some domains undergo en bloc movements, while others show more complex independent internal movements. Interestingly, the substrate-bound apical domains of both the cis (GS(m)-bound ring) and trans (the opposite substrate-free ring) show counterclockwise rotations, in the same direction (though not as dramatic) as those documented for the ATP-GroEL-induced structure changes. The structural changes from the allosteric substrate protein-induced negative cooperativity between the GroEL rings involves upward concerted movements of both cis and trans equatorial domains toward the GS(m)-bound ring, while the inter-ring distances between the heptamer contact residues are maintained. Furthermore, the NMFF analysis identifies the secondary structural elements that are involved in the observed approximately 5 angstroms reduction in the diameter of the cavity opening in the unbound trans ring. Understanding the molecular basis of these substrate protein-induced structural changes across the heptamer rings provides insight into the origins of the allosteric negative cooperative effects that are transmitted over long distances (approximately 140 angstroms).     

Folding, stability and polymerization properties of FtsZ chimeras with inserted tubulin loops involved in the interaction with the cytosolic chaperonin CCT and in microtubule formation

To attain its native conformation, the cytoskeletal protein tubulin needs the concourse of several molecular chaperones, among others the cytosolic chaperonin CCT. It has been previously described that denatured tubulin interacts with CCT in a quasi-folded conformation using several loops located throughout its sequence. These loops are also involved in microtubule formation and are absent in its prokaryote homologue FtsZ, which in vitro folds by itself and does not interact with CCT. Several FtsZ/tubulin chimeric proteins were generated by inserting consecutively one, two or three of the CCT-binding domains of tubulin into the corresponding sequence of FtsZ from Methanococccus jannaschii. The insertion of any of the CCT-binding loops generates in the FtsZ/tubulin chimeras the ability to interact with CCT. The accumulation of CCT-binding loops induces in the FtsZ/tubulin chimeras unfolding and refolding properties that are more similar to tubulin than to its prokaryote counterpart. Finally, the insertion of some of these loops generates in the FtsZ/tubulin chimeras more complex polymeric structures than those found for FtsZ. These results reinforce the notion that CCT has coevolved with tubulin to deal with the folding problems encountered by the eukaryotic protein with the appearance of the new sequences involved in microtubule formation.     

Identification of GATC- and CCGG-recognizing Type II REases and their putative specificity-determining positions using Scan2S--a novel motif scan algorithm with optional secondary structure constraints

Restriction endonucleases (REases) are DNA-cleaving enzymes that have become indispensable tools in molecular biology. Type II REases are highly divergent in sequence despite their common structural core, function and, in some cases, common specificities towards DNA sequences. This makes it difficult to identify and classify them functionally based on sequence, and has hampered the efforts of specificity-engineering. Here, we define novel REase sequence motifs, which extend beyond the PD-(D/E)XK hallmark, and incorporate secondary structure information. The automated search using these motifs is carried out with a newly developed fast regular expression matching algorithm that accommodates long patterns with optional secondary structure constraints. Using this new tool, named Scan2S, motifs derived from REases with specificity towards GATC- and CGGG-containing DNA sequences successfully identify REases of the same specificity. Notably, some of these sequences are not identified by standard sequence detection tools. The new motifs highlight potential specificity-determining positions that do not fully overlap for the GATC- and the CCGG-recognizing REases and are candidates for specificity re-engineering.     

Unfolding study of a trimeric membrane protein AcrB

The folding of a multi‐domain trimeric α‐helical membrane protein, Escherichia coli inner membrane protein AcrB, was investigated. AcrB contains both a transmembrane domain and a large periplasmic domain. Protein unfolding in sodium dodecyl sulfate (SDS) and urea was monitored using the intrinsic fluorescence and circular dichroism spectroscopy. The SDS denaturation curve displayed a sigmoidal profile, which could be fitted with a two‐state unfolding model. To investigate the unfolding of separate domains, a triple mutant was created, in which all three Trp residues in the transmembrane domain were replaced with Phe. The SDS unfolding profile of the mutant was comparable to that of the wild type AcrB, suggesting that the observed signal change was largely originated from the unfolding of the soluble domain. Strengthening of trimer association through the introduction of an inter‐subunit disulfide bond had little effect on the unfolding profile, suggesting that trimer dissociation was not the rate‐limiting step in unfolding monitored by fluorescence emission. Under our experimental condition, AcrB unfolding was not reversible. Furthermore, we experimented with the refolding of a monomeric mutant, AcrB_Δloop, from the SDS unfolded state. The CD spectrum of the refolded AcrB_Δloop superimposed well onto the spectra of the original folded protein, while the fluorescence spectrum was not fully recovered. In summary, our results suggested that the unfolding of the trimeric AcrB started with a local structural rearrangement. While the refolding of secondary structure in individual monomers could be achieved, the re‐association of the trimer might be the limiting factor to obtain folded wild‐type AcrB.     

Inter-ring communication allows the GroEL chaperonin complex to distinguish between different substrates

The productive folding of substrate proteins by the GroEL complex of Escherichia coli requires the activity of both the chaperonin rings. These heptameric rings were shown to regulate the chaperonins' affinity for substrates and co-chaperonin via inter-ring communications; however, the molecular details of the interactions are not well understood. We have investigated the effect of substrate binding on inter-ring communications of the chaperonin complex, both the double-ring GroEL as well as the single-ring SR1 chaperonin in complex with four different substrates by using mass spectrometry. This approach shows that whereas SR1 is unable to distinguish between Rubisco, gp23, gp5, and MDH, GroEL shows clear differences upon binding these substrates. The most distinctive binding behavior is observed for Rubisco, which only occupies one GroEL ring. Both bacteriophage capsid proteins (gp23 and gp5) as well as MDH are able to bind to the two GroEL rings simultaneously. Our data suggest that inter-ring communication allows the chaperonin complex to differentiate between substrates. Using collision induced dissociation in the gas phase, differences between the chaperonin(substrate) complexes are observed only when both rings are present. The data indicate that the size of the substrate is an important factor that determines the degree of stabilization of the chaperonin complex.     

Essential role of the chaperonin folding compartment in vivo

The GroEL/GroES chaperonin system of Escherichia coli forms a nano-cage allowing single protein molecules to fold in isolation. However, as the chaperonin can also mediate folding independently of substrate encapsulation, it remained unclear whether the folding cage is essential in vivo. To address this question, we replaced wild-type GroEL with mutants of GroEL having either a reduced cage volume or altered charge properties of the cage wall. A stepwise reduction in cage size resulted in a gradual loss of cell viability, although the mutants bound non-native protein efficiently. Strikingly, a mild reduction in cage size increased the yield and the apparent rate of green fluorescent protein folding, consistent with the view that an effect of steric confinement can accelerate folding. As shown in vitro, the observed acceleration of folding was dependent on protein encapsulation by GroES but independent of GroES cycling regulated by the GroEL ATPase. Altering the net-negative charge of the GroEL cage wall also strongly affected chaperonin function. Based on these findings, the GroEL/GroES compartment is essential for protein folding in vivo.     

Predicting the structure of the light-harvesting complex II of Rhodospirillum molischianum.

We attempted to predict through computer modeling the structure of the light-harvesting complex II (LH-II) of Rhodospirillum molischianum, before the impending publication of the structure of a homologous protein solved by means of X-ray diffraction. The protein studied is an integral membrane protein of 16 independent polypeptides, 8 alpha-apoproteins and 8 beta-apoproteins, which aggregate and bind to 24 bacteriochlorophyll-a's and 12 lycopenes. Available diffraction data of a crystal of the protein, which could not be phased due to a lack of heavy metal derivatives, served to test the predicted structure, guiding the search. In order to determine the secondary structure, hydropathy analysis was performed to identify the putative transmembrane segments and multiple sequence alignment propensity analyses were used to pinpoint the exact sites of the 20-residue-long transmembrane segment and the 4-residue-long terminal sequence at both ends, which were independently verified and improved by homology modeling. A consensus assignment for the secondary structure was derived from a combination of all the prediction methods used. Three-dimensional structures for the alpha- and the beta-apoprotein were built by comparative modeling. The resulting tertiary structures are combined, using X-PLOR, into an alpha beta dimer pair with bacteriochlorophyll-a's attached under constraints provided by site-directed mutagenesis and spectral data. The alpha beta dimer pairs were then aggregated into a quaternary structure through further molecular dynamics simulations and energy minimization. The structure of LH-II so determined is an octamer of alpha beta heterodimers forming a ring with a diameter of 70 A.     

Symmetry-free cryo-EM structures of the chaperonin TRiC along its ATPase-driven conformational cycle

The eukaryotic group II chaperonin TRiC/CCT is a 16-subunit complex with eight distinct but similar subunits arranged in two stacked rings. Substrate folding inside the central chamber is triggered by ATP hydrolysis. We present five cryo-EM structures of TRiC in apo and nucleotide-induced states without imposing symmetry during the 3D reconstruction. These structures reveal the intra- and inter-ring subunit interaction pattern changes during the ATPase cycle. In the apo state, the subunit arrangement in each ring is highly asymmetric, whereas all nucleotide-containing states tend to be more symmetrical. We identify and structurally characterize an one-ring closed intermediate induced by ATP hydrolysis wherein the closed TRiC ring exhibits an observable chamber expansion. This likely represents the physiological substrate folding state. Our structural results suggest mechanisms for inter-ring-negative cooperativity, intra-ring-positive cooperativity, and protein-folding chamber closure of TRiC. Intriguingly, these mechanisms are different from other group I and II chaperonins despite their similar architecture.     

The N terminus of the head protein of T4 bacteriophage directs proteins to the GroEL chaperonin

The head protein of T4 bacteriophage requires the GroEL chaperonin for its insertion into a growing T4 head. Hundreds of thousands of copies of this protein must pass through the chaperonin in a limited time later in infection, indicating that the protein must use GroEL very efficiently and may contain sequences that bind tightly to GroEL. We show that green fluorescent protein (GFP) fused to the N terminus of the head protein can fold at temperatures higher than those at which the GFP protein can fold well by itself. We present evidence that this folding is promoted by the strong binding of N-terminal head protein sequences to GroEL. This binding is so strong that some fusion proteins can apparently deplete the cell of the GroEL needed for other cellular functions, altering the cellular membranes and slowing growth.     

Penicillin-binding protein 5 can form a homo-oligomeric complex in the inner membrane of Escherichia coli

Penicillin-binding protein 5 (PBP5) is a DD-carboxypeptidase, which cleaves the terminal D-alanine from the muramyl pentapeptide in the peptidoglycan layer of Escherichia coli and other bacteria. In doing so, it varies the substrates for transpeptidation and plays a key role in maintaining cell shape. In this study, we have analyzed the oligomeric state of PBP5 in detergent and in its native environment, the inner membrane. Both approaches indicate that PBP5 exists as a homo-oligomeric complex, most likely as a homo-dimer. As the crystal structure of the soluble domain of PBP5 (i.e., lacking the membrane anchor) shows a monomer, we used our experimental data to generate a model of the homo-dimer. This model extends our understanding of PBP5 function as it suggests how PBP5 can interact with the peptidoglycan layer. It suggests that the stem domains interact and the catalytic domains have freedom to move from the position observed in the crystal structure. This would allow the catalytic domain to have access to pentapeptides at different distances from the membrane.     

The chaperonin TRiC forms an oligomeric complex in the malaria parasite cytosol

The malaria parasite exports numerous proteins into its host red blood cell (RBC). The trafficking of these exported effectors is complex. Proteins are first routed through the secretory system, into the parasitophorous vacuole (PV), a membranous compartment enclosing the parasite. Proteins are then translocated across the PV membrane in a process requiring ATP and unfolding. Once in the RBC compartment the exported proteins are then refolded and further trafficked to their final localizations. Chaperones are important in the unfolding and refolding processes. Recently, it was suggested that the parasite TRiC chaperonin complex is exported, and that it is involved in trafficking of exported effectors. Using a parasite‐specific antibody and epitope‐tagged transgenic parasites we could observe no export of Plasmodium TRiC into the RBC. We tested the importance of the parasite TRiC by creating a regulatable knockdown line of the TRiC‐θ subunit. Loss of the parasite TRiC‐θ led to a severe growth defect in asexual development, but did not alter protein export into the RBC. These observations indicate that the TRiC proteins play a critical role in parasite biology, though their function, within the parasite, appears unrelated to protein trafficking in the RBC compartment.     

Increased sequence hydrophobicity reduces conformational specificity: A mutational case study of the Arc repressor protein

The amino-acid sequences of soluble, globular proteins must have hydrophobic residues to form a stable core, but excess sequence hydrophobicity can lead to loss of native state conformational specificity and aggregation. Previous studies of polar-to-hydrophobic mutations in the β-sheet of the Arc repressor dimer showed that a single substitution at position 11 (N11L) leads to population of an alternate dimeric fold in which the β-sheet is replaced by helix. Two additional hydrophobic mutations at positions 9 and 13 (Q9V and R13V) lead to population of a differently folded octamer along with both dimeric folds. Here we conduct a comprehensive study of the sequence determinants of this progressive loss of fold specificity. We find that the alternate dimer-fold specifically results from the N11L substitution and is not promoted by other hydrophobic substitutions in the β-sheet. We also find that three highly hydrophobic substitutions at positions 9, 11, and 13 are necessary and sufficient for oligomer formation, but the oligomer size depends on the identity of the hydrophobic residue in question. The hydrophobic substitutions increase thermal stability, illustrating how increased hydrophobicity can increase folding stability even as it degrades conformational specificity. The oligomeric variants are predicted to be aggregation-prone but may be hindered from doing so by proline residues that flank the β-sheet region. Loss of conformational specificity due to increased hydrophobicity can manifest itself at any level of structure, depending upon the specific mutations and the context in which they occur.