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.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

Synthetic genes for glycoprotein design and the elucidation of hydroxyproline-O-glycosylation codes

Design of hydroxyproline (Hyp)-rich glycoproteins (HRGPs) offers an approach for the structural and functional analysis of these wall components, which are broadly implicated in plant growth and development. HRGPs consist of multiple small repetitive “glycomodules” extensively O-glycosylated through the Hyp residues. The patterns of Hyp-O-glycosylation are putatively coded by the primary sequence as described by the Hyp contiguity hypothesis, which predicts contiguous Hyp residues to be attachment sites of small arabinooligosaccharides (1–5 Ara residues/Hyp); while clustered, noncontiguous Hyp residues are sites of arabinogalactan polysaccharide attachment. As a test, we designed two simple HRGPs as fusion proteins with green fluorescent protein. The first was a repetitive Ser-Hyp motif that encoded only clustered noncontiguous Hyp residues, predicted polysaccharide addition sites. The resulting glycoprotein had arabinogalactan polysaccharide O-linked to all Hyp residues. The second construct, based on the consensus sequence of a gum arabic HRGP, contained both arabinogalactan and arabinooligosaccharide addition sites and, as predicted, gave a product that contained both saccharide types. These results identify an O-glycosylation code of plants.

A general method to design dominant negatives to B-HLHZip proteins that abolish DNA binding

We describe a method to design dominant-negative proteins (D-N) to the basic helix–loop–helix–leucine zipper (B-HLHZip) family of sequence-specific DNA binding transcription factors. The D-Ns specifically heterodimerize with the B-HLHZip dimerization domain of the transcription factors and abolish DNA binding in an equimolar competition. Thermal denaturation studies indicate that a heterodimer between a Myc B-HLHZip domain and a D-N consisting of a 12-amino acid sequence appended onto the Max dimerization domain (A-Max) is −6.3 kcal⋅mol−1 more stable than the Myc:Max heterodimer. One molar equivalent of A-Max can totally abolish the DNA binding activity of a Myc:Max heterodimer. This acidic extension also has been appended onto the dimerization domain of the B-HLHZip protein Mitf, a member of the transcription factor enhancer binding subfamily, to produce A-Mitf. The heterodimer between A-Mitf and the B-HLHZip domain of Mitf is −3.7 kcal⋅mol−1 more stable than the Mitf homodimer. Cell culture studies show that A-Mitf can inhibit Mitf-dependent transactivation both in acidic extension and in a dimerization-dependent manner. A-Max can inhibit Myc-dependent foci formation twice as well as the Max dimerization domain (HLHZip). This strategy of producing D-Ns may be applicable to other B-HLHZip or B-HLH proteins because it provides a method to inhibit the DNA binding of these transcription factors in a dimerization-specific manner.

Design, synthesis, and characterization of a photoactivatable flavocytochrome molecular maquette

We report the construction of a synthetic flavo-heme protein that incorporates two major physiological activities of flavoproteins: light activation of flavin analogous to DNA photolyase and rapid intramolecular electron transfer between the flavin and heme cofactors as in several oxidoreductases. The functional tetra-α-helix protein comprises two 62-aa helix-loop-helix subunits. Each subunit contains a single cysteine to which flavin (7-acetyl-10-methylisoalloxazine) is covalently attached and two histidines appropriately positioned for bis-his coordination of heme cofactors. Both flavins and hemes are situated within the hydrophobic core of the protein. Intramolecular electron transfer from flavosemiquinone generated by photoreduction from a sacrificial electron donor in solution was examined between protoporphyrin IX and 1-methyl-2-oxomesoheme XIII. Laser pulse-activated electron transfer from flavin to meso heme occurs on a 100-ns time scale, with a favorable free energy of approximately −100 meV. Electron transfer from flavin to the lower potential protoporphyrin IX, with an unfavorable free energy, can be induced after a lag phase under continuous light illumination. Thus, the supporting peptide matrix provides an excellent framework for the positioning of closely juxtaposed redox groups capable of facilitating intramolecular electron transfer and begins to clarify in a simplified and malleable system the natural engineering of flavoproteins.

Design of three-dimensional domain-swapped dimers and fibrous oligomers

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 Å resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

Rational design of a bacterial transcriptional cascade for amplifying gene expression capacity

Cascade regulatory circuits have been described that control numerous cell processes, and may provide models for the design of artificial circuits with novel properties. Here we describe the design of a transcriptional regulatory cascade to amplify the cell response to a given signal. We used the salicylate-responsive activators of Pseudomonas putida NahR of the naphthalene degradation plasmid NAH7 and XylS2, a mutant regulator of the TOL plasmid for catabolism of m-xylene and their respective cognate promoters Psal and Pm. Control of the expression of xylS2 with the nahR/Psal system permitted either their selective activation with specific effectors for each protein or the simultaneous activation of both of them with salicylate. When cells face the common effector of the two regulators, both the increase in XylS2 concentration and the stimulation of its activity act synergistically on the Pm promoter, amplifying the gene expression capacity by at least one order of magnitude with respect to the individual systems. By changing the hierarchy of regulators, we showed that the specific features of the downstream regulator were crucial for the amplification effect. Directed changes in the effector profile of the regulators allowed the extension of the amplifying system to other molecular signals.

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.

Conversion of a maltose receptor into a zinc biosensor by computational design

We have demonstrated that it is possible to radically change the specificity of maltose binding protein by converting it into a zinc sensor using a rational design approach. In this new molecular sensor, zinc binding is transduced into a readily detected fluorescence signal by use of an engineered conformational coupling mechanism linking ligand binding to reporter group response. An iterative progressive design strategy led to the construction of variants with increased zinc affinity by combining binding sites, optimizing the primary coordination sphere, and exploiting conformational equilibria. Intermediates in the design series show that the adaptive process involves both introduction and optimization of new functions and removal of adverse vestigial interactions. The latter demonstrates the importance of the rational design approach in uncovering cryptic phenomena in protein function, which cannot be revealed by the study of naturally evolved systems.

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.

The design of an agent to bend DNA.

An artificial DNA bending agent has been designed to assess helix flexibility over regions as small as a protein binding site. Bending was obtained by linking a pair of 15-base-long triple helix forming oligonucleotides (TFOs) by an adjustable polymeric linker. By design, DNA bending was introduced into the double helix within a 10-bp spacer region positioned between the two sites of 15-base triple helix formation. The existence of this bend has been confirmed by circular permutation and phase-sensitive electrophoresis, and the directionality of the bend has been determined as a compression of the minor helix groove. The magnitude of the resulting duplex bend was found to be dependent on the length of the polymeric linker in a fashion consistent with a simple geometric model. Data suggested that a 50-70 degrees bend was achieved by binding of the TFO chimera with the shortest linker span (18 rotatable bonds). Equilibrium analysis showed that, relative to a chimera which did not bend the duplex, the stability of the triple helix possessing a 50-70 degrees bend was reduced by less than 1 kcal/mol of that of the unbent complex. Based upon this similarity, it is proposed that duplex DNA may be much more flexible with respect to minor groove compression than previously assumed. It is shown that this unusual flexibility is consistent with recent quantitation of protein-induced minor groove bending.

A test of the "jigsaw puzzle" model for protein folding by multiple methionine substitutions within the core of T4 lysozyme.

To test whether the structure of a protein is determined in a manner akin to the assembly of a jigsaw puzzle, up to 10 adjacent residues within the core of T4 lysozyme were replaced by methionine. Such variants are active and fold cooperatively with progressively reduced stability. The structure of a seven-methionine variant has been shown, crystallographically, to be similar to wild type and to maintain a well ordered core. The interaction between the core residues is, therefore, not strictly comparable with the precise spatial complementarity of the pieces of a jigsaw puzzle. Rather, a certain amount of give and take in forming the core structure is permitted. A simplified hydrophobic core sequence, imposed without genetic selection or computer-based design, is sufficient to retain native properties in a globular protein.

Propagation of an attenuated virus by design: engineering a novel receptor for a noninfectious foot-and-mouth disease virus.

To gain entry into cells, viruses utilize a variety of different cell-surface molecules. Foot-and-mouth disease virus (FMDV) binds to cell-surface integrin molecules via an arginine-glycine-aspartic acid (RGD) sequence in capsid protein VP1. Binding to this particular cell-surface molecule influences FMDV tropism, and virus/receptor interactions appear to be responsible, in part, for selection of antigenic variants. To study early events of virus-cell interaction, we engineered an alternative and novel receptor for FMDV. Specifically, we generated a new receptor by fusing a virus-binding, single-chain antibody (scAb) to intracellular adhesion molecule 1 (ICAM1). Cells that are normally not susceptible to FMDV infection became susceptible after being transfected with DNA encoding the scAb/ICAM1 protein. An escape mutant (B2PD.3), derived with the mAb used to generate the genetically engineered receptor, was restricted for growth on the scAb/ICAM1 cells, but a variant of B2PD.3 selected by propagation on scAb/ICAM1 cells grew well on these cells. This variant partially regained wild-type sequence in the epitope recognized by the mAb and also regained the ability to be neutralize by the mAb. Moreover, RGD-deleted virions that are noninfectious in animals and other cell types grew to high titers and were able to form plaques on scAb/ ICAM1 cells. These studies demonstrate the first production of a totally synthetic cell-surface receptor for a virus. This novel approach will be useful for studying virus reception and for the development of safer vaccines against viral pathogens of animals and humans.

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.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.