1
|
Masso M, Mathe E, Parvez N, Hijazi K, Vaisman II. Modeling the functional consequences of single residue replacements in bacteriophage f1 gene V protein. Protein Eng Des Sel 2009; 22:665-71. [PMID: 19690089 DOI: 10.1093/protein/gzp050] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A computational mutagenesis methodology utilizing a four-body, knowledge-based, statistical contact potential is applied toward globally quantifying relative environmental perturbations (residual scores) in bacteriophage f1 gene V protein (GVP) due to single amino acid substitutions. We show that residual scores correlate well with experimentally measured relative changes in protein function upon mutation. Residual scores also distinguish between GVP amino acid positions grouped according to protein structural or functional roles or based on similarities in physicochemical characteristics. For each mutant, the in silico mutagenesis additionally yields local measures of environmental change (EC scores) occurring at every residue position (residual profile) relative to the native protein. Implementation of the random forest (RF) algorithm, utilizing experimental GVP mutants whose feature vector components include EC scores at the mutated position and at six structurally nearest neighbors, correctly classifies mutants based on function with up to 77% cross-validation accuracy while achieving 0.82 area under the receiver operating characteristic curve. A control experiment highlights the effectiveness of mutant feature vector signals, and a variety of learning curves are generated to analyze the impact of GVP mutant data set size on performance measures. An optimally trained RF model is subsequently used for inferring function for all the remaining unexplored GVP mutants.
Collapse
Affiliation(s)
- Majid Masso
- Laboratory for Structural Bioinformatics, Department of Bioinformatics and Computational Biology, George Mason University, Manassas, VA 20110, USA.
| | | | | | | | | |
Collapse
|
2
|
Mou TC, Sreerama N, Terwilliger TC, Woody RW, Gray DM. Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions. Protein Sci 2002; 11:601-13. [PMID: 11847282 PMCID: PMC2373458 DOI: 10.1110/ps.30002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Abstract
The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.
Collapse
Affiliation(s)
- Tung-Chung Mou
- Department of Molecular and Cell Biology, The University of Texas at Dallas, Richardson, Texas 75083-0688, USA
| | | | | | | | | |
Collapse
|
3
|
Su S, Gao YG, Zhang H, Terwilliger TC, Wang AH. Analyses of the stability and function of three surface mutants (R82C, K69H, and L32R) of the gene V protein from Ff phage by X-ray crystallography. Protein Sci 1997; 6:771-80. [PMID: 9098886 PMCID: PMC2144760 DOI: 10.1002/pro.5560060403] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The high-resolution crystal structure of the gene V protein (GVP) from the Ff filamentous phages (M13, fl, fd) has been solved recently for the wild-type and two surface mutant (Y41F and Y41H) proteins, leading to a plausible model for the polymeric GVP-ssDNA complex (Guan Y, Zhang H, Wang AHJ, 1995, Protein Sci 4:187-197). The model of the complex shows extensive contacts between neighboring dimer GVPs involving electrostatic interactions between the K69 from one and the D79 and R82 from the next dimer. In addition, hydrophobic interactions between the amino acids L32 and L44 from one and G23 from the next dimer also contribute to the dimer-dimer interactions. Mutations at the L32, K69, and R82 amino acid sites generally destabilize the protein and many of these affect the function of the phage. We have studied the structural effects of three mutant proteins involving those sites, i.e., L32R, K69H, and R82C, by X-ray crystallographic analysis at 2.0 A resolution. In L32R GVP, the structural perturbation is localized, whereas in K69H and R82C GVPs, some long-range effects are also detected in addition to the local perturbation. We have interpreted the protein stability and the functional properties associated with those mutations in terms of the observed structural perturbations.
Collapse
Affiliation(s)
- S Su
- Department of Cell & Structural Biology, University of Illinois at Urbana-Champaign 61801, USA
| | | | | | | | | |
Collapse
|
4
|
Affiliation(s)
- J J Hill
- School of Pharmacy, University of Wisconsin-Madison 53706, USA
| | | |
Collapse
|
5
|
Villemain JL, Giedroc DP. Characterization of a cooperativity domain mutant Lys3 --> Ala (K3A) T4 gene 32 protein. J Biol Chem 1996; 271:27623-9. [PMID: 8910351 DOI: 10.1074/jbc.271.44.27623] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
The N-terminal "B" domain of T4 gene 32 protein contains a Lys3-Arg4-Lys5 sequence that has been postulated to provide a major determinant for cooperative binding. In this report, the equilibrium binding properties of a Lys3 --> Ala substitution mutant of gp32 (K3A gp32) and described and compared to a set of substitution mutants of Arg4 previously described (Villemain, J. L., and Giedroc, D. P. (1993) Biochemistry 32, 11235-11246) and further characterized here. K3A gp32 exhibits binding behavior which mirrors that of R4Q gp32. Despite an 6-8-fold decrease in overall binding affinity (Kapp = Kint x omega) at pH 8.1, 0.20 M NaCl, 20 degrees C, the magnitude of the cooperativity parameter is at most 2-3-fold smaller than that of the wild-type protein. The magnitude of omega is independent of salt concentration and type over the range in [NaCl] from 0.125 to 0. 225 M and [NaF] between 0.20 and 0.32 M (log omega = 2.86 +/- 0.19). For comparison, log omega for wild-type gp32 is 2.91 (+/- 0.21) resolved at 0.275 M NaCl and 3.39 +/- 0.11 in [NaF] between 0.40 and 0.45 M. In contrast to omega, the [NaCl] dependence of Kapp is large and markedly nonlinear for both wild-type and K3A gp32s over a [NaCl] range extending from 0.05 M to 0.40 M NaCl. Modeling of the complete salt dependence of Kapp for wild-type, K3A, and R4T gp32s in NaCl and NaF with a simple ion-exchange model suggests that substitutions within the basic Lys3-Arg4-Lys5 sequence do not strongly modulate the net displacement of cations and anions upon poly(A) complex formation by gp32.
Collapse
Affiliation(s)
- J L Villemain
- Department of Biochemistry and Biophysics, Center for Macromolecular Design, Institute of Biosciences and Technology, Texas A&M University, College Station, Texas 77843-2128, USA.
| | | |
Collapse
|
6
|
Konings RN, Folmer RH, Folkers PJ, Nilges M, Hilbers CW. Three-dimensional structure of the single-stranded DNA-binding protein encoded by gene V of the filamentous bacteriophage M13 and a model of its complex with single-stranded DNA. FEMS Microbiol Rev 1995. [DOI: 10.1111/j.1574-6976.1995.tb00188.x] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
|
7
|
Guan Y, Zhang H, Wang AH. Electrostatic potential distribution of the gene V protein from Ff phage facilitates cooperative DNA binding: a model of the GVP-ssDNA complex. Protein Sci 1995; 4:187-97. [PMID: 7757008 PMCID: PMC2143068 DOI: 10.1002/pro.5560040206] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
The crystal structure of the gene V protein (GVP) from the Ff filamentous phages (M13, fl, fd) has been solved for the wild-type and two mutant (Y41F and Y41H) proteins at high resolution. The Y41H mutant crystal structure revealed crystal packing interactions, which suggested a plausible scheme for constructing the polymeric protein shell of the GVP-single-stranded DNA (ssDNA) complex (Guan Y, et al., 1994, Biochemistry 33:7768-7778). The electrostatic potentials of the isolated and the cooperatively formed protein shell have been calculated using the program GRASP and they revealed a highly asymmetric pattern of the electrostatic charge distribution. The inner surface of the putative DNA-binding channel is positively charged, whereas the opposite outer surface is nearly neutral. The electrostatic calculation further demonstrated that the formation of the helical protein shell enhanced the asymmetry of the electrostatic distribution. A model of the GVP-ssDNA complex with the n = 4 DNA-binding mode could be built with only minor conformational perturbation to the GVP protein shell. The model is consistent with existing biochemical and biophysical data and provides clues to the properties of GVP, including the high cooperatively of the protein binding to ssDNA. The two antiparallel ssDNA strands form a helical ribbon with the sugar-phosphate backbones at the middle and the bases pointing away from each other. The bases are stacked and the Phe 73 residue is intercalated between two bases. The optimum binding to a tetranucleotide unit requires the participation of four GVP dimers, which may explain the cooperativity of the GVP binding to DNA.
Collapse
Affiliation(s)
- Y Guan
- Biophysics Division, University of Illinois at Urbana-Champaign 61801, USA
| | | | | |
Collapse
|
8
|
Turner GP, Kneale GG. Site-directed mutagenesis of the M13 gene 5 protein: the role of Arg-21, Tyr-26 and Tyr-41. BIOCHIMICA ET BIOPHYSICA ACTA 1995; 1260:79-84. [PMID: 7999799 DOI: 10.1016/0167-4781(94)00174-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The gene 5 protein of bacteriophage M13 is a single stranded DNA binding protein essential for phage replication. We have generated the mutations R21A, Y26F and Y41A in the gene 5 protein and purified the mutant proteins for functional characterisation in vitro. The complex of Y26F with single-stranded DNA is disrupted at 0.8 M NaCl, the same salt concentration as that required to dissociate the native complex. However, the mutant proteins R21A and Y41A are considerably less stable and dissociate from single-stranded DNA at at 0.4 M NaCl. The fluorescence of the mutant proteins and the DNA-protein complexes they form has been compared with the wild-type protein to allow an assessment of the contribution from individual residues. We conclude that the fluorescence of Tyr-26 is 50% quenched in the complex with DNA, whereas that of Tyr-41 is fully quenched. Fluorescence titrations of the mutant proteins with poly(dT) show that all three mutant proteins can bind DNA but, in the case of Y41A, with a change of stoichiometry suggesting a loss of cooperativity. Gel retardation analysis of Y41A also shows anomalous behaviour in binding to oligonucleotides, consistent with the proposed involvement of Tyr-41 in dimer-dimer contacts in the nucleoprotein complex.
Collapse
Affiliation(s)
- G P Turner
- Biophysics Laboratories, University of Portsmouth, UK
| | | |
Collapse
|
9
|
Stassen AP, Folmer RH, Hilbers CW, Konings RN. Single-stranded DNA binding protein encoded by the filamentous bacteriophage M13: structural and functional characteristics. Mol Biol Rep 1995; 20:109-27. [PMID: 7565651 DOI: 10.1007/bf00990543] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
The single-stranded DNA binding protein, or gene V protein (gVp), encoded by gene V of the filamentous bacteriophage M13 is a multifunctional protein that not only regulates viral DNA replication but also gene expression at the level of mRNA translation. It furthermore is implicated as a scaffolding and/or chaperone protein during the phage assembly process at the hostcell membrane. The protein is 87 amino acids long and its biological functional entity is a homodimer. In this manuscript a short description of the life cycle of filamentous phages is presented and our current knowledge of the major functional and structural properties and characteristics of gene V protein are reviewed. In addition models of the superhelical complexes gVp forms with ssDNA are described and their (possible) biological meaning in the infection process are discussed. Finally it is described that the 'DNA binding loop' of gVp is a recurring motif in many ssDNA binding proteins and that the fold of gVp is shared by a large family of evolutionarily conserved gene regulatory proteins.
Collapse
Affiliation(s)
- A P Stassen
- Department of Molecular Biology, University of Nijmegen, The Netherlands
| | | | | | | |
Collapse
|
10
|
Skinner MM, Zhang H, Leschnitzer DH, Guan Y, Bellamy H, Sweet RM, Gray CW, Konings RN, Wang AH, Terwilliger TC. Structure of the gene V protein of bacteriophage f1 determined by multiwavelength x-ray diffraction on the selenomethionyl protein. Proc Natl Acad Sci U S A 1994; 91:2071-5. [PMID: 8134350 PMCID: PMC43311 DOI: 10.1073/pnas.91.6.2071] [Citation(s) in RCA: 88] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- M M Skinner
- Life Sciences Division, Los Alamos National Laboratory, NM 87545
| | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Folkers PJ, van Duynhoven JP, van Lieshout HT, Harmsen BJ, van Boom JH, Tesser GI, Konings RN, Hilbers CW. Exploring the DNA binding domain of gene V protein encoded by bacteriophage M13 with the aid of spin-labeled oligonucleotides in combination with 1H-NMR. Biochemistry 1993; 32:9407-16. [PMID: 8396429 DOI: 10.1021/bi00087a020] [Citation(s) in RCA: 23] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The DNA binding domain of the single-stranded DNA binding protein gene V protein encoded by the bacteriophage M13 was studied by means of 1H nuclear magnetic resonance, through use of a spin-labeled deoxytrinucleotide. The paramagnetic relaxation effects observed in the 1H-NMR spectrum of M13 GVP upon binding of the spin-labeled ligand were made manifest by means of 2D difference spectroscopy. In this way, a vast data reduction was accomplished which enabled us to check and extend the analysis of the 2D spectra carried out previously as well as to probe the DNA binding domain and its surroundings. The DNA binding domain is principally situated on two beta-loops. The major loop of the two is the so-called DNA binding loop (residues 16-28) of the protein where the residues which constitute one side of the beta-ladder (in particular, residues Ser20, Tyr26, and Leu28) are closest to the DNA spin-label. The other loop is part of the so-called dyad domain of the protein (residues 68-78), and mainly its residues at the tip are affected by the spin-label (in particular, Phe73). In addition, a part of the so-called complex domain of the protein (residues 44-51) which runs contiguous to the DNA binding loop is in close vicinity to the DNA. The NMR data imply that the DNA binding domain is divided over two monomeric units of the GVP dimer in which the DNA binding loop and the tip of the dyad loop are part of opposite monomers. The view of the GVP-ssDNA binding interaction which emerges from our data differs from previous molecular modeling proposals which were based on the GVP crystal structure (Brayer & McPherson, 1984; Hutchinson et al., 1990). These models implicate the involvement of one or two tyrosines (Tyr34, Tyr41) of the complex loop of the protein to participate in complex formation with ssDNA. In the NMR studies with the spin-labeled oligonucleotides, no indication of such interactions has been found. Other differences between the models and our NMR data are related to the structural differences found when solution and crystal structures are compared.
Collapse
Affiliation(s)
- P J Folkers
- NSR Center (Nijmegen Son Research Center), University of Nijmegen, The Netherlands
| | | | | | | | | | | | | | | |
Collapse
|
12
|
O'Donohue M, Scarlett G, Kneale G. Tyr26 and Phe73 are essential for full biological activity of the Fd Gene 5 protein. FEMS Microbiol Lett 1993. [DOI: 10.1111/j.1574-6968.1993.tb06171.x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
|