Characterization of molecular alignment in aqueous suspensions of Pf1 bacteriophage

The phase diagram of Pf1 solutions has been studied indirectly by observation of ²H quadrupole splittings of the solvent signal and measurement of dipolar couplings in solute macromolecules. At low volume fractions of Pf1 and at high ionic strength, alignment of both the phage and the solute depends strongly on the strength of the magnetic field. Both the theoretical and experimentally determined phase diagram of Pf1 show that at low concentrations and high ionic strengths the solution becomes isotropic. However, just below the nematic phase boundary the behavior of the system is paranematic, with cooperative alignment which depends on the strength of the applied magnetic field. Above 16 mg/ml Pf1 is fully nematic up to 600 mM NaCl. Alignment of proteins with a significant electric dipole moment, which tends to be strong in Pf1, can be reduced by either high ionic strength or low phage concentration. Because ionic strength modulates both the orientation and magnitude of the alignment tensor in Pf1 medium, measurement at two ionic strengths can yield linearly independent alignment tensors.     

Ribonucleotide reductases: divergent evolution of an ancient enzyme

Ribonucleotide reductases (RNRs) are uniquely responsible for converting nucleotides to deoxynucleotides in all dividing cells. The three known classes of RNRs operate through a free radical mechanism but differ in the way in which the protein radical is generated. Class I enzymes depend on oxygen for radical generation, class II uses adenosylcobalamin, and the anaerobic class III requires S-adenosylmethionine and an iron–sulfur cluster. Despite their metabolic prominence, the evolutionary origin and relationships between these enzymes remain elusive. This gap in RNR knowledge can, to a major extent, be attributed to the fact that different RNR classes exhibit greatly diverged polypeptide chains, rendering homology assessments inconclusive. Evolutionary studies of RNRs conducted until now have focused on comparison of the amino acid sequence of the proteins, without considering how they fold into space. The present study is an attempt to understand the evolutionary history of RNRs taking into account their three-dimensional structure. We first infer the structural alignment by superposing the equivalent stretches of the three-dimensional structures of representatives of each family. We then use the structural alignment to guide the alignment of all publicly available RNR sequences. Our results support the hypothesis that the three RNR classes diverged from a common ancestor currently represented by the anaerobic class III. Also, lateral transfer appears to have played a significant role in the evolution of this protein family.     

Solution structure determination of the two DNA-binding domains in the Schizosaccharomyces pombe Abp1 protein by a combination of dipolar coupling and diffusion anisotropy restraints

We have solved the solution structure of the N-terminal region of the fission yeast centromere protein, Abp1, bound to a 21-base pair DNA fragment bearing its recognition site (Mw = 30 kDa). Although the two DNA-binding domains in the Abp1 protein were defined well by a conventional NOE-based NMR methodology, the overall structure of the Abp1 protein was poorly defined, due to the lack of interdomain distance restraints. Therefore, we additionally used residual dipolar couplings measured in a weakly aligned state, and rotational diffusion anisotropies. Neither the NH residual dipolar couplings nor the backbone ¹⁵N T ₁/T ₂ data were sufficient to determine the overall structure of the Abp1 protein, due to spectral overlap. We used a combination of these two orientational restraints (residual dipolar coupling and rotational diffusion anisotropy), which significantly improved the convergence of the overall structures. The range of the observed T ₁/T ₂ ratios was wider (20–50 for the secondary structure regions of Abp1) than the previously reported data for several globular proteins, indicating that the overall shape of the Abp1DNA complex is ellipsoid. This extended form would facilitate the recognition of the two separate sites in the relatively long DNA sequence by the DNA-binding domains of Apb1.     

Charged acrylamide copolymer gels as media for weak alignment

The use of mechanically strained acrylamide/acrylate copolymers is reported as a new alignment medium for biomacromolecules. Compared to uncharged, strained polyacrylamide gels, the negative charges of the acrylamide/acrylate copolymer strongly alter the alignment tensor and lead to pronounced electroosmotic swelling. The swelling itself can be used to achieve anisotropic, mechanical strain. The method is demonstrated for the alignment of TipAS, a 17 kDa antibiotic resistance protein, as well as for human ubiquitin, where alignment tensors with an A_ZZ,NH of up to 60 Hz are achieved at a gel concentration of 2% (w/v). The alignment can be modulated by the variation of pH, ionic strength, and gel concentration. The high mechanical stability of the swollen gels makes it possible to obtain alignment at polymer concentrations of less than 1% (w/v).     

The PRINTS database: a resource for identification of protein families

The PRINTS database houses a collection of protein fingerprints, which may be used to assign family and functional attributes to uncharacterised sequences, such as those currently emanating from the various genome-sequencing projects. The April 2002 release includes 1,700 family fingerprints, encoding approximately 10,500 motifs, covering a range of globular and membrane proteins, modular polypeptides and so on. Fingerprints are groups of conserved motifs that, taken together, provide diagnostic protein family signatures. They derive much of their potency from the biological context afforded by matching motif neighbours; this makes them at once more flexible and powerful than single-motif approaches. The technique further departs from other pattern-matching methods by readily allowing the creation of fingerprints at superfamily-, family- and subfamily-specific levels, thereby allowing more fine-grained diagnoses. Here, we provide an overview of the method of protein fingerprinting and how the results of fingerprint analyses are used to build PRINTS and its relational cousin, PRINTS-S.     

ProDom: automated clustering of homologous domains

The ProDom database is a comprehensive set of protein domain families automatically generated from the SWISS-PROT and TrEMBL sequence databases. An associated database, ProDom-CG, has been derived as a restriction of ProDom to completely sequenced genomes. The ProDom construction method is based on iterative PSI-BLAST searches and multiple alignments are generated for each domain family. The ProDom web server provides the user with a set of tools to visualise multiple alignments, phylogenetic trees and domain architectures of proteins, as well as a BLAST-based server to analyse new sequences for homologous domains. The comprehensive nature of ProDom makes it particularly useful to help sustain the growth of InterPro.     

Performance evaluation of a new algorithm for the detection of remote homologs with sequence comparison

A detailed analysis of the performance of hybrid, a new sequence alignment algorithm developed by Yu and coworkers that combines Smith Waterman local dynamic programming with a local version of the maximum-likelihood approach, was made to access the applicability of this algorithm to the detection of distant homologs by sequence comparison. We analyzed the statistics of hybrid with a set of nonhomologous protein sequences from the SCOP database and found that the statistics of the scores from hybrid algorithm follows an Extreme Value Distribution with lambda approximately 1, as previously shown by Yu et al. for the case of artificially generated sequences. Local dynamic programming was compared to the hybrid algorithm by using two different test data sets of distant homologs from the PFAM and COGs protein sequence databases. The studies were made with several score functions in current use including OPTIMA, a new score function originally developed to detect remote homologs with the Smith Waterman algorithm. We found OPTIMA to be the best score function for both both dynamic programming and the hybrid algorithms. The ability of dynamic programming to discriminate between homologs and nonhomologs in the two sets of distantly related sequences is slightly better than that of hybrid algorithm. The advantage of producing accurate score statistics with only a few simulations may overcome the small differences in performance and make this new algorithm suitable for detection of homologs in conjunction with a wide range of score functions and gap penalties.     

M13 endopeptidases: New conserved motifs correlated with structure,and simultaneous phylogenetic occurrence of PHEX and the bony fish

M13 endopeptidase alignments have focused mainly on mammalian sequences and on the active site region defining the catalytic sequence signatures. Aligning all available M13 from bacteria to human on a full-length basis, we have performed a sequence analysis. This enabled us to highlight the origin and function of the M13 PHEX subtype family endopeptidase (phosphate regulating gene with homologies to endopeptidases on the X chromosome). New evolutionary conserved regions in both prokaryotes and eukaryotes have been detected and eukaryotic-specific regions clearly delineated. Using the recently solved neprilysin structure, we have observed that all new motifs, except one, localize in the spatial vicinity of the previously reported catalytic signatures. Interestingly, a highly hydrophobic pocket containing three newly reported motifs is centered by the C-terminal tryptophan residue. Extensive M13 searches in complete and in progress higher eukaryotic genomes have lead to the identification of Danio rerio as the simplest organism having PHEX. Finally, the human PHEX substrate, the parathyroid hormone-related peptide, PTHrP(107-139), is absent in bony fish: this suggests the existence of further PHEX substrates common to both bony fishes and higher vertebrates.     

Use of conserved key amino acid positions to morph protein folds

By using three-dimensional (3D) structure alignments and a previously published method to determine Conserved Key Amino Acid Positions (CKAAPs) we propose a theoretical method to design mutations that can be used to morph the protein folds. The original Paracelsus challenge, met by several groups, called for the engineering of a stable but different structure by modifying less than 50% of the amino acid residues. We have used the sequences from the Protein Data Bank (PDB) identifiers 1ROP, and 2CRO, which were previously used in the Paracelsus challenge by those groups, and suggest mutation to CKAAPs to morph the protein fold. The total number of mutations suggested is less than 40% of the starting sequence theoretically improving the challenge results. From secondary structure prediction experiments of the proposed mutant sequence structures, we observe that each of the suggested mutant protein sequences likely folds to a different, non-native potentially stable target structure. These results are an early indicator that analyses using structure alignments leading to CKAAPs of a given structure are of value in protein engineering experiments.     

Subunit E of mitochondrial ATP synthase: a bioinformatic analysis reveals a phosphopeptide binding motif supporting a multifunctional regulatory role and identifies a related human brain protein with the same motif

The mitochondrial adenosine triphosphate (ATP) synthase is located in the inner membrane and consists of at least 16 subunit types in animals, one of which is subunit e, the function of which is not clearly defined. A highly homologous protein is located in the nucleus and named progesterone receptor binding protein (RBF), to designate its role in this organelle. In addition, the expression level of subunit e in mammalian cells fluctuates greatly and is induced by certain carcinogens and elevated in liver cancers. Because these previous observations suggested to us that subunit e may play multifunctional regulatory roles, we employed a bioinformatic approach to test this view. First, from sequence alignment studies, secondary structure analyses, and basic local alignment search tool (BLAST) searches, we concluded that mitochondrial subunit e and the homologous nuclear protein RBF are most likely the same protein. Second, we examined the known sequence and structure of one of the most common multifunctional cell regulatory proteins, the 14-3-3 protein, involved in phosphopeptide binding, and deduced that it has an apparent binding motif (-KX(6)R---RY-). Third, from careful examination of the conserved residues within all subunit e sequences in the database, we discovered that this protein has a comparable binding motif (-RY---KX(6)R-). Finally, in a BLAST search for additional homologs of subunit e, we found a human brain protein, KIAA1578, the C-terminal 30 amino acids of which are identical to those of human subunit e. This protein also contains a potential phosphopeptide binding motif. In summary, these studies provide support for the view that subunit e is a multifunctional cell regulator involved in cell signaling, and implicate the involvement of the KIAA1578 protein in cell signaling as well. These studies suggest also that, while functioning as a subunit of mitochondrial ATP synthases, subunit e may help regulate these complexes by binding to phosphopeptides within one or more of the other subunit types.