Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Crystal structures of bovine milk xanthine dehydrogenase and xanthine oxidase: Structure-based mechanism of conversion

Mammalian xanthine oxidoreductases, which catalyze the last two steps in the formation of urate, are synthesized as the dehydrogenase form xanthine dehydrogenase (XDH) but can be readily converted to the oxidase form xanthine oxidase (XO) by oxidation of sulfhydryl residues or by proteolysis. Here, we present the crystal structure of the dimeric (Mr, 290,000) bovine milk XDH at 2.1-Å resolution and XO at 2.5-Å resolution and describe the major changes that occur on the proteolytic transformation of XDH to the XO form. Each molecule is composed of an N-terminal 20-kDa domain containing two iron sulfur centers, a central 40-kDa flavin adenine dinucleotide domain, and a C-terminal 85-kDa molybdopterin-binding domain with the four redox centers aligned in an almost linear fashion. Cleavage of surface-exposed loops of XDH causes major structural rearrangement of another loop close to the flavin ring (Gln 423—Lys 433). This movement partially blocks access of the NAD substrate to the flavin adenine dinucleotide cofactor and changes the electrostatic environment of the active site, reflecting the switch of substrate specificity observed for the two forms of this enzyme.

Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

Structure-based design of selective and potent G quadruplex-mediated telomerase inhibitors

The telomerase enzyme is a potential therapeutic target in many human cancers. A series of potent inhibitors has been designed by computer modeling, which exploit the unique structural features of quadruplex DNA. These 3,6,9-trisubstituted acridine inhibitors are predicted to interact selectively with the human DNA quadruplex structure, as a means of specifically inhibiting the action of human telomerase in extending the length of single-stranded telomeric DNA. The anilino substituent at the 9-position of the acridine chromophore is predicted to lie in a third groove of the quadruplex. Calculated relative binding energies predict enhanced selectivity compared with earlier 3,6-disubstituted compounds, as a result of this substituent. The ranking order of energies is in accord with equilibrium binding constants for quadruplex measured by surface plasmon resonance techniques, which also show reduced duplex binding compared with the disubstituted compounds. The 3,6,9-trisubstututed acridines have potent in vitro inhibitory activity against human telomerase, with EC50 values of up to 60 nM.

A structure-based catalytic mechanism for the xanthine oxidase family of molybdenum enzymes.

The crystal structure of the xanthine oxidase-related molybdenum-iron protein aldehyde oxido-reductase from the sulfate reducing anaerobic Gram-negative bacterium Desulfovibrio gigas (Mop) was analyzed in its desulfo-, sulfo-, oxidized, reduced, and alcohol-bound forms at 1.8-A resolution. In the sulfo-form the molybdenum molybdopterin cytosine dinucleotide cofactor has a dithiolene-bound fac-[Mo, = O, = S, ---(OH2)] substructure. Bound inhibitory isopropanol in the inner compartment of the substrate binding tunnel is a model for the Michaelis complex of the reaction with aldehydes (H-C = O,-R). The reaction is proposed to proceed by transfer of the molybdenum-bound water molecule as OH- after proton transfer to Glu-869 to the carbonyl carbon of the substrate in concert with hydride transfer to the sulfido group to generate [MoIV, = O, -SH, ---(O-C = O, -R)). Dissociation of the carboxylic acid product may be facilitated by transient binding of Glu-869 to the molybdenum. The metal-bound water is replenished from a chain of internal water molecules. A second alcohol binding site in the spacious outer compartment may cause the strong substrate inhibition observed. This compartment is the putative binding site of large inhibitors of xanthine oxidase.

Analysis of homeodomain function by structure-based design of a transcription factor.

The homeodomain is a 60-amino acid module which mediates critical protein-DNA and protein-protein interactions for a large family of regulatory proteins. We have used structure-based design to analyze the ability of the Oct-1 homeodomain to nucleate an enhancer complex. The Oct-1 protein regulates herpes simplex virus (HSV) gene expression by participating in the formation of a multiprotein complex (C1 complex) which regulates alpha (immediate early) genes. We recently described the design of ZFHD1, a chimeric transcription factor containing zinc fingers 1 and 2 of Zif268, a four-residue linker, and the Oct-1 homeodomain. In the presence of alpha-transinduction factor and C1 factor, ZFHD1 efficiently nucleates formation of the C1 complex in vitro and specifically activates gene expression in vivo. The sequence specificity of ZFHD1 recruits C1 complex formation to an enhancer element which is not efficiently recognized by Oct-1. ZFHD1 function depends on the recognition of the Oct-1 homeodomain surface. These results prove that the Oct-1 homeodomain mediates all the protein-protein interactions that are required to efficiently recruit alpha-transinduction factor and C1 factor into a C1 complex. The structure-based design of transcription factors should provide valuable tools for dissecting the interactions of DNA-bound domains in other regulatory circuits.

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

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

PALI—a database of Phylogeny and ALIgnment of homologous protein structures

PALI (release 1.2) contains three-dimensional (3-D) structure-dependent sequence alignments as well as structure-based phylogenetic trees of homologous protein domains in various families. The data set of homologous protein structures has been derived by consulting the SCOP database (release 1.50) and the data set comprises 604 families of homologous proteins involving 2739 protein domain structures with each family made up of at least two members. Each member in a family has been structurally aligned with every other member in the same family (pairwise alignment) and all the members in the family are also aligned using simultaneous super­position (multiple alignment). The structural alignments are performed largely automatically, with manual interventions especially in the cases of distantly related proteins, using the program STAMP (version 4.2). Every family is also associated with two dendrograms, calculated using PHYLIP (version 3.5), one based on a structural dissimilarity metric defined for every pairwise alignment and the other based on similarity of topologically equivalent residues. These dendrograms enable easy comparison of sequence and structure-based relationships among the members in a family. Structure-based alignments with the details of structural and sequence similarities, superposed coordinate sets and dendrograms can be accessed conveniently using a web interface. The database can be queried for protein pairs with sequence or structural similarities falling within a specified range. Thus PALI forms a useful resource to help in analysing the relationship between sequence and structure variation at a given level of sequence similarity. PALI also contains over 653 ‘orphans’ (single member families). Using the web interface involving PSI_BLAST and PHYLIP it is possible to associate the sequence of a new protein with one of the families in PALI and generate a phylogenetic tree combining the query sequence and proteins of known 3-D structure. The database with the web interfaced search and dendrogram generation tools can be accessed at http://pa uling.mbu.iisc.ernet.in/~pali.

Crystal structure of a thermostable type B DNA polymerase from Thermococcus gorgonarius

Most known archaeal DNA polymerases belong to the type B family, which also includes the DNA replication polymerases of eukaryotes, but maintain high fidelity at extreme conditions. We describe here the 2.5 Å resolution crystal structure of a DNA polymerase from the Archaea Thermococcus gorgonarius and identify structural features of the fold and the active site that are likely responsible for its thermostable function. Comparison with the mesophilic B type DNA polymerase gp43 of the bacteriophage RB69 highlights thermophilic adaptations, which include the presence of two disulfide bonds and an enhanced electrostatic complementarity at the DNA–protein interface. In contrast to gp43, several loops in the exonuclease and thumb domains are more closely packed; this apparently blocks primer binding to the exonuclease active site. A physiological role of this “closed” conformation is unknown but may represent a polymerase mode, in contrast to an editing mode with an open exonuclease site. This archaeal B DNA polymerase structure provides a starting point for structure-based design of polymerases or ligands with applications in biotechnology and the development of antiviral or anticancer agents.

Structure of the HIV-1 integrase catalytic domain complexed with an inhibitor: A platform for antiviral drug design

HIV integrase, the enzyme that inserts the viral DNA into the host chromosome, has no mammalian counterpart, making it an attractive target for antiviral drug design. As one of the three enzymes produced by HIV, it can be expected that inhibitors of this enzyme will complement the therapeutic use of HIV protease and reverse transcriptase inhibitors. We have determined the structure of a complex of the HIV-1 integrase core domain with a novel inhibitor, 5ClTEP, 1-(5-chloroindol-3-yl)-3-hydroxy-3-(2H-tetrazol-5-yl)-propenone, to 2.1-Å resolution. The inhibitor binds centrally in the active site of the integrase and makes a number of close contacts with the protein. Only minor changes in the protein accompany inhibitor binding. This inhibitor complex will provide a platform for structure-based design of an additional class of inhibitors for antiviral therapy.

Antigen presentation subverted: Structure of the human cytomegalovirus protein US2 bound to the class I molecule HLA-A2

Many persistent viruses have evolved the ability to subvert MHC class I antigen presentation. Indeed, human cytomegalovirus (HCMV) encodes at least four proteins that down-regulate cell-surface expression of class I. The HCMV unique short (US)2 glycoprotein binds newly synthesized class I molecules within the endoplasmic reticulum (ER) and subsequently targets them for proteasomal degradation. We report the crystal structure of US2 bound to the HLA-A2/Tax peptide complex. US2 associates with HLA-A2 at the junction of the peptide-binding region and the α3 domain, a novel binding surface on class I that allows US2 to bind independently of peptide sequence. Mutation of class I heavy chains confirms the importance of this binding site in vivo. Available data on class I-ER chaperone interactions indicate that chaperones would not impede US2 binding. Unexpectedly, the US2 ER-luminal domain forms an Ig-like fold. A US2 structure-based sequence alignment reveals that seven HCMV proteins, at least three of which function in immune evasion, share the same fold as US2. The structure allows design of further experiments to determine how US2 targets class I molecules for degradation.

(4-Aminomethyl)phenylguanidine derivatives as nonpeptidic highly selective inhibitors of human urokinase

Increased expression of the serine protease urokinase-type plasminogen activator (uPA) in tumor tissues is highly correlated with tumor cell migration, invasion, proliferation, progression, and metastasis. Thus inhibition of uPA activity represents a promising target for antimetastatic therapy. So far, only the x-ray crystal structure of uPA inactivated by H-Glu-Gly-Arg-chloromethylketone has been reported, thus limited data are available for a rational structure-based design of uPA inhibitors. Taking into account the trypsin-like arginine specificity of uPA, (4-aminomethyl)phenylguanidine was selected as a potential P1 residue and iterative derivatization of its amino group with various hydrophobic residues, and structure–activity relationship-based optimization of the spacer in terms of hydrogen bond acceptor/donor properties led to N-(1-adamantyl)-N′-(4-guanidinobenzyl)urea as a highly selective nonpeptidic uPA inhibitor. The x-ray crystal structure of the uPA B-chain complexed with this inhibitor revealed a surprising binding mode consisting of the expected insertion of the phenylguanidine moiety into the S1 pocket, but with the adamantyl residue protruding toward the hydrophobic S1′ enzyme subsite, thus exposing the ureido group to hydrogen-bonding interactions. Although in this enzyme-bound state the inhibitor is crossing the active site, interactions with the catalytic residues Ser-195 and His-57 are not observed, but their side chains are spatially displaced for steric reasons. Compared with other trypsin-like serine proteases, the S2 and S3/S4 pockets of uPA are reduced in size because of the 99-insertion loop. Therefore, the peculiar binding mode of the new type of uPA inhibitors offers the possibility of exploiting optimized interactions at the S1′/S2′ subsites to further enhance selectivity and potency. Because crystals of the uPA/benzamidine complex allow inhibitor exchange by soaking procedures, the structure-based design of new generations of uPA inhibitors can rely on the assistance of x-ray analysis.

Estimating the probability for a protein to have a new fold: A statistical computational model

Structural genomics aims to solve a large number of protein structures that represent the protein space. Currently an exhaustive solution for all structures seems prohibitively expensive, so the challenge is to define a relatively small set of proteins with new, currently unknown folds. This paper presents a method that assigns each protein with a probability of having an unsolved fold. The method makes extensive use of protomap, a sequence-based classification, and scop, a structure-based classification. According to protomap, the protein space encodes the relationship among proteins as a graph whose vertices correspond to 13,354 clusters of proteins. A representative fold for a cluster with at least one solved protein is determined after superposition of all scop (release 1.37) folds onto protomap clusters. Distances within the protomap graph are computed from each representative fold to the neighboring folds. The distribution of these distances is used to create a statistical model for distances among those folds that are already known and those that have yet to be discovered. The distribution of distances for solved/unsolved proteins is significantly different. This difference makes it possible to use Bayes' rule to derive a statistical estimate that any protein has a yet undetermined fold. Proteins that score the highest probability to represent a new fold constitute the target list for structural determination. Our predicted probabilities for unsolved proteins correlate very well with the proportion of new folds among recently solved structures (new scop 1.39 records) that are disjoint from our original training set.

The design, synthesis, and evaluation of molecules that enable or enhance cellular uptake: Peptoid molecular transporters

Certain proteins contain subunits that enable their active translocation across the plasma membrane into cells. In the specific case of HIV-1, this subunit is the basic domain Tat49–57 (RKKRRQRRR). To establish the optimal structural requirements for this translocation process, and thereby to develop improved molecular transporters that could deliver agents into cells, a series of analogues of Tat49–57 were prepared and their cellular uptake into Jurkat cells was determined by flow cytometry. All truncated and alanine-substituted analogues exhibited diminished cellular uptake, suggesting that the cationic residues of Tat49–57 play a principal role in its uptake. Charge alone, however, is insufficient for transport as oligomers of several cationic amino acids (histidine, lysine, and ornithine) are less effective than Tat49–57 in cellular uptake. In contrast, a 9-mer of l-arginine (R9) was 20-fold more efficient than Tat49–57 at cellular uptake as determined by Michaelis–Menton kinetic analysis. The d-arginine oligomer (r9) exhibited an even greater uptake rate enhancement (>100-fold). Collectively, these studies suggest that the guanidinium groups of Tat49–57 play a greater role in facilitating cellular uptake than either charge or backbone structure. Based on this analysis, we designed and synthesized a class of polyguanidine peptoid derivatives. Remarkably, the subset of peptoid analogues containing a six-methylene spacer between the guanidine head group and backbone (N-hxg), exhibited significantly enhanced cellular uptake compared to Tat49–57 and even to r9. Overall, a transporter has been developed that is superior to Tat49–57, protease resistent, and more readily and economically prepared.

CKAAPs DB: a conserved key amino acid positions database

The Conserved Key Amino Acid Positions DataBase (CKAAPs DB) provides access to an analysis of structurally similar proteins with dissimilar sequences where key residues within a common fold are identified. The derivation and significance of CKAAPs starting from pairwise structure alignments is described fully in Reddy et al. [Reddy,B.V.B., Li,W.W., Shindyalov,I.N. and Bourne,P.E. (2000) Proteins, in press]. The CKAAPs identified from this theoretical analysis are provided to experimentalists and theoreticians for potential use in protein engineering and modeling. It has been suggested that CKAAPs may be crucial features for protein folding, structural stability and function. Over 170 substructures, as defined by the Combinatorial Extension (CE) database, which are found in approximately 3000 representative polypeptide chains have been analyzed and are available in the CKAAPs DB. CKAAPs DB also provides CKAAPs of the representative set of proteins derived from the CE and FSSP databases. Thus the database contains over 5000 representative poly­peptide chains, covering all known structures in the PDB. A web interface to a relational database permits fast retrieval of structure-sequence alignments, CKAAPs and associated statistics. Users may query by PDB ID, protein name, function and Enzyme Classification number. Users may also submit protein alignments of their own to obtain CKAAPs. An interface to display CKAAPs on each structure from a web browser is also being implemented. CKAAPs DB is maintained by the San Diego Supercomputer Center and accessible at the URL http://ckaaps.sdsc.edu.