Influence of interior packing and hydrophobicity on the stability of a protein.

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.

Rapid protein-folding assay using green fluorescent protein.

Formation of the chromophore of green fluorescent protein (GFP) depends on the correct folding of the protein. We constructed a "folding reporter" vector, in which a test protein is expressed as an N-terminal fusion with GFP. Using a test panel of 20 proteins, we demonstrated that the fluorescence of Escherichia coli cells expressing such GFP fusions is related to the productive folding of the upstream protein domains expressed alone. We used this fluorescent indicator of protein folding to evolve proteins that are normally prone to aggregation during expression in E. coli into closely related proteins that fold robustly and are fully soluble and functional. This approach to improving protein folding does not require functional assays for the protein of interest and provides a simple route to improving protein folding and expression by directed evolution.

Structure of translation initiation factor 5A from Pyrobaculum aerophilum at 1.75 A resolution.

BACKGROUND: Translation initiation factor 5A (IF-5A) is reported to be involved in the first step of peptide bond formation in translation, to be involved in cell-cycle regulation and to be a cofactor for the Rev and Rex transactivator proteins of human immunodeficiency virus-1 and T-cell leukemia virus I, respectively. IF-5A contains an unusual amino acid, hypusine (N-epsilon-(4-aminobutyl-2-hydroxy)lysine), that is required for its function. The first step in the post-translational modification of lysine to hypusine is catalyzed by the enzyme deoxyhypusine synthase, the structure of which has been published recently. RESULTS: IF-5A from the archebacterium Pyrobaculum aerophilum has been heterologously expressed in Escherichia coli with selenomethionine substitution. The crystal structure of IF-5A has been determined by multiwavelength anomalous diffraction and refined to 1.75 A. Unmodified P. aerophilum IF-5A is found to be a beta structure with two domains and three separate hydrophobic cores. CONCLUSIONS: The lysine (Lys42) that is post-translationally modified by deoxyhypusine synthase is found at one end of the IF-5A molecule in an turn between beta strands beta4 and beta5; this lysine residue is freely solvent accessible. The C-terminal domain is found to be homologous to the cold-shock protein CspA of E. coli, which has a well characterized RNA-binding fold, suggesting that IF-5A is involved in RNA binding.

Osmotic water permeability of human red cells.

The osmotic water permeability of human red cells has been reexamined with a stopped-flow device and a new perturbation technique. Small osmotic gradients are used to minimize the systematic error caused by nonlinearities in the relationship between cell volume and light scattering. Corrections are then made for residual systematic error. Our results show that the hydraulic conductivity, Lp, is essentially independent of the direction of water flow and of osmolality in the range 184-365 mosM. the mean value of Lp obtained obtained was 1.8 +/- 0.1 (SEM) X 10-11 cm3 dyne -1 s-1.

Isolation and in vitro characterization of temperature-sensitive mutants of the bacteriophage f1 gene V protein.

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.

In vivo characterization of mutants of the bacteriophage f1 gene V protein isolated by saturation mutagenesis.

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).

Context dependence of mutational effects in a protein: the crystal structures of the V35I, I47V and V35I/I47V gene V protein core mutants.

The basis for the context dependence of the effects of core mutations on protein stability was investigated by comparing the structures of three gene V protein mutants with that of the wild-type protein. We previously examined a "swapped" mutant in which core residues Val35 and Ile47 were simply reversed so that the mutant had no hydrophobicity change from the native protein. The swapped mutant was destabilized by 3 kcal/mol per gene V protein dimer relative to the wild-type protein, demonstrating that factors other than hydrophobicity must make substantial contributions to the effects of mutations on the stability of the protein. Here we have determined the structure of this swapped mutant (V35I/I47V) as well as those of the two constituent mutants (V35I and I47V). We find that the structures of the mutant proteins are very similar to that of the wild-type protein except for the necessary addition or deletion of methylene groups and for slight positional shifts of atoms around each mutated residue. The structure of the double mutant is a composite of the structures of the two single mutants. In the mutant structures, the V35I mutation fills a cavity that exists in the wild-type protein and the I47V mutation creates a new cavity. The structures of the mutants indicate further that the reason the V35I and I47V mutations do not have opposite effects on stability is that the cavity in the wild-type protein filled by the V35I mutation is not optimally shaped for accommodating the additional methylene group of the isoleucine. These results support the concepts that the details of core packing have substantial influence on the effects of core mutations on protein stability and that these packing effects are major determinants of the context dependence of core mutation effects on stability.

Scission of DNA at a preselected sequence using a single-strand-specific chemical nuclease.

BACKGROUND: We were interested in developing a protocol for cleaving large DNAs specifically. Previous attempts to develop such methods have failed to work because of high levels of nonspecific background scission. RESULTS: R-loop formation was chosen for sequence-specific targeting, a method of hybridization whereby an RNA displaces a DNA strand of identical sequence in 70% formamide using Watson-Crick base-pairing, leading to a three-stranded structure. R-loops are stabilized in aqueous solution by modifying the bases with chemical reagents. The R-loop was cleaved using a novel nuclease prepared from the Thr48-->Cys mutant of the single-strand-specific M-13 gene V protein (GVP), which was alkylated with 5-(iodoacetamido-beta-alanyl)1,10-phenanthroline. The cleavage products of the pGEM plasmid were cloned in to the pCR 2.1-TOPO vector. Adenovirus 2 DNA (35.8 kb; tenfold larger than the pGEM plasmid) was also cleaved quantitatively at a preselected sequence. CONCLUSIONS: A new method for cleaving duplex DNA at any preselected sequence was developed. The cleavage method relies on the chemical conversion of M-13 GVP into a nuclease, reflecting GVP's specificity for single-stranded DNA. The GVP chimera is the first example of a semisynthetic secondary structure specific nuclease. The chemical nuclease activity of 1,10-phenanthroline-copper is uniquely suited to this technique because it oxidizes the deoxyribose moiety without generating diffusible intermediates, providing clonable DNA fragments. The protocol could be useful in generating large DNA fragments for mapping the contiguity of probes or defining the exon-intron structure of transcription units.

Repacking protein interiors.

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.

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.

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.

Class-directed structure determination: foundation for a protein structure initiative.

The recent sequencing of many complete genomes, combined with the development of methods that allow rapid structure determination for many proteins, has changed the way in which protein structure determinations can be approached. One-by-one determinations of individual protein structures will soon be augmented by class-directed structure analyses in which a group of proteins is targeted and structures of representative members are determined and used to represent the entire group. Such a shift in approach would be the foundation for a broad protein structure initiative targeting classes of proteins important for biotechnology and for a fundamental understanding of protein function.

Simple and highly efficient site-specific mutagenesis, by ligation of an oligodeoxyribonucleotide into gapped heteroduplex DNA in which the template strand contains deoxyuridine.

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.

Construction of a synthetic variant of the bacteriophage f1 gene V by assembling oligodeoxynucleotides corresponding to only one strand of DNA.

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.

The TB structural genomics consortium: a resource for Mycobacterium tuberculosis biology.

The TB Structural Genomics Consortium is an organization devoted to encouraging, coordinating, and facilitating the determination and analysis of structures of proteins from Mycobacterium tuberculosis. The Consortium members hope to work together with other M. tuberculosis researchers to identify M. tuberculosis proteins for which structural information could provide important biological information, to analyze and interpret structures of M. tuberculosis proteins, and to work collaboratively to test ideas about M. tuberculosis protein function that are suggested by structure or related to structural information. This review describes the TB Structural Genomics Consortium and some of the proteins for which the Consortium is in the progress of determining three-dimensional structures.

High throughput crystallography of TB drug targets.

Tuberculosis (TB) infects one-third of the world population. Despite 50 years of available drug treatments, TB continues to increase at a significant rate. The failure to control TB stems in part from the expense of delivering treatment to infected individuals and from complex treatment regimens. Incomplete treatment has fueled the emergence of multi-drug resistant (MDR) strains of Mycobacterium tuberculosis (Mtb). Reducing non-compliance by reducing the duration of chemotherapy will have a great impact on TB control. The development of new drugs that either kill persisting organisms, inhibit bacilli from entering the persistent phase, or convert the persistent bacilli into actively growing cells susceptible to our current drugs will have a positive effect. We are taking a multidisciplinary approach that will identify and characterize new drug targets that are essential for persistent Mtb. Targets are exposed to a battery of analyses including microarray experiments, bioinformatics, and genetic techniques to prioritize potential drug targets from Mtb for structural analysis. Our core structural genomics pipeline works with the individual laboratories to produce diffraction quality crystals of targeted proteins, and structural analysis will be completed by the individual laboratories. We also have capabilities for functional analysis and the virtual ligand screening to identify novel inhibitors for target validation. Our overarching goals are to increase the knowledge of Mtb pathogenesis using the TB research community to drive structural genomics, particularly related to persistence, develop a central repository for TB research reagents, and discover chemical inhibitors of drug targets for future development of lead compounds.

Maximum-likelihood density modification using pattern recognition of structural motifs.

The likelihood-based approach to density modification [Terwilliger (2000), Acta Cryst. D56, 965-972] is extended to include the recognition of patterns of electron density. Once a region of electron density in a map is recognized as corresponding to a known structural element, the likelihood of the map is reformulated to include a term that reflects how closely the map agrees with the expected density for that structural element. This likelihood is combined with other aspects of the likelihood of the map, including the presence of a flat solvent region and the electron-density distribution in the protein region. This likelihood-based pattern-recognition approach was tested using the recognition of helical segments in a largely helical protein. The pattern-recognition method yields a substantial phase improvement over both conventional and likelihood-based solvent-flattening and histogram-matching methods. The method can potentially be used to recognize any common structural motif and incorporate prior knowledge about that motif into density modification.

Map-likelihood phasing.

The recently developed technique of maximum-likelihood density modification [Terwilliger (2000), Acta Cryst. D56, 965-972] allows a calculation of phase probabilities based on the likelihood of the electron-density map to be carried out separately from the calculation of any prior phase probabilities. Here, it is shown that phase-probability distributions calculated from the map-likelihood function alone can be highly accurate and that they show minimal bias towards the phases used to initiate the calculation. Map-likelihood phase probabilities depend upon expected characteristics of the electron-density map, such as a defined solvent region and expected electron-density distributions within the solvent region and the region occupied by a macromolecule. In the simplest case, map-likelihood phase-probability distributions are largely based on the flatness of the solvent region. Though map-likelihood phases can be calculated without prior phase information, they are greatly enhanced by high-quality starting phases. This leads to the technique of prime-and-switch phasing for removing model bias. In prime-and-switch phasing, biased phases such as those from a model are used to prime or initiate map-likelihood phasing, then final phases are obtained from map-likelihood phasing alone. Map-likelihood phasing can be applied in cases with solvent content as low as 30%. Potential applications of map-likelihood phasing include unbiased phase calculation from molecular-replacement models, iterative model building, unbiased electron-density maps for cases where 2F(o) - F(c) or sigma(A)-weighted maps would currently be used, structure validation and ab initio phase determination from solvent masks, non-crystallographic symmetry or other knowledge about expected electron density.

Maximum-likelihood density modification.

A likelihood-based approach to density modification is developed that can be applied to a wide variety of cases where some information about the electron density at various points in the unit cell is available. The key to the approach consists of developing likelihood functions that represent the probability that a particular value of electron density is consistent with prior expectations for the electron density at that point in the unit cell. These likelihood functions are then combined with likelihood functions based on experimental observations and with others containing any prior knowledge about structure factors to form a combined likelihood function for each structure factor. A simple and general approach to maximizing the combined likelihood function is developed. It is found that this likelihood-based approach yields greater phase improvement in model and real test cases than either conventional solvent flattening and histogram matching or a recent reciprocal-space solvent-flattening procedure [Terwilliger (1999), Acta Cryst. D55, 1863-1871].