Introduction of a [4Fe-4S (S-cys)4]+1,+2 iron-sulfur center into a four-alpha helix protein using design parameters from the domain of the Fx cluster in the Photosystem I reaction center.

We describe the insertion of an iron-sulfur center into a designed four alpha-helix model protein. The model protein was re-engineered by introducing four cysteine ligands required for the coordination of the mulinucleate cluster into positions in the main-chain directly analogous to the domain predicted to ligand the interpeptide [4Fe-4S (S-cys)4] cluster, Fx, from PsaA and PsaB of the Photosystem I reaction center. This was achieved by inserting the sequence, CDGPGRGGTC, which is conserved in PsaA and PsaB, into interhelical loops 1 and 3 of the four alpha-helix model. The holoprotein was characterized spectroscopically after insertion of the iron-sulfur center in vitro. EPR spectra confirmed the cluster is a [4Fe-4S] type, indicating that the cysteine thiolate ligands were positioned as designed. The midpoint potential of the iron-sulfur center in the model holoprotein was determined via redox titration and shown to be -422 mV (pH 8.3, n = 1). The results support proposals advanced for the structure of the domain of the [4Fe-4S] Fx cluster in Photosystem I based upon sequence predictions and molecular modeling. We suggest that the lower potential of the Fx cluster is most likely due to factors in the protein environment of Fx rather than the identity of the residues proximal to the coordinating ligands.     

Cross-strand side-chain interactions versus turn conformation in beta-hairpins.

A series of designed peptides has been analyzed by 1H-NMR spectroscopy in order to investigate the influence of cross-strand side-chain interactions in beta-hairpin formation. The peptides differ in the N-terminal residues of a previously designed linear decapeptide that folds in aqueous solution into two interconverting beta-hairpin conformations, one with a type I turn (beta-hairpin 4:4) and the other with a type I + G1 beta-bulge turn (beta-hairpin 3:5). Analysis of the conformational behavior of the peptides studied here demonstrates three favorable and two unfavorable cross-strand side-chain interactions for beta-hairpin formation. These results are in agreement with statistical data on side-chain interactions in protein beta-sheets. All the peptides in this study form significant populations of the beta-hairpin 3:5, but only some of them also adopt the beta-hairpin 4:4. The formation of beta-hairpin 4:4 requires the presence of at least two favorable cross-strand interactions, whereas beta-hairpin 3:5 seems to be less susceptible to side-chain interactions. A protein database analysis of beta-hairpins 3:5 and beta-hairpins 4:4 indicates that the former occur more frequently than the latter. In both peptides and proteins, beta-hairpins 3:5 have a larger right-handed twist than beta-hairpins 4:4, so that a factor contributing to the higher stability of beta-hairpin 3:5 relative to beta-hairpin 4:4 is due to an appropriate backbone conformation of the type I + G1 beta-bulge turn toward the right-handed twist usually observed in protein beta-sheets. In contrast, as suggested previously, backbone geometry of the type I turn is not adequate for the right-handed twist. Because analysis of buried hydrophobic surface areas on protein beta-hairpins reveals that beta-hairpins 3:5 bury more hydrophobic surface area than beta-hairpins 4:4, we suggest that the right-handed twist observed in beta-hairpin 3:5 allows a better packing of side chains and that this may also contribute to its higher intrinsic stability.     

De novo design of the hydrophobic core of ubiquitin.

We have previously reported the development and evaluation of a computational program to assist in the design of hydrophobic cores of proteins. In an effort to investigate the role of core packing in protein structure, we have used this program, referred to as Repacking of Cores (ROC), to design several variants of the protein ubiquitin. Nine ubiquitin variants containing from three to eight hydrophobic core mutations were constructed, purified, and characterized in terms of their stability and their ability to adopt a uniquely folded native-like conformation. In general, designed ubiquitin variants are more stable than control variants in which the hydrophobic core was chosen randomly. However, in contrast to previous results with 434 cro, all designs are destabilized relative to the wild-type (WT) protein. This raises the possibility that beta-sheet structures have more stringent packing requirements than alpha-helical proteins. A more striking observation is that all variants, including random controls, adopt fairly well-defined conformations, regardless of their stability. This result supports conclusions from the cro studies that non-core residues contribute significantly to the conformational uniqueness of these proteins while core packing largely affects protein stability and has less impact on the nature or uniqueness of the fold. Concurrent with the above work, we used stability data on the nine ubiquitin variants to evaluate and improve the predictive ability of our core packing algorithm. Additional versions of the program were generated that differ in potential function parameters and sampling of side chain conformers. Reasonable correlations between experimental and predicted stabilities suggest the program will be useful in future studies to design variants with stabilities closer to that of the native protein. Taken together, the present study provides further clarification of the role of specific packing interactions in protein structure and stability, and demonstrates the benefit of using systematic computational methods to predict core packing arrangements for the design of proteins.     

Crystal structure of a non-toxic mutant of heat-labile enterotoxin, which is a potent mucosal adjuvant.

Two closely related bacterial toxins, heat-labile enterotoxin (LT-I) and cholera toxin (CT), not only invoke a toxic activity that affects many victims worldwide but also contain a beneficial mucosal adjuvant activity that significantly enhances the potency of vaccines in general. For the purpose of vaccine design it is most interesting that the undesirable toxic activity of these toxins can be eliminated by the single-site mutation Ser63Lys in the A subunit while the mucosal adjuvant activity is still present. The crystal structure of the Ser63Lys mutant of LT-I is determined at 2.0 A resolution. Its structure appears to be essentially the same as the wild-type LT-I structure. The substitution Ser63Lys was designed, based on the wild-type LT-I crystal structure, to decrease toxicity by interfering with NAD binding and/or catalysis. In the mutant crystal structure, the newly introduced lysine side chain is indeed positioned such that it could potentially obstruct the productive binding mode of the substrate NAD while at the same time its positive charge could possibly interfere with the critical function of nearby charged groups in the active site of LT-I. The fact that the Ser63Lys mutant of LT-I does not disrupt the wild-type LT-I structure makes the non-toxic mutant potentially suitable, from a structural point of view, to be used as a vaccine to prevent enterotoxigenic E. coli infections. The structural similarity of mutant and wild-type toxin might also be the reason why the inactive Ser63Lys variant retains its adjuvant activity.     

Characterization of the allosteric binding pocket of human liver fructose-1,6-bisphosphatase by protein crystallography and inhibitor activity studies.

The structures of three complexes of human fructose-1,6-bisphosphatase (FB) with the allosteric inhibitor AMP and two AMP analogues have been determined and all fully refined. The data used for structure determination were collected at cryogenic temperature (110 K), and with the use of synchrotron radiation. The structures reveal a common mode of binding for AMP and formycine monophosphate (FMP). 5-Amino-4-carboxamido-1 beta-D-5-phosphate-ribofuranosyl-1H-imidazole (AICAR-P) shows an unexpected mode of binding to FB, different from that of the other two ligands. The imidazole ring of AICAR-P is rotated 180 degrees compared to the AMP and FMP bases. This rotation results in a slightly different hydrogen bonding pattern and minor changes in the water structure in the binding pocket. Common features of binding are seen for the ribose and phosphate moieties of all three compounds. Although binding in a different mode, AICAR-P is still capable of making all the important interactions with the residues building the allosteric binding pocket. The IC50 values of AMP, FMP, and AICAR-P were determined to be 1.7, 1.4, and 20.9 microM, respectively. Thus, the approximately 10 times lower potency of AICAR-P is difficult to explain solely from the variations observed in the binding pocket. Only one water molecule in the allosteric binding pocket was found to be conserved in all four subunits in all three structures. This water molecule coordinates to a phosphate oxygen atom and the N7 atom of the AMP molecule, and to similarly situated atoms in the FMP and AICAR-P complexes. This implies an important role of the conserved water molecule in binding of the ligand.     

Multiple modes of ligand recognition: Crystal structures of cyclin-dependent protein kinase 2 in complex with ATP and two inhibitors,olomoucine and isopentenyladenine

Cyclin-dependent kinases (CDKs) are conserved regulators of the eukaryotic cell cycle with different isoforms controlling specific phases of the cell cycle. Mitogenic or growth inhibitory signals are mediated, respectively, by activation or inhibition of CDKs which phosphorylate proteins associated with the cell cycle. The central role of CDKs in cell cycle regulation makes them a potential new target for inhibitory molecules with anti-proliferative and/or anti-neoplastic effects. We describe the crystal structures of the complexes of CDK2 with a weakly specific CDK inhibitor, N6-(δ2-isopentenyl)adenine, and a strongly specific inhibitor, olomoucine. Both inhibitors are adenine derivatives and bind in the adenine binding pocket of CDK2, but in an unexpected and different orientation from the adenine of the authentic ligand ATP. The N6-benzyl substituent in olomoucine binds outside the conserved binding pocket and is most likely responsible for its specificity. The structural information from the CDK2-olomoucine complex will be useful in directing the search for the next generation inhibitors with improved properties. © 1995 Wiley-Liss, Inc.     

Atomic and residue hydrophilicity in the context of folded protein structures

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes. © 1995 Wiley-Liss, Inc.     

Fermenter production of an artificial fab fragment,rationally designed for the antigen cystatin,and its optimized crystallization through constant domain shuffling

The synthetic antibody model “M41” was rationally designed with a binding site complementary to chicken egg white cystatin as the prescribed antigen. In order to permit comparison between the computer model and an experimental three-dimensional structure of the artificial protein, its X-ray crystallographic analysis was pursued. For this purpose, M41 was expressed as a recombinant Fab fragment in E. coli by medium cell density fermentation employing the tightly regulated tetracycline promoter. The Fab fragment was efficiently purified via a His-6 tail fused to its heavy chain and immobilized metal affinity chromatography. To raise the chances for the productive formation of crystal packing contacts, three versions of the Fab fragment were generated with differing constant domains. One of these, the variant with murine C_K and C_H 1_γ1 domains, was successfully crystallized by microseeding in a sitting drop. The orthorhombic crystals exhibited symmetry of the space group P2₁2₁2₁ with unit cell dimensions a = 104.7 Å, b = 113.9 Å, c = 98.8 Å and diffracted X-rays to a nominal resolution of 2.5 Å. © 1995 Wiley-Liss, Inc.     

A multivariate approach to the selection of biological reserves

Multivariate analysis provides an effective context for the examination of some significant aspects of biodiversity and conservation. The framework is a multidimensional space that integrates sample sites, taxa and environments. This approach enables terms such as representativeness, complementarity and irreplaceability to be integrated within an intuitive and practical framework for reserve design. Cluster analysis is proposed to determine what is there by defining a set of complementary clusters. These clusters are sampled in a representative manner; from the core outward. The degree of irreplaceability of a site is defined as the multivariate distance of each potential reserve site to its nearest neighbour.