Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Designed protein tetramer zipped together with a hydrophobic Alzheimer homology: A structural clue to amyloid assembly

Limited solubility and precipitation of amyloidogenic sequences such as the Alzheimer peptide (β-AP) are major obstacles to a molecular understanding of protein fibrillation and deposition processes. Here we have circumvented the solubility problem by stepwise engineering a β-AP homology into a soluble scaffold, the monomeric protein S6. The S6 construct with the highest β-AP homology crystallizes as a tetramer that is linked by the β-AP residues forming intermolecular antiparallel β-sheets. This construct also shows increased coil aggregation during refolding, and a 14-mer peptide encompassing the engineered sequence forms fibrils. Mutational analysis shows that intermolecular association is linked to the overall hydrophobicity of the sticky sequence and implies the existence of “structural gatekeepers” in the wild-type protein, that is, charged side chains that prevent aggregation by interrupting contiguous stretches of hydrophobic residues in the primary sequence.

Rapid generation of incremental truncation libraries for protein engineering using α-phosphothioate nucleotides

Incremental truncation for the creation of hybrid enzymes (ITCHY) is a novel tool for the generation of combinatorial libraries of hybrid proteins independent of DNA sequence homology. We herein report a fundamentally different methodology for creating incremental truncation libraries using nucleotide triphosphate analogs. Central to the method is the polymerase catalyzed, low frequency, random incorporation of α-phosphothioate dNTPs into the region of DNA targeted for truncation. The resulting phosphothioate internucleotide linkages are resistant to 3′→5′ exonuclease hydrolysis, rendering the target DNA resistant to degradation in a subsequent exonuclease III treatment. From an experimental perspective the protocol reported here to create incremental truncation libraries is simpler and less time consuming than previous approaches by combining the two gene fragments in a single vector and eliminating additional purification steps. As proof of principle, an incremental truncation library of fusions between the N-terminal fragment of Escherichia coli glycinamide ribonucleotide formyltransferase (PurN) and the C-terminal fragment of human glycinamide ribonucleotide formyltransferase (hGART) was prepared and successfully tested for functional hybrids in an auxotrophic E.coli host strain. Multiple active hybrid enzymes were identified, including ones fused in regions of low sequence homology.

Topological testing of the mechanism of homology search promoted by RecA protein

To initiate homologous recombination, sequence similarity between two DNA molecules must be searched for and homology recognized. How the search for and recognition of homology occurs remains unproven. We have examined the influences of DNA topology and the polarity of RecA–single-stranded (ss)DNA filaments on the formation of synaptic complexes promoted by RecA. Using two complementary methods and various ssDNA and duplex DNA molecules as substrates, we demonstrate that topological constraints on a small circular RecA–ssDNA filament prevent it from interwinding with its duplex DNA target at the homologous region. We were unable to detect homologous pairing between a circular RecA–ssDNA filament and its relaxed or supercoiled circular duplex DNA targets. However, the formation of synaptic complexes between an invading linear RecA–ssDNA filament and covalently closed circular duplex DNAs is promoted by supercoiling of the duplex DNA. The results imply that a triplex structure formed by non-Watson–Crick hydrogen bonding is unlikely to be an intermediate in homology searching promoted by RecA. Rather, a model in which RecA-mediated homology searching requires unwinding of the duplex DNA coupled with local strand exchange is the likely mechanism. Furthermore, we show that polarity of the invading RecA–ssDNA does not affect its ability to pair and interwind with its circular target duplex DNA.

The SBASE protein domain library, release 8.0: a collection of annotated protein sequence segments

SBASE 8.0 is the eighth release of the SBASE library of protein domain sequences that contains 294 898 annotated structural, functional, ligand-binding and topogenic segments of proteins, cross-referenced to most major sequence databases and sequence pattern collections. The entries are clustered into over 2005 statistically validated domain groups (SBASE-A) and 595 non-validated groups (SBASE-B), provided with several WWW-based search and browsing facilities for online use. A domain-search facility was developed, based on non-parametric pattern recognition methods, including artificial neural networks. SBASE 8.0 is freely available by anonymous ‘ftp’ file transfer from ftp.icgeb.trieste.it. Automated searching of SBASE can be carried out with the WWW servers http://www.icgeb.trieste.it/sbase/ and http://sbase.abc.hu/sbase/.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

Mapping the early steps in the pH-induced conformational conversion of the prion protein

Under certain conditions, the prion protein (PrP) undergoes a conformational change from the normal cellular isoform, PrPC, to PrPSc, an infectious isoform capable of causing neurodegenerative diseases in many mammals. Conversion can be triggered by low pH, and in vivo this appears to take place in an endocytic pathway and/or caveolae-like domains. It has thus far been impossible to characterize the conformational change at high resolution by experimental methods. Therefore, to investigate the effect of acidic pH on PrP conformation, we have performed 10-ns molecular dynamics simulations of PrPC in water at neutral and low pH. The core of the protein is well maintained at neutral pH. At low pH, however, the protein is more dynamic, and the sheet-like structure increases both by lengthening of the native β-sheet and by addition of a portion of the N terminus to widen the sheet by another two strands. The side chain of Met-129, a polymorphic codon in humans associated with variant Creutzfeldt–Jakob disease, pulls the N terminus into the sheet. Neutralization of Asp-178 at low pH removes interactions that inhibit conversion, which is consistent with the Asp-178–Asn mutation causing human prion diseases.

Identification of protein oligomerization states by analysis of interface conservation

The discrimination of true oligomeric protein–protein contacts from nonspecific crystal contacts remains problematic. Criteria that have been used previously base the assignment of oligomeric state on consideration of the area of the interface and/or the results of scoring functions based on statistical potentials. Both techniques have a high success rate but fail in more than 10% of cases. More importantly, the oligomeric states of several proteins are incorrectly assigned by both methods. Here we test the hypothesis that true oligomeric contacts should be identifiable on the basis of an increased degree of conservation of the residues involved in the interface. By quantifying the degree of conservation of the interface and comparing it with that of the remainder of the protein surface, we develop a new criterion that provides a highly effective complement to existing methods.

Attomole level protein sequencing by Edman degradation coupled with accelerator mass spectrometry

Edman degradation remains the primary method for determining the sequence of proteins. In this study, accelerator mass spectrometry was used to determine the N-terminal sequence of glutathione S-transferase at the attomole level with zeptomole precision using a tracer of 14C. The transgenic transferase was labeled by growing transformed Escherichia coli on [14C]glucose and purified by microaffinity chromatography. An internal standard of peptides on a solid phase synthesized to release approximately equal amounts of all known amino acids with each cycle were found to increase yield of gas phase sequencing reactions and subsequent semimicrobore HPLC as did a lactoglobulin carrier. This method is applicable to the sequencing of proteins from cell culture and illustrates a path to more general methods for determining N-terminal sequences with high sensitivity.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Copper-catalyzed oxidation of the recombinant SHa(29–231) prion protein

Metal-catalyzed oxidation may result in structural damage to proteins and has been implicated in aging and disease, including neurological disorders such as Alzheimer's disease and amyotrophic lateral sclerosis. The selective modification of specific amino acid residues with high metal ion affinity leads to subtle structural changes that are not easy to detect but may have dramatic consequences on physical and functional properties of the oxidized protein molecules. PrP contains a histidine-rich octarepeat domain that binds copper. Because copper-binding histidine residues are particularly prone to metal-catalyzed oxidation, we investigated the effect of this reaction on the recombinant prion protein SHaPrP(29–231). Using Cu2+/ascorbate, we oxidized SHaPrP(29–231) in vitro. Oxidation was demonstrated by liquid chromatography/mass spectrometry, which showed the appearance of protein species of higher mass, including increases in multiples of 16, characteristic of oxygen incorporation. Digestion studies using Lys C indicate that the 29–101 region, which includes the histidine-containing octarepeats, is particularly affected by oxidation. Oxidation was time- and copper concentration-dependent and was evident with copper concentrations as low as 1 μM. Concomitant with oxidation, SHaPrP(29–231) suffered aggregation and precipitation, which was nearly complete after 15 min, when the prion protein was incubated at 37°C with a 6-fold molar excess of Cu2+. These findings indicate that PrP, a copper-binding protein, may be particularly susceptible to metal-catalyzed oxidation and that oxidation triggers an extensive structural transition leading to aggregation.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

Hepatitis B virus transactivator protein, HBx, associates with the components of TFIIH and stimulates the DNA helicase activity of TFIIH.

Human hepatitis B virus genome encodes a protein, termed HBx, that is widely recognized as a transcriptional transactivator. While HBx does not directly bind cis-acting transcriptional control elements, it has been shown to associate with cellular proteins that bind DNA. Because HBx transactivated a large number of viral/cellular transcriptional control elements, we looked for its targets within the components of the basal transcriptional machinery. This search led to the identification of its interactions with TFIIH. Here, we show that HBx interacts with yeast and mammalian TFIIH complexes both in vitro and in vivo. These interactions between HBx and the components of TFIIH are supported by several lines of evidence including results from immunoprocedures and direct methods of measuring interactions. We have identified ERCC3 and ERCC2 DNA helicase subunits of holoenzyme TFIIH as targets of HBx interactions. Furthermore, the DNA helicase activity of purified TFIIH from rat liver and, individually, the ERCC2 component of TFIIH is stimulated in the presence of HBx. These observations suggest a role for HBx in transcription and DNA repair.

Phosducin is a ubiquitous G-protein regulator.

Phosducin is a 33-kDa cytosolic regulator of G-protein-mediated signaling that has previously been thought to be specific for retina and pineal gland. In this study, we show widespread tissue distribution of phosducin by the amplification of its cDNA and the detection of two different transcripts in Northern analyses in liver, lung, heart, brain, and retina. On the protein level, phosducin could be detected in 12 bovine tissues by immune precipitation and subsequent Western analysis using anti-phosducin antibodies generated in two different species. Masking of phosducin in direct Western blots appears to explain the failure to detect phosducin in earlier studies. The concentration of phosducin in bovine brain was calculated in the range of 10 pmol/mg total cytosolic protein (approximately 1 microM), whereas in the other tissues, it was slightly less. In these concentrations, phosducin inhibited receptor-stimulated adenylyl cyclase activity in cell membranes by about 50%. Taken together, our results indicate that phosducin is a ubiquitous regulator of G-protein function.

Reverse two-hybrid and one-hybrid systems to detect dissociation of protein-protein and DNA-protein interactions.

Macromolecular interactions define many biological phenomena. Although genetic methods are available to identify novel protein-protein and DNA-protein interactions, no genetic system has thus far been described to identify molecules or mutations that dissociate known interactions. Herein, we describe genetic systems that detect such events in the yeast Saccharomyces cerevisiae. We have engineered yeast strains in which the interaction of two proteins expressed in the context of the two-hybrid system or the interaction between a DNA-binding protein and its binding site in the context of the one-hybrid system is deleterious to growth. Under these conditions, dissociation of the interaction provides a selective growth advantage, thereby facilitating detection. These methods referred to as the "reverse two-hybrid system" and "reverse one-hybrid system" facilitate the study of the structure-function relationships and regulation of protein-protein and DNA-protein interactions. They should also facilitate the selection of dissociator molecules that could be used as therapeutic agents.