Showing papers in "Proteins in 1992"

Stereochemical quality of protein structure coordinates.

Abstract: Methods have been developed to assess the stereochemical quality of any protein structure both globally and locally using various criteria. Several parameters can be derived from the coordinates of a given structure. Global parameters include the distribution of phi, psi and chi 1 torsion angles, and hydrogen bond energies. There are clear correlations between these parameters and resolution; as the resolution improves, the distribution of the parameters becomes more clustered. These features show a broad distribution about ideal values derived from high-resolution structures. Some structures have tightly clustered distributions even at relatively low resolutions, while others show abnormal scatter though the data go to high resolution. Additional indicators of local irregularity include proline phi angles, peptide bond planarities, disulfide bond lengths, and their chi 3 torsion angles. These stereochemical parameters have been used to generate measures of stereochemical quality which provide a simple guide as to the reliability of a structure, in addition to the most important measures, resolution and R-factor. The parameters used in this evaluation are not novel, and are easily calculated from structure coordinates. A program suite is currently being developed which will quickly check a given structure, highlighting unusual stereochemistry and possible errors.

Multiple protein sequence alignment from tertiary structure comparison: assignment of global and residue confidence levels.

Abstract: An algorithm is presented for the accurate and rapid generation of multiple protein sequence alignments from tertiary structure comparisons. A preliminary multiple sequence alignment is performed using sequence information, which then determines an initial superposition of the structures. A structure comparison algorithm is applied to all pairs of proteins in the superimposed set and a similarity tree calculated. Multiple sequence alignments are then generated by following the tree from the branches to the root. At each branchpoint of the tree, a structure-based sequence alignment and coordinate transformations are output, with the multiple alignment of all structures output at the root. The algorithm encoded in STAMP (STructural Alignment of Multiple Proteins) is shown to give alignments in good agreement with published structural accounts within the dehydrogenase fold domains, globins, and serine proteinases. In order to reduce the need for visual verification, two similarity indices are introduced to determine the quality of each generated structural alignment. Sc quantifies the global structural similarity between pairs or groups of proteins, whereas Pij' provides a normalized measure of the confidence in the alignment of each residue. STAMP alignments have the quality of each alignment characterized by Sc and Pij' values and thus provide a reproducible resource for studies of residue conservation within structural motifs.

The structure of the complex between influenza virus neuraminidase and sialic acid, the viral receptor

Abstract: Crystallographic studies of neuraminidase-sialic acid complexes indicate that sialic acid is distorted on binding the enzyme. Three arginine residues on the enzyme interact with the carboxylate group of the sugar which is observed to be equatorial to the saccharide ring as a consequence of its distorted geometry. The glycosidic oxygen is positioned within hydrogen-bonding distance of Asp-151, implicating this residue in catalysis.

Identification of novel peptide antagonists for GPIIb/IIIa from a conformationally constrained phage peptide library.

Abstract: Methods have recently been developed to present vast libraries of random peptides on the surface of filamentous phage. To introduce a degree of conformational constraint into random peptides, a library of hexapeptides flanked by cysteine residues (capable of forming cyclic disulfides) was constructed. This library was screened using the platelet glycoprotein, IIb/IIIa, which mediates the aggregation of platelets through binding of fibrinogen. A variety of peptides containing the sequence Arg-Gly-Asp or Lys-Gly-Asp were discovered and synthesized. The cyclic, disulfide bonded forms of the peptides bound IIb/IIIa with dissociation constants in the nanomolar range, while reduced forms or an analogue in which Ser replaced the Cys residues bound considerably less tightly. These results demonstrate the feasibility for introducing conformational constraints into random peptide libraries and also demonstrates the potential for using phage peptide libraries to discover pharmacologically active lead compounds. © 1992 Wiley-Liss, Inc.

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

Abstract: We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base andthe associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. Withone exception we find that for all globin sequences one of the known globinfolds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structurlresearch and future development of our approach. © 1992 Wiley-Liss, Inc.

Subunit assembly and active site location in the structure of glutamate dehydrogenase.

Abstract: The three-dimensional crystal structure of the NAD(+)-linked glutamate dehydrogenase from Clostridium symbiosum has been solved to 1.96 A resolution by a combination of isomorphous replacement and molecular averaging and refined to a conventional crystallographic R factor of 0.227. Each subunit in this multimeric enzyme is organised into two domains separated by a deep cleft. One domain directs the self-assembly of the molecule into a hexameric oligomer with 32 symmetry. The other domain is structurally similar to the classical dinucleotide binding fold but with the direction of one of the strands reversed. Difference Fourier analysis on the binary complex of the enzyme with NAD+ shows that the dinucleotide is bound in an extended conformation with the nicotinamide moiety deep in the cleft between the two domains. Hydrogen bonds between the carboxyamide group of the nicotinamide ring and the side chains of T209 and N240, residues conserved in all hexameric GDH sequences, provide a positive selection for the syn conformer of this ring. This results in a molecular arrangement in which the A face of the nicotinamide ring is buried against the enzyme surface and the B face is exposed, adjacent to a striking cluster of conserved residues including K89, K113, and K125. Modeling studies, correlated with chemical modification data, have implicated this region as the glutamate/2-oxoglutarate binding site and provide an explanation at the molecular level for the B type stereospecificity of the hydride transfer of GDH during the catalytic cycle.

Investigation of the function of mutated cellulose-binding domains of Trichoderma reesei cellobiohydrolase I.

Abstract: The function of the cellulose-binding domain (CBD) of the cellobiohydrolase I of Trichoderma reesei was studied by site-directed mutagenesis of two amino acid residues identified by analyzing the 3D structure of this domain. The mutant enzymes were produced in yeast and tested for binding and activity on crystalline cellulose. Mutagenesis of the tyrosine residue (Y492) located at the tip of the wedge-shaped domain to alanine or aspartate reduced the binding and activity on crystalline cellulose to the level of the core protein lacking the CBD. However, there was no effect on the activity toward small oligosaccharide (4-methylumbelliferyl beta-D-lactoside). The mutation tyrosine to histidine (Y492H) lowered but did not destroy the cellulose binding, suggesting that the interaction of the pyranose ring of the substrate with an aromatic side chain is important. However, the catalytic activity of this mutant on crystalline cellulose was identical to the other two mutants. The mutation P477R on the edge of the other face of the domain reduces both binding and activity of CBHI. These results support the hypothesis that both surfaces of the CBD are involved in the interaction of the binding domain with crystalline cellulose.

A multiple-start Monte Carlo docking method

Abstract: We present a method to search for possible binding modes of molecular fragments at a specific site of a potential drug target of known structure. Our method is based on a Monte Carlo (MC) algorithm applied to the translational and rotational degrees of freedom of the probe fragment. Starting from a randomly generated initial configuration, favorable binding modes are generated using a two-step process. An MC run is first performed in which the energy in the Metropolis algorithm is substituted by a score function that measures the average distance of the probe to the target surface. This has the effect of making buried probes move toward the target surface and also allows enhanced sampling of deep pockets. In a second MC run, a pairwise atom potential function is used, and the temperature parameter is slowly lowered during the run (Simulated Annealing). We repeat this procedure starting from a large number of different randomly generated initial configurations in order to find all energetically favorable docking modes in a specified region around the target. We test this method using two inhibitor-receptor systems: Streptomyces griseus proteinase B in complex with the third domain of the ovomucoid inhibitor from turkey, and dihydrofolate reductase from E. coli in complex with methotrexate. The method could consistently reproduce the complex found in the crystal structure searching from random initial positions in cubes ranging from 25 A to 50 A about the binding site. In the case of SGPB, we were also successful in docking to the native structure. In addition, we were successful in docking small probes in a search that included the entire protein surface.

A metal-mediated hydride shift mechanism for xylose isomerase based on the 1.6 A Streptomyces rubiginosus structures with xylitol and D-xylose.

Abstract: The crystal structure of recombinant Streptomyces rubiginousD-xylose isomerase (D-xylose keto-isomerase, EC 5.3.1.5) solved by the multiple isomorphous replacement technique has been refined to R = 0.16 at 1.64 A resolution. As observed in an earlier study at 4.0 A (Carrell et al., J. Biol. Chem. 259: 3230–3236, 1984), xylose isomerase is a tetramer composed of four identical subunits. The monomer consists of an eight-stranded parallel β-barrel surrounded by eight helices with an extended C-terminal tail that provides extensive contacts with a neighboring monomer. The active site pocket is defined by an opening in the barrel whose entrance is lined with hydrophobic residues while the bottom of the pocket consists mainly of glutamate, aspartate, and histidine residues coordinated to two manganese ions. The structures of the enzyme in the presence of MnCl2, the inhibitor xylitol, and the substrate D-xylose in the presence and absence of MnCl2 have also been refined to R = 0.14 at 1.60 A, R = 0.15 at 1.71 A, R = 0.15 at 1.60 A, and R = 0.14 at 1.60 A, respectively. Both the ring oxygen of the cyclic α-D-xylose and its C1 hydroxyl are within hydrogen bonding distance of NE2 of His-54 in the structure crystallized in the presence of D-xylose. Both the inhibitor, xylitol, and the extended form of the substrate, D-xylose, bind such that the C2 and C4 OH groups interact with one of the two divalent cations found in the active site and the C1 OH with the other cation. The remainder of the OH groups hydrogen bond with neighboring amino acid side chains. A detailed mechanism for D-xylose isomerase is proposed. Upon binding of cyclic α-D-xylose to xylose isomerase, His-54 acts as the catalytic base in a ring opening reaction. The ring opening step is followed by binding of D-xylose, in volving two divalent cations, in an extended conformation. The isomerization of D-xylose to D-xylulose involves a metal-mediated 1,2-hydride shift. The final step in the mechanism is a ring closure to produce α-D-xylulose. The ring closing is the reverse of the ring opening step. This mechanism accounts for the majority of xylose iomerase's biochemical properties, in cluding (1) the lack of solvent exchange between the 2-position of D-xylose and the 1-pro-R position of D-xylulose, (2) the chemical modification of histidine and lysine, (3) the pH vs. activity profile, and (4) the requirement for two divalent cations in the mechanism.

Hydrogen exchange in native and denatured states of hen egg-white lysozyme.

Abstract: The hydrogen exchange kinetics of 68 individual amide protons in the native state of hen lysozyme have been measured at pH 7.5 and 30 degrees C by 2D NMR methods. These constitute the most protected subset of amides, with exchange half lives some 10(5)-10(7) times longer than anticipated from studies of small model peptides. The observed distribution of rates under these conditions can be rationalized to a large extent in terms of the hydrogen bonding of individual amides and their burial from bulk solvent. Exchange rates have also been measured in a reversibly denatured state of lysozyme; this was made possible under very mild conditions, pH 2.0 35 degrees C, by lowering the stability of the native state through selective cleavage of the Cys-6-Cys-127 disulfide cross-link (CM6-127 lysozyme). In this state the exchange rates for the majority of amides approach, within a factor of 5, the values anticipated from small model peptides. For a few amides, however, there is evidence for significant retardation (up to nearly 20-fold) relative to the predicted rates. The pattern of protection observed under these conditions does not reflect the behavior of the protein under strongly native conditions, suggesting that regions of native-like structure do not persist significantly in the denatured state of CM6-127 lysozyme. The pattern of exchange rates from the native protein at high temperature, pH 3.8 69 degrees C, resembles that of the acid-denatured state, suggesting that under these conditions the exchange kinetics are dominated by transient global unfolding. The rates of folding and unfolding under these conditions were determined independently by magnetization transfer NMR methods, enabling the intrinsic exchange rates from the denatured state to be deduced on the basis of this model, under conditions where the predominant equilibrium species is the native state. Again, in the case of most amides these rates showed only limited deviation from those predicted by a simple random coil model. This reinforces the view that these denatured states of lysozyme have little persistent residual order and contrasts with the behavior found for compact partially folded states of proteins, including an intermediate detected transiently during the refolding of hen lysozyme.

Fast and simple Monte Carlo algorithm for side chain optimization in proteins: application to model building by homology.

Abstract: An unknown protein structure can be predicted with fair accuracy once an evolutionary connection at the sequence level has been made to a protein of known 3-D structure. In model building by homology, one typically starts with a backbone framework, rebuilds new loop regions, and replaces nonconserved side chains. Here, we use an extremely efficient Monte Carlo algorithm in rotamer space with simulated annealing and simple potential energy functions to optimize the packing of side chains on given backbone models. Optimized models are generated within minutes on a workstation, with reasonable accuracy (average of 81% side chain chi 1 dihedral angles correct in the cores of proteins determined at better than 2.5 A resolution). As expected, the quality of the models decreases with decreasing accuracy of backbone coordinates. If the back-bone was taken from a homologous rather than the same protein, about 70% side chain chi 1 angles were modeled correctly in the core in a case of strong homology and about 60% in a case of medium homology. The algorithm can be used in automated, fast, and reproducible model building by homology.

Structural principles for the propeller assembly of β‐sheets: The preference for seven‐fold symmetry

Abstract: Twisted beta-sheets, packed face to face, may be arranged in circular formation like blades of a propeller or turbine. This beta-propeller fold has been found in three proteins: that in neuraminidase consists of six beta-sheets while those in methylamine dehydrogenase and galactose oxidase are composed of seven beta-sheets. A model for multisheet packing in the beta-propeller fold is proposed. This model gives both geometrical parameters of the beta-propellers composed of different numbers of sheets and patterns of residue packing at their sheet-to-sheet interfaces. All the known beta-propeller structures have been analyzed, and the observed geometries and residue packing are found to be in good agreement with those predicted by models. It is shown that unusual seven-fold symmetry is preferable to six- or eight-fold symmetry for propeller-like multi-sheet assembly. According to the model, a six-beta-sheet propeller has to have predominantly small residues in the beta-strands closed to its six-fold axis, but no strong sequence constraints are necessary for a seven-fold beta-propeller.

Fast structure alignment for protein databank searching.

Abstract: A fast method is described for searching and analyzing the protein structure databank. It uses secondary structure followed by residue matching to compare protein structures and is developed from a previous structural alignment method based on dynamic programming.Linear representations of secondary structures are derived and their features compared to identify equivalent elements in two proteins. The secondary structure alignment then constrains the residue alignment, which compares only residues within aligned secondary structures and with similar buried areas and torsional angles. The initial secondary structure alignment improves accuracy and provides a means of filtering out unrelated proteins before the slower residue alignment stage.It is possible to search or sort the protein structure databank very quickly using just secondary structure comparisons. A search through 720 structures with a probe protein of 10 secondary structures required 1.7 CPU hours on a Sun 4/280. Alternatively, combined secondary structure and residue alignments, with a cutoff on the secondary structure score to remove pairs of unrelated proteins from further analysis, took 10.1 CPU hours. The method was applied in searches on different classes of proteins and to cluster a subset of the databank into structurally related groups. Relationships were consistent with known families of protein structure.

CLIX: a search algorithm for finding novel ligands capable of binding proteins of known three-dimensional structure.

Abstract: A computer algorithm, CLIX, capable of searching a crystallographic data-base of small molecules for candidates which have both steric and chemical likelihood of binding a protein of known three-dimensional structure is presented. The algorithm is a significant advance over previous strategies which consider solely steric or chemical requirements for binding. The algorithm is shown to be capable of predicting the correct binding geometry of sialic acid to a mutant influenza-virus hemagglutinin and of proposing a number of potential new ligands to this protein.

Redesigning the DNA-binding specificity of a zinc finger protein: a data base-guided approach.

Abstract: A peptide corresponding to the three zinc finger domains of the human transcription factor Sp1 has been expressed and found to bind a consensus Sp1 binding site with the sequence 5'-GGGGCGGGG-3'. Examination of the amino acid distributions within a large zinc finger sequence data base and chemical arguments suggested that a particular Arg to Gln sequence change might convert binding specificity to 5'-GGGGCAGGG-3'. Experimental tests of this hypothesis revealed that such a change could be induced only when two other sequence changes, deduced from examination of sequence correlations, were made as well. These results provide the most direct information to date about how zinc finger proteins might recognize adenine-containing binding sites and bear on the existence and nature of any code between zinc finger protein and binding site sequences.

Selection of a representative set of structures from brookhaven protein data-bank

Abstract: Reliable structural and statistical analyses of three dimensional protein structures should be based on unbiased data. The Protein Data Bank is highly redundant, containing several entries for identical or very similar sequences. A technique was developed for clustering the known structures based on their sequences and contents of alpha- and beta-structures. First, sequences were aligned pairwise. A representative sample of sequences was then obtained by grouping similar sequences together, and selecting a typical representative from each group. The similarity significance threshold needed in the clustering method was found by analyzing similarities of random sequences. Because three dimensional structures for proteins of same structural class are generally more conserved than their sequences, the proteins were clustered also according to their contents of secondary structural elements. The results of these clusterings indicate conservation of alpha- and beta-structures even when sequence similarity is relatively low. An unbiased sample of 103 high resolution structures, representing a wide variety of proteins, was chosen based on the suggestions made by the clustering algorithm. The proteins were divided into structural classes according to their contents and ratios of secondary structural elements. Previous classifications have suffered from subjectice view of secondary structures, whereas here the classification was based on backbone geometry. The concise view lead to reclassification of some structures. The representative set of structures facilitates unbiased analyses of relationships between protein sequence, function, and structure as well as of structural characteristics. (Less)

Protease pro region required for folding is a potent inhibitor of the mature enzyme.

Abstract: alpha-Lytic protease, an extracellular bacterial serine protease, is synthesized with a large pro region that is required in vivo for the proper folding of the protease domain. To allow detailed mechanistic study, we have reconstituted pro region-dependent folding in vitro. The pro region promotes folding of the protease domain in the absence of other protein factors or exogenous energy sources. Surprisingly, we find that the pro region is a high affinity inhibitor of the mature protease. The pro region also inhibits the closely related Streptomyces griseus protease B, but not the more distantly related, yet structurally similar protease, elastase. Based on these data, we suggest a mechanism in which pro region binding reduces the free energy of a late folding transition state having native-like structure.

Deconvolution of the circular dichroism spectra of proteins: The circular dichroism spectra of the antiparallel β-sheet in proteins

Abstract: A recently developed algorithm, called Convex Constraint Analysis (CCA), was successfully applied to determine the circular dichroism (CD) spectra of the pure β-pleated sheet in globular proteins. On the basis of X-ray diffraction determined secondary structures, the original data set used (Perczel, A., Hollosi, M., Tusnady, G. Fasman, G.D. Convex constraint analysis: A natural deconvolution of circular dichroism curves of proteins, Prot. Eng., 4:-669–679, 1991), was improved by the addition of proteins with high β-pleated sheet content. The analysis yielded CD curves of the pure components of the main secondary structural elements (α-helix, antiparallel β-pleated sheet, β-turns, and unordered conformation), as well as a curve attributed to the “aromatic contribution” in the wavelength range of 195–240 nm. Upon deconvolution the curves obtained were assigned to various secondary structures. The calculated weights (percentages determining the contributions of each pure component curve in the measured CD spectra of a given protein) were correlated with the X-ray diffraction determined percentages in an assignment procedure and were evaluated. The Pearson product correlation coefficients (R) are significant for all five components. The new pure component curves, which were obtained through deconvolution of the protein CD spectra alone, are promising candidates for determining the percentages of the secondary structural components in globular proteins without the necessity of adopting an X-ray database. The CD spectrum of the CheY protein was interesting because it has the characteristic shape associated with the α-helical structure, but upon analysis yielded a considerable amount of β-sheet in agreement with the X-ray structure. © 1992 Wiley-Liss, Inc.

Modeling the anti-CEA antibody combining site by homology and conformational search.

Abstract: A model for an antibody specific for the carcinoembryonic antigen (CEA) has been constructed using a method which combines the concept of canonical structures with conformational search. A conformational search technique is introduced which couples random generation of backbone loop conformations to a simulated annealing method for assigning side chain conformations. This technique was used both to verify conformations selected from the set of known canonical structures and to explore conformations available to the H3 loop in CEA ab initio. Canonical structures are not available for H3 due to its variability in length, sequence, and observed conformation in known antibody structures. Analysis of the results of conformational search resulted in three equally probable conformations for H3 loop in CEA. Force field energies, solvation free energies, exposure of charged residues and burial of hydrophobic residues, and packing of hydrophobic residues at the base of the loop were used as selection criteria. The existence of three equally plausible structures may reflect the high degree of flexibility expected for an exposed loop of this length. The nature of the combining site and features which could be important to interaction with antigen are discussed.

The molecular structure of UDP‐galactose 4‐epimerase from Escherichia coli determined at 2.5 Å resolution

Abstract: UDP-galactose 4-epimerase catalyzes the conversion of UDP-galactose to UDP-glucose during normal galactose metabolism. The molecular structure of UDP-galactose 4-epimerase from Escherichia coli has now been solved to a nominal resolution of 2.5 A. As isolated from E. coli, the molecule is a dimer of chemically identical subunits with a total molecular weight of 79,000. Crystals of the enzyme used for this investigation were grown as a complex with the substrate analogue, UDP-benzene, and belonged to the space group P2(1)2(1)2(1) with unit cell dimensions of a = 76.3 A, b = 83.1 A, c = 132.1 A, and one dimer per asymmetric unit. An interpretable electron density map calculated to 2.5 A resolution was obtained by a combination of multiple isomorphous replacement with six heavy atom derivatives, molecular averaging, and solvent flattening. Each subunit of epimerase is divided into two domains. The larger N-terminal domain, composed of amino acid residues 1-180, shows a classic NAD+ binding motif with seven strands of parallel beta-pleated sheet flanked on either side of alpha-helices. The seventh strand of the beta-pleated sheet is contributed by amino acid residues from the smaller domain. In addition, this smaller C-terminal domain, consisting of amino acid residues 181-338, contains three strands of beta-pleated sheet, two major alpha-helices and one helical turn. The substrate analogue, UDP-benzene, binds in the cleft located between the two domains with its phenyl ring in close proximity to the nicotinamide ring of NAD+. Contrary to the extensive biochemical literature suggesting that epimerase binds only one NAD+ per functional dimer, the map clearly shows electron density for two nicotinamide cofactors binding in symmetry-related positions in the dimer. Likewise, each subunit in the dimer also binds one substrate analogue.

The refined crystal structure of Pseudomonas putida lipoamide dehydrogenase complexed with NAD+ at 2.45 A resolution.

Abstract: The three-dimensional structure of one of the three lipoamide dehydrogenases occurring in Pseudomonas putida, LipDH Val, has been determined at 2.45 A resolution. The orthorhombic crystals, grown in the presence of 20 mM NAD+, contain 458 residues per asymmetric unit. A crystallographic 2-fold axis generates the dimer which is observed in solution. The final crystallographic R-factor is 21.8% for 18,216 unique reflections and a model consisting of 3,452 protein atoms, 189 solvent molecules and 44 NAD+ atoms, while the overall B-factor is unusually high: 47 angstrom2. The structure of LipDH Val reveals the conformation of the C-terminal residues which fold "back" into the putative lipoamide binding region. The C-terminus has been proven to be important for activity by site-directed mutagenesis. However, the distance of the C-terminus to the catalytically essential residues is surprisingly large, over 6 angstrom, and the precise role of the C-terminus still needs to be elucidated. In this crystal form LipDH Val contains one NAD+ molecule per subunit. Its adenine-ribose moiety occupies an analogous position as in the structure of glutathione reductase. However, the nicotinamide-ribose moiety is far removed from its expected position near the isoalloxazine ring and points into solution. Comparison of LipDH Val with Azotobacter vinelandii lipoamide dehydrogenase yields an rms difference of 1.6 angstrom for 440 well defined C(alpha) atoms per subunit. Comparing LipDH Val with glutathione reductase shows large differences in the tertiary and quaternary structure of the two enzymes. For instance, the two subunits in the dimer are shifted by 6 angstrom with respect to each other. So, LipDH Val confirms the surprising differences in molecular architecture between glutathione reductase and lipoamide dehydrogenase, which were already observed in Azotobacter vinelandii LipDH. This is the more remarkable since the active sites are located at the subunit interface and are virtually identical in all three enzymes.

A model of the platelet factor 4 complex with heparin.

Abstract: A model of heparin bound to bovine platelet factor 4 (BPF4) was completed using a graphically designed heparin molecule and the crystallographic coordinates of the native bovine platelet factor 4 tetramer. The oligosaccharides had a chain length of at least eight disaccharide units with the major repeating disaccharide unit consisting of (I→4)-O-(α-L-idopyranosyluronic acid 2-sulfate)-(l→4)-(2-deoxy-2-sulfamino-2-D-glucopyra-nosyl 6-sulfate). Each disaccharide unit carried a −4.0 charge. The structure of BPF4 was solved to 2.6 A resolution with R = 0.237. Each monomer of BPF4 contains an α-helix lying across 3 strands of antiparallel β-sheet. Each helix has four lysines, which have been implicated in heparin binding. These lysine residues are predominantly on one side of the helix and are solvent accessible. Electrostatic calculations performed on the BPF4 tetramer show a ring of strong, positive charge which runs perpendicularly across the helices. Included in this ring of density is His-38, which has been shown by NMR to have a large pKa shift when heparin binds to BPF4. Our model of heparin bound to PF4 has the anionic polysaccharide perpendicular to the α-helices, wrapped about the tetramer along the ring of positive charge, and salt linked to all four lysines on the helix of each monomer. © 1992 Wiley-Liss, Inc.

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction

Abstract: The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.

Crystal structure of the binary complex of pig muscle phosphoglycerate kinase and its substrate 3-phospho-D-glycerate.

Abstract: Pig muscle phosphoglycerate kinase has been crystallized from polyethyleneglycol in the presence of its substrate 3-phospho-D-glycerate (3-PG) and the structure has been determined at 2.0 A resolution. The structure was solved using the known structure of the substrate-free horse muscle enzyme and has been refined to a crystallographic R-factor of 21.5%. 3-Phospho-D-glycerate is bound to the N-domain of the enzyme through a network of hydrogen bonds to a cluster of basic amino acid residues and by electrostatic interactions between the negatively charged phosphate and these basic protein side chains. This binding site is in good agreement with earlier proposals [Banks et al., Nature (London) 279:773-777, 1979]. The phosphate oxygen atoms are hydrogen bonded to His-62, Arg-65, Arg-122, and Arg-170. The 2-hydroxyl group, which defines the D-isomer of 3PG, is hydrogen bonded to Asp-23 and Asn-25. The carboxyl group of 3-PG points away from the N-domain towards the C-domain and is hydrogen bonded via a water molecule to main chain nitrogen atoms of helix-14. The present structure of the 3-PG-bound pig muscle enzyme is compared with the structure of the substrate-free horse enzyme. Major changes include an ordering of helix-13 and a domain movement, which brings the N-domain closer to the ATP-binding C-domain. This domain movement consists of a 7.7 degree rotation, which is less than previously estimated for the ternary complex. Local changes close to the 3-PG binding site include an ordering of Arg-65 and a shift of helix-5.

Domain flexibility in aspartic proteinases

Abstract: Comparison of the three-dimensional structures of native endothiapepsin (EC 3.4.23.6) and 15 endothiapepsin oligopeptide inhibitor complexes defined at high resolution by X-ray crystallography shows that endothiapepsin exists in two forms differing in the relative orientation of a domain comprising residues 190-302. There are relatively few interactions between the two parts of the enzyme; consequently, they can move as separate rigid bodies. A translational, librational, and screw analysis of the thermal parameters of endothiapepsin also supports a model in which the two parts can move relative to each other. In the comparison of different aspartic proteinases, the rms values are reduced by up to 47% when the two parts of the structure are superposed independently. This justifies description of the differences, including those between pepsinogen and pepsin (EC 3.4.34.1), as a rigid movement of one part relative to another although considerable distortions within the domains also occur. The consequence of the rigid body movement is a change in the shape of the active site cleft that is largest around the S3 pocket. This is associated with a different position and conformation of the inhibitors that are bound to the two endothiapepsin forms. The relevance of these observations to a model of the hydrolysis by aspartic proteinases is briefly discussed.

Comparative analysis of the sequences and structures of HIV-1 and HIV-2 proteases.

Abstract: The different isolates available for HIV-1 and HIV-2 were compared for the region of the protease (PR) sequence, and the variations in amino acids were analyzed with respect to the crystal structure of HIV-1 PR with inhibitor. Based on the extensive homology (39 identical out of 99 residues), models were built of the HIV-2 PR complexed with two different aspartic protease inhibitors, acetylpepstatin and a renin inhibitor, H-261. Comparison of the HIV-1 PR crystal structure and the HIV-2 PR model structure and the analysis of the changes found in different isolates showed that correlated substitutions occur in the hydrophobic interior of the molecule and at surface residues involved in ionic or hydrogen bond interactions. The substrate binding residues of HIV-1 and HIV-2 PRs show conservative substitutions of four residues. The difference in affinity of HIV-1 and HIV-2 PRs for the two inhibitors appears to be due in part to the change of Val 32 in HIV-1 PR to Ile in HIV-2 PR.

Folding protein α-carbon chains into compact forms by Monte Carlo methods

Abstract: A method is presented for generating folded chains of specific amino acid sequences on a simple cubic lattice. Monte Carlo simulations are used to transform extended geometries of simplified alpha-carbon chains for eight small monomeric globular proteins into folded states. Permitted chain transitions are limited to a few types of moves, all restricted to occur on the lattice. Crude residue-residue potentials derived from statistical structure data are used to describe the energies for each conformer. The low resolution structures obtained by this procedure contain many of the correct gross features of the native folded architectures with respect to average residue energy per nonbonded contact, segment density, and location of surface loops and disulfide pairs. Rms deviations between these and the native X-ray structures and percentage of native long-range contacts found in these final folded structures are 7.6 +/- 0.7 A and 48 +/- 3%, respectively. This procedure can be useful for predicting approximate tertiary interactions from amino acid sequence.

Monte Carlo docking of oligopeptides to proteins.

Abstract: A new two-step procedure has been developed for the docking of flexible oligopeptide chains of unkown conformation to static proteins ofknown structure. In the first step positions and conformations are sampled and the association energy, minimized starting from an approximate preselected docking position. The resulting conformations are further optimized in the second step by a Metropolis Monte Carlo minimization, which optimizes each of these structures. The method has been tested on the HIV-1 aspartic proteinase complex with an inhibitor, whose crystallographic structure is known at 2.3 A resolution. Furthermore, the application of this method to the docking of the hendecapeptide 58–68 of the influenza A virus matrix protein to the HLA-A2 molecule produced results which are in agreement with experimental observations in identifying side chains critical for T cell recognition and residues responsible of MHC protein binding. © 1992 Wiley-Liss, Inc.

The high-resolution crystal structure of porcine pepsinogen.

Abstract: The structure of porcine pepsinogen at pH 6.1 has been refined to an R-factor of 0.173 for data extending to 1.65 A. The final model contains 180 solvent molecules and lacks density for residues 157-161. The structure of this aspartic proteinase zymogen possesses many of the characteristics of pepsin, the mature enzyme. The secondary structure of the zymogen consists predominantly of beta-sheet, with an approximate 2-fold axis of symmetry. The activation peptide packs into the active site cleft, and the N-terminus (1P-9P) occupies the position of the mature N-terminus (1-9). Thus changes upon activation include excision of the activation peptide and proper relocation of the mature N-terminus. The activation peptide or residues of the displaced mature N-terminus make specific interactions with the substrate binding subsites. The active site of pepsinogen is intact; thus the lack of activity of pepsinogen is not due to a deformation of the active site. Nine ion pairs in pepsinogen may be important in the advent of activation and involve the activation peptide or regions of the mature N-terminus which are relocated in the mature enzyme. The activation peptide-pepsin junction, 44P-1, is characterized by high thermal parameters and weak density, indicating a flexible structure which would be accessible to cleavage. Pepsinogen is an appropriate model for the structures of other zymogens in the aspartic proteinase family.

Mutating the charged residues in the binding pocket of cellular retinoic acid-binding protein simultaneously reduces its binding affinity to retinoic acid and increases its thermostability

Abstract: Three-dimensional modeling of the complex between retinoic acid-binding protein (CRABP) and retinoic acid suggests thatbinding of the ligand is mediated by interaction between the carboxyl groupof retinoic acid and two charged amino acids (Arg-111 and Arg-131) whose side chains project into the barrel of the protein. To assess the contribution of these amino acids to protein–ligand interaction, amino acid substitutions were made by oligonucleotide-directed, site-specific mutagenesis. The wild-type and mutant proteins were expressed in E. coli and subsequently purified. Like wild-type CRABP, the mutant proteins are composed mainly of β-strands as determined by circular dichroism in the presence and absence of ligand, and thus presumably are folded into the same compact barrel structure as the wild-type protein. Mutants in which Arg-111 and Arg-131 are replaced by glutamine bind retinoic acid with significantly lower affinity than the wild-type protein, arguing that these two residues indeed interact with the ligand. The mutant proteins are more resistant to thermal denaturation than wild-type CRABP in the absence of retinoic acid, but they are not as thermostable as the CRABP–retinoic acid complex. These data suggest a model for CRABP–retinoic acid interaction in which the repulsive forces between the positively-charged arginine residues provide conformational flexibility to the native protein for retinoic acid to enter the binding pocket. Elimination of the positively-charged pair of amino acids produces a protein that is more thermostable than wild-type CRABP but less effective at ligand-binding. © 1992 Wiley-Liss, Inc.