Computer-assisted dissection of rolling circle DNA replication

A comparative analysis of the proteins involved in initiation and termination of rolling circle replication (RCR) was performed using computer-assisted methods of data based screening, motif search and multiple amino acid sequence alignment. Two vast classes of such proteins were delineated, one of these being associated with RCR proper, and the other with mobilization (conjugal transfer) of plasmid DNA. The common denominator of the two classes was found to be a conserved amino acid motif that consists of the sequence HisUHisUUU (U--bulky hydrophobic residue; hereafter HUH motif). Based on analogies with metalloenzymes, it is hypothesized that the two conserved His residues this motif may be involved in metal ion coordination required for the activity of the RCR and mobilization proteins. The proteins of the replication (Rep) class contained two additional conserved motifs, with the motif around the Tyr residue(s) forming the covalent link with nicked DNA being located C-proximally of the HUH motif. This class further split into two large superfamilies and several smaller families, with the proteins belonging to a single but not to different (super)families demonstrating statistically significant similarity to each other. Superfamily I, prototyped by the gene A proteins of small isometric single-stranded (ss) DNA bacteriophages, included also Rep proteins of P2-related double-stranded (ds) DNA bacteriophages, the small phage-plasmid hybrid phasyl, and several cyanobacterial and archaebacterial plasmids. These proteins contained two invariant Tyr residues separated by three partially conserved amino acids, suggesting that they all may share the cleavage-ligation mechanism proposed for phi X174 A protein and involving alternate covalent binding of both tyrosines to DNA (Van Mansfeld, A.D., Van Teeffelen, H.A., Baas, P.D., Jansz, H.S., 1986. Nucl. Acids Res. 14, 4229-4238). Superfamily II included Rep proteins of a number of ssDNA plasmids replicating mainly in gram-positive bacteria that unexpectedly were shown to be related to the Rep proteins of plant geminiviruses. Conservation of the "HUH" motif and a motif around the putative DNA-linking Tyr residue was observed also in the Rep proteins of animal parvoviruses containing linear ssDNA with a terminal hairpin and replicating via the rolling hairpin mechanism. The class of plasmid mobilization (Mob) proteins was characterized by the opposite orientation of the conserved motifs, with the (putative) DNA-linking Tyr being located N-proximally of the "HUH" motif.(ABSTRACT TRUNCATED AT 400 WORDS)     

Central nervous system infection due to Herpes simplex virus in AIDS

The secondary structure and global fold of the AVR9 elicitor protein of Cladosporium fulvum has been determined by 2D NMR and distance-geometry protocols. The protein consists of three anti-parallel strands forming a rigid region of beta-sheet. On the basis of the NMR-derived parameters and distance geometry calculations, it is evident that the AVR9 protein is structurally very homologuous to carboxy peptidase inhibitor (CPI) of which the X-ray structure is known. The AVR9 protein reveals the presence of a cystine knot, which consists of a ring formed by two disulfide bridges and the interconnecting backbone through which the third disulfide bridge penetrates. This structural motif is found in several small proteins such as proteinase inhibitors, ion channel blockers and growth factors. The implications of the structural relationship between AVR9 and other biologically active proteins are discussed.     

Solution structure of human intestinal fatty acid binding protein: implications for ligand entry and exit

The human intestinal fatty acid binding protein (I-FABP) is a small (131 amino acids) protein which binds dietary long-chain fatty acids in the cytosol of enterocytes. Recently, an alanine to threonine substitution at position 54 in I-FABP has been identified which affects fatty acid binding and transport, and is associated with the development of insulin resistance in several populations including Mexican-Americans and Pima Indians. To investigate the molecular basis of the binding properties of I-FABP, the 3D solution structure of the more common form of human I-FABP (Ala54) was studied by multidimensional NMR spectroscopy. Recombinant I-FABP was expressed from E. coli in the presence and absence of 15N-enriched media. The sequential assignments for non-delipidated I-FABP were completed by using 2D homonuclear spectra (COSY, TOCSY and NOESY) and 3D heteronuclear spectra (NOESY-HMQC and TOCSY-HMQC). The tertiary structure of human I-FABP was calculated by using the distance geometry program DIANA based on 2519 distance constraints obtained from the NMR data. Subsequent energy minimization was carried out by using the program SYBYL in the presence of distance constraints. The conformation of human I-FABP consists of 10 antiparallel beta-strands which form two nearly orthogonal beta-sheets of five strands each, and two short alpha-helices that connect the beta-strands A and B. The interior of the protein consists of a water-filled cavity between the two beta-sheets. The NMR solution structure of human I-FABP is similar to the crystal structure of rat I-FABP. The NMR results show significant conformational variability of certain backbone segments around the postulated portal region for the entry and exit of fatty acid ligand.     

Recognition of spatial motifs in protein structures

As the structural database continues to expand, new methods are required to analyse and compare protein structures. Whereas the recognition, comparison, and classification of folds is now more or less a solved problem, tools for the study of constellations of small numbers of residues are few and far between. In this paper, two programs are described for the analysis of spatial motifs in protein structures. The first, SPASM, can be used to find the occurrence of a motif consisting of arbitrary main-chain and/or side-chains in a database of protein structures. The program also has a unique capability to carry out "fuzzy pattern matching" with relaxed requirements on the types of some or all of the matching residues. The second program, RIGOR, scans a single protein structure for the occurrence of any of a set of pre-defined motifs from a database. In one application, spatial motif recognition combined with profile analysis enabled the assignment of the structural and functional class of an uncharacterised hypothetical protein in the sequence database. In another application, the occurrence of short left-handed helical segments in protein structures was investigated, and such segments were found to be fairly common. Potential applications of the techniques presented here lie in the analysis of (newly determined) structures, in comparative structural analysis, in the design and engineering of novel functional sites, and in the prediction of structure and function of uncharacterised proteins.     

Insect prothoracicotropic hormone: a new member of the vertebrate growth factor superfamily

Prothoracicotropic hormone (PTTH) is a brain neurosecretory protein that controls insect development. PTTH of the silkmoth Bombyx mori is a homodimeric protein, the subunit of which consists of 109 amino acids. Clear-cut sequence similarity to any other proteins has not been observed. By disulfide-bond pattern analysis and modeling of the PTTH structure based on the known three-dimensional (3D) structures of growth factor family with cystine-knot motif, we propose that the PTTH protomer adopts the fold unique to the structural superfamily of the growth factors, beta-nerve growth factor (beta-NGF), transforming growth factor-beta 2 (TGF-beta 2), and platelet-derived growth factor-BB (PDGF-BB). The insect neurohormone PTTH appears to be a member of the growth factor superfamily, sharing a common ancestral gene with the three vertebrate growth factors, beta-NGF, TGF-beta 2 and PDGF-BB.     

A graph-theoretic algorithm for comparative modeling of protein structure

The interconnected nature of interactions in protein structures appears to be the major hurdle in preventing the construction of accurate comparative models. We present an algorithm that uses graph theory to handle this problem. Each possible conformation of a residue in an amino acid sequence is represented using the notion of a node in a graph. Each node is given a weight based on the degree of the interaction between its side-chain atoms and the local main-chain atoms. Edges are then drawn between pairs of residue conformations/nodes that are consistent with each other (i.e. clash-free and satisfying geometrical constraints). The edges are weighted based on the interactions between the atoms of the two nodes. Once the entire graph is constructed, all the maximal sets of completely connected nodes (cliques) are found using a clique-finding algorithm. The cliques with the best weights represent the optimal combinations of the various main-chain and side-chain possibilities, taking the respective environments into account. The algorithm is used in a comparative modeling scenario to build side-chains, regions of main chain, and mix and match between different homologs in a context-sensitive manner. The predictive power of this method is assessed by applying it to cases where the experimental structure is not known in advance.     

Molecular shape comparisons in searches for active sites and functional similarity

Here we examine the reliability of surface comparisons in searches for active sites in proteins. Detection of a patch of surface on one protein which is similar to an active site in another, may suggest similarities in enzymatic mechanisms, in enzyme functions and implicate a potential target for ligand/inhibitor design. Specifically, we compare the efficacy of molecular surface comparisons with comparisons of surface atoms and of C(alpha) backbone atoms. We further investigate comparisons of specific atoms, belonging to a predefined pattern of catalytic residues versus comparisons of molecular surfaces and, separately, of surface atoms. This aspect is particularly relevant, as catalytic residues may be (partially) buried. We also explore active site comparisons versus comparisons in which the entire molecular surfaces are scanned. While here we focus on the geometrical aspect of the problem, we also investigate the effect of adding residue labels in these comparisons. Our extensive studies cover the serine proteases, containing the highly conserved triad motif, and the chorismate mutases. Since such active site comparisons entail comparisons between unconnected points in 3D space, an order-independent comparison technique is necessary. The geometric hashing algorithm is ideally suited to handling such a task. It can perform both global shape matching for the whole surfaces of large protein molecules and searching for local shape similarities for small surface motifs. Our results show that molecular surface comparisons work best when the similarity is high. As the similarity deteriorates, the number of potential solutions increases rapidly, making their ranking difficult, particularly when scanning entire molecular surfaces. Utilizing atomic coordinates directly appears more adequate under such circumstances.     

Structure and distribution of pentapeptide repeats in bacteria

We report the discovery of a novel family of proteins, each member contains tandem pentapeptide (five residue) repeats, described by the motif A(D/N)LXX. Members of this family are both membrane bound and cytoplasmic. The function of these repeats is uncertain, but they may have a targeting or structural function rather than enzymatic activity. This family is most common in cyanobacteria, suggesting a function related to cyanobacterial-specific metabolism. Although no experimental information is available for the structure of this family, it is predicted that the tandem pentapeptide repeats will form a right-handed beta-helical structure. A structural model of the pentapeptide repeats is presented.     

Oscillin, an extracellular, Ca2+-binding glycoprotein essential for the gliding motility of cyanobacteria

Electron microscopic studies have demonstrated that various gliding filamentous cyanobacteria have trichome surfaces with a common structural organization. They contain an S-layer attached to the outer membrane and an array of parallel fibrils on top of the S-layer. In all species studied, the helical arrangement of these fibrils corresponds to the sense of rotation of the organism during the gliding movement. We have investigated the surface fibrils of Phormidium uncinatum using electron microscopic, spectroscopic and biochemical techniques. The fibrils consist of a single rod-shaped protein, which we refer to as oscillin. Oscillin is a 646 amino acid residue protein (Mr 65807; pI 3.63) and appears to be glycosylated. Sequence analysis reveals a two-domain structure: a 554 residue domain contains 46 repeats of a Ca2+-binding motif; it is followed by a 92 residue C-terminal domain, which might mediate its export. Filaments that do not express oscillin lose their ability to move. Homology studies suggest that similar proteins play comparable roles in other motile cyanobacteria. The structure of oscillin appears to favour a passive role in gliding.     

Mutational analysis of the hepatitis C virus RNA helicase

The carboxyl-terminal three-fourths of the hepatitis C virus (HCV) NS3 protein has been shown to possess an RNA helicase activity, typical of members of the DEAD box family of RNA helicases. In addition, the NS3 protein contains four amino acid motifs conserved in DEAD box proteins. In order to inspect the roles of individual amino acid residues in the four conserved motifs (AXXXXGKS, DECH, TAT, and QRRGRTGR) of the NS3 protein, mutational analysis was used in this study. Thirteen mutant proteins were constructed, and their biochemical activities were examined. Lys1235 in the AXXXXGKS motif was important for basal nucleoside triphosphatase (NTPase) activity in the absence of polynucleotide cofactor. A serine in the X position of the DEXH motif disrupted the NTPase and RNA helicase activities. Alanine substitution at His1318 of the DEXH motif made the protein possess high NTPase activity. In addition, we now report inhibition of NTPase activity of NS3 by polynucleotide cofactor. Gln1486 was indispensable for the enzyme activity, and this residue represents a distinguishing feature between DEAD box and DEXH proteins. There are four Arg residues in the QRRGRTGR motif of the HCV NS3 protein, and the second, Arg1488, was important for RNA binding and enzyme activity, even though it is less well conserved than other Arg residues. Arg1490 and Arg1493 were essential for the enzymatic activity. As the various enzymatic activities were altered by mutation, the enzyme characteristics were also changed.     

Evidence supporting a role for histidine-235 in cation binding to human 3-hydroxy-3-methyglutaryl-CoA lyase

Histidine-235 of human 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) lyase is the second basic residue in a conserved HXH motif. This residue is solvent accessible, readily reacting with the group specific reagent diethyl pyrocarbonate. Site-directed mutagenesis has been employed to substitute alanine or aspartate for H235. Characterization of the isolated H235A and H235D lyase mutants indicates that their tertiary structure is substantially intact. The mutant proteins, like the wild-type enzyme, are stoichiometrically modified by the affinity label, 2-butynoyl-CoA. Catalytic activity of the mutants is diminished by 15-fold and Km for HMG-CoA elevated approximately 4-fold in comparison with the values for wild-type enzyme. The function of H235 is suggested by investigation of the interaction of these enzymes with the dissociable divalent cation (e.g. Mg2+ or Mn2+) that is required for activity. ESR experiments show that wild-type enzyme forms a stable binary E*M complex. In contrast, H235A and H235D proteins do not efficiently form a binary complex. Significant interaction with cation (Mn2+) only occurs in the presence of the substrate analog, 3-hydroxyglutaryl-CoA. Similarly, when cation interaction is estimated in the presence of substrate using steady-state kinetic approaches, activator constants (Ka) and divalent cation Km values are measurable but are elevated by 15-90-fold over comparable estimates for the wild-type enzyme. The data confirm our earlier suggestion that both protein and substrate contribute ligands to HMG-CoA lyase's divalent cation activator. More specifically, the current observations suggest that H235 has an important function in cation binding.     

In vitro assembly of the mouse U14 snoRNP core complex and identification of a 65-kDa box C/D-binding protein

The eukaryotic nucleolus contains a diverse population of small nucleolar RNAs (snoRNAs) that have been categorized into two major families based on evolutionarily conserved sequence elements. U14 snoRNA is a member of the larger, box C/D snoRNA family and possesses nucleotide box C and D consensus sequences. In previous studies, we have defined a U14 box C/D core motif that is essential for intronic U14 snoRNA processing. These studies also revealed that nuclear proteins that recognize boxes C/D are required. We have now established an in vitro U14 snoRNP assembly system to characterize protein binding. Electrophoretic mobility-shift analysis demonstrated that all the sequences and structures of the box C/D core motif required for U14 processing are also necessary for protein binding and snoRNP assembly. These required elements include a base paired 5',3' terminal stem and the phylogenetically conserved nucleotides of boxes C and D. The ability of other box C/D snoRNAs to compete for protein binding demonstrated that the box C/D core motif-binding proteins are common to this family of snoRNAs. UV crosslinking of nuclear proteins bound to the U14 core motif identified a 65-kDa mouse snoRNP protein that requires boxes C and D for binding. Two additional core motif proteins of 55 and 50 kDa were also identified by biochemical fractionation of the in vitro-assembled U14 snoRNP complex. Thus, the U14 snoRNP core complex is a multiprotein particle whose assembly requires nucleotide boxes C and D.     

Envelope formation is blocked by mutation of a sequence related to the HKD phospholipid metabolism motif in the vaccinia virus F13L protein

The outer envelope of the extracellular form of vaccinia virus is derived from Golgi membranes that have been modified by the insertion of specific viral proteins, of which the major component is the 37-kDa, palmitylated, nonglycosylated product of the F13L gene. The F13L protein contains a variant of the HKD (His-Lys-Asp) motif, which is conserved in numerous enzymes of phospholipid metabolism. Vaccinia virus mutants with a conservative substitution of either the K (K314R) or the D (D319E) residue of the F13L protein formed only tiny plaques similar to those produced by an F13L deletion mutant, were unable to produce extracellular enveloped virions, and failed to mediate low-pH-induced fusion of infected cells. Membrane-wrapped forms of intracellular virus were rarely detected in electron microscopic images of cells infected with either of the mutants. Western blotting and pulse-chase experiments demonstrated that the D319E protein was less stable than either the K314R or wild-type F13L protein. Most striking, however, was the failure of either of the two mutated proteins to concentrate in the Golgi compartment. Palmitylation, oleation, and partitioning of the F13L protein in Triton X-114 detergent were unaffected by the K314R substitution. These results indicated that the F13L protein must retain the K314 and D319 for it to localize in the Golgi compartment and function in membrane envelopment of vaccinia virus.     

TESS: a geometric hashing algorithm for deriving 3D coordinate templates for searching structural databases. Application to enzyme active sites

It is well established that sequence templates such as those in the PROSITE and PRINTS databases are powerful tools for predicting the biological function and tertiary structure for newly derived protein sequences. The number of X-ray and NMR protein structures is increasing rapidly and it is apparent that a 3D equivalent of the sequence templates is needed. Here, we describe an algorithm called TESS that automatically derives 3D templates from structures deposited in the Brookhaven Protein Data Bank. While a new sequence can be searched for sequence patterns, a new structure can be scanned against these 3D templates to identify functional sites. As examples, 3D templates are derived for enzymes with an O-His-O "catalytic triad" and for the ribonucleases and lysozymes. When these 3D templates are applied to a large data set of nonidentical proteins, several interesting hits are located. This suggests that the development of a 3D template database may help to identify the function of new protein structures, if unknown, as well as to design proteins with specific functions.     

Uncoupling hydrophobicity and helicity in transmembrane segments. Alpha-helical propensities of the amino acids in non-polar environments

Although the chains of amino acids in proteins that span the membrane are demonstrably helical and hydrophobic, little attention has been paid toward addressing the range of helical propensities of individual amino acids in the non-polar environment of membranes. Because it is inappropriate to apply soluble protein-based structure prediction algorithms to membrane proteins, we have used de novo designed peptides (KKAAAXAAAAAXAAWAAXAAAKKKK-amide, where X indicates one of the 20 commonly occurring amino acids) that mimic a protein membrane-spanning domain to determine the alpha-helical proclivity of each residue in the isotropic non-polar environment of n-butanol. Peptide helicities measured by circular dichroism spectroscopy were found to range from theta222 = -17,000 degrees (Pro) to -38,800 degrees (Ile) in n-butanol. The relative helicity of each amino acid is shown to be well correlated with its occurrence frequency in natural transmembrane segments, indicating that the helical propensity of individual residues in concert with their hydrophobicity may be a key determinant of the conformations of protein segments in membranes.     

Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker's yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.     

Mutations in the human immunodeficiency virus type 1 integrase D,D(35)E motif do not eliminate provirus formation

The core domain of human immunodeficiency virus type 1 (HIV-1) integrase (IN) contains a D,D(35)E motif, named for the phylogenetically conserved glutamic acid and aspartic acid residues and the invariant 35 amino acid spacing between the second and third acidic residues. Each acidic residue of the D,D(35)E motif is independently essential for the 3'-processing and strand transfer activities of purified HIV-1 IN protein. Using a replication-defective viral genome with a hygromycin selectable marker, we recently reported that a mutation at any of the three residues of the D,D(35)E motif produces a 10(3)- to 10(4)-fold reduction in infectious titer compared with virus encoding wild-type IN (A. D. Leavitt et al., J. Virol. 70:721-728. 1996). The infectious titer, as measured by the number of hygromycin-resistant colonies formed following infection of cells in culture, was less than a few hundred colonies per microg of p24. To understand the mechanism by which the mutant virions conferred hygromycin resistance, we characterized the integrated viral DNA in cells infected with virus encoding mutations at each of the three residues of the D,D(35)E motif. We found the integrated viral DNA to be colinear with the incoming viral genome. DNA sequencing of the junctions between integrated viral DNA and host DNA showed that (i) the characteristic 5-bp direct repeat of host DNA flanking the HIV-1 provirus was not maintained, (ii) integration often produced a deletion of host DNA, (iii) integration sometimes occurred without the viral DNA first undergoing 3'-processing, (iv) integration sites showed a strong bias for a G residue immediately adjacent to the conserved viral CA dinucleotide, and (v) mutations at each of the residues of the D,D(35)E motif produced essentially identical phenotypes. We conclude that mutations at any of the three acidic residues of the conserved D,D(35)E motif so severely impair IN activity that most, if not all, integration events by virus encoding such mutations are not IN mediated. IN-independent provirus formation may have implications for anti-IN therapeutic agents that target the IN active site.     

AT-hook motifs identified in a wide variety of DNA-binding proteins

The AT-hook is a small DNA-binding protein motif which was first described in the high mobility group non-histone chromosomal protein HMG-I(Y). Since its discovery, this motif has been observed in other DNA-binding proteins from a wide range of organisms. Using pattern searches and position-dependent matrices, we have extracted the AT-hook motifs present in a non-redundant protein sequence database. We have classified these motifs into three types according to their sequence similarity and have found that they are prevalent in many eukaryotic nuclear proteins in single or multiple copies. Furthermore, AT-hook motifs are frequently associated with known functional domains seen in chromatin proteins and in DNA-binding proteins (e.g. histone folds, homeodomains and zinc fingers). In general, it appears that the AT-hook motif is an auxiliary protein motif cooperating with other DNA-binding activities and facilitating changes in the structure of the DNA either as a polypeptide on its own [e.g. HMG-I(Y)] or as part of a multidomain protein [e.g. Swi2p in Saccharomyces cerevisiae or HRX (ALL-1) in Homo sapiens]. It is most interesting that this motif seems to be quite specific to known or predicted chromosomal/DNA-binding proteins, suggesting that it may act as a versatile minor groove tether.     

Conformational and dynamic changes of Yersinia protein tyrosine phosphatase induced by ligand binding and active site mutation and revealed by H/D exchange and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry

Protein tyrosine phosphatases (PTPase) play important roles in the intracellular signal transduction pathways that regulate cell transformation, growth, and proliferation. Here, solvent accessibility is determined for backbone amide protons from various segments of wild-type Yersinia PTPase in the presence or absence of 220 microM vanadate, a competitive inhibitor, as well as an active site mutant in which the essential cysteine 403 has been replaced by serine (C403S). The method consists of solution-phase H/D exchange, followed by pepsin digestion, high-performance liquid chromatography, and electrospray ionization high-field (9.4 T) Fourier transform ion cyclotron resonance mass spectrometry. Proteolytic segments spanning approximately 93.5% of the primary sequence are analyzed. Binding of vanadate reduces the H/D exchange rate throughout the protein, both for the WpD loop and for numerous other residues that are shielded when that loop is pulled down over the active site on binding of the inhibitor. The single active site C403S mutation reduces solvent access to the WpD loop itself, but opens up the structure in several other segments. Although the 3D structure of the ligand-bound C403S mutant is similar to that of the wild-type PTPase, and the C403S mutant and the wild-type enzyme display similar affinities for vanadate, the thermodynamics for binding of vanadate is different for the two proteins. Collectively, these results establish the flexibility of the WpD loop (previously inferred by comparing PTPase X-ray single-cyrstal diffraction structures in the presence and absence of a tungstate inhibitor), as well as several other signficant changes in segment exposure and/or flexibility that are not evident from X-ray structures.