Two-dimensional 1H NMR studies of histidine-containing protein from Escherichia coli. 3. Secondary and tertiary structure as determined by NMR

Sequence-specific resonance assignments of the 1H NMR spectrum of the 85-residue histidine-containing phosphocarrier protein (HPr) are complete [Klevit, R. E., Drobny, G. P., & Waygood, E. B. (1986) Biochemistry (first paper of three in this issue)]. Additional side-chain assignments have been made with long-range coherence transfer experiments [Klevit, R. E., & Drobny, G. P. (1986) Biochemistry (second paper of three in this issue)]. In this paper, the NMR assignments were used to determine the secondary structure and the tertiary folding of HPr in solution. The secondary structural elements of the protein were determined by visual inspection of the pattern of nearest-neighbor nuclear Overhauser effects (NOEs) and the presence of persistent amide resonances. Escherichia coli HPr consists of four beta-strands, three alpha-helices, four reverse turns, and several regions of extended backbone structure. Long-range NOEs, especially among side-chain protons, were used to determine the tertiary structure of the protein by use of the secondary structural components. The four beta-strands form a single antiparallel beta-pleated sheet. The hydrophobic faces of the alpha-helices interact to form a hydrophobic core and sit above the hydrophobic face of the beta-sheet, forming an open-face beta-sheet sandwich structure. The active site histidine, His-15, is on a short kinked segment of backbone that is accessible to the solvent. The positively charged phosphorylation site (His-15 and Arg-17) interacts with the negatively charged carboxyl terminus of the protein (Glu-85).(ABSTRACT TRUNCATED AT 250 WORDS)     

Structural plasticity associated with the beta-propeller architecture

The beta-propeller architecture observed in protein tertiary structure and classified into the five different types according to number of 'blades' (or beta-sheets) and a sixth type classified according to the secondary structure composition of the blades (the beta beta alpha beta-molecular unit) is characterized by variations (or plasticity) in the structure. These correspond to the number of beta-strands associated with the blade, the number of amino acid residues associated with equivalent beta-strands in the different blades and the presence of alpha-helices and twisted beta-strands. We have generated a beta-sheet associated beta-strand pattern that may be important for protein structure prediction and modeling. Analysis of the beta-propellers extracted primarily from the SCOP database revealed there are 179 beta-propellers. The examination of the secondary structure corresponding to the beta-propeller using PDBsum that was useful to define the beta-sheet associated beta-strand pattern, combined with visualization on graphics display revealed structural plasticity associated with the beta-propeller architecture. Particularly, the type 6- and 7-bladed beta-propellers known to be associated with sequence and functional diversity are more common and associated with relatively more structural variations compared to the other beta-propeller types.     

Multinuclear magnetic resonance studies of Escherichia coli adenylate kinase in free and bound forms. Resonance assignment, secondary structure and ligand binding.

The crystal structure of Escherichia coli adenylate kinase (AKe) revealed three main components: a CORE domain, composed of a five-stranded parallel beta-sheet surrounded by alpha-helices, and two peripheral domains involved in covering the ATP in the active site (LID) and binding of the AMP (NMPbind). We initiated a long-term NMR study aiming to characterize the solution structure, binding mechanism and internal dynamics of the various domains. Using single (15N) and double-labeled (13C and 15N) samples and double- and triple-resonance NMR experiments we assigned 97% of the 1H, 13C and 15N backbone resonances, and proton and 13Cbeta resonances for more than 40% of the side chains in the free protein. Analysis of a 15N-labeled enzyme in complex with the bi-substrate analogue [P1,P5-bis(5'-adenosine)-pentaphosphate] (Ap5A) resulted in the assignment of 90% of the backbone 1H and 15N resonances and 42% of the side chain resonances. Based on short-range NOEs and 1H and 13C secondary chemical shifts, we identified the elements of secondary structure and the topology of the beta-strands in the unliganded form. The alpha-helices and the beta-strands of the parallel beta-sheet in solution have the same limits (+/- 1 residue) as those observed in the crystal. The first helix (alpha1) appears to have a frayed N-terminal side. Significant differences relative to the crystal were noticed in the LID domain, which in solution exhibits four antiparallel beta-strands. The secondary structure of the nucleoside-bound form, as deduced from intramolecular NOEs and the 1Halpha chemical shifts, is similar to that of the free enzyme. The largest chemical shift differences allowed us to map the regions of protein-ligand contacts. 1H/2H exchange experiments performed on free and Ap5A-bound enzymes showed a general decrease of the structural flexibility in the complex which is accompanied by a local increased flexibility on the N-side of the parallel beta-sheet.     

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.     

RNA recognition motifs: boring? Not quite

The RNA recognition motif (RRM) is one of the most abundant protein domains in eukaryotes. While the structure of this domain is well characterized by the packing of two alpha-helices on a four-stranded beta-sheet, the mode of protein and RNA recognition by RRMs is not clear owing to the high variability of these interactions. Here we report recent structural data on RRM-RNA and RRM-protein interactions showing the ability of this domain to modulate its binding affinity and specificity using each of its constitutive elements (beta-strands, loops, alpha-helices). The extreme structural versatility of the RRM interactions explains why RRM-containing proteins have so diverse biological functions.     

Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.     

Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses

For a very long time, Type II restriction enzymes (REases) have been a paradigm of ORFans: proteins with no detectable similarity to each other and to any other protein in the database, despite common cellular and biochemical function. Crystallographic analyses published until January 2008 provided high-resolution structures for only 28 of 1637 Type II REase sequences available in the Restriction Enzyme database (REBASE). Among these structures, all but two possess catalytic domains with the common PD-(D/E)XK nuclease fold. Two structures are unrelated to the others: R.BfiI exhibits the phospholipase D (PLD) fold, while R.PabI has a new fold termed 'half-pipe'. Thus far, bioinformatic studies supported by site-directed mutagenesis have extended the number of tentatively assigned REase folds to five (now including also GIY-YIG and HNH folds identified earlier in homing endonucleases) and provided structural predictions for dozens of REase sequences without experimentally solved structures. Here, we present a comprehensive study of all Type II REase sequences available in REBASE together with their homologs detectable in the nonredundant and environmental samples databases at the NCBI. We present the summary and critical evaluation of structural assignments and predictions reported earlier, new classification of all REase sequences into families, domain architecture analysis and new predictions of three-dimensional folds. Among 289 experimentally characterized (not putative) Type II REases, whose apparently full-length sequences are available in REBASE, we assign 199 (69%) to contain the PD-(D/E)XK domain. The HNH domain is the second most common, with 24 (8%) members. When putative REases are taken into account, the fraction of PD-(D/E)XK and HNH folds changes to 48% and 30%, respectively. Fifty-six characterized (and 521 predicted) REases remain unassigned to any of the five REase folds identified so far, and may exhibit new architectures. These enzymes are proposed as the most interesting targets for structure determination by high-resolution experimental methods. Our analysis provides the first comprehensive map of sequence-structure relationships among Type II REases and will help to focus the efforts of structural and functional genomics of this large and biotechnologically important class of enzymes.     

Structural trees for proteins containing φ-motifs

In the present study, a novel structural motif of proteins referred to as the phi-motif is considered, and two novel structural trees in which the phi-motif is taken as the root structure have been constructed. The simplest phi-motif is formed by three adjacent beta-strands connected by loops and packed in one beta-sheet so that its overall fold resembles the Greek letter phi. Construction of the structural trees and modeling of folding pathways have shown that all structures of the protein superfamilies can be obtained by stepwise addition of alpha-helices and/or beta-strands to the root phi-motif taking into account a restricted set of rules inferred from known principles of protein structure. The structural trees are a good tool for structure comparison, structural classification of proteins, as well as for searching for all possible protein folds and folding pathways.     

Solution structure of hypothetical Nudix hydrolase DR0079 from extremely radiation-resistant Deinococcus radiodurans bacterium

Using nuclear magnetic resonance (NMR) based methods, including residual dipolar coupling restraints, we have determined the solution structure of the hypothetical Deinococcus radiodurans Nudix protein DR0079 (171 residues, MW = 19.3 kDa). The protein contains eight beta-strands and three alpha-helices organized into three subdomains: an N-terminal beta-sheet (1-34), a central Nudix core (35-140), and a C-terminal helix-turn-helix (141-171). The Nudix core and the C-terminal helix-turn-helix form the fundamental fold common to the Nudix family, a large mixed beta-sheet sandwiched between alpha-helices. The residues that compose the signature Nudix sequence, GX5EX7REUXEEXGU (where U = I, L, or V and X = any amino acid), are contained in a turn-helix-turn motif on the face of the mixed beta-sheet. Chemical shift mapping experiments suggest that DR0079 binds Mg2+. Experiments designed to determine the biological function of the protein indicate that it is not a type I isopentenyl-diphosphate delta-isomerase and that it does not bind alpha,beta-methyleneadenosine 5'-triphosphate (AMPCPP) or guanosine 5'-[beta,gamma-imido]triphosphate (GMPPNP). In this article, the structure of DR0079 is compared to other known Nudix protein structures, a potential substrate-binding surface is proposed, and its possible biological function is discussed.     

A new protein folding algorithm based on hydrophobic compactness: Rigid Unconnected Secondary Structure Iterative Assembly (RUSSIA). II: Applications

The RUSSIA procedure (Rigid Unconnected Secondary Structure Iterative Assembly) produces structural models of cores of small- and medium-sized proteins. Loops are omitted from this treatment and regular secondary structures are reduced to points, the centers of their hydrophobic faces. This methodology relies on the maximum compactness of the hydrophobic residues, as described in detail in Part I. Starting data are the sequence and the predicted limits and natures of regular secondary structures (alpha or beta). Helices are treated as rigid cylinders, whereas beta-strands are collectively taken into account within beta-sheets modeled by helicoid surfaces. Strands are allowed to shift along their mean axis to allow some flexibility and the alpha-helices can be placed on either side of beta-sheets. Numerous initial conformations are produced by discrete rotations of the helices and sheets around the direction going from the center of their hydrophobic face to the global center of the protein. Selection of proposed models is based upon a criterion lying on the minimization of distances separating hydrophobic residues belonging to different regular secondary structures. The procedure is rapid and appears to be robust relative to the quality of starting data (nature and length of regular secondary structures). This dependence of the quality of the model on secondary structure prediction and in particular the beta-sheet topology, is one of the limits of the present algorithm. We present here some results for a set of 12 proteins (alpha, beta and alpha/beta classes) of lengths 40-166 amino acids. The r.m.s. deviations for core models with respect to the native proteins are in the range 1.4-3.7 A.     

Structure of pvu II DNA-(cytosine N4) methyltransferase, an example of domain permutation and protein fold assignment.

We have determined the structure of Pvu II methyltransferase (M. Pvu II) complexed with S -adenosyl-L-methionine (AdoMet) by multiwavelength anomalous diffraction, using a crystal of the selenomethionine-substituted protein. M. Pvu II catalyzes transfer of the methyl group from AdoMet to the exocyclic amino (N4) nitrogen of the central cytosine in its recognition sequence 5'-CAGCTG-3'. The protein is dominated by an open alpha/beta-sheet structure with a prominent V-shaped cleft: AdoMet and catalytic amino acids are located at the bottom of this cleft. The size and the basic nature of the cleft are consistent with duplex DNA binding. The target (methylatable) cytosine, if flipped out of the double helical DNA as seen for DNA methyltransferases that generate 5-methylcytosine, would fit into the concave active site next to the AdoMet. This M. Pvu IIalpha/beta-sheet structure is very similar to those of M. Hha I (a cytosine C5 methyltransferase) and M. Taq I (an adenine N6 methyltransferase), consistent with a model predicting that DNA methyltransferases share a common structural fold while having the major functional regions permuted into three distinct linear orders. The main feature of the common fold is a seven-stranded beta-sheet (6 7 5 4 1 2 3) formed by five parallel beta-strands and an antiparallel beta-hairpin. The beta-sheet is flanked by six parallel alpha-helices, three on each side. The AdoMet binding site is located at the C-terminal ends of strands beta1 and beta2 and the active site is at the C-terminal ends of strands beta4 and beta5 and the N-terminal end of strand beta7. The AdoMet-protein interactions are almost identical among M. Pvu II, M. Hha I and M. Taq I, as well as in an RNA methyltransferase and at least one small molecule methyltransferase. The structural similarity among the active sites of M. Pvu II, M. Taq I and M. Hha I reveals that catalytic amino acids essential for cytosine N4 and adenine N6 methylation coincide spatially with those for cytosine C5 methylation, suggesting a mechanism for amino methylation.     

The effect of long-range interactions on the secondary structure formation of proteins

The influence of long-range residue interactions on defining secondary structure in a protein has long been discussed and is often cited as the current limitation to accurate secondary structure prediction. There are several experimental examples where a local sequence alone is not sufficient to determine its secondary structure, but a comprehensive survey on a large data set has not yet been done. Interestingly, some earlier studies denied the negative effect of long-range interactions on secondary structure prediction accuracy. Here, we have introduced the residue contact order (RCO), which directly indicates the separation of contacting residues in terms of the position in the sequence, and examined the relationship between the RCO and the prediction accuracy. A large data set of 2777 nonhomologous proteins was used in our analysis. Unlike previous studies, we do find that prediction accuracy drops as residues have contacts with more distant residues. Moreover, this negative correlation between the RCO and the prediction accuracy was found not only for beta-strands, but also for alpha-helices. The prediction accuracy of beta-strands is lower if residues have a high RCO or a low RCO, which corresponds to the situation that a beta-sheet is formed by beta-strands from different chains in a protein complex. The reason why the current study draws the opposite conclusion from the previous studies is examined. The implication for protein folding is also discussed.     

Is tensegrity a unifying concept of protein folds?

We suggest that the three-dimensional architecture of globular proteins can be described in terms of tensegrets, i.e. structural elements that are held together through attractive and repulsive forces. Hard elements of tensegrets are represented by secondary structure elements, i.e. alpha-helices and beta-strands, while the role of elastic elements is played by attractive and repulsive atomic forces. Characteristics of tensegrets is that they can auto-assemble and that they respond to changes of tension in some part of the entire object through a deformation in another part, thus partially preserving their structure, despite their deformation. This latter property well explains both the folding process and the behavior of globular proteins under mild denaturing conditions, as revealed by the molten globule state.     

An algorithm for predicting protein-protein interaction sites: Abnormally exposed amino acid residues and secondary structure elements

Multiprotein systems mediate most regulatory processes in living organisms. Although the structures of the individual proteins are often defined, less is known of the structures of multiprotein systems. Computational methods for predicting interfaces, using evolutionary conservation and/or physicochemical data, have been developed. Here we consider the use of solvent accessibility, residue propensity, and hydrophobicity, in conjunction with secondary structure data, as prediction parameters. We analyze the influence of residue type and secondary structure on solvent accessibility and define a measure of "relative exposedness." Clustering abnormally high scoring residues provides a basis for predicting interaction sites. The analysis is extended to investigate abnormally exposed secondary structure elements, particularly beta-sheet strands. We show that surface-exposed beta-strands lacking protective features are more likely to be found at protein-protein interfaces, allowing us to create an algorithm with approximately 68% and approximately 75% accuracy in differentiating between interacting and edge strands in isolated beta-strands and beta-sheet strands, respectively. These methods of identifying abnormally exposed surface regions are combined in an algorithm, which, on a data set of 77 unbound and disjoint (single chain extracted from complex) structures, predicts 79% of the protein-protein interfaces correctly. If enzyme-inhibitor complexes, where the inhibitor mimics a nonprotein substrate, are excluded, the accuracy increases to 85%.     

Novel protein structural motifs containing two-turn and longer 3(10)-helices.

The 3(10)-helix constitutes a small but significant fraction of secondary structural elements in proteins. Protein data base surveys have shown these helices to be present as alpha-helical extensions, in loops and as connectors between beta-strands. The present work focuses on two-turn and longer 3(10)-helices where we establish that two-turn and longer 3(10) helices, unlike the more abundant single-turn 3(10)-helices, frequently occur independent of any other contiguous secondary structural elements. More importantly, a large fraction of these independent two-turn and longer 3(10)-helices, along with alpha-helices and beta-strands, are found to form novel super-secondary structural motifs in several proteins with possible implications for protein folding, local conformational relaxation and biological functions.     

Type I restriction enzymes and their relatives

Type I restriction enzymes (REases) are large pentameric proteins with separate restriction (R), methylation (M) and DNA sequence-recognition (S) subunits. They were the first REases to be discovered and purified, but unlike the enormously useful Type II REases, they have yet to find a place in the enzymatic toolbox of molecular biologists. Type I enzymes have been difficult to characterize, but this is changing as genome analysis reveals their genes, and methylome analysis reveals their recognition sequences. Several Type I REases have been studied in detail and what has been learned about them invites greater attention. In this article, we discuss aspects of the biochemistry, biology and regulation of Type I REases, and of the mechanisms that bacteriophages and plasmids have evolved to evade them. Type I REases have a remarkable ability to change sequence specificity by domain shuffling and rearrangements. We summarize the classic experiments and observations that led to this discovery, and we discuss how this ability depends on the modular organizations of the enzymes and of their S subunits. Finally, we describe examples of Type II restriction–modification systems that have features in common with Type I enzymes, with emphasis on the varied Type IIG enzymes.     

Structural study of the zinc and cadmium complexes of a type 2 plant (Quercus suber) metallothionein: insights by vibrational spectroscopy

Zn- and Cd-complexes of Quercus suber metallothionein (QsMT) were obtained by in vivo-synthesis, in order to obtain physiologically representative aggregates, and characterized by spectrometric and spectroscopic methods. The secondary structure elements and the coordination environments of the metal binding sites of the two aggregates were determined, as well as the main metal-containing species formed. The results obtained from the analysis of the Raman and IR spectra reveal that these metal-MT complexes predominantly contain beta-sheet elements (about 60%), whereas they lack alpha-helices. These structural features slightly depend on the divalent metal bound. In particular, Cd(II) binding to QsMT induces a slight increase of the beta-sheet percentage, as well as a decrease in beta-turn elements with respect to Zn(II) binding. Conversely, the in vivo capability of QsMT to inglobe metal and sulfide ions is metal-depending. Spectroscopic vibrational data also confirm the presence of sulfide ligands in the metal clusters of both Zn- and Cd-QsMT, while the participation of the spacer His residue in metal coordination was only found in Cd-QsMT, in agreement with the CD results. Overall data suggest different coordination environments for Zn(II) and Cd(II) ions in QsMT.     

Type II restriction endonucleases—a historical perspective and more

This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss ‘Type II’ REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.     

Structure- and function-based characterization of a new phosphoglycolate phosphatase from Thermoplasma acidophilum

The protein TA0175 has a large number of sequence homologues, most of which are annotated as unknown and a few as belonging to the haloacid dehalogenase superfamily, but has no known biological function. Using a combination of amino acid sequence analysis, three-dimensional crystal structure information, and kinetic analysis, we have characterized TA0175 as phosphoglycolate phosphatase from Thermoplasma acidophilum. The crystal structure of TA0175 revealed two distinct domains, a larger core domain and a smaller cap domain. The large domain is composed of a centrally located five-stranded parallel beta-sheet with strand order S10, S9, S8, S1, S2 and a small beta-hairpin, strands S3 and S4. This central sheet is flanked by a set of three alpha-helices on one side and two helices on the other. The smaller domain is composed of an open faced beta-sandwich represented by three antiparallel beta-strands, S5, S6, and S7, flanked by two oppositely oriented alpha-helices, H3 and H4. The topology of the large domain is conserved; however, structural variation is observed in the smaller domain among the different functional classes of the haloacid dehalogenase superfamily. Enzymatic assays on TA0175 revealed that this enzyme catalyzed the dephosphorylation of phosphoglycolate in vitro with similar kinetic properties seen for eukaryotic phosphoglycolate phosphatase. Activation by divalent cations, especially Mg2+, and competitive inhibition behavior with Cl- ions are similar between TA0175 and phosphoglycolate phosphatase. The experimental evidence presented for TA0175 is indicative of phosphoglycolate phosphatase.