Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis

Proteins are intrinsically flexible molecules. The role of internal motions in a protein''s designated function is widely debated. The role of protein structure in enzyme catalysis is well established, and conservation of structural features provides vital clues to their role in function. Recently, it has been proposed that the protein function may involve multiple conformations: the observed deviations are not random thermodynamic fluctuations; rather, flexibility may be closely linked to protein function, including enzyme catalysis. We hypothesize that the argument of conservation of important structural features can also be extended to identification of protein flexibility in interconnection with enzyme function. Three classes of enzymes (prolyl-peptidyl isomerase, oxidoreductase, and nuclease) that catalyze diverse chemical reactions have been examined using detailed computational modeling. For each class, the identification and characterization of the internal protein motions coupled to the chemical step in enzyme mechanisms in multiple species show identical enzyme conformational fluctuations. In addition to the active-site residues, motions of protein surface loop regions (>10 Å away) are observed to be identical across species, and networks of conserved interactions/residues connect these highly flexible surface regions to the active-site residues that make direct contact with substrates. More interestingly, examination of reaction-coupled motions in non-homologous enzyme systems (with no structural or sequence similarity) that catalyze the same biochemical reaction shows motions that induce remarkably similar changes in the enzyme–substrate interactions during catalysis. The results indicate that the reaction-coupled flexibility is a conserved aspect of the enzyme molecular architecture. Protein motions in distal areas of homologous and non-homologous enzyme systems mediate similar changes in the active-site enzyme–substrate interactions, thereby impacting the mechanism of catalyzed chemistry. These results have implications for understanding the mechanism of allostery, and for protein engineering and drug design.

Author''s Summary

Enzymes are nature''s molecular machines that catalyze biochemical reactions with remarkable efficiency. Recent evidence suggests that enzyme function may involve not only direct structural interactions between the enzyme and its substrate, but also internal motions of the enzyme itself. Here, we describe a computational investigation of three classes of enzymes that catalyze completely different biochemical reactions. Remarkably, the mobile enzyme regions and the nature of these motions are the same across species ranging from single-celled organisms to complex life-forms. Also surprisingly, non-homologous enzymes that catalyze the same chemical reaction but do not share sequence or structural similarity reveal a similar impact of enzyme motions on their reaction mechanisms. Flexible enzyme regions are found to be connected by conserved networks of coupled interactions that connect surface regions to active-site residues. These networks may provide a mechanism for the solvent on an enzyme''s surface to couple to the reaction catalyzed by the enzyme. These results have implications for understanding the mechanism of allostery (long-range effects), and for protein engineering and drug design.     

Bi-Factor Analysis Based on Noise-Reduction (BIFANR): A New Algorithm for Detecting Coevolving Amino Acid Sites in Proteins

Previous statistical analyses have shown that amino acid sites in a protein evolve in a correlated way instead of independently. Even though located distantly in the linear sequence, the coevolved amino acids could be spatially adjacent in the tertiary structure, and constitute specific protein sectors. Moreover, these protein sectors are independent of one another in structure, function, and even evolution. Thus, systematic studies on protein sectors inside a protein will contribute to the clarification of protein function. In this paper, we propose a new algorithm BIFANR (Bi-factor Analysis Based on Noise-reduction) for detecting protein sectors in amino acid sequences. After applying BIFANR on S1A family and PDZ family, we carried out internal correlation test, statistical independence test, evolutionary rate analysis, evolutionary independence analysis, and function analysis to assess the prediction. The results showed that the amino acids in certain predicted protein sector are closely correlated in structure, function, and evolution, while protein sectors are nearly statistically independent. The results also indicated that the protein sectors have distinct evolutionary directions. In addition, compared with other algorithms, BIFANR has higher accuracy and robustness under the influence of noise sites.     

Glycosylation Site Analysis of Human Platelets by Electrostatic Repulsion Hydrophilic Interaction Chromatography

Introduction

Glycosylations range among the most common posttranslational modifications with an estimated 50% of all proteins supposed to be glycosylated. These modifications are required for essential cellular processes including cell–cell recognition, protein structure and activity, e.g., of surface receptors, as well as subcellular localization of proteins. Beside the elucidation of the carbohydrate structures, the annotation of glycosylation sites is of primary interest as a basis for subsequent functional characterization. Although mass spectrometry is the method of choice for large-scale analysis of glycosylation sites, it requires initial enrichment of glycopeptides prior mass spectrometric detection in most cases.

Materials and Methods

In this paper, we present a novel approach for glycopeptide enrichment by electrostatic repulsion hydrophilic interaction chromatography (ERLIC). Glycopeptides were separated from the bulk of non-modified peptides and gradually eluted from the stationary phase with potential for isoform resolution. Applied to human platelets, 125 glycosylation sites on 66 proteins were identified including major platelet glycoproteins responsible for cellular function.

Conclusion

These sites add a major contribution to the now more than 250 glycosylation sites annotated for platelets, which enable the clinically relevant design of quantification assays for platelet glycoproteins.     

Structure and increased thermostability of Rhodococcus sp. naphthalene 1,2-dioxygenase

Rieske nonheme iron oxygenases form a large class of aromatic ring-hydroxylating dioxygenases found in microorganisms. These enzymes enable microorganisms to tolerate and even exclusively utilize aromatic compounds for growth, making them good candidates for use in synthesis of chiral intermediates and bioremediation. Studies of the chemical stability and thermostability of these enzymes thus become important. We report here the structure of free and substrate (indole)-bound forms of naphthalene dioxygenase from Rhodococcus sp. strain NCIMB12038. The structure of the Rhodococcus enzyme reveals that, despite a approximately 30% sequence identity between these naphthalene dioxygenases, their overall structures superpose very well with a root mean square deviation of less than 1.6 A. The differences in the active site of the two enzymes are pronounced near the entrance; however, indole binds to the Rhodococcus enzyme in the same orientation as in the Pseudomonas enzyme. Circular dichroism spectroscopy experiments show that the Rhodococcus enzyme has higher thermostability than the naphthalene dioxygenase from Pseudomonas species. The Pseudomonas enzyme has an apparent melting temperature of 55 degrees C while the Rhodococcus enzyme does not completely unfold even at 95 degrees C. Both enzymes, however, show similar unfolding behavior in urea, and the Rhodococcus enzyme is only slightly more tolerant to unfolding by guanidine hydrochloride. Structure analysis suggests that the higher thermostability of the Rhodococcus enzyme may be attributed to a larger buried surface area and extra salt bridge networks between the alpha and beta subunits in the Rhodococcus enzyme.     

Structural insights into the U-box,a domain associated with multi-ubiquitination

The structure of the U-box in the essential Saccharomyces cerevisiae pre-mRNA splicing factor Prp19p has been determined by NMR. The conserved zinc-binding sites supporting the cross-brace arrangement in RING-finger domains are replaced by hydrogen-bonding networks in the U-box. These hydrogen-bonding networks are necessary for the structural stabilization and activity of the U-box. A conservative Val-->Ile point mutation in the Prp19p U-box domain leads to pre-mRNA splicing defects in vivo. NMR analysis of this mutant shows that the substitution disrupts structural integrity of the U-box domain. Furthermore, comparison of the Prp19p U-box domain with known RING-E2 complex structures demonstrates that both U-box and RING-fingers contain a conserved interaction surface. Mutagenesis of residues at this interface, while not perturbing the structure of the U-box, abrogates Prp19p function in vivo. These comparative structural and functional analyses imply that the U-box and its associated ubiquitin ligase activity are critical for Prp19p function in vivo.     

Changes in the relative occupancy of metal-binding sites in the profile structure of the sarcoplasmic reticulum membrane induced by phosphorylation of the Ca2+ATPase enzyme in the presence of terbium: a time-resolved, resonance x-ray diffraction study.

Time-resolved, terbium resonance x-ray diffraction experiments have provided the locations of three different high-affinity metal-binding/transport sites on the Ca2+ATPase enzyme in the profile structure of the sarcoplasmic reticulum (SR) membrane. By considering these results in conjunction with the known, moderate-resolution profile structure of the SR membrane (derived from nonresonance x-ray and neutron diffraction studies), it was determined that the three metal-binding sites are located at the "headpiece/stalk" junction in the Ca2+ATPase profile structure, in the "transbilayer" portion of the enzyme profile near the center of the membrane phospholipid bilayer, and at the intravesicular surface of the membrane profile. All three metal-binding sites so identified are simultaneously occupied in the unphosphorylated enzyme conformation. Phosphorylation of the ATPase causes a redistribution of metal density among the sites, resulting in a net movement of metal density toward the intravesicular side of the membrane, i.e., in the direction of calcium active transport. We propose that this redistribution of metal density is caused by changes in the relative binding affinities of the three sites, mediated by local structural changes at the sites resulting from the large-scale (i.e., long-range) changes in the profile structure of the Ca2+ATPase induced by phosphorylation, as reported in an accompanying paper. The implications of these results for the mechanism of calcium active transport by the SR Ca2+ATPase are discussed briefly.     

Residue centrality, functionally important residues, and active site shape: analysis of enzyme and non-enzyme families

The representation of protein structures as small-world networks facilitates the search for topological determinants, which may relate to functionally important residues. Here, we aimed to investigate the performance of residue centrality, viewed as a family fold characteristic, in identifying functionally important residues in protein families. Our study is based on 46 families, including 29 enzyme and 17 non-enzyme families. A total of 80% of these central positions corresponded to active site residues or residues in direct contact with these sites. For enzyme families, this percentage increased to 91%, while for non-enzyme families the percentage decreased substantially to 48%. A total of 70% of these central positions are located in catalytic sites in the enzyme families, 64% are in hetero-atom binding sites in those families binding hetero-atoms, and only 16% belong to protein-protein interfaces in families with protein-protein interaction data. These differences reflect the active site shape: enzyme active sites locate in surface clefts, hetero-atom binding residues are in deep cavities, while protein-protein interactions involve a more planar configuration. On the other hand, not all surface cavities or clefts are comprised of central residues. Thus, closeness centrality identifies functionally important residues in enzymes. While here we focus on binding sites, we expect to identify key residues for the integration and transmission of the information to the rest of the protein, reflecting the relationship between fold and function. Residue centrality is more conserved than the protein sequence, emphasizing the robustness of protein structures.     

Structures at the proteolytic processing region of cathepsin D

The amino acid sequences at the "proteolytic processing regions" of cathepsin Ds have been determined for the enzymes from cows, pigs, and rats in order to deduce the sites of cleavage as well as the function of the proteolytic processing of cathepsin D. For bovine cathepsin D, the "processing region" sequence was determined from a peptide isolated from the single-chain enzyme. The COOH-terminal sequence of the light chain and the NH2-terminal sequence of the heavy chain were also determined. The processing region sequence of porcine cathepsin D was determined from its cDNA structure, and the same structure from rat cathepsin D was determined from the peptide sequence of the single-chain rat enzyme. From sequence homology to other aspartic proteases whose x-ray crystallographic structures are known, such as pepsinogen and penicillopepsin, it is clear that the processing regions are insertions to form an extended beta-hairpin loop between residues 91 and 92 (porcine pepsin numbers). However, the sizes of the processing regions of cathepsin Ds from different species are considerably different. For the enzymes from rats, cows, pigs, and human, the sizes of the processing regions are 6, 9, 9, and 11 amino acid residues, respectively. The amino acid sequences within the processing regions are considerably different. In addition, the proteolytic processing sites were found to be completely different in the bovine and porcine cathepsin Ds. While in the porcine enzyme, an Asn-Ser bond and a Gly-Val bond are cleaved to release 5 residues as a consequence of the processing; in the bovine enzyme, two Ser-Ser bonds are cleaved to release 2 serine residues. These findings would argue that the in vivo proteolytic processing of the cathepsin D single chain is probably not carried out by a specific "processing protease." Model building of the cathepsin D processing region conformation was conducted utilizing the homology between procathepsin D and porcine pepsinogen. The beta-hairpin structure of the processing region was found to (i) interact with the activation peptide of the procathepsin D in a beta-structure and (ii) place the Cys residue in the processing region within disulfide linkage distance to Cys-27 of cathepsin D light chain. These observations support the view that the processing region of cathepsin D may function to stabilize the conformation of procathepsin D and may play a role in its activation.     

Similar binding sites and different partners: implications to shared proteins in cellular pathways

We studied a data set of structurally similar interfaces that bind to proteins with different binding-site structures and different functions. Our multipartner protein interface clusters enable us to address questions like: What makes a given site bind different proteins? How similar/different are the interactions? And, what drives the apparently less-specific association? We find that proteins with common binding-site motifs preferentially use conserved interactions at similar interface locations, despite the different partners. Helices are major vehicles for binding different partners, allowing alternate ways to achieve favorable association. The binding sites are characterized by imperfect packing, planar architectures, bridging water molecules, and, on average, smaller size. Interestingly, analysis of the connectivity of these proteins illustrates that they have more interactions with other proteins. These findings are important in predicting "date hubs," if we assume that "date hubs" are shared proteins with binding sites capable of transient binding to multipartners, linking higher-order networks.     

Organization of enzymes on subcellular structures: relay at the surface

There is a large body of evidence that soluble cytoplasmic enzymes of eukaryotic cells, e.g., glycolytic enzymes and proteins of the translational machinery, are organized in some way in space and in time. The following features of such organization emerge from the experimental data: (1) metabolites are transferred between enzymes directly "from hand to hand" in short-living enzyme-enzyme complexes rather than by diffusion in aqueous media; (2) enzymes show a tendency to be absorbed on surfaces of subcellular structures, such as membranes, cytoskeleton and polyribosomes; (3) enzymes are desorbed from a surface of a subcellular structure after binding specific metabolites, i.e., substrates and/or products of the reactions catalyzed by these enzymes. These features are suggestive of a relay mechanism for the enzyme systems functioning in a cell; an enzyme adsorbed on a surface of a subcellular structure is desorbed after binding its substrate or in the course of the catalytic act. Within a complex with its product the enzyme diffuses into the environment, until it reaches the next enzyme adsorbed on the same surface; then a short-living enzyme-enzyme complex is formed, and a direct "from hand to hand" transfer of the metabolite takes place. As a result, the overall metabolic process appears to be localized near the surface. We termed this mechanism as a "relay at the surface".     

Glutathione reductase and thioredoxin reductase at the crossroad: the structure of Schistosoma mansoni thioredoxin glutathione reductase

Thioredoxin glutathione reductase (TGR) is a key flavoenzyme expressed by schistosomes that bridges two detoxification pathways crucial for the parasite survival in the host's organism. In this article we report the crystal structure (at 2.2 A resolution) of TGR from Schistosoma mansoni (SmTGR), deleted in the last two residues. The structure reveals the peculiar architecture of this chimeric enzyme: the small Glutaredoxin (Grx) domain at the N-terminus is joined to the large thioredoxin reductase (TR) one via an extended complementary surface, involving residues not conserved in the Grx superfamily; the TR domain interacts with an identical partner via its C-terminal domain, forming a dimer with a twisted "W" shape. Although lacking the penultimate Selenocysteine residue (Sec), the enzyme is still able to reduce oxidized glutathione. These data update the interpretation of the interdomain communication in TGR enzymes. The possible function of this enzyme in pathogenic parasites is discussed.     

zapERtrap: A light-regulated ER release system reveals unexpected neuronal trafficking pathways

Here we introduce zapalog-mediated endoplasmic reticulum trap (zapERtrap), which allows one to use light to precisely trigger forward trafficking of diverse integral membrane proteins from internal secretory organelles to the cell surface with single cell and subcellular spatial resolution. To demonstrate its utility, we use zapERtrap in neurons to dissect where synaptic proteins emerge at the cell surface when processed through central (cell body) or remote (dendrites) secretory pathways. We reveal rapid and direct long-range trafficking of centrally processed proteins deep into the dendritic arbor to synaptic sites. Select proteins were also trafficked to the plasma membrane of the axon initial segment, revealing a novel surface trafficking hotspot. Proteins locally processed through dendritic secretory networks were widely dispersed before surface insertion, challenging assumptions for precise trafficking at remote sites. These experiments provide new insights into compartmentalized secretory trafficking and showcase the tunability and spatiotemporal control of zapERtrap, which will have broad applications for regulating cell signaling and function.     

Calcium modulates endopeptidase 24.15 (EC 3.4.24.15) membrane association, secondary structure and substrate specificity

The metalloendopeptidase 24.15 (EP24.15) is ubiquitously present in the extracellular environment as a secreted protein. Outside the cell, this enzyme degrades several neuropeptides containing from 5 to 17 amino acids (e.g. gonadotropin releasing hormone, bradykinin, opioids and neurotensin). The constitutive secretion of EP24.15 from glioma C6 cells was demonstrated to be stimulated linearly by reduced concentrations of extracellular calcium. In the present report we demonstrate that extracellular calcium concentration has no effect on the total amount of the extracellular (cell associated + medium) enzyme. Indeed, immuno-cytochemical analyses by confocal and electron microscopy suggested that the absence of calcium favors the enzyme shedding from the plasma membrane into the medium. Two putative calcium-binding sites on EP24.15 (D93 and D159) were altered by site-directed mutagenesis to investigate their possible contribution to binding of the enzyme at the cell surface. These mutated recombinant proteins behave similarly to the wild-type enzyme regarding enzymatic activity, secondary structure, calcium sensitivity and immunoreactivity. However, immunocytochemical analyses by confocal microscopy consistently show a reduced ability of the D93A mutant to associate with the plasma membrane of glioma C6 cells when compared with the wild-type enzyme. These data and the model of the enzyme's structure as determined by X-ray diffraction suggest that D93 is located at the enzyme surface and is consistent with membrane association of EP24.15. Moreover, calcium was also observed to induce a major change in the EP24.15 cleavage site on distinctive fluorogenic substrates. These data suggest that calcium may be an important modulator of ep24.15 cell function.     

DNA abasic lesions in a different light: solution structure of an endogenous topoisomerase II poison

Topoisomerase II is the target for several anticancer drugs that "poison" the enzyme and convert it to a cellular toxin by increasing topoisomerase II-mediated DNA cleavage. In addition to these "exogenous topoisomerase II poisons," DNA lesions such as abasic sites act as "endogenous poisons" of the enzyme. Drugs and lesions are believed to stimulate DNA scission by altering the structure of the double helix within the cleavage site of the enzyme. However, the structural alterations that enhance cleavage are unknown. Since abasic sites are an intrinsic part of the genetic material, they represent an attractive model to assess DNA distortions that lead to altered topoisomerase II function. Therefore, the structure of a double-stranded dodecamer containing a tetrahydrofuran apurinic lesion at the +2 position of a topoisomerase II DNA cleavage site was determined by NMR spectroscopy. Three major features distinguished the apurinic structure ( = 0.095) from that of wild-type ( = 0.077). First, loss of base stacking at the lesion collapsed the major groove and reduced the distance between the two scissile phosphodiester bonds. Second, the apurinic lesion induced a bend that was centered about the topoisomerase II cleavage site. Third, the base immediately opposite the lesion was extrahelical and relocated to the minor groove. All of these structural alterations have the potential to influence interactions between topoisomerase II and its DNA substrate.     

Folding on the chaperone: yield enhancement through loose binding

A variety of small cageless chaperones have been discovered that can assist protein folding without the consumption of ATP. These include mini-chaperones (catalytically active fragments of larger chaperones), as well as small proteins such as alpha-casein and detergents acting as "artificial chaperones." These chaperones all possess exposed hydrophobic patches on their surface that act as recognition sites for misfolded proteins. They lack the complexity of chaperonins (that encapsulate proteins in their inner rings) and their study can offer insight into the minimal requirements for chaperone function. We use molecular dynamics simulations to investigate how a cageless chaperone, modeled as a sphere of tunable hydrophobicity, can assist folding of a substrate protein. We find that under steady-state (non-stress) conditions, cageless chaperones that bind to a single substrate protein increase folding yields by reducing the time the substrate spends in an aggregation-prone state in a dual manner: (a) by competing for aggregation-prone hydrophobic sites on the surface of a protein, hence reducing the time the protein spends unprotected in the bulk and (b) by accelerating folding rates of the protein. In both cases, the chaperone must bind to and hold the protein loosely enough to allow the protein to change its conformation and fold while bound. Loose binding may enable small cageless chaperones to help proteins fold and avoid aggregation under steady-state conditions, even at low concentrations, without the consumption of ATP.