Engineering the stability and function of gene V protein.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to T C Terwilliger.
Explore the source record for details and available documents.
Gene V protein (GVP) encoded by the filamentous phage Ff (M13, fl, fd) is a homodimeric protein of 87 amino acids that binds to single-stranded DNA (ssDNA) nonspecifically and cooperatively. The structure (monoclinic C2 form) of the wild-type protein has been determined and refined at 1.8-A resolution [Skinner et al. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 2071-2075]. The monomer structure consists of a somewhat distorted five-stranded beta-barrel core with three prominent loops: a DNA-binding loop, a dyad loop, and a dimer contact loop. The amino acid residue at position 41 plays an important role in the dimer-dimer interactions of the protein-ssDNA complex. Two Y41 mutant structures have been studied by X-ray crystallography. The Y41F GVP structure has been refined to an R-factor of 0.180 at 2.2-A resolution and is very similar to the wild-type (wt) structure (rmsd of all C alpha atoms = 0.30 A). In contrast, Y41H GVP forms a new crystal lattice in the space group P2(1)2(1)2(1) with a = 77.18 A, b = 84.17 A, and c = 28.62 A. Its structure has been solved by the molecular replacement method and refined to an R-factor of 0.170 at 2.5-A resolution. The two monomers of Y41H are crystallographically independent, and their structures remain similar to wt-GVP but with significant differences, particularly in the DNA-binding hairpin region. In both crystals, the loop (residues 36-43) that contains the Y41 residue is involved in the crystal dimer packings but in a different manner. The dimer-dimer contacts found in the wt-GVP crystal may be important for GVP aggregation in the absence of DNA. In the presence of DNA, the dimer-dimer contacts may switch to the type found in the Y41H crystal, allowing the GVP-ssDNA complex to form cooperatively. A model of the complex, consistent with existing biochemical and biophysical data, has been constructed from those crystal packing data.
The crystal structure of the dimeric gene V protein of bacteriophage f1 was determined using multiwavelength anomalous diffraction on the selenomethionine-containing wild-type and isoleucine-47-->methionine mutant proteins with x-ray diffraction data phased to 2.5 A resolution. The structure of the wild-type protein has been refined to an R factor of 19.2% using native data to 1.8 A resolution. The structure of the gene V protein was used to obtain a model for the protein portion of the gene V protein-single-stranded DNA complex.
The gene V protein of bacteriophage f1 binds to single-stranded nucleic acids and is essential for propagation of phage f1. We tested the function of gene V protein mutants with single amino acid substitutions in two ways: by the ability of the mutant proteins to support phage growth, and by the ability of the mutant proteins, when expressed at high levels, to inhibit the growth of Escherichia coli. The results of the tests were used to identify sites in the protein that are relatively tolerant or intolerant to substitution, where tolerant sites are defined as those where most substitutions do not affect the function of the protein. The two assays generally yielded similar results for the tolerance of sites to substitution. Many sites that are less than 10% exposed to the solvent are relatively intolerant of substitution, with even very conservative substitutions leading to loss of function in some cases such as Ile to Leu at residue 6. Some buried sites such as Ile47 are more tolerant, with even a substitution of Ile to Thr leading to a functional protein based on the ability of the proteins to inhibit the growth of E. coli. Some surface sites in the protein (> 10% exposure to solvent) that are thought to be near the location of bound oligonucleotides, such as Arg16, Val19, Ser20, Arg21, Tyr26, Lys46 and Arg80, are sensitive to substitution. Other side-chains thought to be close to bound oligonucleotides, including Leu28, can be replaced with a number of amino acids with little loss of function based on either assay. Most non-Gly/Pro surface residues thought to be distant from the locations of bound oligonucleotides are relatively tolerant of substitution, except for two small residues (Ala11 and Thr14), two aromatic residues (Tyr34 and Tyr56), two residues that are only partially exposed to solvent (Asn29 and Val70), and three residues that have been proposed to be at the dimer-dimer interface formed when gene V protein binds to nucleic acids (Glu40, Tyr41 and Arg82).
A Bayesian approach is applied to the calculation of Patterson functions and cross-Fourier maps in the analysis of multi-wavelength anomalous-diffraction (MAD) data. This procedure explicitly incorporates information available a priori on the likely magnitudes of partial structure factors (F(A)) corresponding to the anomalously scattering atoms, uses weighted-average estimates of F(A), and incorporates estimates of errors in the data that are not represented in the instrumental uncertainties. The method is demonstrated by application to MAD data collected on selenomethionine-containing gene V protein.
A framework analogous to that used for the analysis of data in the method of isomorphous replacement is applied to the analysis of multi-wavelength anomalous-diffraction (MAD) data. The present method is advantageous in that refinement of parameters describing the anomalously scattering atoms in the structure can be carried out using nearly all the data, and in that phase information can be readily combined. The procedure described here involves an approximation in which it is assumed that the magnitude of the structure factor corresponding to anomalously scattering atoms in the unit cell is small compared to that corresponding to all other atoms in the structure. A model calculation is applied to a protein crystal with 682 non-H atoms in the asymmetric unit and two Se atoms as the anomalous scatterers. It is shown using this model calculation that the approximation used in this analysis does not substantially affect the accuracy of phase calculations for this MAD data. The method is demonstrated by application to MAD data collected on gene V protein.
A method for simultaneously engineering multiple properties of a protein, based on the observed additivity of effects of individual mutations, is presented. We show that, for the gene V protein of bacteriophage f1, effects of double mutations on both protein stability and DNA binding affinity are approximately equal to the sums of the effects of the constituent single mutations. This additivity of effects implies that it is possible to deliberately construct mutant proteins optimized for multiple properties by combination of appropriate single mutations chosen from a characterized library.
The gene V protein of bacteriophage f1 is a single-stranded DNA and RNA-binding protein composed of two identical subunits. We have constructed single-chain variants of the protein using short peptide linkers of five or six amino acids to connect the carboxyl terminus of one monomer to the amino terminus of the second monomer. The resulting subunit-fusion gene V proteins were found to bind single-stranded DNA nearly as tightly as the wild-type protein. Denaturation measurements show that the subunit-fusion gene V proteins are 5 kcal/mol (1 kcal = 4.18 kJ) more stable than the wild-type protein at a protein concentration of 10 microM. The rate of unfolding of the protein is essentially unaffected by the fusion of monomeric subunits, whereas the rate of folding is greatly enhanced. Our results suggest a simple way of obtaining a substantial thermodynamic stabilization for some oligomeric proteins.
The relative activities of 313 mutants of the gene V protein of bacteriophage f1, assayed in vivo, have been used to evaluate two approaches to predicting the effects of single amino acid substitutions on the function of a protein. First, we tested methods that only depend on the properties of the wild-type and substituting amino acids. None of the properties or measures of the functional equivalence of amino acids we tested, including the frequency of exchange of amino acids among homologous proteins as well as changes in side-chain size, hydrophobicity, and charge, were found to be more than weakly correlated with the activities of mutants. The principal reason for this poor correlation was found to be that the effect of a particular substitution varies considerably from site to site. We then tested an approach using the activities of several mutants with substitutions at a site to predict the activity of another mutant, and we find that this is a relatively good indicator of whether the other mutant at that site will be functional. A predictive scheme was developed that combines the weak information from the models depending on the properties of the wild-type and substituting amino acids with the stronger information from the tolerance of a site to substitution. Although this scheme requires no knowledge of the structure of a mutant protein, it is useful in predicting the activities of mutants.
In vivo selections were used to isolate 43 temperature-sensitive gene V mutants of the bacteriophage f1 from a collection of mutants constructed by saturation mutagenesis of the gene. The sites of temperature-sensitive substitutions are found in both the beta-sheets and the turns of the protein, and some sites are exposed to the solvent while others are not. Thirteen of the variant proteins were purified and characterized to evaluate their free energy changes upon unfolding and their affinities for single-stranded DNA, and eight were tested for their tendencies to aggregate at 42 degrees C. Each of the three temperature-sensitive mutants at buried sites and six of ten at surface sites had free energy changes of unfolding substantially lower (less stabilizing) than the wild-type at 25 degrees C. A seventh mutant at a surface site had a substantially altered unfolding transition and its free energy of unfolding was not estimated. The affinities of the mutant proteins for single-stranded DNA varied considerably, but two mutants at a surface site, Lys69, had much weaker binding to single-stranded DNA than any of the other mutants, while two mutants at another surface site, Glu30, had the highest DNA-binding affinities. The wild-type gene V protein is stable at 42 degrees C, but six of the eight mutants tested aggregated within a few minutes and the remaining two aggregated within 30 minutes at this temperature. Overall, each of the temperature-sensitive proteins tested had a tendency to aggregate at 42 degrees C, and most also had either a low free energy of unfolding (at 25 degrees C), or weak DNA binding. We suggest that any of these properties can lead to a temperature-sensitive gene V phenotype.
The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.
To test whether interactions in the hydrophobic core of a protein can be adequately modeled based on the properties of a liquid hydrocarbon, we measured the unfolding free energies of the wild-type bacteriophage f1 gene V protein and 29 mutants with apolar substitutions at positions 35 and 47. Stability changes arising from identical mutations at these two buried sites are quite different, suggesting that one site is more rigid than the other. Reversals of residues at positions 35 and 47 confirm that their environments are distinct. Mutants containing weakly polar residues at these two sites suggest that the protein interior is more polar than a liquid hydrocarbon. Interactions between residues at the two sites appear to be minimal. These observations are compatible with a view of protein interiors that incorporates properties of liquid hydrocarbons but also includes polar interactions and a site-dependent "packing energy" associated with changes in internal structure.
Several goals of protein engineering may be achieved through redesign and repacking of protein interiors. The effects of interior apolar substitutions on protein stability depend strongly on the site of the substitution. One reason for this is that protein interiors have properties both of apolar liquids and of crystalline solids. Substitutions at interior sites affect the stability of a protein by changing the hydrophobicity, but each site in a protein has a characteristic energy associated with introducing packing changes, and the net stability depends on both of these factors.
Protein interiors contain many tightly packed apolar atoms in a nearly crystalline state. Both shielding of apolar atoms from solvent and efficient interior packing arrangements affect protein stability, but their relative importance is unclear. To separate these effects, the stabilities of wild-type and mutant gene V proteins from bacteriophage fl were studied by measuring resistance to denaturation. The effects of subtle interior packing changes, both separate from and combined with changes in buried side chain hydrophobicity, were measured. For the interior apolar-to-apolar substitutions studied, the two effects were of the same magnitude and alteration of packing without accompanying hydrophobicity changes substantially destabilized the protein.
A simple and widely applicable procedure for constructing synthetic variants of a gene, involving the synthesis of only one strand of DNA, has been developed. The method is suited for cases in which a cloned DNA with a sequence related to the gene to be constructed is available. First, a heteroduplex DNA which is single-stranded throughout the region of interest is made. This single-stranded region is then used as a template to correctly align and allow ligation of synthetic oligos corresponding to the entire gene. To favor the replication of the strand encoding the synthetic gene, a template strand containing some substitutions of deoxyuridine for deoxythymidine is used. This procedure was used to construct a synthetic bacteriophage f1 gene V which differs from the wild-type (wt) gene at 45 positions out of 298. The synthetic gene was designed to include nine restriction sites without altering the sequence of the encoded DNA-binding protein. The gene construction was found to be very efficient, and about 40% of the resulting plasmids contained the desired synthetic gene. The synthetic gene was found to be fully active and could substitute for the wt gene in bacteriophage f1.
A simple and highly efficient procedure for oligodeoxynucleotide (oligo)-directed mutagenesis has been developed. In this procedure, a gapped heteroduplex DNA is first constructed and purified. The gapped heteroduplex consists of a circular 'template' strand of DNA, which contains some misincorporated deoxyuridine nucleotides, and a complementary strand which does not contain deoxyuridine, and which lacks a defined segment. Making a specific change in the sequence of the DNA within the gapped region then only requires ligation and transformation. An oligo, exactly the same length as the gap, and with the desired sequence, is synthesized, purified, and ligated directly into the gap in the heteroduplex. When this DNA is used to transform wt (ung+) Escherichia coli, about 80% of the resulting plasmids contain the sequence determined by the synthetic oligo. One gapped heteroduplex preparation can be used for many mutagenesis experiments, so that this procedure is well-suited for producing a series of defined mutations within a defined target region flanked by sites for restriction enzyme cleavage. As the method does not require a polymerase, the effects of primer displacement and polymerase infidelity are avoided.
Complementary negative and positive genetic selections based on the activity of a plasmid-encoded bacteriophage f1 gene V are developed. The negative selection is based on an activity of the gene V protein in E. coli cells which markedly reduces the infection of those cells by f1-related viruses. In order to select against cells expressing active gene V protein, the cells are infected with the p'age R386, a derivative of f1 which confers resistance to chloramphenicol, and are plated in the presence of the antibiotic. Those cells which contain gene V protein are infrequently infected with the virus and are unable to grow in the presence of chloramphenicol; those which do not contain the gene V protein are readily infected and can grow in the presence of the antibiotic. The positive genetic selection consists of excising the gene V sequences from the plasmids and using them to replace the gene V of a bacteriophage f1 derivative containing an amber mutation in gene V. Only those genes which encode an active gene V protein can support phage growth and yield plaques. The two genetic selections can be combined in order to yield a substantial enrichment for genes encoding temperature-sensitive gene V proteins.
A method for determining the extent of methyl esterification of each of the four potential sites on the aspartate receptors involved in chemotaxis in Escherichia coli and Salmonella typhimurium is presented. In this procedure, radioactive methyl esters are incorporated into the receptors, the receptors are cleaved by trypsin and the V8 protease from Staphylococcus aureus, and the four fragments containing sites of methylation are separated by high performance liquid chromatography. Using this technique, we find that the rate of methyl esterification increases at all four sites after stimulation with the "attractant" aspartate, suggesting that all four sites of modification are involved in adaptation to aspartate. We also find that the rate of methyl esterification at each site is correlated with the homology between the protein sequence at that site and the "consensus" sequence, Glu-Glu-X-X-Ala-Thr/Ser.