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.

Efficient generation of recombinant adenoviruses using adenovirus DNA-terminal protein complex and a cosmid bearing the full-length virus genome.

An efficient method of constructing recombinant adenoviruses (Ads) has been established. The expression unit to be introduced into recombinant Ad was first inserted into the unique Swa I site of the full-length Ad genome cloned in a cassette cosmid. The cassette bearing the expression unit was then cotransfected into human embryonic kidney 293 cells together with the Ad DNA-terminal protein complex digested at several sites with Eco T22I or Ase I/EcoRI. The use of the parent Ad DNA-terminal protein complex instead of the deproteinized Ad genome DNA allowed very efficient recovery of the desired recombinant Ad, and the above restriction digestion drastically reduced regeneration of the parent virus. Several hundred virus clones were readily obtained in each experiment, and about 70% of the clones were the desired recombinant viruses. Furthermore, because the cassette contained the full-length Ad genome, any position of the genome could be easily modified to develop a new vector design. We established construction systems for two types of Ad vectors, the E1-substitution type and the E4-insertion type. This method may greatly facilitate the application of recombinant Ads and should be useful for further improvement of Ad vectors.

Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution: implications for catalytic proton transfer and inhibitor design.

The three-dimensional structure of murine mitochondrial carbonic anhydrase V has been determined and refined at 2.45-A resolution (crystallographic R factor = 0.187). Significant structural differences unique to the active site of carbonic anhydrase V are responsible for differences in the mechanism of catalytic proton transfer as compared with other carbonic anhydrase isozymes. In the prototypical isozyme, carbonic anhydrase II, catalytic proton transfer occurs via the shuttle group His-64; carbonic anhydrase V has Tyr-64, which is not an efficient proton shuttle due in part to the bulky adjacent side chain of Phe-65. Based on analysis of the structure of carbonic anhydrase V, we speculate that Tyr-131 may participate in proton transfer due to its proximity to zinc-bound solvent, its solvent accessibility, and its electrostatic environment in the protein structure. Finally, the design of isozyme-specific inhibitors is discussed in view of the complex between carbonic anhydrase V and acetazolamide, a transition-state analogue. Such inhibitors may be physiologically important in the regulation of blood glucose levels.

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.

Apoptosis induced by a human milk protein.

To the breast-fed infant, human milk is more than a source of nutrients; it furnishes a wide array of molecules that restrict microbes, such as antibodies, bactericidins, and inhibitors of bacterial adherence. However, it has rarely been considered that human milk may also contain substances bioactive toward host cells. While investigating the effect of human milk on bacterial adherence to a human lung cancer cell line, we were surprised to discover that the milk killed the cells. Analysis of this effect revealed that a component of milk in a particular physical state--multimeric alpha-lact-albumin--is a potent Ca(2+)-elevating and apoptosis-inducing agent with broad, yet selective, cytotoxic activity. Multimeric alpha-lactalbumin killed all transformed, embryonic, and lymphoid cells tested but spared mature epithelial elements. These findings raise the possibility that milk contributes to mucosal immunity not only by furnishing antimicrobial molecules but also by policing the function of lymphocytes and epithelium. Finally, analysis of the mechanism by which multimeric alpha-lactalbumin induces apoptosis in transformed epithelial cells could lead to the design of antitumor agents.