Showing papers in "Proteins in 1990"

Automated docking of substrates to proteins by simulated annealing.

Abstract: The Metropolis technique of conformation searching is combined with rapid energy evaluation using molecular affinity potentials to give an efficient procedure for docking substrates to macromolecules of known structure. The procedure works well on a number of crystallographic test systems, functionally reproducing the observed binding modes of several substrates.

An expectation maximization (EM) algorithm for the identification and characterization of common sites in unaligned biopolymer sequences.

Abstract: Statistical methodology for the identification and characterization of protein binding sites in a set of unaligned DNA fragments is presented. Each sequence must contain at least one common site. No alignment of the sites is required. Instead, the uncertainty in the location of the sites is handled by employing the missing information principle to develop an "expectation maximization" (EM) algorithm. This approach allows for the simultaneous identification of the sites and characterization of the binding motifs. The reliability of the algorithm increases with the number of fragments, but the computations increase only linearly. The method is illustrated with an example, using known cyclic adenosine monophosphate receptor protein (CRP) binding sites. The final motif is utilized in a search for undiscovered CRP binding sites.

Hormone phage: an enrichment method for variant proteins with altered binding properties.

Abstract: Human growth hormone (hGH), a 191 residue protein containing two disulfide bonds, was fused to the carboxyl-terminal domain of the gene III protein, a minor coat protein exposed at one end of the filamentous phage M13. The gene fusion was cloned into a plasmid containing origins of replication for Escherichia coli and filamentous phage and was packaged into phagemid particles upon infection by an M13KO7 helper phage. Transcription of the hGH-gene III fusion was controlled so that usually no more than one copy of the fusion protein was displayed along with the four copies of the wild-type gene III protein. The hGH-gene III fusion protein was properly folded, as judged by reactivity with six hGH monoclonal antibodies whose epitopes are sensitive to the folded conformation of hGH. Moreover, the hGH-gene III phagemid particles were enriched over 5000-fold from non-hGH phage, and 8-fold from a mutant hGH phagemid following a single hGH-specific elution step from hGH receptor-coated beads. The hGH phagemid should be useful for isolating new receptor binding mutants of hGH. More generally, this expression system may allow other large proteins with discontinuous binding epitopes to be displayed, and binding selections applied to their mutated gene III fusions on filamentous phage.

Comparative modeling methods: Application to the family of the mammalian serine proteases

Abstract: Comparative modeling methods are described that can be used to construct a three-dimensional model structure of a new protein from knowledge of its sequence and of the experimental structures and sequences of other members of its homology family. The methods are illustrated with the mammalian serine protease family, for which seven experimental structures have been reported in the literature, and the sequences for over 35 different protein members of the family are available. The strategy for modeling these proteins is presented, and criteria are developed for determining and assigning the reliability of the modeled structure. Criteria are described that are specially designed to help detect cases in which it is likely that the local structure diverges significantly from the usual conformation of the family.

Crystallographic refinement of human serum retinol binding protein at 2A resolution.

Abstract: Human serum retinol binding protein (RBP) in complex with retinol has been crystallographically refined to an R-factor of 181% with 2A resolution data The protein topology results in an anti-parallel beta-barrel that encapsulates the retinol ligand A detailed description of the protein and the binding site is provided Our structural work has helped to define a family of proteins, many of which are carrier proteins for smaller ligand molecules We describe the structural basis for the conservation of sequence within the family

On the nature of antibody combining sites: unusual structural features that may confer on these sites an enhanced capacity for binding ligands

Abstract: A detailed analysis of the structural aspects of antibody-antigen interactions has been made possible by the availability of X-ray structures for three complexes of antilysozyme Fabs to lysozyme (reviewed by Davies et al.: J. Biol. Chem. 263:10541-10544, 1988.) Examination of the antigen-contacting residues in the three antilysozyme Fabs reveals the occurrence of a large number of aromatics, particularly tyrosines, and the absence of apolar, aliphatic residues. Calculation of the frequency of occurrence of the various amino acid types reveals that tyrosines are three times, and histidines and asparagines eight times, more likely to be found in the complementarity-determining regions than in the framework of the variable domains. Analysis of the solvent accessibility of the residues in Fvs (the modules containing variable domains of the light and heavy chains) of known three-dimensional structure indicates that tyrosines and tryptophans are more exposed when they occur in the complementarity-determining regions than when in the framework. Furthermore, many more of the asparagines in the complementarity-determining regions than in the framework are buried. These asparagines appear to have a structural role in that they hydrogen-bond through their side chains to other side chains and, even more so, to the protein backbone. The stabilizing effect of the asparagines, plus the rigidity of the framework, may serve to allow the greater exposure of the aromatic residues to solvent. In view of the greater potential contribution of aromatic side chains to the total binding energy, these results suggest that antibody combining sites have structural features that make them especially suited for interacting with ligands.

Crystallographic studies of inhibitor binding sites in human carbonic anhydrase II: a pentacoordinated binding of the SCN- ion to the zinc at high pH.

Abstract: The binding of four inhibitors--mercuric ion, 3-acetoxymercuri-4-aminobenzenesulfonamide (AMS), acetazolamide (Diamox), and thiocyanate ion--to human carbonic anhydrase II (HCA II) has been studied with X-ray crystallography The binding of mercury to HCA II at pH 70 has been investigated at 31 A resolution Mercuric ions are observed at both nitrogens in the His-64 ring One of these sites is pointing toward the zinc ion The only other binding site for mercury is at Cys-206 The binding of the two sulfonamide inhibitors AMS and Diamox, has been reinvestigated at 20 and 30 A, respectively Only the nitrogen of the sulfonamide group binds to the zinc ion replacing the hydroxyl ion The sulfonamide oxygen closest to the zinc ion is 31 A away Thus the tetrahedral geometry of the zinc is retained, refuting earlier models of a pentacoordinated zinc The structure of the thiocyanate complex has been investigated at pH 85 and the structure has been refined at 19 A resolution using the least-squares refinement program PROLSQ The crystallographic R factor is 176% The zinc ion is pentacoordinated with the anion as well as a water molecule bound in addition to the three histidine residues The nitrogen atom of the SCN- ion is 19 A from the zinc ion but shifted 13 A with respect to the hydroxyl ion in the native structure and at van der Waals' distance from the O gamma l atom of Thr-199 This is due to the inability of the O gamma l atom of Thr-199 to serve as a hydrogen bond donor, thus repelling the nonprotonated nitrogen The SCN- molecule reaches into the deep end of the active site cavity where the sulfur atom has displaced the so-called "deep" water molecule of the native enzyme The zinc-bound water molecule is 22 A from the zinc ion and 24 A from the SCN- nitrogen In addition, this water is hydrogen bonded to the O gamma l atom of Thr-199 and to another water molecule We have observed that solvent and inhibitor molecules have three possible binding sites on the zinc ion and their significance for the catalysis and inhibition of HCA II will be discussed All available crystallographic data are consistent with a proposed catalytic mechanism in which both the OH moiety and one oxygen of the substrate HCO3- ion are ligated to the zinc ion

Beta-lactamase of Bacillus licheniformis 749/C at 2 A resolution.

Abstract: Two crystal forms (A and B) of the 29,500 Da Class A beta-lactamase (penicillinase) from Bacillus licheniformis 749/C have been examined crystallographically. The structure of B-form crystals has been solved to 2 A resolution, the starting model for which was a 3.5 A structure obtained from A-form crystals. The beta-lactamase has an alpha + beta structure with 11 helices and 5 beta-strands seen also in a penicillin target DD-peptidase of Streptomyces R61. Atomic parameters of the two molecules in the asymmetric unit were refined by simulated annealing at 2.0 A resolution. The R factor is 0.208 for the 27,330 data greater than 3 sigma (F), with water molecules excluded from the model. The catalytic Ser-70 is at the N-terminus of a helix and is within hydrogen bonding distance of conserved Lys-73. Also interacting with the Lys-73 are Asn-132 and the conserved Glu-166, which is on a potentially flexible helix-containing loop. The structure suggests the binding of beta-lactam substrates is facilitated by interactions with Lys-234, Thr-235, and Ala-237 in a conserved beta-strand peptide, which is antiparallel to the beta-lactam's acylamido linkage; an exposed cavity near Asn-170 exists for acylamido substituents. The reactive double bond of clavulanate-type inhibitors may interact with Arg-244 on the fourth beta-strand. A very similar binding site architecture is seen in the DD-peptidase.

Hydrophobicity of amino acid subgroups in proteins

Abstract: Protein folding studies often utilize areas and volumes to assess the hydrophobic contribution to conformational free energy (Richards, F.M. Annu. Rev. Biophys. Bioeng. 6:151-176, 1977). We have calculated the mean area buried upon folding for every chemical group in each residue within a set of X-ray elucidated proteins. These measurements, together with a standard state cavity size for each group, are documented in a table. It is observed that, on average, each type of group buries a constant fraction of its standard state area. The mean area buried by most, though not all, groups can be closely approximated by summing contributions from three characteristic parameters corresponding to three atom types: (1) carbon or sulfur, which turn out to be 86% buried, on average; (2) neutral oxygen or nitrogen, which are 40% buried, on average; and (3) charged oxygen or nitrogen, which are 32% buried, on average.

Refinement of the NMR structures for acyl carrier protein with scalar coupling data.

Abstract: Structure determination of small proteins using NMR data is most commonly pursued by combining NOE derived distance constraints with inherent constraints based on chemical bonding. Ideally, one would make use of a variety of experimental observations, not just distance constraints. Here, coupling constant constraints have been added to molecular mechanics and molecular dynamics protocols for structure determination in the form of a psuedoenergy function that is minimized in a search for an optimum molecular conformation. Application is made to refinement of a structure for a 77 amino acid protein involved in fatty acid synthesis, Escherichia coli acyl carrier protein (ACP). 54 3JHN alpha coupling constants, 12 coupling constants for stereospecifically assigned side chain protons, and 450 NOE distance constraints were used to calculate the 3-D structure of ACP. A three-step protocol for a molecular dynamics calculation is described, in analogy to the protocol previously used in molecular mechanics calculations. The structures calculated with the molecular mechanics approach and the molecular dynamics approach using a rigid model for the protein show similar molecular energies and similar agreement with experimental distance and coupling constant constraints. The molecular dynamics approach shows some advantage in overcoming local minimum problems, but only when a two-state averaging model for the protein was used, did molecular energies drop significantly.

Plastic adaptation toward mutations in proteins: Structural comparison of thymidylate synthases

Abstract: The structure of thymidylate synthase (TS) from Escherichia coli was solved from cubic crystals with a = 133 A grown under reducing conditions at pH 7.0, and refined to R = 22% at 2.1 A resolution. The structure is compared with that from Lactobacillus casei solved to R = 21% at 2.3 A resolution. The structures are compared using a difference distance matrix, which identifies a common core of residues that retains the same relationship to one another in both species. After subtraction of the effects of a 50 amino acid insert present in Lactobacillus casei, differences in position of atoms correlate with temperature factors and with distance from the nearest substituted residue. The dependence of structural difference on thermal factor is parameterized and reflects both errors in coordinates that correlate with thermal factor, and the increased width of the energy well in which atoms of high thermal factor lie. The dependence of structural difference on distance from the nearest substitution also depends on thermal factors and shows an exponential dependence with half maximal effect at 3.0 A from the substitution. This represents the plastic accommodation of the protein which is parameterized in terms of thermal B factor and distance from a mutational change.

The three‐dimensional structure of recombinant bovine chymosin at 2.3 Å resolution

Abstract: The crystal structure of recombinant bovine chymosin (EC 3.4.23.4; renin), which was cloned and expressed in Escherichia coli, has been determined using X-ray data extending to 2.3 A resolution. The crystals of the enzyme used in this study belong to the space group I222 with unit cell dimensions alpha = 72.7 A, b = 80.3 A, and c = 114.8 A. The structure was solved by the molecular replacement method and was refined by a restrained least-squares procedure. The crystallographic R factor is 0.165 and the deviation of bond distances from ideality is 0.020 A. The resulting model includes all 323 amino acid residues, as well as 297 water molecules. The enzyme has an irregular shape with approximate maximum dimensions of 40 x 50 x 65 A. The secondary structure consists primarily of parallel and antiparallel beta-strands with a few short alpha-helices. The enzyme can be subdivided into N- and C-terminal domains which are separated by a deep cleft containing the active aspartate residues Asp-34 and Asp-216. The amino acid residues and waters at the active site form an extensive hydrogen-bonded network which maintains the pseudo 2-fold symmetry of the entire structure. A comparison of recombinant chymosin with other acid proteinases reveals the high degree of structural similarity with other members of this family of proteins as well as the subtle differences which make chymosin unique. In particular, Tyr-77 of the flap region of chymosin does not hydrogen bond to Trp-42 but protrudes out in the P1 pocket forming hydrophobic interactions with Phe-119 and Leu-32. This may have important implications concerning the mechanism of substrate binding and substrate specificity.

Revised 2.3 A structure of porcine pepsin: evidence for a flexible subdomain

Abstract: A revised three-dimensional crystal structure of ethanol-inhibited porcine pepsin refined to an R-factor of 0.171 at 2.3 A resolution is presented and compared to the refined structures of the fungal aspartic proteinases: penicillopepsin, rhizopuspepsin, and endothiapepsin. Pepsin is composed of two nearly equal N and C domains related by an intra dyad. The overall polypeptide fold and active site structures are homologous for pepsin and the fungal enzymes. The weak inhibition of pepsin by ethanol can be explained by the presence of one or more ethanol molecules, in the vicinity of the active site carboxylates, which slightly alter the hydrogen-bonding network and which may compete with substrate binding in the active site. Structural superposition analysis showed that the N domains aligned better than the C-domains for pepsin and the fungal aspartic proteinases: 107-140 C alpha pairs aligned to 0.72-0.85 A rms for the N domains; 64-95 C alpha pairs aligned to 0.78-1.03 A rms for the C domains. The major structural difference between pepsin and the fungal enzymes concerns a newly described subdomain whose conformation varies markedly among these enzyme structures. The subdomain in pepsin comprises nearly 100 residues and is composed of two contiguous segments within the C domain (residues 192-212 and 223-299). the subdomain is connected, or "hinged," to a mixed beta-sheet that forms one of the structurally invariant, active site psi-loops. Relative subdomain displacements as large as a 21.0 degrees rotation and a 5.9 A translation were observed among the different enzymes. There is some suggestion in pepsin that the subdomain may be flexible and perhaps plays a structural role in mediating substrate binding, determining the substrate specificity, or in the activation of the zymogen.

Identification of protein folds: matching hydrophobicity patterns of sequence sets with solvent accessibility patterns of known structures.

Abstract: Hydrophobic side chains often are buried in the interior of a protein, and evolutionarily related proteins usually maintain the hydrophobic character of buried positions. In this paper we show that a pattern of hydrophobicity values derived from a set of related protein sequences is well correlated with the linear pattern of side-chain solvent accessibility values, calculated from a known protein structure representative of the sequences. In several cases, information from aligned sequences can be used to select the correct tertiary fold from a large data base of protein structures.

Understanding structural relationships in proteins of unsolved three-dimensional structure.

Abstract: The locations of functionally important sequences and general structural motifs have been assigned to Ile-tRNA synthetase. However, a function has not been established for some segments of the protein (e.g., CP1). The method of structural modeling described here cannot establish the details of a 3 A crystal structure, and, in contrast to a crystal structure, the precision of the model varies according to the extent of a sequence similarity or the functional importance of a region. In Ile-tRNA synthetase, the signature sequence and the flanking regions are likely to be similar in structure to the proteins on which the model is based. For other regions, it may be possible to build a three-dimensional model by connecting well defined regions and refining the positions of the connecting elements by energy minimization. Structural modelling of this kind must be done cautiously, because the order and orientation of the elements of a structural motif can change in subtle ways. In the case of Tyr-tRNA synthetase, the beta-strand nearest the N-terminus is the outermost strand of the nucleotide binding fold; in Met-tRNA synthetase, the same strand is innermost. Furthermore, the orientation of this strand may be antiparallel (Tyr-tRNA synthetase) or parallel (Met-tRNA synthetase). Because multiple structures that differ in their orientations of structural elements are possible, the structural analogies between proteins should not be naively extrapolated without independent experimental support. As described above, some regions of proteins tolerate internal deletions and insertions. This provides further experimental support for the practice of allowing for gaps in computer-generated sequence alignments. Nevertheless, because some regions are more tolerant of insertions and deletions than others, the structural and functional significance of a region of broken alignment must be assessed carefully. All gaps in sequence alignments cannot be treated equally, and each must be evaluated within its own context. In the synthetases of known structure, structural analogy can be used to identify important functional elements. For example, the amino acid binding site of Met-tRNA synthetase might be formed, at least in part, by a peptide that encompasses Ala50; this amino acid aligns with Gly94 of the Ile-tRNA synthetase. This is an example in which results on a protein of unknown structure (Ile-tRNA synthetases) can lead to identification of a potential substrate binding site in a protein of known structure (Met-tRNA synthetase).

Solution structure of a zinc finger domain of yeast ADR1

Abstract: The "zinc finger" is a 30-residue repeating motif that has been identified in a variety of eukaryotic transcription factors. Each domain is capable of binding a Zn2+ ion through invariant Cys and His residues. We have determined the three-dimensional structure of a synthetic peptide that corresponds to one of the two zinc finger domains in the yeast transcription factor ADR1, using two-dimensional nuclear magnetic resonance spectroscopy. The Zn2(+)-bound structure of the peptide consists of a loop containing the two Cys residues, a "fingertip," a 12- to 13-residue alpha-helix containing the two His residues, and a C-terminal tail. A majority of the interresidue contacts observed involve the seven conserved residues of the prototypic zinc finger (i.e., the four zinc ligands and the three hydrophobic residues), indicating that these residues are largely responsible for the three-dimensional structure of the domain and that all the zinc finger domains of the TFIIIA class will have similar structures. Potential DNA-binding residues are found throughout the structure, with the highest concentration of such residues on the external face of the alpha-helix.

Crystal structure of myoglobin from a synthetic gene.

Abstract: Crystals have been grown of myoglobin produced in Escherichia coli from a synthetic gene, and the structure has been solved to 1.9 A resolution. The space group of the crystals is P6, which is different from previously solved myoglobin crystal forms. The synthetic myoglobin is essentially identical to myoglobin isolated from sperm whale tissue, except for the retention of the initiator methionine at the N-terminus and the substitution of asparagine for aspartic acid at position 122. Superposition of the coordinates of native and synthetic sperm whale myoglobins reveals only minor changes in the positions of main chain atoms and reorientation of some surface side chains. Crystals of variants of the "synthetic" myoglobin have also been grown for structural analysis of the role of key amino acid residues in ligand binding and specificity.

Polyproline, β-turn helices. Novel secondary structures proposed for the tandem repeats within rhodopsin, synaptophysin, synexin, gliadin, RNA polymerase II, hordein, and gluten

Abstract: Seven proteins each contain 8 to 52 tandem repeats of a unique class of oligopeptide. The consensus peptide for each is rhodopsin Tyr Pro Pro Gln Gly synaptophysin Tyr Gly Pro Gln Gly synexin Tyr Pro Pro Pro Pro Gly gliadin Tyr Pro Pro Pro Gln Pro RNA polymerase II Tyr Ser Pro Thr Ser Pro Ser hordein Phe Pro Gln Gln Pro Gln Gln Pro gluten Tyr Pro Thr Ser Pro Gln Gln Gly Tyr Although there is obvious variation of sequence and of length, the penta- to nonapeptides share an initial Tyr (or Phe) and have high Pro contents and abundant Gly, Gln, and Ser. We have evaluated helical models that both recognize the uniqueness of these sequence repeats and accommodate variations on the basic theme. We have developed a group of related helical models for these proteins with about three oligopeptide repeats per turn of 10-20 A. These models share several common features: Most of the phi dihedral angles are -54 degrees, to accommodate Pro at all positions except the first (Tyr). Except for the beta-turns, most psi dihedral angles are near +140 degrees as found in polyproline. Each oligopeptide has at least one beta-turn; several have two. Some contain a cis-Tyr, Pro peptide bond; a few have a cis-bond plus one beta-turn. Tyr side chains vary from totally exposed to buried within the helices and could move to accommodate either external hydrophobic interactions or phosphorylation. The several related structures seem to be readily interconverted without major change in the overall helical parameters, and therein may lie the key to their functions.

HERA—A program to draw schematic diagrams of protein secondary structures

Abstract: A program is described which generates hydrogen bonding diagrams of protein structures and optionally helical wheels and helical nets. The program can also be used simply to calculate the connectivities of β-strands and to automatically extract simple structural motifs such as hairpins or Greek keys. The program greatly reduces the effort required to produce these diagrams and offers considerable flexibility in the information which can be represented. The usefulness of the program is illustrated by several examples including comparing homologous families, correlating protein structure with attributes of individual residues, and extracting all examples of the ψ-loop motif from the Brookhaven Data Bank.

Comparison of the structures of globins and phycocyanins: evidence for evolutionary relationship.

Abstract: Globins and phycocyanins are two classes of proteins with different function, different ligands, and no substantial sequence similarity, yet the conformations of their polypeptide chains show very similar folding patterns. Does this arise from a genuine, albeit very distant, evolutionary relationship, or does it represent a common solution of a structural problem? We address this question by a very detailed comparison of the structures of the two protein families. An analysis of the helices and their interactions shows many features common to globins and phycocyanins, including some exceptional features of the globins such as a 3–10 C helix and the unusual “crossed-ridge” packing pattern at the B/E helix interfaces. We conclude that the evidence supports the hypothesis of distant evolutionary relationship between globins and phycocyanins.

Protein stability and electrostatic interactions between solvent exposed charged side chains

Abstract: To investigate the contribution to protein stability of electrostatic interactions between charged surface residues, we have studied the effect of substituting three negatively charged solvent exposed residues with their side-chain amide analogs in bovine calbindin D9k--a small (Mr 8,500) globular protein of the calmodulin superfamily. The free energy of urea-induced unfolding for the wild-type and seven mutant proteins has been measured. The mutant proteins have increased stability towards unfolding relative to the wild-type. The experimental results correlate reasonably well with theoretically calculated relative free energies of unfolding and show that electrostatic interactions between charges on the surface of a protein can have significant effects on protein stability.

Normal mode analysis of human lysozyme: study of the relative motion of the two domains and characterization of the harmonic motion.

Abstract: A normal mode analysis of human lysozyme has been carried out at room temperature. Human lysozyme is an enzyme constituted of two domains separated by an active site cleft, the motion of which is thought to be relevant for biological function. This motion has been described as a hinge bending motion. McCammon et al. have determined the characteristics of the hinge bending motion but they assumed a prior knowledge of the hinge axis. In this work we propose a method which is free from this assumption and determines the hinge axis and root mean square (rms) rotation angle which give the best agreement with the pattern of changes in all the distances between nonhydrogen atoms in the two domains, obtained by the normal mode analysis. The hinge axis we found is notably different from the one previously determined and goes, roughly, through the C alpha 55 and C alpha 76, i.e., it is located at the base of the beta-sheet of the second domain. The rms value for the rotation angle is also twice as large as the previous one: 3.37 degrees. It is shown that this hinge bending motion provides a fairly good approximation of the dynamics of human lysozyme and that the normal mode with the lowest frequency has a dominating contribution to this hinge bending motion. A study of the accessible surface area of the residues within the cleft reveals that the motion does not result in a better exposure to the solvent of these residues. A characterization of the thermally excited state (under the hypothesis of the harmonicity of the potential energy surface) has been done using the concept of topology of atom packing. Under this hypothesis the thermal fluctuations result only in a small change of the topology of atom packing, leading therefore to nearly elastic deformations of the protein.

Comparative molecular model building of two serine proteinases from cytotoxic T lymphocytes

Abstract: Two genes that are expressed when precursor cytotoxic T lymphocytes are transformed to T killer cells have been cloned and sequenced. The derived amino acid sequences, coding for cytotoxic cell protease 1 (CCP1) and Hannuka factor (HF) are highly homologous to members of the serine proteinase family. Comparative molecular model building using the known three-dimensional structures and the derived amino acid sequences of the lymphocyte enzymes has provided useful structural information, especially in predicting the conformations of the substrate binding sites. In applying this modelling procedure, we used the X-ray structures of four serine proteinases to provide a structurally based sequence alignment: alpha-chymotrypsin (CHT), bovine trypsin (BT), Streptomyces griseus trypsin (SGT), and rat mast cell protease 2 (RMCP2). The root mean square differences in alpha-carbon atom positions among these four structures when compared in a pairwise fashion range from 0.79 to 0.97 A for structurally equivalent residues. The sequences of the two lymphocyte enzymes were then aligned to these proteinases using chemical criteria and the superimposed X-ray structures as guides. The alignment showed that the sequence of CCP1 was most similar to RMCP2, whereas HF has regions of homology with both RMCP2 and BT. With RMCP2 as a template for CCP1 and the two enzymes RMCP2 and BT as templates for HF, the molecular models were constructed. Intramolecular steric clashes that resulted from the replacement of amino acid side chains of the templates by the aligned residues of CCP1 and HF were relieved by adjustment of the side chain conformational angles in an interactive computer graphics device. This process was followed by energy minimization of the enzyme model to optimize the stereochemical geometry and to relieve any remaining unacceptably close nonbonded contacts. The resulting model of CCP1 has an arginine residue at position 226 in the specificity pocket, thereby predicting a substrate preference for P1 aspartate or glutamate residues. The model also predicts favorable binding for a small hydrophobic residue at the P2 position of the substrate. The primary specificity pocket of HF resembles that of BT and therefore predicts a lysine or arginine preference for the P1 residue. The arginine at position 99 in the model of HF suggests a preference for aspartate or glutamate side chains in the P2 position of the substrate. Both CCP1 and HF have a free cysteine in the segment of polypeptide 88 to 93.(ABSTRACT TRUNCATED AT 400 WORDS)

Functional interaction among catalytic residues in subtilisin BPN

Abstract: Variants of the serine protease, subtilisin BPN', in which the catalytic triad residues (Ser-221, His-64, and Asp-32) are replaced singly or in combination by alanine retain activities with the substrate N-succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide (sAAPF-pna) that are at least 10(3) to 10(4) above the non-enzymatic rate [Carter, P., Wells, J.A. Nature (London) 322:564-568, 1988]. A possible source of the residual activity was the hydrogen bond with the N delta 2 of Asn-155 that helps to stabilize the oxyanion generated in the tetrahedral transition state during amide bond hydrolysis by the wild-type enzyme. Replacing Asn-155 by Gly (N155G) lowers the turnover number (kcat) for sAAPF-pna by 150-fold with virtually no change in the Michaelis constant (KM). However, upon combining the N155G and S221A mutations to give N155G:S221A, kcat is actually 5-fold greater than for the S221A enzyme. Thus, the catalytic role of Asn-155 is dependent upon the presence of Ser-221. The residual activity of the N155G:S221A enzyme (approximately 10(4)-fold above the uncatalyzed rate) is not an artifact because it can be completely inhibited by the third domain of the turkey ovomucoid inhibitor (OMTKY3), which forms a strong 1:1 complex with the active site. The mutations N155G and S221A individually weaken the interaction between subtilisin and OMTKY3 by 1.8 and 2.0 kcal/mol, respectively, and in combination by 2.1 kcal/mol. This is consistent with disruption of stabilizing interactions around the reactive site carbonyl of the OMTKY3 inhibitor.(ABSTRACT TRUNCATED AT 250 WORDS)

A mutant T4 lysozyme (Val 131 → Ala) designed to increase thermostability by the reduction of strain within an α-helix

Abstract: An attempt has been made to identify residues in T4 phage lyoszyme that may have strained conformations and, by appropriate site-directed replacements, to reduce this strain and thus increase the thermostability of the protein. Valine 131, within α-helix 126–134, was identified as a potential candidate. Its side-chain rotational as a potential candidate. Its side-chain rotation angle, χ1, differs by approximately 18° from the low-energy trans configuration. In addition, it is largely solvent exposed, yet is held in a rigid conformation. The mutant protein with Val 131 replaced by alanine temperature 0.9°C higher than that of wildtype lyoszyme at pH 2.8. As a control, The mutant Val 131 → Thr was also constructed and its melting temperature was found to be marginally lower than wild type. Higher-resolution crystal structure determination of the mutant lysozymes show that their structure are virtually identical with that of wild-type lyoszyme, except for the Val → Ala or Val → Thr replacement. Analysis of the different structures suggests that the design of the Val→Ala substitution was, in principle, successful, although the apparent gain in stability caused by reduction in strain is modest and is somewhat offset by the loss of hydrophobic interactions and by entrophic effects. The results also help to provide a structural retionalization that alanine has a higher helix propensity than valine or theronine.

Functionally acceptable substitutions in two α-helical regions of λ repressor

Abstract: A method of targeted random mutagenesis has been used to investigate the informational content of 25 residue positions in two alpha-helical regions of the N-terminal domain of lambda repressor. Examination of the functionally allowed sequences indicates that there is a wide range in tolerance to amino acid substitution at these positions. At positions that are buried in the structure, there are severe limitations on the number and type of residues allowed. At most surface positions, many different residues and residue types are tolerated. However, at several surface positions there is a strong preference for hydrophilic amino acids, and at one surface position proline is absolutely conserved. The results reveal the high level of degeneracy in the information that specifies a particular protein fold.

Accommodation of single amino acid insertions by the native state of staphylococcal nuclease

Abstract: Single alanine and glycine insertions were introduced at 20 randomly selected positions in staphylococcal nuclease. The resulting changes in catalytic activity and in stability to guanidine hydrochloride denaturation indicate that the native state structure is frequently able to accommodate the extra residue without great difficulty, even insertions within secondary structural elements such as alpha helices and beta sheets. On average, an inserted residue reduces the free energy of denaturation (delta GH2O) by an amount roughly comparable to an alanine or glycine substitution for one of the residues flanking the site of insertion. Several positions outside of the enzyme active site were found where insertions, but not substitutions, lead to structural changes that modify catalytic activity and the circular dichroism spectrum. Amino acid insertions represent a virtually unexplored class of genetic mutation that may prove complementary to amino acid substitutions for engineering proteins with altered functional and structural properties.