Patterned library analysis: A method for the quantitative assessment of hypotheses concerning the determinants of protein structure

Site-directed mutagenesis and combinatorial libraries are powerful tools for providing information about the relationship between protein sequence and structure. Here we report two extensions that expand the utility of combinatorial mutagenesis for the quantitative assessment of hypotheses about the determinants of protein structure. First, we show that resin-splitting technology, which allows the construction of arbitrarily complex libraries of degenerate oligonucleotides, can be used to construct more complex protein libraries for hypothesis testing than can be constructed from oligonucleotides limited to degenerate codons. Second, using eglin c as a model protein, we show that regression analysis of activity scores from library data can be used to assess the relative contributions to the specific activity of the amino acids that were varied in the library. The regression parameters derived from the analysis of a 455-member sample from a library wherein four solvent-exposed sites in an α-helix can contain any of nine different amino acids are highly correlated (P < 0.0001, R2 = 0.97) to the relative helix propensities for those amino acids, as estimated by a variety of biophysical and computational techniques.

The use of site-directed fluorophore labeling and donor–donor energy migration to investigate solution structure and dynamics in proteins

The use of molecular genetics for introducing fluorescent molecules enables the use of donor–donor energy migration to determine intramolecular distances in a variety of proteins. This approach can be applied to examine the overall molecular dimensions of proteins and to investigate structural changes upon interactions with specific target molecules. In this report, the donor–donor energy migration method is demonstrated by experiments with the latent form of plasminogen activator inhibitor type 1. Based on the known x-ray structure of plasminogen activator inhibitor type 1, three positions forming the corners of a triangle were chosen. Double Cys substitution mutants (V106C-H185C, H185C-M266C, and M266C-V106C) and corresponding single substitution mutants (V106C, H185C, and M266C) were created and labeled with a sulfhydryl specific derivative of BODIPY (=the D molecule). The side lengths of this triangle were obtained from analyses of the experimental data. The analyses account for the local anisotropic order and rotational motions of the D molecules, as well as for the influence of a partial DD-labeling. The distances, as determined from x-ray diffraction, between the Cα-atoms of the positions V106C–H185C, H185C–M266C, and M266C–V106C were 60.9, 30.8, and 55.1 Å, respectively. These are in good agreement with the distances of 54 ± 4, 38 ± 3, and 55 ± 3 Å, as determined between the BODIPY groups attached via linkers to the same residues. Although the positions of the D-molecules and the Cα-atoms physically cannot coincide, there is a reasonable agreement between the methods.

The structure of the acto-myosin subfragment 1 complex: Results of searches using data from electron microscopy and x-ray crystallography

Surmises of how myosin subfragment 1 (S1) interacts with actin filaments in muscle contraction rest upon knowing the relative arrangement of the two proteins. Although there exist crystallographic structures for both S1 and actin, as well as electron microscopy data for the acto–S1 complex (AS1), modeling of this arrangement has so far only been done “by eye.” Here we report fitted AS1 structures obtained using a quantitative method that is both more objective and makes more complete use of the data. Using undistorted crystallographic results, the best-fit AS1 structure shows significant differences from that obtained by visual fitting. The best fit is produced using the F-actin model of Holmes et al. [Holmes, K. C., Popp, D., Gebhard, W. & Kabsch, W. (1990) Nature (London) 347, 44–49]. S1 residues at the AS1 interface are now found at a higher radius as well as being translated axially and rotated azimuthally. Fits using S1 plus loops missing from the crystal structure were achieved using a homology search method to predict loop structures. These improved fits favor an arrangement in which the loop at the 50- to 20-kDa domain junction of S1 is located near the N terminus of actin. Rigid-body movements of the lower 50-kDa domain, which further improve the fit, produce closure of the large 50-kDa domain cleft and bring conserved residues in the lower 50-kDa domain into an apparently appropriate orientation for close interaction with actin. This finding supports the idea that binding of ATP to AS1 at the end of the ATPase cycle disrupts the actin binding site by changing the conformation of the 50-kDa cleft of S1.

Large-scale protein structure modeling of the Saccharomyces cerevisiae genome

The function of a protein generally is determined by its three-dimensional (3D) structure. Thus, it would be useful to know the 3D structure of the thousands of protein sequences that are emerging from the many genome projects. To this end, fold assignment, comparative protein structure modeling, and model evaluation were automated completely. As an illustration, the method was applied to the proteins in the Saccharomyces cerevisiae (baker’s yeast) genome. It resulted in all-atom 3D models for substantial segments of 1,071 (17%) of the yeast proteins, only 40 of which have had their 3D structure determined experimentally. Of the 1,071 modeled yeast proteins, 236 were related clearly to a protein of known structure for the first time; 41 of these previously have not been characterized at all.

A method that allows the assembly of kinetochore components onto chromosomes condensed in clarified Xenopus egg extracts

Kinetochores are complex macromolecular structures that link mitotic chromosomes to spindle microtubules. Although a small number of kinetochore components have been identified, including the kinesins CENP-E and XKCM1 as well as cytoplasmic dynein, neither how these and other proteins are organized to produce a kinetochore nor their exact functions within this structure are understood. For this reason, we have developed an assay that allows kinetochore components to assemble onto discrete foci on in vitro-condensed chromosomes. The source of the kinetochore components is a clarified cell extract from Xenopus eggs that can be fractionated or immunodepleted of individual proteins. Kinetochore assembly in these clarified extracts requires preincubating the substrate sperm nuclei in an extract under low ATP conditions. Immunodepletion of XKCM1 from the extracts prevents the localization of kinetochore-associated XKCM1 without affecting the targeting of CENP-E and cytoplasmic dynein or the binding of monomeric tubulin to the kinetochore. Extension of this assay for the analysis of other components should help to dissect the protein–protein interactions involved in kinetochore assembly and function.

Structure and function in rhodopsin: Kinetic studies of retinal binding to purified opsin mutants in defined phospholipid–detergent mixtures serve as probes of the retinal binding pocket

In the current standard procedure for preparation of mammalian rhodopsin mutants, transfected COS-1 cells expressing the mutant opsin genes are treated with 5 μM 11-cis-retinal before detergent solubilization for purification. We found that binding of 11-cis-retinal to opsin mutants with single amino acid changes at Trp-265 (W265F,Y,A) and a retinitis pigmentosa mutant (A164V) was far from complete and required much higher concentrations of 11-cis-retinal. By isolation of the expressed opsins in a stable form, kinetic studies of retinal binding to the opsins in vitro have been carried out by using defined phospholipid–detergent mixtures. The results show wide variation in the rates of 11-cis-retinal binding. Thus, the in vitro reconstitution procedure serves as a probe of the retinal-binding pocket in the opsins. Further, a method is described for purification and characterization of the rhodopsin mutants after retinal binding to the opsins in vitro.

Statistical significance of protein structure prediction by threading

In this study, we estimate the statistical significance of structure prediction by threading. We introduce a single parameter ɛ that serves as a universal measure determining the probability that the best alignment is indeed a native-like analog. Parameter ɛ takes into account both length and composition of the query sequence and the number of decoys in threading simulation. It can be computed directly from the query sequence and potential of interactions, eliminating the need for sequence reshuffling and realignment. Although our theoretical analysis is general, here we compare its predictions with the results of gapless threading. Finally we estimate the number of decoys from which the native structure can be found by existing potentials of interactions. We discuss how this analysis can be extended to determine the optimal gap penalties for any sequence-structure alignment (threading) method, thus optimizing it to maximum possible performance.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

Trends in lipoplex physical properties dependent on cationic lipid structure, vehicle and complexation procedure do not correlate with biological activity

Using a group of structurally related cytofectins, the effects of different vehicle constituents and mixing techniques on the physical properties and biological activity of lipoplexes were systematically examined. Physical properties were examined using a combination of dye accessibility assays, centrifugation, gel electrophoresis and dynamic light scattering. Biological activity was examined using in vitro transfection. Lipoplexes were formulated using two injection vehicles commonly used for in vivo delivery (PBS pH 7.2 and 0.9% saline), and a sodium phosphate vehicle previously shown to enhance the biological activity of naked pDNA and lipoplex formulations. Phosphate was found to be unique in its effect on lipoplexes. Specifically, the accessible pDNA in lipoplexes formulated with cytofectins containing a γ-amine substitution in the headgroup was dependent on alkyl side chain length and sodium phosphate concentration, but the same effects were not observed when using cytofectins containing a β-OH headgroup substitution. The physicochemical features of the phosphate anion, which give rise to this effect in γ-amine cytofectins, were deduced using a series of phosphate analogs. The effects of the formulation vehicle on transfection were found to be cell type-dependent; however, of the formulation variables examined, the liposome/pDNA mixing method had the greatest effect on transgene expression in vitro. Thus, though predictive physical structure relationships involving the vehicle and cytofectin components of the lipoplex were uncovered, they did not extrapolate to trends in biological activity.

Prediction of protein deamidation rates from primary and three-dimensional structure

A method for the quantitative estimation of instability with respect to deamidation of the asparaginyl (Asn) residues in proteins is described. The procedure involves the observation of several simple aspects of the three-dimensional environment of each Asn residue in the protein and a calculation that includes these observations, the primary amino acid residue sequence, and the previously reported complete set of sequence-dependent rates of deamidation for Asn pentapeptides. This method is demonstrated and evaluated for 23 proteins in which 31 unstable and 167 stable Asn residues have been reported and for 7 unstable and 63 stable Asn residues that have been reported in 61 human hemoglobin variants. The relative importance of primary structure and three-dimensional structure in Asn deamidation is estimated.

Galaxies and large scale structure at high redshifts

It is now straightforward to assemble large samples of very high redshift (z ∼ 3) field galaxies selected by their pronounced spectral discontinuity at the rest frame Lyman limit of hydrogen (at 912 Å). This makes possible both statistical analyses of the properties of the galaxies and the first direct glimpse of the progression of the growth of their large-scale distribution at such an early epoch. Here I present a summary of the progress made in these areas to date and some preliminary results of and future plans for a targeted redshift survey at z = 2.7–3.4. Also discussed is how the same discovery method may be used to obtain a “census” of star formation in the high redshift Universe, and the current implications for the history of galaxy formation as a function of cosmic epoch.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

Structure and function in rhodopsin: high level expression of a synthetic bovine opsin gene and its mutants in stable mammalian cell lines.

Stable mammalian cell lines harboring a synthetic bovine opsin gene have been derived from the suspension-adapted HEK293 cell line. The opsin gene is under the control of the immediate-early cytomegalovirus promoter/enhancer in an expression vector that also contains a selectable marker (Neo) governed by a relatively weak promoter. The cell lines expressing the opsin gene at high levels are selected by growth in the presence of high concentrations of the antibiotic geneticin. Under the conditions used for cell growth in suspension, opsin is produced at saturated culture levels of more than 2 mg/liter. After reconstitution with 11-cis-retinal, rhodopsin is purified to homogeneity in a single step by immunoaffinity column chromatography. Rhodopsin thus prepared (> 90% recovery at concentrations of up to 15 microM) is indistinguishable from rhodopsin purified from bovine rod outer segments by the following criteria: (i) UV/Vis absorption spectra in the dark and after photobleaching and the rate of metarhodopsin II decay, (ii) initial rates of transducin activation, and (iii) the rate of phosphorylation by rhodopsin kinase. Although mammalian cell opsin migrates slower than rod outer segment opsin on SDS/polyacrylamide gels, presumably due to a different N-glycosylation pattern, their mobilities after deglycosylation are identical. This method has enabled the preparation of several site-specific mutants of bovine opsin in comparable amounts.

Structure of the ripple phase in lecithin bilayers.

The phases of the x-ray form factors are derived for the ripple (Pbeta') thermodynamic phase in the lecithin bilayer system. By combining these phases with experimental intensity data, the electron density map of the ripple phase of dimyristoyl-phosphatidylcholine is constructed. The phases are derived by fitting the intensity data to two-dimensional electron density models, which are created by convolving an asymmetric triangular ripple profile with a transbilayer electron density profile. The robustness of the model method is indicated by the result that many different models of the transbilayer profile yield essentially the same phases, except for the weaker, purely ripple (0,k) peaks. Even with this residual ambiguity, the ripple profile is well determined, resulting in 19 angstroms for the ripple amplitude and 10 degrees and 26 degrees for the slopes of the major and the minor sides, respectively. Estimates for the bilayer head-head spacings show that the major side of the ripple is consistent with gel-like structure, and the minor side appears to be thinner with lower electron density.