Computational design of a Zn2+ receptor that controls bacterial gene expression

Engineered Metal-Binding Sites to Probe Protein Folding Transition States: Psi Analysis

The formation of the transition state ensemble (TSE) represents the rate-limiting step in protein folding. The TSE is the least populated state on the pathway, and its characterization remains a challenge. Properties of the TSE can be inferred from the effects on folding and unfolding rates for various perturbations. A difficulty remains on how to translate these kinetic effects to structural properties of the TSE. Several factors can obscure the translation of point mutations in the frequently used method, "mutational Phi analysis." We take a complementary approach in "Psi analysis," employing rationally inserted metal binding sites designed to probe pairwise contacts in the TSE. These contacts can be confidently identified and used to construct structural models of the TSE. The method has been applied to multiple proteins and consistently produces a considerably more structured and native-like TSE than Phi analysis. This difference has significant implications to our understanding of protein folding mechanisms. Here we describe the application of the method and discuss how it can be used to study other conformational transitions such as binding.

129Xe NMR-Protein Sensor Reveals Cellular Ribose Concentration

Dysregulation of cellular ribose uptake can be indicative of metabolic abnormalities or tumorigenesis. However, analytical methods are currently limited for quantifying ribose concentration in complex biological samples. Here, we utilize the highly specific recognition of ribose by ribose-binding protein (RBP) to develop a single-protein ribose sensor detectable via a sensitive NMR technique known as hyperpolarized 129Xe chemical exchange saturation transfer (hyper-CEST). We demonstrate that RBP, with a tunable ribose-binding site and further engineered to bind xenon, enables the quantitation of ribose over a wide concentration range (nM to mM). Ribose binding induces the RBP "closed" conformation, which slows Xe exchange to a rate detectable by hyper-CEST. Such detection is remarkably specific for ribose, with the minimal background signal from endogenous sugars of similar size and structure, for example, glucose or ribose-6-phosphate. Ribose concentration was measured for mammalian cell lysate and serum, which led to estimates of low-mM ribose in a HeLa cell line. This highlights the potential for using genetically encoded periplasmic binding proteins such as RBP to measure metabolites in different biological fluids, tissues, and physiologic states.

Probing the minimal determinants of zinc binding with computational protein design

Structure-based protein design tests our understanding of the minimal determinants of protein structure and function. Previous studies have demonstrated that placing zinc binding amino acids (His, Glu, Asp or Cys) near each other in a folded protein in an arrangement predicted to be tetrahedral is often sufficient to achieve binding to zinc. However, few designs have been characterized with high-resolution structures. Here, we use X-ray crystallography, binding studies and mutation analysis to evaluate three alternative strategies for designing zinc binding sites with the molecular modeling program Rosetta. While several of the designs were observed to bind zinc, crystal structures of two designs reveal binding configurations that differ from the design model. In both cases, the modeling did not accurately capture the presence or absence of second-shell hydrogen bonds critical in determining binding-site structure. Efforts to more explicitly design second-shell hydrogen bonds were largely unsuccessful as evidenced by mutation analysis and low expression of proteins engineered with extensive primary and secondary networks. Our results suggest that improved methods for designing interaction networks will be needed for creating metal binding sites with high accuracy.

Designing pH induced fold switch in proteins

This work investigates the computational design of a pH induced protein fold switch based on a self-consistent mean-field approach by identifying the ensemble averaged characteristics of sequences that encode a fold switch. The primary challenge to balance the alternative sets of interactions present in both target structures is overcome by simultaneously optimizing two foldability criteria corresponding to two target structures. The change in pH is modeled by altering the residual charge on the amino acids. The energy landscape of the fold switch protein is found to be double funneled. The fold switch sequences stabilize the interactions of the sites with similar relative surface accessibility in both target structures. Fold switch sequences have low sequence complexity and hence lower sequence entropy. The pH induced fold switch is mediated by attractive electrostatic interactions rather than hydrophobic-hydrophobic contacts. This study may provide valuable insights to the design of fold switch proteins.

Escherichia [corrected] coli ribose binding protein based bioreporters revisited

Bioreporter bacteria, i.e., strains engineered to respond to chemical exposure by production of reporter proteins, have attracted wide interest because of their potential to offer cheap and simple alternative analytics for specified compounds or conditions. Bioreporter construction has mostly exploited the natural variation of sensory proteins, but it has been proposed that computational design of new substrate binding properties could lead to completely novel detection specificities at very low affinities. Here we reconstruct a bioreporter system based on the native Escherichia coli ribose binding protein RbsB and one of its computationally designed variants, reported to be capable of binding 2,4,6-trinitrotoluene (TNT). Our results show in vivo reporter induction at 50 nM ribose, and a 125 nM affinity constant for in vitro ribose binding to RbsB. In contrast, the purified published TNT-binding variant did not bind TNT nor did TNT cause induction of the E. coli reporter system.

Design of catalytically amplified sensors for small molecules

Catalytically amplified sensors link an allosteric analyte binding site with a reactive site to catalytically convert substrate into colored or fluorescent product that can be easily measured. Such an arrangement greatly improves a sensor’s detection limit as illustrated by successful application of ELISA-based approaches. The ability to engineer synthetic catalytic sites into non-enzymatic proteins expands the repertoire of analytes as well as readout reactions. Here we review recent examples of small molecule sensors based on allosterically controlled enzymes and organometallic catalysts. The focus of this paper is on biocompatible, switchable enzymes regulated by small molecules to track analytes both in vivo and in the environment.

Designing hydrolytic zinc metalloenzymes

Zinc is an essential element required for the function of more than 300 enzymes spanning all classes. Despite years of dedicated study, questions regarding the connections between primary and secondary metal ligands and protein structure and function remain unanswered, despite numerous mechanistic, structural, biochemical, and synthetic model studies. Protein design is a powerful strategy for reproducing native metal sites that may be applied to answering some of these questions and subsequently generating novel zinc enzymes. From examination of the earliest design studies introducing simple Zn(II)-binding sites into de novo and natural protein scaffolds to current studies involving the preparation of efficient hydrolytic zinc sites, it is increasingly likely that protein design will achieve reaction rates previously thought possible only for native enzymes. This Current Topic will review the design and redesign of Zn(II)-binding sites in de novo-designed proteins and native protein scaffolds toward the preparation of catalytic hydrolytic sites. After discussing the preparation of Zn(II)-binding sites in various scaffolds, we will describe relevant examples for reengineering existing zinc sites to generate new or altered catalytic activities. Then, we will describe our work on the preparation of a de novo-designed hydrolytic zinc site in detail and present comparisons to related designed zinc sites. Collectively, these studies demonstrate the significant progress being made toward building zinc metalloenzymes from the bottom up.

Crosstalk between endogenous and synthetic components--synthetic signaling meets endogenous components

Synthetic biology uses biological components to engineer new functionality in living organisms. We have used the tools of synthetic biology to engineer detector plants that can sense man-made chemicals, such as the explosive trinitrotoluene, and induce a response detectable by eye or instrumentation. A goal of this type of work is to make the designed system orthogonal, that is, able to function independently of systems in the host. In this review, the design and function of two partially synthetic signaling pathways for use in plants is discussed. We describe observed interactions (crosstalk) with endogenous signaling components. This crosstalk can be beneficial, allowing the creation of hybrid synthetic/endogenous signaling pathways, or detrimental, resulting in system noise and/or false positives. Current approaches in the field of synthetic biology applicable to the design of orthogonal signaling systems, including the design of synthetic components, partially synthetic systems that utilize crosstalk to signal through endogenous components, computational redesign of proteins, and the use of heterologous components, are discussed.

Engineering a zinc binding site into the de novo designed protein DS119 with a βαβ structure

Functional proteins designed de novo have potential application in chemical engineering, agriculture and healthcare. Metal binding sites are commonly used to incorporate functions. Based on a de novo designed protein DS119 with a βαβ structure, we have computationally engineered zinc binding sites into it using a home-made searching program. Seven out of the eight designed sequences tested were shown to bind Zn(2+) with micromolar affinity, and one of them bound Zn(2+) with 1:1 stoichiometry. This is the first time that metalloproteins with an α, β mixed structure have been designed from scratch.

Tunable signal processing in synthetic MAP kinase cascades

The flexibility of MAPK cascade responses enables regulation of a vast array of cell fate decisions, but elucidating the mechanisms underlying this plasticity is difficult in endogenous signaling networks. We constructed insulated mammalian MAPK cascades in yeast to explore how intrinsic and extrinsic perturbations affect the flexibility of these synthetic signaling modules. Contrary to biphasic dependence on scaffold concentration, we observe monotonic decreases in signal strength as scaffold concentration increases. We find that augmenting the concentration of sequential kinases can enhance ultrasensitivity and lower the activation threshold. Further, integrating negative regulation and concentration variation can decouple ultrasensitivity and threshold from the strength of the response. Computational analyses show that cascading can generate ultrasensitivity and that natural cascades with different kinase concentrations are innately biased toward their distinct activation profiles. This work demonstrates that tunable signal processing is inherent to minimal MAPK modules and elucidates principles for rational design of synthetic signaling systems.

Developing a synthetic signal transduction system in plants

One area of focus in the emerging field of plant synthetic biology is the manipulation of systems involved in sensing and response to environmental signals. Sensing and responding to signals, including ligands, typically involves biological signal transduction. Plants use a wide variety of signaling systems to sense and respond to their environment. One of these systems, a histidine kinase (HK) based signaling system, lends itself to manipulation using the tools of synthetic biology. Both plants and bacteria use HKs to relay signals, which in bacteria can involve as few as two proteins (two-component systems or TCS). HK proteins are evolutionarily conserved between plants and bacteria and plant HK components have been shown to be functional in bacteria. We found that this conservation also applies to bacterial HK components which can function in plants. This conservation of function led us to hypothesize that synthetic HK signaling components can be designed and rapidly tested in bacteria. These novel HK signaling components form the foundation for a synthetic signaling system in plants, but typically require modifications such as codon optimization and proper targeting to allow optimal function. We describe the process and methodology of producing a synthetic signal transduction system in plants. We discovered that the bacterial response regulator (RR) PhoB shows HK-dependent nuclear translocation in planta. Using this discovery, we engineered a partial synthetic pathway in which a synthetic promoter (PlantPho) is activated using a plant-adapted PhoB (PhoB-VP64) and the endogenous HK-based cytokinin signaling pathway. Building on this work, we adapted an input or sensing system based on bacterial chemotactic binding proteins and HKs, resulting in a complete eukaryotic signal transduction system. Input to our eukaryotic signal transduction system is provided by a periplasmic binding protein (PBP), ribose-binding protein (RBP). RBP interacts with the membrane-localized chemotactic receptor Trg. PBPs like RBP have been computationally redesigned to bind small ligands, such as the explosive 2,4,6-trinitrotoluene (TNT). A fusion between the chemotactic receptor Trg and the HK, PhoR, enables signal transduction via PhoB, which undergoes nuclear translocation in response to phosphorylation, resulting in transcriptional activation of an output gene under control of a synthetic plant promoter. Collectively, these components produce a novel ligand-responsive signal transduction system in plants and provide a means to engineer a eukaryotic synthetic signaling system.

Design strategies of fluorescent biosensors based on biological macromolecular receptors

Fluorescent biosensors to detect the bona fide events of biologically important molecules in living cells are increasingly demanded in the field of molecular cell biology. Recent advances in the development of fluorescent biosensors have made an outstanding contribution to elucidating not only the roles of individual biomolecules, but also the dynamic intracellular relationships between these molecules. However, rational design strategies of fluorescent biosensors are not as mature as they look. An insatiable request for the establishment of a more universal and versatile strategy continues to provide an attractive alternative, so-called modular strategy, which permits facile preparation of biosensors with tailored characteristics by a simple combination of a receptor and a signal transducer. This review describes an overview of the progress in design strategies of fluorescent biosensors, such as auto-fluorescent protein-based biosensors, protein-based biosensors covalently modified with synthetic fluorophores, and signaling aptamers, and highlights the insight into how a given receptor is converted to a fluorescent biosensor. Furthermore, we will demonstrate a significance of the modular strategy for the sensor design.

OptGraft: A computational procedure for transferring a binding site onto an existing protein scaffold

One of the many challenging tasks of protein design is the introduction of a completely new function into an existing protein scaffold. In this study, we introduce a new computational procedure OptGraft for placing a novel binding pocket onto a protein structure so as its geometry is minimally perturbed. This is accomplished by introducing a two-level procedure where we first identify where are the most appropriate locations to graft the new binding pocket into the protein fold by minimizing the departure from a set of geometric restraints using mixed-integer linear optimization. On identifying the suitable locations that can accommodate the new binding pocket, CHARMM energy calculations are employed to identify what mutations in the neighboring residues, if any, are needed to ensure that the minimum energy conformation of the binding pocket conserves the desired geometry. This computational framework is benchmarked against the results available in the literature for engineering a copper binding site into thioredoxin protein. Subsequently, OptGraft is used to guide the transfer of a calcium-binding pocket from thermitase protein (PDB: 1thm) into the first domain of CD2 protein (PDB:1hng). Experimental characterization of three de novo redesigned proteins with grafted calcium-binding centers demonstrated that they all exhibit high affinities for terbium (Kd) approximately 22, 38, and 55 microM) and can selectively bind calcium over magnesium.

Kinetic buffering of cross talk between bacterial two-component sensors

Two-component systems are a class of sensors that enable bacteria to respond to environmental and cell-state signals. The canonical system consists of a membrane-bound sensor histidine kinase that autophosphorylates in response to a signal and transfers the phosphate to an intracellular response regulator. Bacteria typically have dozens of two-component systems. The key questions are whether these systems are linear and, if they are, how cross talk between systems is buffered. In this work, we studied the EnvZ/OmpR and CpxA/CpxR systems from Escherichia coli, which have been shown previously to exhibit slow cross talk in vitro. Using in vitro radiolabeling and a rapid quenched-flow apparatus, we experimentally measured 10 biochemical parameters capturing the cognate and non-cognate phosphotransfer reactions between the systems. These data were used to parameterize a mathematical model that was used to predict how cross talk is affected as different genes are knocked out. It was predicted that significant cross talk between EnvZ and CpxR only occurs for the triple mutant DeltaompR DeltacpxA DeltaactA-pta. All seven combinations of these knockouts were made to test this prediction and only the triple mutant demonstrated significant cross talk, where the cpxP promoter was induced 280-fold upon the activation of EnvZ. Furthermore, the behavior of the other knockouts agrees with the model predictions. These results support a kinetic model of buffering where both the cognate bifunctional phosphatase activity and the competition between regulator proteins for phosphate prevent cross talk in vivo.

Catalytic mechanism and performance of computationally designed enzymes for Kemp elimination

A series of enzymes for Kemp elimination of 5-nitrobenzisoxazole has been recently designed and tested. In conjunction with the design process, extensive computational analyses were carried out to evaluate the potential performance of four of the designs, as presented here. The enzyme-catalyzed reactions were modeled using mixed quantum and molecular mechanics (QM/MM) calculations in the context of Monte Carlo (MC) statistical mechanics simulations. Free-energy perturbation (FEP) calculations were used to characterize the free-energy surfaces for the catalyzed reactions as well as for reference processes in water. The simulations yielded detailed information about the catalytic mechanisms, activation barriers, and structural evolution of the active sites over the course of the reactions. The catalytic mechanism for the designed enzymes KE07, KE10(V131N), and KE15 was found to be concerted with proton transfer, generally more advanced in the transition state than breaking of the isoxazolyl N-O bond. On the basis of the free-energy results, all three enzymes were anticipated to be active. Ideas for further improvement of the enzyme designs also emerged. On the technical side, the synergy of parallel QM/MM and experimental efforts in the design of artificial enzymes is well illustrated.

Origin of the activity drop with the E50D variant of catalytic antibody 34E4 for Kemp elimination

In enzymes, multiple structural effects cooperatively lead to the high catalytic activity, while individually these effects can be small. The design of artificial enzymes requires the understanding and ability to manipulate such subtle effects. The 34E4 catalytic antibody, catalyzing Kemp elimination of 5-nitrobenzisoxazole, and its Glu50Asp (E50D) variant are the subject of the present investigation. This removal of only a methylene group yields an approximately 30-fold reduction in the rate for the catalyzed Kemp elimination. Here, the aim is to understand this difference in the catalytic performance. The mechanism of Kemp elimination catalyzed by 34E4 and the E50D mutant is elucidated using QM/MM Monte Carlo simulations and free energy perturbation theory. In both proteins, the reaction is shown to follow a single-step, concerted mechanism. In the mutant, the activation barrier rises by 2.4 kcal/mol, which corresponds to a 62-fold rate deceleration, which is in good agreement with the experimental data. The positions and functionality of the residues in the active site are monitored throughout the reaction. It is concluded that the looser contact with the base, shorter base-Asn58 contact, less favorable pi-stacking with Trp91 in the transition state of the reaction, and different solvation pattern all contribute to the reduction of the reaction rate in the E50D variant of 34E4.

Peptide and protein building blocks for synthetic biology: from programming biomolecules to self-organized biomolecular systems

There are several approaches to creating synthetic-biological systems. Here, we describe a molecular-design approach. First, we lay out a possible synthetic-biology space, which we define with a plot of complexity of components versus divergence from nature. In this scheme, there are basic units, which range from natural amino acids to totally synthetic small molecules. These are linked together to form programmable tectons, for example, amphipathic alpha-helices. In turn, tectons can interact to give self-assembled units, which can combine and organize further to produce functional assemblies and systems. To illustrate one path through this vast landscape, we focus on protein engineering and design. We describe how, for certain protein-folding motifs, polypeptide chains can be instructed to fold. These folds can be combined to give structured complexes, and function can be incorporated through computational design. Finally, we describe how protein-based systems may be encapsulated to control and investigate their functions.

Biology by design: reduction and synthesis of cellular components and behaviour

Biological research is experiencing an increasing focus on the application of knowledge rather than on its generation. Thanks to the increased understanding of cellular systems and technological advances, biologists are more frequently asking not only 'how can I understand the structure and behaviour of this biological system?', but also 'how can I apply that knowledge to generate novel functions in different biological systems or in other contexts?' Active pursuit of the latter has nurtured the emergence of synthetic biology. Here, we discuss the motivation behind, and foundational technologies enabling, the development of this nascent field. We examine some early successes and applications while highlighting the challenges involved. Finally, we consider future directions and mention non-scientific considerations that can influence the field's growth.

Structure-based design of robust glucose biosensors using a Thermotoga maritima periplasmic glucose-binding protein

We report the design and engineering of a robust, reagentless fluorescent glucose biosensor based on the periplasmic glucose-binding protein obtained from Thermotoga maritima (tmGBP). The gene for this protein was cloned from genomic DNA and overexpressed in Escherichia coli, the identity of its cognate sugar was confirmed, ligand binding was studied, and the structure of its glucose complex was solved to 1.7 Angstrom resolution by X-ray crystallography. TmGBP is specific for glucose and exhibits high thermostability (midpoint of thermal denaturation is 119 +/- 1 degrees C and 144 +/- 2 degrees C in the absence and presence of 1 mM glucose, respectively). A series of fluorescent conjugates was constructed by coupling single, environmentally sensitive fluorophores to unique cysteines introduced by site-specific mutagenesis at positions predicted to be responsive to ligand-induced conformational changes based on the structure. These conjugates were screened to identify engineered tmGBPs that function as reagentless fluorescent glucose biosensors. The Y13C*Cy5 conjugate is bright, gives a large response to glucose over concentration ranges appropriate for in vivo monitoring of blood glucose levels (1-30 mM), and can be immobilized in an orientation-specific manner in microtiter plates to give a reversible response to glucose. The immobilized protein retains its response after long-term storage at room temperature.

As easy as flipping a switch?

Proteins that behave as switches help to establish the complex molecular logic that is central to biological systems. Aspiring to be nature's equal, researchers have successfully created protein switches of their own design; in particular, numerous and varied zinc-triggered switches have been made. Recent studies in which such switches have been readily identified from combinatorial protein libraries support the notion that proteins are primed to show allosteric behavior and that newly created ligand-binding sites will often be functionally coupled to the original activity of the protein. If true, this notion suggests that switch engineering might be more tractable than previously thought, boding well for the basic science, sensing and biomedical applications for which protein switches hold much promise.

Protein design: reengineering cellular retinoic acid binding protein II into a rhodopsin protein mimic

Rational redesign of the binding pocket of Cellular Retinoic Acid Binding Protein II (CRABPII) has provided a mutant that can bind retinal as a protonated Schiff base, mimicking the binding observed in rhodopsin. The reengineering was accomplished through a series of choreographed manipulations to ultimately orient the reactive species (the epsilon-amino group of Lys132 and the carbonyl of retinal) in the proper geometry for imine formation. The guiding principle was to achieve the appropriate Bürgi-Dunitz trajectory for the reaction to ensue. Through crystallographic analysis of protein mutants incapable of forming the requisite Schiff base, a highly ordered water molecule was identified as a key culprit in orienting retinal in a nonconstructive manner. Removal of the ordered water, along with placing reinforcing mutations to favor the desired orientation of retinal, led to a triple mutant CRABPII protein capable of nanomolar binding of retinal as a protonated Schiff base. The high-resolution crystal structure of all-trans-retinal bound to the CRABPII triple mutant (1.2 A resolution) unequivocally illustrates the imine formed between retinal and the protein.

The Use of Quantum Molecular Calculations to Guide a Genetic Algorithm: A Way to Search for New Chemistry

The process of gene-based molecular evolution has been simulated in silico by using massively parallel density functional theory quantum calculations, coupled with a genetic algorithm, to test for fitness with respect to a target chemical reaction in populations of genetically encoded molecules. The goal of this study was the identification of transition-metal complexes capable of mediating a known reaction, namely the cleavage of N(2) to give the metal nitride. Each complex within the search space was uniquely specified by a nanogene consisting of an eight-digit number. Propagation of an individual nanogene into successive generations was determined by the fitness of its phenotypic molecule to perform the target reaction and new generations were created by recombination and mutation of surviving nanogenes. In its simplest implementation, the quantum-directed genetic algorithm (QDGA) quickly located a local minimum on the evolutionary fitness hypersurface, but proved incapable of progressing towards the global minimum. A strategy for progressing beyond local minima consistent with the Darwinian paradigm by the use of environmental variations coupled with mass extinctions was therefore developed. This allowed for the identification of nitriding complexes that are very closely related to known examples from the chemical literature. Examples of mutations that appear to be beneficial at the genetic level but prove to be harmful at the phenotypic level are described. As well as revealing fundamental aspects of molecular evolution, QDGA appears to be a powerful tool for the identification of lead compounds capable of carrying out a target chemical reaction.

Local encoding of computationally designed enzyme activity

One aim of computational protein design is to introduce novel enzyme activity into proteins of known structure by predicting mutations that stabilize transition states. Previously, we showed that it is possible to introduce triose phosphate isomerase activity into the ribose-binding protein of Escherichia coli by constructing 17 mutations in the first two layers of residues that surround the wild-type ligand-binding site. Here, we report that these mutations can be "transplanted" into a homologous ribose-binding protein, isolated from the hyperthermophilic bacterium Thermoanaerobacter tengcongensis, with retention of catalytic activity, substrate affinity, and reaction pH dependence. The observed 10(5)-10(6)-fold rate enhancement corresponds to 70% of the maximally known transition-state binding energy. The wild-type sequences in these two homologues are almost perfectly conserved in the vicinity of their ribose-binding sites, but diverge significantly at increasing distance from these sites. The results demonstrate that the computationally designed mutations are sufficient to encode the observed enzyme activity, that all the observed activity is encoded locally within the layer of residues directly in contact with the substrate and that, in this case, at least 70% of transition state stabilization energy can be achieved using straightforward considerations of stereochemical complementarity between enzyme and reactants.

Role of Na+and K+in Enzyme Function

Metal complexation is a key mediator or modifier of enzyme structure and function. In addition to divalent and polyvalent metals, group IA metals Na+and K+play important and specific roles that assist function of biological macromolecules. We examine the diversity of monovalent cation (M+)-activated enzymes by first comparing coordination in small molecules followed by a discussion of theoretical and practical aspects. Select examples of enzymes that utilize M+as a cofactor (type I) or allosteric effector (type II) illustrate the structural basis of activation by Na+and K+, along with unexpected connections with ion transporters. Kinetic expressions are derived for the analysis of type I and type II activation. In conclusion, we address evolutionary implications of Na+binding in the trypsin-like proteases of vertebrate blood coagulation. From this analysis, M+complexation has the potential to be an efficient regulator of enzyme catalysis and stability and offers novel strategies for protein engineering to improve enzyme function.

Genetic parts to program bacteria

Genetic engineering is entering a new era, where microorganisms can be programmed using synthetic constructs of DNA encoding logic and operational commands. A toolbox of modular genetic parts is being developed, comprised of cell-based environmental sensors and genetic circuits. Systems have already been designed to be interconnected with each other and interfaced with the control of cellular processes. Engineering theory will provide a predictive framework to design operational multicomponent systems. On the basis of these developments, increasingly complex cellular machines are being constructed to build specialty chemicals, weave biomaterials, and to deliver therapeutics.

Guanidinium derivatives bind preferentially and trigger long-distance conformational changes in an engineered T4 lysozyme

The binding of guanidinium ion has been shown to promote a large-scale translation of a tandemly duplicated helix in an engineered mutant of T4 lysozyme. The guanidinium ion acts as a surrogate for the guanidino group of an arginine side chain. Here we determine whether methyl- and ethylguanidinium provide better mimics. The results show that addition of the hydrophobic moieties to the ligand enhances the binding affinity concomitant with reduction in ligand solubility. Crystallographic analysis confirms that binding of the alternative ligands to the engineered site still drives the large-scale conformational change. Thermal analysis and NMR data show, in comparison to guanidinium, an increase in protein stability and in ligand affinity. This is presumably due to the successive increase in hydrophobicity in going from guanidinium to ethylguanidinium. A fluorescence-based optical method was developed to sense the ligand-triggered helix translation in solution. The results are a first step in the de novo design of a molecular switch that is not related to the normal function of the protein.

Knowledge-based potentials in protein design

Knowledge-based potentials are statistical parameters derived from databases of known protein properties that empirically capture aspects of the physical chemistry of protein structure and function. These potentials play a key role in protein design by improving the accuracy of physics-based models of interatomic interactions and enhancing the computational efficiency of the design process by limiting the complexity of searching sequence space. Recently, knowledge-based potentials (in isolation or in combination with physics-based potentials) have been applied to the modification of existing protein function, the redesign of natural protein folds and the complete design of a non-natural protein fold. In addition, knowledge-based potentials appear to be providing important information about the global topology of amino acid interactions in natural proteins. A detailed study of the methods and products of these protein design efforts promises to greatly expand our understanding of proteins and the evolutionary process that created them.

Identification of cognate ligands for the Escherichia coli phnD protein product and engineering of a reagentless fluorescent biosensor for phosphonates

The Escherichia coli phnD gene is hypothesized to code for the periplasmic binding component of a phosphonate uptake system. Here we report the characterization of the phosphonate-binding properties of the phnD protein product. We find that PhnD exhibits high affinity for 2-aminoethylphosphonate (5 nM), the most commonly occurring natural phosphonate produced by lower eukaryotes, but also binds several other phosphonates with micromolar affinities. A significant number of man-made phosphonates, such as insecticides and chemical warfare agents, are chemical threats and environmental pollutants. Consequently, there is an interest in developing methods for the detection and bioremediation of phosphonates. Bacterial periplasmic-binding proteins have been utilized for developing reagentless biosensors that report analytes by coupling ligand-binding events to changes in the emission properties of a covalently conjugated environmentally-sensitive fluorophore. Several PhnD conjugates described here show large changes in fluorescence upon binding to methylphosphonate (MP), with two conjugates exhibiting up to 50% decrease in emission intensity. Since MP is the final degradation product of many nerve agents, these PhnD conjugates can function as components in a biosensor system for chemical warfare agents.

IPRO: an iterative computational protein library redesign and optimization procedure

A number of computational approaches have been developed to reengineer promising chimeric proteins one at a time through targeted point mutations. In this article, we introduce the computational procedure IPRO (iterative protein redesign and optimization procedure) for the redesign of an entire combinatorial protein library in one step using energy-based scoring functions. IPRO relies on identifying mutations in the parental sequences, which when propagated downstream in the combinatorial library, improve the average quality of the library (e.g., stability, binding affinity, specific activity, etc.). Residue and rotamer design choices are driven by a globally convergent mixed-integer linear programming formulation. Unlike many of the available computational approaches, the procedure allows for backbone movement as well as redocking of the associated ligands after a prespecified number of design iterations. IPRO can also be used, as a limiting case, for the redesign of a single or handful of individual sequences. The application of IPRO is highlighted through the redesign of a 16-member library of Escherichia coli/Bacillus subtilis dihydrofolate reductase hybrids, both individually and through upstream parental sequence redesign, for improving the average binding energy. Computational results demonstrate that it is indeed feasible to improve the overall library quality as exemplified by binding energy scores through targeted mutations in the parental sequences.

Molego-based definition of the architecture and specificity of metal-binding sites

Decomposing proteins into "molegos," building blocks that are conserved in sequence and 3D-structure, can identify functional elements. To demonstrate the specificity of the decomposition method, the PCPMer program suite was used to numerically define physical chemical property motifs corresponding to the molegos that make up the metal-containing active sites of three distinct enzyme families, from the dimetallic phosphatases, DNase 1 related nucleases/phosphatases, and dioxygenases. All three superfamilies bind metal ions in a beta-strand core region but differ in the number and type of ions needed for activity. The motifs were then used to automatically identify proteins in the ASTRAL40 database that contained similar motifs. The proteins with the highest PCPMer score in the database were primarily metal-binding enzymes that were related in function to those in the alignment used to generate the PCPMer motif lists. The proteins that contained motifs similar to the dioxygenases differed from those found with PCP-motifs for phosphatases and nucleases. Relatively few metal-binding enzymes were detected when the search was done with PCP-motifs defined for interleukin-1 related proteins, which have a beta-strand core but do not bind metal ions. While the box architecture was constant in each superfamily, the specificity for the metal ion preferred for enzymatic activity is determined by the pattern of carbonyl, hydroxyl or imadazole groups in key positions in the molegos. These results have implications for the design of metal-binding enzymes, and illustrate the ability of the PCPMer approach to distinguish, at the sequence level, structural and functional elements.

Intelligent design: the de novo engineering of proteins with specified functions

One of the principal successes of de novo protein design has been the creation of small, robust protein-cofactor complexes which can serve as simplified models, or maquettes, of more complicated multicofactor protein complexes commonly found in nature. Different maquettes, generated by us and others, recreate a variety of aspects, or functional elements, recognized as parts of natural enzyme function. The current challenge is to both expand the palette of functional elements and combine and/or integrate them in recreating familiar enzyme activities or generating novel catalysis in the simplest protein scaffolds.

ZiCo: A Peptide Designed to Switch Folded State upon Binding Zinc

We describe a novel approach to the design of a metal-triggered conformational switch. Specifically, two distinct protein-folding motifs were merged into one polypeptide sequence. The target structures were an alpha-helical coiled-coil trimer and zinc-bound monomer. Solution-phase spectroscopic, sedimentation, and binding studies confirmed the key aspects of the design. Both forms of the peptide were cooperatively folded, and the switch between them was reversible. This design process potentially presents a novel route to peptide-based biosensors.

Studying zinc biology with fluorescence: ain’t we got fun?

Zinc has emerged as a metal ion of substantial interest in biology and medicine, especially in neuroscience, gene transcription, the immune response, and mammalian reproduction. Fueling these advances in understanding has been the development of new fluorescence-based indicator systems for zinc with unprecedented sensitivity and selectivity. This review summarizes recent progress in the development of fluorescence-based sensors and biosensors for zinc, with a view to evaluating their suitability for use with biologically derived specimens, especially in vivo and in situ.

Analysis of allosteric signal transduction mechanisms in an engineered fluorescent maltose biosensor

We previously reported the construction of a family of reagentless fluorescent biosensor proteins by the structure-based design of conjugation sites for a single, environmentally sensitive small molecule dye, thus providing a mechanism for the transduction of ligand-induced conformational changes into a macroscopic fluorescence observable. Here we investigate the microscopic mechanisms that may be responsible for the macroscopic fluorescent changes in such Fluorescent Allosteric Signal Transduction (FAST) proteins. As case studies, we selected three individual cysteine mutations (F92C, D95C, and S233C) of Escherichia coli maltose binding protein (MBP) covalently labeled with a single small molecule fluorescent probe, N-((2-iodoacetoxy)ethyl)-N-methyl)amino-7-nitrobenz-2-oxa-1,3-diazole (NBD), each giving rise to a robust FAST protein with a distinct maltose-dependent fluorescence response. The fluorescence emission intensity, anisotropy, lifetime, and iodide-dependent fluorescence quenching were determined for each conjugate in the presence and absence of maltose. Structure-derived solvent accessible surface areas of the three FAST proteins are consistent with experimentally observed quenching data. The D95C protein exhibits the largest fluorescence change upon maltose binding. This mutant was selected for further characterization, and residues surrounding the fluorophore coupling site were mutagenized. Analysis of the resulting mutant FAST proteins suggests that specific hydrogen-bonding interactions between the fluorophore molecule and two tyrosine side-chains, Tyr171 and Tyr176, in the open state but not the closed, are responsible for the dramatic fluorescence response of this construct. Taken together these results provide insights that can be used in future design cycles to construct fluorescent biosensors that optimize signaling by engineering specific hydrogen bonds between a fluorophore and protein.

Engineering of non-natural receptors

Rational design, usually guided by computational prediction, and selection from libraries of variants of natural proteins have been used with success in the engineering of novel non-natural receptors. Many of these engineered protein binders will find use in biotechnological, diagnostic and medical applications, sometimes in the place of natural antibodies.

Periplasmic binding proteins: a versatile superfamily for protein engineering

The diversity of biological function, ligand binding, conformational changes and structural adaptability of the periplasmic binding protein superfamily have been exploited to engineer biosensors, allosteric control elements, biologically active receptors and enzymes using a combination of techniques, including computational design. Extensively redesigned periplasmic binding proteins have been re-introduced into bacteria to function in synthetic signal transduction pathways that respond to extracellular ligands and as biologically active enzymes.

Use of sequence duplication to engineer a ligand-triggered, long-distance molecular switch in T4 lysozyme

We have designed a molecular switch in a T4 lysozyme construct that controls a large-scale translation of a duplicated helix. As shown by crystal structures of the construct with the switch on and off, the conformational change is triggered by the binding of a ligand (guanidinium ion) to a site that in the wild-type protein was occupied by the guanidino head group of an Arg. In the design template, a duplicated helix is flanked by two loop regions of different stabilities. In the "on" state, the N-terminal loop is weakly structured, whereas the C-terminal loop has a well defined conformation that is stabilized by means of nonbonded interactions with the Arg head group. The truncation of the Arg to Ala destabilizes this loop and switches the protein to the "off" state, in which the duplicated helix is translocated approximately 20 A. Guanidinium binding restores the key interactions, restabilizes the C-terminal loop, and restores the "on" state. Thus, the presence of an external ligand, which is unrelated to the catalytic activity of the enzyme, triggers the inserted helix to translate 20 A away from the binding site. The results illustrate a proposed mechanism for protein evolution in which sequence duplication followed by point mutation can lead to the establishment of new function.

Computational Design of a Biologically Active Enzyme

Rational design of enzymes is a stringent test of our understanding of protein chemistry and has numerous potential applications. Here, we present and experimentally validate the computational design of enzyme activity in proteins of known structure. We have predicted mutations that introduce triose phosphate isomerase activity into ribose-binding protein, a receptor that normally lacks enzyme activity. The resulting designs contain 18 to 22 mutations, exhibit 10(5)- to 10(6)-fold rate enhancements over the uncatalyzed reaction, and are biologically active, in that they support the growth of Escherichia coli under gluconeogenic conditions. The inherent generality of the design method suggests that many enzymes can be designed by this approach.

Computational design of receptors for an organophosphate surrogate of the nerve agent soman

We report the computational design of soluble protein receptors for pinacolyl methyl phosphonic acid (PMPA), the predominant hydrolytic product of the nerve agent soman. Using recently developed computational protein design techniques, the ligand-binding pockets of two periplasmic binding proteins, glucose-binding protein and ribose-binding protein, were converted to bind PMPA instead of their cognate sugars. The designs introduce 9-12 mutations in the parent proteins. Twelve of 20 designs tested exhibited PMPA-dependent changes in emission intensity of a fluorescent reporter with affinities between 45 nM and 10 microM. The contributions to ligand binding by individual residues were determined in two designs by alanine-scanning mutagenesis, and are consistent with the molecular models. These results demonstrate that designed receptors with radically altered binding specificities and affinities that rival or exceed those of the parent proteins can be successfully predicted. The designs vary in parent scaffold, sequence diversity, and orientation of docked ligand, suggesting that the number of possible solutions to the design problem is large and degenerate. This observation has implications for the genesis of biological function by random processes. The designed receptors reported here may have utility in the development of fluorescent biosensors for monitoring nerve agents.

Cobalt(II) and zinc(II) binding to a ferredoxin maquette

We have examined the Co(II) and Zn(II) affinity of the prototype ferredoxin maquette ligand, NH(2)-KLCEGG.CIACGAC.GGW-CONH2 (IAA), which was originally designed to bind a [4Fe-4S] cluster. UV-Vis spectroscopy demonstrates tight 1:1 complex formation between Co(II) and IAA. The intensity of the S-->Co(II) charge transfer bands at 304 and 340 nm and the ligand field bands between 630 and 728 nm indicate Co(II) coordination by the four cysteine thiolates of IAA in a pseudo-tetrahedral geometry. A dissociation constant value of 5.3 microM was determined for the Co(II)-IAA complex at pH 6.5. Zn(II) readily displaces Co(II) from IAA as evinced by loss of the Co(II) spectral features. The dissociation constant for Zn(II), 20 pM at pH 6.5, was determined be competition experiments with Co(II)-IAA. These results demonstrate that the ferredoxin maquette ligand is an excellent ligand for Zn(II).