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.

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.

The protein kinases of Caenorhabditis elegans: A model for signal transduction in multicellular organisms

Caenorhabditis elegans should soon be the first multicellular organism whose complete genomic sequence has been determined. This achievement provides a unique opportunity for a comprehensive assessment of the signal transduction molecules required for the existence of a multicellular animal. Although the worm C. elegans may not much resemble humans, the molecules that regulate signal transduction in these two organisms prove to be quite similar. We focus here on the content and diversity of protein kinases present in worms, together with an assessment of other classes of proteins that regulate protein phosphorylation. By systematic analysis of the 19,099 predicted C. elegans proteins, and thorough analysis of the finished and unfinished genomic sequences, we have identified 411 full length protein kinases and 21 partial kinase fragments. We also describe 82 additional proteins that are predicted to be structurally similar to conventional protein kinases even though they share minimal primary sequence identity. Finally, the richness of phosphorylation-dependent signaling pathways in worms is further supported with the identification of 185 protein phosphatases and 128 phosphoprotein-binding domains (SH2, PTB, STYX, SBF, 14-3-3, FHA, and WW) in the worm genome.

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.

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.

A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

Structural Analysis and Molecular Model of a Self-Incompatibility RNase from Wild Tomato1

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Crystal structure at 2.6-A resolution of human macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) was the first cytokine to be described, but for 30 years its role in the immune response remained enigmatic. In recent studies, MIF has been found to be a novel pituitary hormone and the first protein identified to be released from immune cells on glucocorticoid stimulation. Once secreted, MIF counterregulates the immunosuppressive effects of steroids and thus acts as a critical component of the immune system to control both local and systemic immune responses. We report herein the x-ray crystal structure of human MIF to 2.6 angstrom resolution. The protein is a trimer of identical subunits. Each monomer contains two antiparallel alpha-helices that pack against a four-stranded beta-sheet. The monomer has an additional two beta-strands that interact with the beta-sheets of adjacent subunits to form the interface between monomers. The three beta-sheets are arranged to form a barrel containing a solvent-accessible channel that runs through the center of the protein along a molecular 3-fold axis. Electrostatic potential maps reveal that the channel has a positive potential, suggesting that it binds negatively charged molecules. The elucidated structure for MIF is unique among cytokines or hormonal mediators, and suggests that this counterregulator of glucocorticoid action participates in novel ligand-receptor interactions.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.

On the reaction mechanism of l-lactate oxidase: Quantitative structure-activity analysis of the reaction with para-substituted l-mandelates

The rate constants for reduction of the flavoenzyme, l-lactate oxidase, and a mutant (in which alanine 95 is replaced by glycine), by a series of para-substituted mandelates, in both the 2-1H- and 2-2H- forms, have been measured by rapid reaction spectrophotometry. In all cases, significant isotope effects (1H/2H = 3–7) on the rate constants of flavin reduction were found, indicating that flavin reduction is a direct measure of α-C-H bond breakage. The rate constants show only a small influence of the electronic characteristics of the substituents, but show a good correlation when combined with some substituent volume parameters. A surprisingly good correlation is found with the molecular mass of the substrate. The results are compatible with any mechanism in which there is little development of charge in the transition state. This could be a transfer of hydride to the flavin N(5) position or a synchronous mechanism in which the α-C-H is formally abstracted as a H+ while the resulting charge is simultaneously neutralized by another event.

A link between protein structure and enzyme catalyzed hydrogen tunneling

We present evidence that the size of an active site side chain may modulate the degree of hydrogen tunneling in an enzyme-catalyzed reaction. Primary and secondary kH/kT and kD/kT kinetic isotope effects have been measured for the oxidation of benzyl alcohol catalyzed by horse liver alcohol dehydrogenase at 25°C. As reported in earlier studies, the relationship between secondary kH/kT and kD/kT isotope effects provides a sensitive probe for deviations from classical behavior. In the present work, catalytic efficiency and the extent of hydrogen tunneling have been correlated for the alcohol dehydrogenase-catalyzed hydride transfer among a group of site-directed mutants at position 203. Val-203 interacts with the opposite face of the cofactor NAD+ from the alcohol substrate. The reduction in size of this residue is correlated with diminished tunneling and a two orders of magnitude decrease in catalytic efficiency. Comparison of the x-ray crystal structures of a ternary complex of a high-tunneling (Phe-93 → Trp) and a low-tunneling (Val-203 → Ala) mutant provides a structural basis for the observed effects, demonstrating an increase in the hydrogen transfer distance for the low-tunneling mutant. The Val-203 → Ala ternary complex crystal structure also shows a hyperclosed interdomain geometry relative to the wild-type and the Phe-93 → Trp mutant ternary complex structures. This demonstrates a flexibility in interdomain movement that could potentially narrow the distance between the donor and acceptor carbons in the native enzyme and may enhance the role of tunneling in the hydride transfer reaction.

Cerebral protein synthesis in a genetic mouse model of phenylketonuria

Local rates of cerebral protein synthesis (lCPSleu) were measured with the quantitative autoradiographic [1-14C]leucine method in a genetic mouse model (Pahenu2) of phenylketonuria. As in the human disease, Pahenu2 mice have a mutation in the gene for phenylalanine hydroxylase. We compared adult homozygous (HMZ) and heterozygous (HTZ) Pahenu2 mice with the background strain (BTBR). Arterial plasma concentrations of phenylalanine (Phe) were elevated in both HMZ and HTZ mutants by 21 times and 38%, respectively. In the total acid-soluble pool in brain concentrations of Phe were higher and other neutral amino acids lower in HMZ mice compared with either HTZ or BTBR mice indicating a partial saturation of the l-amino acid carrier at the blood brain barrier by the elevated plasma Phe concentrations. In a series of steady-state experiments, the contribution of leucine from the arterial plasma to the tRNA-bound pool in brain was found to be statistically significantly reduced in HMZ mice compared with the other groups, indicating that a greater fraction of leucine in the precursor pool for protein synthesis is derived from protein degradation. We found reductions in lCPSleu of about 20% throughout the brain in the HMZ mice compared with the other two groups, but no reductions in brain concentrations of tRNA-bound neutral amino acids. Our results in the mouse model suggest that in untreated phenylketonuria in adults, the partial saturation of the l-amino acid transporter at the blood–brain barrier may not underlie a reduction in cerebral protein synthesis.

Mutant G protein α subunit activated by Gβγ: A model for receptor activation?

How receptors catalyze exchange of GTP for GDP bound to the Gα subunit of trimeric G proteins is not known. One proposal is that the receptor uses the G protein's βγ heterodimer as a lever, tilting it to pull open the guanine nucleotide binding pocket of Gα. To test this possibility, we designed a mutant Gα that would bind to βγ in the tilted conformation. To do so, we excised a helical turn (four residues) from the N-terminal region of αs, the α subunit of GS, the stimulatory regulator of adenylyl cyclase. In the presence, but not in the absence, of transiently expressed β1 and γ2, this mutant (αsΔ), markedly stimulated cAMP accumulation. This effect depended on the ability of the coexpressed β protein to interact normally with the lip of the nucleotide binding pocket of αsΔ. We substituted alanine for an aspartate in β1 that binds to a lysine (K206) in the lip of the α subunit's nucleotide binding pocket. Coexpressed with αsΔ and γ2, this mutant, β1-D228A, elevated cAMP much less than did β1-wild type; it did bind to αsΔ normally, however, as indicated by its unimpaired ability to target αsΔ to the plasma membrane. We conclude that βγ can activate αs and that this effect probably involves both a tilt of βγ relative to αs and interaction of β with the lip of the nucleotide binding pocket. We speculate that receptors use a similar mechanism to activate trimeric G proteins.

Cloning of a Tobacco Apoplasmic Invertase Inhibitor1 : Proof of Function of the Recombinant Protein and Expression Analysis during Plant Development

Higher plants express several isoforms of vacuolar and cell wall invertases (CWI), some of which are inactivated by inhibitory proteins at certain stages of plant development. We have purified an apoplasmic inhibitor (INH) of tobacco (Nicotiana tabacum) CWI to homogeneity. Based on sequences from tryptic fragments, we have isolated a full-length INH-encoding cDNA clone (Nt-inh1) via a reverse transcriptase-polymerase chain reaction. Southern-blot analysis revealed that INH is encoded by a single- or low-copy gene. Comparison with expressed sequence tag clones from Arabidopsis thaliana and Citrus unshiu indicated the presence of Nt-inh1-related proteins in other plants. The recombinant Nt-inh1-encoded protein inhibits CWI from tobacco and Chenopodium rubrum suspension-cultured cells and vacuolar invertase from tomato (Lycopersicon esculentum) fruit, whereas yeast invertase is not affected. However, only in the homologous system is the inhibition modulated by the concentration of Suc as previously shown for INH isolated from tobacco cells. Highly specific binding of INH to CWI could be shown by affinity chromatography of a total cell wall protein fraction on immobilized recombinant Nt-inh1 protein. RNA-blot analysis of relative transcript ratios for Nt-inh1 and CWI in different parts of adult tobacco plants revealed that the expression of both proteins is not always coordinate.