Design principles of protein switches

Protein switches perform essential roles in many biological processes and are exciting targets for de novo protein design, which aims to produce proteins of arbitrary shape and functionality. However, the biophysical requirements for switch function — multiple conformational states, fine-tuned energetics, and stimuli-responsiveness — pose a formidable challenge for design by computation (or intuition). A variety of methods have been developed toward tackling this challenge, usually taking inspiration from the wealth of sequence and structural information available for naturally occurring protein switches. More recently, modular switches have been designed computationally, and new methods have emerged for sampling unexplored structure space, providing promising new avenues toward the generation of purpose-built switches and de novo signaling systems for cellular engineering.     

Designing redox potential-controlled protein switches based on mutually exclusive proteins

Synthetic/artificial protein switches provide an efficient means of controlling protein functions using chemical signals and stimuli. Mutually exclusive proteins, in which only the host or guest domain can remain folded at a given time owing to conformational strain, have been used to engineer novel protein switches that can switch enzymatic functions on and off in response to ligand binding. To further explore the potential of mutually exclusive proteins as protein switches and sensors, we report here a new redox-based approach to engineer a mutually exclusive folding-based protein switch. By introducing a disulfide bond into the host domain of a mutually exclusive protein, we demonstrate that it is feasible to use redox potential to switch the host domain between its folded and unfolded conformations via the mutually exclusive folding mechanism, and thus switching the functionality of the host domain on and off. Our study opens a new and potentially general avenue that uses mutually exclusive proteins to design novel switches able to control the function of a variety of proteins.     

A divalent metal-mediated switch controlling protein-induced DNA bending

Architectural proteins that reconfigure the paths of DNA segments are required for the establishment of functional interfaces in many genomic transactions. A single-chain derivative of the DNA architectural protein integration host factor was found to adopt two stable conformational states in complex with a specific DNA target. In the so-called open state, the degree of protein-induced DNA bending is reduced significantly compared with the closed state. The conformational switch between these states is controlled by divalent metal binding in two electronegative zones arising from the lysine-to-glutamate substitution in the protein body proximal to the phosphate backbone of one DNA arm. We show that this switch can be employed to control the efficiency of site-specific recombination catalyzed by lambda integrase. Introduction of acidic residues at the protein-DNA interface holds potential for the design of metal-mediated switches for the investigation of functional relationships.     

The stability of a stochastic CaMKII switch: dependence on the number of enzyme molecules and protein turnover

Molecular switches have been implicated in the storage of information in biological systems. For small structures such as synapses, these switches are composed of only a few molecules and stochastic fluctuations are therefore of importance. Such fluctuations could potentially lead to spontaneous switch reset that would limit the lifetime of information storage. We have analyzed a model of the calcium/calmodulin-dependent protein kinase II (CaMKII) switch implicated in long-term memory in the nervous system. The bistability of this switch arises from autocatalytic autophosphorylation of CaMKII, a reaction that is countered by a saturable phosphatase-1-mediated dephosphorylation. We sought to understand the factors that control switch stability and to determine the functional relationship between stability and the number of molecules involved. Using Monte Carlo simulations, we found that the lifetime of states of the switch increase exponentially with the number of CaMKII holoenzymes. Switch stability requires a balance between the kinase and phosphatase rates, and the kinase rate must remain high relative to the rate of protein turnover. Thus, a critical limit on switch stability is set by the observed turnover rate (one per 30 h on average). Our computational results show that, depending on the timescale of fluctuations in enzyme numbers, for a switch composed of about 15 CaMKII holoenzymes, the stable persistent activation can span from a few years to a human lifetime.     

A nucleic acid switch triggered by the HIV-1 nucleocapsid protein

A unimolecular oligonucleotide switch, termed here an AlloSwitch, binds the mature HIV-1 nucleocapsid protein, NCp7. This switch can be used as an indicator for the presence of free NCp7 and NC domains in precursor and fusion proteins. It is thermodynamically stable in two conformations, H and O. A FRET pair is covalently attached to the strands to report on the molecular state of the switch. The results show that NC has an affinity for O 170 times higher than its affinity for H and that in the absence of NC the equilibrium ratio K1 = [O]/[H] = 0.10 +/- 0.03 for the switch sequence reported here. The change between the two states happens on a rapid kinetic time scale. A framework is introduced to aid in the design of AlloSwitches aimed at other targets. A high-affinity probe segment must be available to bind the target in the O-form, while a cover segment hides the probe in H. A key is adjusting the cover sequence to favor the H-form by a factor of 10-1000. This affords a robust response to small changes in target concentration, while saturation produces more than 90% of the maximal change in fluorescence. When a competitor displaces the switch from the NC-O complex, the released switch reverts to the H-form. This is the basis for a mix-and-read strategy for high-throughput screening of anti-nucleocapsid drug candidates that is much simpler to execute than traditional assays that require immobilization and washing steps.     

Design of multispecific protein sequences using probabilistic graphical modeling

In nature, proteins partake in numerous protein– protein interactions that mediate their functions. Moreover, proteins have been shown to be physically stable in multiple structures, induced by cellular conditions, small ligands, or covalent modifications. Understanding how protein sequences achieve this structural promiscuity at the atomic level is a fundamental step in the drug design pipeline and a critical question in protein physics. One way to investigate this subject is to computationally predict protein sequences that are compatible with multiple states, i.e., multiple target structures or binding to distinct partners. The goal of engineering such proteins has been termed multispecific protein design. We develop a novel computational framework to efficiently and accurately perform multispecific protein design. This framework utilizes recent advances in probabilistic graphical modeling to predict sequences with low energies in multiple target states. Furthermore, it is also geared to specifically yield positional amino acid probability profiles compatible with these target states. Such profiles can be used as input to randomly bias high‐throughput experimental sequence screening techniques, such as phage display, thus providing an alternative avenue for elucidating the multispecificity of natural proteins and the synthesis of novel proteins with specific functionalities. We prove the utility of such multispecific design techniques in better recovering amino acid sequence diversities similar to those resulting from millions of years of evolution. We then compare the approaches of prediction of low energy ensembles and of amino acid profiles and demonstrate their complementarity in providing more robust predictions for protein design. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

A fast method for predicting amino acid mutations that lead to unfolding

Amino acid mutation(s) that cause(s) partial or total unfolding of a protein can lead to disease states and failure to produce mutants. It is therefore very useful to be able to predict which mutations can retain the conformation of a wild-type protein and which mutations will lead to local or global unfolding of the protein. We have developed a fast and reasonably accurate method based on a backbone-dependent side-chain rotamer library to predict the (folded or unfolded) conformation of a protein upon mutation. This method has been tested on proteins whose wild-type 3D structures are known and whose mutant conformations have been experimentally characterized to be folded or unfolded. Furthermore, for the cases studied here, the predicted partially folded or denatured mutant conformation correlate with a decrease in the stability of the mutant relative to the wild-type protein. The key advantage of our method is that it is very fast and predicts locally or globally unfolded states fairly accurately. Hence, it may prove to be useful in designing site-directed mutagenesis, X-ray crystallography and drug design experiments as well as in free energy simulations by helping to ascertain whether a mutation will alter or retain the wild-type conformation.     

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.     

Modeling genetic switches with positive feedback loops

In this paper, we develop a new methodology to design synthetic genetic switch networks with multiple genes and time delays, by using monotone dynamical systems. We show that the networks with only positive feedback loops have no stable oscillation but stable equilibria whose stability is independent of the time delays. In other words, such systems have ideal properties for switch networks and can be designed without consideration of time delays, because the systems can be reduced from functional spaces to Euclidian spaces. Therefore, we can ensure that the designed switches function correctly even with uncertain delays. We first prove the basic properties of the genetic networks composed of only positive feedback loops, and then propose a procedure to design the switches, which drastically simplifies analysis of the switches and makes theoretical analysis and design tractable even for large-scaled systems. Finally, to demonstrate our theoretical results, we show biologically plausible examples by designing a synthetic genetic switch with experimentally well investigated lacI, tetR, and cI genes for numerical simulation.     

Protein folding: a problem with multiple solutions

There is continued interest in predicting the structure of proteins either at the simplest level of identifying their fold class or persevering all the way to an atomic resolution structure. Protein folding methods have become very sophisticated and many successes have been recorded with claims to have solved the native structure of the protein. But for any given protein, there may be more than one solution. Many proteins can exist in one of the other two (or more) different forms and some populate multiple metastable states. Here, the two-state case is considered and the key structural changes that take place when the protein switches from one state to the other are identified. Analysis of these results show that hydrogen bonding patterns and hydrophobic contacts vary considerably between different conformers. Contrary to what has often been assumed previously, these two types of interaction operate essentially independently of one another. Core packing is critical for proper protein structure and function and it is shown that there are considerable changes in internal cavity volumes in many cases. The way in which these switches are made is fold dependent. Considerations such as these need to be taken into account in protein structure prediction.     

Off the beaten paths: alternative and crosstalk regulation of Rho GTPases

Rho proteins are small GTPases of the Ras superfamily that regulate a wide variety of biological processes, ranging from gene expression to cell migration. Mechanistically, the major Rho GTPases function as molecular switches cycling between an inactive GDP-bound and an active GTP-bound conformation, although several Rho proteins spontaneously exchange nucleotides or are simply devoid of GTPase activity. For over a decade, RhoGEFs and RhoGAPs have been established as the mainstream regulators of Rho proteins, respectively flipping the switch on or off. However, regulation by GEFs and GAPs leaves several fundamental questions on the operation of the Rho switch unanswered, indicating that the regulation of Rho proteins does not rely exclusively on RhoGEFs and RhoGAPs. Recent evidence indeed suggests that Rho GTPases are finely tuned by multiple alternative regulatory mechanisms, including post-translational modifications and protein degradation, as well as crosstalk mechanisms between Rho proteins. Here we review these alternative mechanisms and discuss how they alter Rho protein function and signaling. We also envision how the classic binary Rho switch may indeed function more like a switchboard with multiple switches and dials that can all contribute to the regulation of Rho protein function.     

Model‐guided design of ligand‐regulated RNAi for programmable control of gene expression

Progress in constructing biological networks will rely on the development of more advanced components that can be predictably modified to yield optimal system performance. We have engineered an RNA‐based platform, which we call an shRNA switch, that provides for integrated ligand control of RNA interference (RNAi) by modular coupling of an aptamer, competing strand, and small hairpin (sh)RNA stem into a single component that links ligand concentration and target gene expression levels. A combined experimental and mathematical modelling approach identified multiple tuning strategies and moves towards a predictable framework for the forward design of shRNA switches. The utility of our platform is highlighted by the demonstration of fine‐tuning, multi‐input control, and model‐guided design of shRNA switches with an optimized dynamic range. Thus, shRNA switches can serve as an advanced component for the construction of complex biological systems and offer a controlled means of activating RNAi in disease therapeutics.     

Modular enzyme design: regulation by mutually exclusive protein folding

A regulatory mechanism is introduced whereupon the catalytic activity of a given enzyme is controlled by ligand binding to a receptor domain of choice. A small enzyme (barnase) and a ligand-binding polypeptide (GCN4) are fused so that a simple topological constraint prevents them from existing simultaneously in their folded states. The two domains consequently engage in a thermodynamic tug-of-war in which the more stable domain forces the less stable domain to unfold. In the absence of ligand, the barnase domain is more stable and is therefore folded and active; the GCN4 domain is substantially unstructured. DNA binding induces folding of GCN4, forcibly unfolding and inactivating the barnase domain. Barnase-GCN4 is thus a "natively unfolded" protein that uses ligand binding to switch between partially folded forms. The key characteristics of each parent protein (catalytic efficiency of barnase, DNA binding affinity and sequence specificity of GCN4) are retained in the chimera. Barnase-GCN4 thus defines a modular approach for assembling enzymes with novel sensor capabilities from a variety of catalytic and ligand binding domains.     

A generic program for multistate protein design

Some protein design tasks cannot be modeled by the traditional single state design strategy of finding a sequence that is optimal for a single fixed backbone. Such cases require multistate design, where a single sequence is threaded onto multiple backbones (states) and evaluated for its strengths and weaknesses on each backbone. For example, to design a protein that can switch between two specific conformations, it is necessary to to find a sequence that is compatible with both backbone conformations. We present in this paper a generic implementation of multistate design that is suited for a wide range of protein design tasks and demonstrate in silico its capabilities at two design tasks: one of redesigning an obligate homodimer into an obligate heterodimer such that the new monomers would not homodimerize, and one of redesigning a promiscuous interface to bind to only a single partner and to no longer bind the rest of its partners. Both tasks contained negative design in that multistate design was asked to find sequences that would produce high energies for several of the states being modeled. Success at negative design was assessed by computationally redocking the undesired protein-pair interactions; we found that multistate design's accuracy improved as the diversity of conformations for the undesired protein-pair interactions increased. The paper concludes with a discussion of the pitfalls of negative design, which has proven considerably more challenging than positive design.     

Rapid DNA detection by beacon-assisted detection amplification

This protocol describes a new and rapid isothermal reaction process designed to amplify and detect a specific DNA sequence in purified DNA extracted from cultured cells. The protocol uses a DNA nanomachine that comprises two molecular switches that function in concert to isothermally amplify and detect a DNA target. First, a molecular beacon detection switch is 'activated' only if a DNA target sequence is present. A DNA primer and DNA polymerase are used to lock the beacon in an activated conformation. Second, an amplification and signal-transduction switch is initiated following successful activation. A nicking endonuclease and the DNA polymerase are used to replicate the DNA target. Both switches operate simultaneously at 40 °C in a single reaction to rapidly generate multiple copies of the DNA target in a cyclic polymerization reaction. This protocol enables femtomole amounts of a DNA target to be reproducibly amplified and detected in <40 min. We demonstrate the successful use of this protocol in assays containing synthetic DNA components and purified DNA extracted from biological samples.     

How to make a biological switch

Some biological regulatory systems must "remember" a state for long periods of time. A simple type of system that can accomplish this task is one in which two regulatory elements negatively regulate one another. For example, two repressor proteins might control one another's synthesis. Qualitative reasoning suggests that such a system will have two stable states, one in which the first element is "on" and the second "off", and another in which these states are reversed. Quantitative analysis shows that the existence of two stable steady states depends on the details of the system. Among other things, the shapes of functions describing the effect of one regulatory element on the other must meet certain criteria in order for two steady states to exist. Many biologically reasonable functions do not meet these criteria. In particular, repression that is well described by a Michaelis-Menten-type equation cannot lead to a working switch. However, functions describing positive cooperativity of binding, non-additive effects of multiple operator sites, or depletion of free repressor can lead to working switches.     

Design and construction of a double inversion recombination switch for heritable sequential genetic memory

Background

Inversion recombination elements present unique opportunities for computing and information encoding in biological systems. They provide distinct binary states that are encoded into the DNA sequence itself, allowing us to overcome limitations posed by other biological memory or logic gate systems. Further, it is in theory possible to create complex sequential logics by careful positioning of recombinase recognition sites in the sequence.

Methodology/Principal Findings

In this work, we describe the design and synthesis of an inversion switch using the fim and hin inversion recombination systems to create a heritable sequential memory switch. We have integrated the two inversion systems in an overlapping manner, creating a switch that can have multiple states. The switch is capable of transitioning from state to state in a manner analogous to a finite state machine, while encoding the state information into DNA. This switch does not require protein expression to maintain its state, and “remembers” its state even upon cell death. We were able to demonstrate transition into three out of the five possible states showing the feasibility of such a switch.

Conclusions/Significance

We demonstrate that a heritable memory system that encodes its state into DNA is possible, and that inversion recombination system could be a starting point for more complex memory circuits. Although the circuit did not fully behave as expected, we showed that a multi-state, temporal memory is achievable.     

Characterizing multistationarity regimes in biochemical reaction networks

Switch like responses appear as common strategies in the regulation of cellular systems. Here we present a method to characterize bistable regimes in biochemical reaction networks that can be of use to both direct and reverse engineering of biological switches. In the design of a synthetic biological switch, it is important to study the capability for bistability of the underlying biochemical network structure. Chemical Reaction Network Theory (CRNT) may help at this level to decide whether a given network has the capacity for multiple positive equilibria, based on their structural properties. However, in order to build a working switch, we also need to ensure that the bistability property is robust, by studying the conditions leading to the existence of two different steady states. In the reverse engineering of biological switches, knowledge collected about the bistable regimes of the underlying potential model structures can contribute at the model identification stage to a drastic reduction of the feasible region in the parameter space of search. In this work, we make use and extend previous results of the CRNT, aiming not only to discriminate whether a biochemical reaction network can exhibit multiple steady states, but also to determine the regions within the whole space of parameters capable of producing multistationarity. To that purpose we present and justify a condition on the parameters of biochemical networks for the appearance of multistationarity, and propose an efficient and reliable computational method to check its satisfaction through the parameter space.