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.

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.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having α,β-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme’s catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

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.

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.

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.

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.

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.

(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.

Minimizing a binding domain from protein A.

We present a systematic approach to minimizing the Z-domain of protein A, a three-helix bundle (59 residues total) that binds tightly (Kd = 10 nM) to the Fc portion of an immunoglobin IgG1. Despite the fact that all the contacts seen in the x-ray structure of the complex with the IgG are derived from residues in the first two helices, when helix 3 is deleted, binding affinity is reduced > 10(5)-fold (Kd > 1 mM). By using structure-based design and phage display methods, we have iteratively improved the stability and binding affinity for a two-helix derivative, 33 residues in length, such that it binds IgG1, with a Kd of 43 nM. This was accomplished by stepwise selection of random mutations from three regions of the truncated Z-peptide: the 4 hydrophobic residues from helix 1 and helix 2 that contacted helix 3 (the exoface), followed by 5 residues between helix 1 and helix 2 (the intraface), and lastly by 19 residues at or near the interface that interacts with Fc (the interface). As selected mutations from each region were compiled (12 in total), they led to progressive increases in affinity for IgG, and concomitant increases in alpha-helical content reflecting increased stabilization of the two-helix scaffold. Thus, by sequential increases in the stability of the structure and improvements in the quality of the intermolecular contacts, one can reduce larger binding domains to smaller ones. Such mini-protein binding domains are more amenable to synthetic chemistry and thus may be useful starting points for the design of smaller organic mimics. Smaller binding motifs also provide simplified and more tractable models for understanding determinants of protein function and stability.

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.

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.

Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.